Treatment of radioactive decay in pharmacokinetic modeling:: Influence on parameter estimation in cardiac 13N-PET

Radioactive decay during measurement can be accounted for by either a (d) under bar ecay correction of the measured data (b) under bar efore (m) under bar odeling (DbM) or by direct implementation of (d) under bar ecay (i) under bar nto the pharmacokinetic (m) under bar odel (DiM). The purpose of this study was to quantify the influence of the type of decay correction on the calculated parameters for the example of a three-compartment model used for the calculation of myocardial perfusion with N-13 ammonia and positron emission tomography (PET). For a given input function [C-a(t)proportional to t exp(-kt), k = 1.72/min] the tissue uptake for two parameter sets of K1, k2, k3, TBV were calculated for 20 frames (12 X 10 s, 3 X 30 s, 3 X 120 s, 1 X 300 s). These values were mathematically deteriorated by various noise levels according to Poisson statistics and fitted by a Levenberg-Marquardt algorithm. Estimated parameter means and coefficients of variation of the fitted parameters were calculated for the DbM and DiM case. The estimated parameter means for both decay correction methods were of comparable quality. The important measure for a single fit is the relative variability of the fitted parameters. This value is up to a factor 1.15 smaller for K-1 obtained with DiM and a reasonable noise level of 10%. Therefore, decay correction should be taken into account during modeling to reduce the variability in the fitted parameters. (C) 1999 American Association of Physicists in Medicine. [S0094-2405(99)00204-7].