Simulated Annealing in pharmacokinetic modeling of PET neuroreceptor studies: Accuracy and precision compared with other optimization algorithms

Maqsood Yaqub, Ronald Boellaard, Marc A. Kropholler, Mark Lubberink, Adriaan A. Lammertsma

Research output: Contribution to conferencePaperAcademic

Abstract

Fitting of PET pharmacokinetic parameters may suffer from bias or unrealistic outcomes, especially for noisy data. There are many readily available optimization algorithms, but each has different characteristics when fitting PET pharmacokinetic models. The purpose of this study is to evaluate the performance of four different types of optimization algorithms, including a modified simulated annealing method, in terms of precision and accuracy of PET pharmacokinetic parameters. The simulated annealing algorithm (SA), called Basin-hoping, was modified for present application. Input data, taken from [ 11C]-PK11195 neuroreceptor ligand studies, was used to simulate time activity curves at various noise levels. Also the influence of incorrect weighting factors on algorithm performance was studied. Surprisingly, effects of using incorrect but reasonable weighting factors on bias and precision were negligible. Except when extreme and unrealistic weighting factors were used, an increase in bias and decrease in precision was observed. In general the modified SA provided smallest weighted squared residual error and was able to find the global minimum without the need for a proper start parameter selection. However, occasionally better fits were found with the interior-reflective Newton method, but only when implemented using a range of start parameters, centered on the expected value.

Original languageEnglish
Pages3222-3225
Number of pages4
Publication statusPublished - 1 Dec 2004
Event2004 Nuclear Science Symposium, Medical Imaging Conference, Symposium on Nuclear Power Systems and the 14th International Workshop on Room Temperature Semiconductor X- and Gamma- Ray Detectors - Rome, Italy
Duration: 16 Oct 200422 Oct 2004

Conference

Conference2004 Nuclear Science Symposium, Medical Imaging Conference, Symposium on Nuclear Power Systems and the 14th International Workshop on Room Temperature Semiconductor X- and Gamma- Ray Detectors
CountryItaly
CityRome
Period16/10/200422/10/2004

Cite this

Yaqub, M., Boellaard, R., Kropholler, M. A., Lubberink, M., & Lammertsma, A. A. (2004). Simulated Annealing in pharmacokinetic modeling of PET neuroreceptor studies: Accuracy and precision compared with other optimization algorithms. 3222-3225. Paper presented at 2004 Nuclear Science Symposium, Medical Imaging Conference, Symposium on Nuclear Power Systems and the 14th International Workshop on Room Temperature Semiconductor X- and Gamma- Ray Detectors, Rome, Italy.
Yaqub, Maqsood ; Boellaard, Ronald ; Kropholler, Marc A. ; Lubberink, Mark ; Lammertsma, Adriaan A. / Simulated Annealing in pharmacokinetic modeling of PET neuroreceptor studies : Accuracy and precision compared with other optimization algorithms. Paper presented at 2004 Nuclear Science Symposium, Medical Imaging Conference, Symposium on Nuclear Power Systems and the 14th International Workshop on Room Temperature Semiconductor X- and Gamma- Ray Detectors, Rome, Italy.4 p.
@conference{649c6d91764c48fb92f4bd318698d185,
title = "Simulated Annealing in pharmacokinetic modeling of PET neuroreceptor studies: Accuracy and precision compared with other optimization algorithms",
abstract = "Fitting of PET pharmacokinetic parameters may suffer from bias or unrealistic outcomes, especially for noisy data. There are many readily available optimization algorithms, but each has different characteristics when fitting PET pharmacokinetic models. The purpose of this study is to evaluate the performance of four different types of optimization algorithms, including a modified simulated annealing method, in terms of precision and accuracy of PET pharmacokinetic parameters. The simulated annealing algorithm (SA), called Basin-hoping, was modified for present application. Input data, taken from [ 11C]-PK11195 neuroreceptor ligand studies, was used to simulate time activity curves at various noise levels. Also the influence of incorrect weighting factors on algorithm performance was studied. Surprisingly, effects of using incorrect but reasonable weighting factors on bias and precision were negligible. Except when extreme and unrealistic weighting factors were used, an increase in bias and decrease in precision was observed. In general the modified SA provided smallest weighted squared residual error and was able to find the global minimum without the need for a proper start parameter selection. However, occasionally better fits were found with the interior-reflective Newton method, but only when implemented using a range of start parameters, centered on the expected value.",
author = "Maqsood Yaqub and Ronald Boellaard and Kropholler, {Marc A.} and Mark Lubberink and Lammertsma, {Adriaan A.}",
year = "2004",
month = "12",
day = "1",
language = "English",
pages = "3222--3225",
note = "2004 Nuclear Science Symposium, Medical Imaging Conference, Symposium on Nuclear Power Systems and the 14th International Workshop on Room Temperature Semiconductor X- and Gamma- Ray Detectors ; Conference date: 16-10-2004 Through 22-10-2004",

}

Yaqub, M, Boellaard, R, Kropholler, MA, Lubberink, M & Lammertsma, AA 2004, 'Simulated Annealing in pharmacokinetic modeling of PET neuroreceptor studies: Accuracy and precision compared with other optimization algorithms' Paper presented at 2004 Nuclear Science Symposium, Medical Imaging Conference, Symposium on Nuclear Power Systems and the 14th International Workshop on Room Temperature Semiconductor X- and Gamma- Ray Detectors, Rome, Italy, 16/10/2004 - 22/10/2004, pp. 3222-3225.

Simulated Annealing in pharmacokinetic modeling of PET neuroreceptor studies : Accuracy and precision compared with other optimization algorithms. / Yaqub, Maqsood; Boellaard, Ronald; Kropholler, Marc A.; Lubberink, Mark; Lammertsma, Adriaan A.

2004. 3222-3225 Paper presented at 2004 Nuclear Science Symposium, Medical Imaging Conference, Symposium on Nuclear Power Systems and the 14th International Workshop on Room Temperature Semiconductor X- and Gamma- Ray Detectors, Rome, Italy.

Research output: Contribution to conferencePaperAcademic

TY - CONF

T1 - Simulated Annealing in pharmacokinetic modeling of PET neuroreceptor studies

T2 - Accuracy and precision compared with other optimization algorithms

AU - Yaqub, Maqsood

AU - Boellaard, Ronald

AU - Kropholler, Marc A.

AU - Lubberink, Mark

AU - Lammertsma, Adriaan A.

PY - 2004/12/1

Y1 - 2004/12/1

N2 - Fitting of PET pharmacokinetic parameters may suffer from bias or unrealistic outcomes, especially for noisy data. There are many readily available optimization algorithms, but each has different characteristics when fitting PET pharmacokinetic models. The purpose of this study is to evaluate the performance of four different types of optimization algorithms, including a modified simulated annealing method, in terms of precision and accuracy of PET pharmacokinetic parameters. The simulated annealing algorithm (SA), called Basin-hoping, was modified for present application. Input data, taken from [ 11C]-PK11195 neuroreceptor ligand studies, was used to simulate time activity curves at various noise levels. Also the influence of incorrect weighting factors on algorithm performance was studied. Surprisingly, effects of using incorrect but reasonable weighting factors on bias and precision were negligible. Except when extreme and unrealistic weighting factors were used, an increase in bias and decrease in precision was observed. In general the modified SA provided smallest weighted squared residual error and was able to find the global minimum without the need for a proper start parameter selection. However, occasionally better fits were found with the interior-reflective Newton method, but only when implemented using a range of start parameters, centered on the expected value.

AB - Fitting of PET pharmacokinetic parameters may suffer from bias or unrealistic outcomes, especially for noisy data. There are many readily available optimization algorithms, but each has different characteristics when fitting PET pharmacokinetic models. The purpose of this study is to evaluate the performance of four different types of optimization algorithms, including a modified simulated annealing method, in terms of precision and accuracy of PET pharmacokinetic parameters. The simulated annealing algorithm (SA), called Basin-hoping, was modified for present application. Input data, taken from [ 11C]-PK11195 neuroreceptor ligand studies, was used to simulate time activity curves at various noise levels. Also the influence of incorrect weighting factors on algorithm performance was studied. Surprisingly, effects of using incorrect but reasonable weighting factors on bias and precision were negligible. Except when extreme and unrealistic weighting factors were used, an increase in bias and decrease in precision was observed. In general the modified SA provided smallest weighted squared residual error and was able to find the global minimum without the need for a proper start parameter selection. However, occasionally better fits were found with the interior-reflective Newton method, but only when implemented using a range of start parameters, centered on the expected value.

UR - http://www.scopus.com/inward/record.url?scp=23844488222&partnerID=8YFLogxK

M3 - Paper

SP - 3222

EP - 3225

ER -

Yaqub M, Boellaard R, Kropholler MA, Lubberink M, Lammertsma AA. Simulated Annealing in pharmacokinetic modeling of PET neuroreceptor studies: Accuracy and precision compared with other optimization algorithms. 2004. Paper presented at 2004 Nuclear Science Symposium, Medical Imaging Conference, Symposium on Nuclear Power Systems and the 14th International Workshop on Room Temperature Semiconductor X- and Gamma- Ray Detectors, Rome, Italy.