Structured Total Least Squares for Volumes of Distribution Estimation in PET Imaging

Hongbin Guo,Kewei Chen,Rosemary Renaut,Eric M. Reiman
DOI: https://doi.org/10.1016/j.jalz.2010.05.1451
2010-01-01
Abstract:The volume distribution (DV), an index for amyloid deposition, is routinely estimated by the Graphical analysis (GA) for PiB (Pittsburgh Compound-B)- Positron emission tomography (PET) studies of Alzheimer's disease (AD). We found that the well-known bias of the GA estimated DV is primarily due to model error and noise in the PET data. We recently developed an alternative model to minimize the bias due to model error. In this study, we propose a procedure to minimize the bias due to the measurement noise. A commonly used method to account of the noise is the simple linear least squares (LLS) which assume the noise is non-correlated and identically distributed. The noise in the GA Logan's equation is actually structured and non-identical. Therefore, we propose to use structured total least squares (STLS) instead of the LLS to estimate DV. Ignoring the model error temporarily in this study, a efficient and effective numerical algorithm was derived for DV estimation. Simulated PiB PET data from AD related brain regions in patients and normal controls to examine the proposed algorithm. 1000 random samples for each of 11 regions of interest (ROI) are tested. The simulation results show that the proposed method significantly reduces the DV errors from -15.7% ± 4.4% to -3.9% ± 10.5 as compared with LLS. We also note that STLS is much better than the simple total least squares (TLS) and MA1, which includes huge outliers. Statistical descriptions for three representative ROIs are listed in the attached table while histograms for for normalized error (in percentage) for all eleven ROIs and four methods, LGA, MA1, TLS and STLS are presented in the attached figure. We conclude that the proposed STLS method is appropriate for DV estimation in AD PiB PET studies. Additional studied are under-going to incorporate procedures for bias correction related to model error and noise structure.
What problem does this paper attempt to address?