PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Nucl Med. Author manuscript; available in PMC Mar 1, 2012.
Published in final edited form as:
PMCID: PMC3290866
NIHMSID: NIHMS358702
An Approach for Balancing Diagnostic Image Quality with Cancer Risk: Application to Pediatric Diagnostic Imaging of 99mTc-Dimercaptosuccinic Acid
George Sgouros,1 Eric C. Frey,1 Wesley E. Bolch,2 Michael B. Wayson,2 Andres F. Abadia,2 and S. Ted Treves3
1Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University, School of Medicine, Baltimore, Maryland
2Nuclear Engineering Department, University of Florida, Gainesville, Florida
3Boston Children’s Hospital, Harvard University, Boston, Massachusetts
For correspondence or reprints contact: George Sgouros, Johns Hopkins University, School of Medicine CRB II, 4M.61 1550 Orleans St., Baltimore, MD 21231. gsgouros/at/jhmi.edu
A recent survey of pediatric hospitals showed a large variability in the activity administered for diagnostic nuclear medicine imaging of children. Imaging guidelines, especially for pediatric patients, must balance the risks associated with radiation exposure with the need to obtain the high-quality images necessary to derive the benefits of an accurate clinical diagnosis.
Methods
Pharmacokinetic modeling and a pediatric series of nonuniform rational B-spline–based phantoms have been used to simulate 99mTc-dimercaptosuccinic acid SPECT images. Images were generated for several different administered activities and for several lesions with different target-to-background activity concentration ratios; the phantoms were also used to calculate organ S values for 99mTc. Channelized Hotelling observer methodology was used in a receiver-operating-characteristic analysis of the diagnostic quality of images with different modeled administered activities (i.e., count densities) for anthropomorphic reference phantoms representing two 10-y-old girls with equal weights but different body morphometry. S value–based dosimetry was used to calculate the mean organ-absorbed doses to the 2 pediatric patients. Using BEIR VII age- and sex-specific risk factors, we converted absorbed doses to excess risk of cancer incidence and used them to directly assess the risk of the procedure.
Results
Combined, these data provided information about the tradeoff between cancer risk and diagnostic image quality for 2 phantoms having the same weight but different body morphometry. The tradeoff was different for the 2 phantoms, illustrating that weight alone may not be sufficient for optimally scaling administered activity in pediatric patients.
Conclusion
The study illustrates implementation of a rigorous approach for balancing the benefits of adequate image quality against the radiation risks and also demonstrates that weight-based adjustment to the administered activity is suboptimal. Extension of this methodology to other radiopharmaceuticals would yield the data required to generate objective and well-founded administered activity guidelines for pediatric and other patients.
Keywords: pediatric imaging, dose reduction, image simulation, dosimetry, DMSA, optimization
The goal of every nuclear medicine study is to obtain the maximum diagnostic information as cost-effectively and in as short a period of time as possible and with the lowest patient radiation exposure (1). Children are generally more sensitive to radiation and have their entire lifetimes to manifest radiation effects, so these considerations are even more pressing in pediatric populations. Consensus pediatric dosing guidelines for diagnostic studies in nuclear medicine have recently been published in the United States (2). These update previously described guidelines (1,39). Administered activities (AAs) in pediatric nuclear medicine have been developed empirically, taking into account the whole-body radiation absorbed dose, type of examination, available photon flux, instrumentation, and examination time: high AAs, which do not directly result in improved diagnostic sensitivity or accuracy, or low AAs, which do not permit adequate examination, result in unnecessary radiation exposures. AA estimates for pediatric patients based on adult AAs corrected for body weight or body surface area are generally used for children over 1 y of age. Premature infants and newborns require special consideration, and the concept of minimum total AA should be considered. Minimum total AA can be defined as that AA of the radiopharmaceutical below which the study will be inadequate regardless of the patient’s body weight or surface area.
In this work, we used modeling-derived biodistribution data to simulate images generated using different AAs or imaging durations (i.e., different image count densities) for 2 simulated 10-y-old girls having the same weight but different body morphometry. The resulting images were evaluated for task-based diagnostic quality and, combined with the corresponding calculated absorbed and effective doses, were used to define the tradeoff between diagnostic image quality and risk to the patient of radiogenic cancer.
Pharmacokinetic Modeling and Dosimetry
The pharmacokinetics of dimercaptosuccinic acid (DMSA) in children of various ages have already been characterized (10,11). The whole-body clearance rate, λeff, was not explicitly listed in these references; it was calculated from the urine data provided by fitting a monoexponential function to the complement of the urine data (100% excreted) at each time point. These parameters (Table 1) were used with the appropriate exponential expressions to generate the kinetic data used in the imaging simulations and also in the dosimetry calculations of the present study (Fig. 1).
TABLE 1
TABLE 1
99mTc-DMSA Pharmacokinetic Parameters for Pediatric Patients
FIGURE 1
FIGURE 1
99mTc-DMSA time–activity curves for 10-y-old. Vertical dotted line corresponds to time after injection at which DMSA imaging is performed. WB = whole body.
The analytic expressions were integrated from zero to infinity to give the number of nuclear transformations (NTs) in individual organs for the absorbed dose calculations. NTs were also provided for lungs and heart with contents. These were obtained by apportioning a calculated blood NT concentration to each organ according to its blood volume. Blood NT was obtained from the work of Treves (1) wherein clearance of DMSA from blood is described as following a single exponential function with a 56-min half-life. Absorbed doses were calculated according to the MIRD Committee methodology as described in MIRD pamphlet 21 (12). The working expression for absorbed dose is:
equation M1
Eq. 1
The radionuclide and anatomic model-dependent portion (i.e., the term to the right) of Equation 1 is designated the S value and is indexed by source (rS) and target region (rT) (12,13). The S value is given as:
equation M2
Eq. 2
wherein a sum over all radionuclide emissions is included. The numerator in this expression is the product of the energy emitted per disintegration of the radionuclide for emission i, Δi, and the fraction, ϕ(rTrS, Ei), of this emitted energy that is absorbed by the target tissue rT, when emissions with energy Ei originate in source tissue rS. The denominator is the mass of the target region, M(rT).
The time-integrated activity Ã(rS) (or number of NTs) in source tissue rS is given by:
equation M3
Eq. 3
Organ-absorbed doses obtained in this way were used to calculate both the tissue-factor-weighted organ-equivalent dose, H(rT), and the effective dose, E:
equation M4
Eq. 4
where the first weighting parameter, wR, is related to the biologic effectiveness of different types of radiation particles for stochastic effects and is equal to 1 for 99mTc emissions. The second weighting parameter—the tissue-weighting factor wT—is defined by International Commission on Radiological Protection Committee (ICRP) 1 and represents the proportion of total radiation detriment attributed to individual tissue rT in the body. Radiation detriment is defined in ICRP publication 60 (14) and more recently publication 103 (15) and includes a variety of stochastic endpoints such as cancer mortality, cancer morbidity, and years of life lost.
Cancer Risk Estimation
To balance the risk of an imaging procedure with the image quality needed for accurate diagnosis, we examined several alternatives to evaluating risk from the absorbed and effective dose calculations. As noted in ICRP 103, effective dose is appropriate in assessing radiation detriment, which includes cancer incidence risk for workers and the general population (15). Risk assessment for medical diagnosis is “best evaluated using appropriate risk values for the individual tissues at risk and for the age and sex distribution of the individuals undergoing the medical procedures.” Consistent with this recommendation and to further optimize AA recommendations for pediatric diagnostic imaging, we used sex- and age-specific cancer incidence risk estimates, expressed as lifetime attributable risk (LAR), as tabulated in the National Academy of Science’s BEIR VII report (16). Although these risk estimates have high uncertainties for individual patients, they can be used collectively in optimizing values of AA and resulting organ-absorbed dose to classes of pediatric nuclear medicine patients. In the case illustrated in this work, the product of organ-absorbed doses and corresponding LAR for a 10-y-old girl was summed over all organs to give the overall risk for cancer incidence. This approach was compared with a risk estimate based on sex- and age-averaged LAR values and with an overall radiation detriment derived from the effective dose using ICRP 103 tissue-weighting factors (15). The latter comparison was made using a detriment-adjusted nominal risk coefficient for cancer incidence for the whole population that is reported in ICRP publication 103 as 5.5 × 10−2 Sv−1. This value was multiplied by the effective dose to give an effective dose–derived overall risk for radiation detriment, which includes cancer incidence. The organ LAR and ICRP 103 weighting factors used in the calculations are shown in Table 2.
TABLE 2
TABLE 2
Organ wT and Lifetime Attributable Risk (LAR) Used in Risk Calculations
Pediatric Phantom
The implications of administering the same weight-based activity to a tall, thin child and a short, heavy child both of the same weight were examined using 2 versions of the 10-y-old female phantom. The Advanced Laboratory for Radiation Dosimetry Studies (ALRADS) UF NHANES–based phantom series (17,18) was used to select a 90th weight-percentile (fifth height bin) child and a 10th weight-percentile (seventh height bin) child. The 2 phantoms are shown in Figure 2. To calculate S values, the phantoms were voxelized at 1.6- and 1.5-mm isotropic resolutions, respectively; 1.2-mm voxelization was used for the image simulations. An in-house MATLAB (The MathWorks, Inc.) code, Voxelizer 4 described previously by Lee et al. (19), was used to perform the voxelization. The resulting hybrid voxel phantom was imported into the radiation transport code MCNPX 2.5 (20). Photon and electron emission spectra from 99mTc were then abstracted from a recently updated decay scheme monograph published by the MIRD committee (21). Two Monte Carlo simulations were thus run—one for photon emissions (x and γ rays) and another for electron emissions (β, IC, and Auger electrons)—each with 5 × 106 particle histories yielding this number of histories per source–target region pair. Photon and electron components of the S value were then combined to yield total spectral S values as listed in Table 3. Source organs correspond to those in the biokinetic model of 99mTc-labeled DMSA. Target tissues were those assigned values of wT in the ICRP 103 definition of effective dose, which includes the collective target region remainder tissues. This approach—direct sampling of radionuclide emission spectrum in MCNPX—yields higher S value convergence with fewer Monte Carlo histories than does the traditional approach, in which monoenergetic values of specific absorbed fractions are assembled and then spectrum-weighted.
FIGURE 2
FIGURE 2
Ten-year-old female phantoms weighing 32 kg: anatomic depiction (A) and coronal slice through corresponding attenuation map (B), generated from anatomic description of A. %ile = percentile; H = height.
TABLE 3
TABLE 3
wT-Weighted Equivalent Organ Doses, Absorbed Doses, and LAR Values for Two 32-kg Pediatric Patients for 99mTc-DMSA Renal Imaging Studies
Imaging and Observer Simulations
The 1.2-mm voxel representation of the two 10-y female UF phantoms was used as input to an analytic projection code that models attenuation, the spatially varying collimator–detector response (22), and object-dependent scatter (23). Organ uptake was based on values from the pharmacokinetic modeling, with imaging at 3 h after injection. From the models, we estimated the percentage injected activities in the left and right kidneys, spleen, liver, and body. In each kidney we distributed the activity to the cortex, medulla, and pelvis on the basis of observed image intensity ratios in pediatric SPECT images. The total activity in the phantom was based on the 32-kg weight of the phantom times a dose of 1.85 MBq per kg, the standard dosing at Boston Children’s Hospital. Projection data were simulated using an analytic projection code that models attenuation, detector, and scatter. This code has been extensively validated for imaging of a variety of isotopes by comparison to Monte Carlo and experimental projection data. We modeled a low-energy, high-resolution collimator, 128 × 128 projection images, a 360° body-contouring orbit, 120 projection views, and a 0.22-cm projection bin size. Organ uptake was varied randomly with a uniform distribution over a range of ± 10% of the modeled values. We created defects with uptake reduced by uniformly distributed random amounts in the range of 15%–25% relative to the activity concentration of surrounding normal kidney, giving defect–to–background activity ratios of 0.15–0.25. In each kidney, the defects were placed in the inferior, lateral, and superior positions, with a size resulting from the intersection of a 2.2-cm-radius sphere centered on the defect location within the cortex of the kidney. We created 50 random uptake variations for each defect location and images with and without the defect, for a total of 600 phantom projections. The projections were scaled to represent the count level that would be obtained using a low-energy, high-resolution collimator (sensitivity of 85.6 counts/s/MBq) and a 30-min acquisition using a dual-camera system. The images were then scaled to represent 25%, 50%, 75%, 100%, 125%, and 150% of this count level (corresponding to a change in the AA by the same proportion), and Poisson noise was additionally simulated. There was an average of approximately 31 kilo counts in the projection data for one 0.22-cm-thick slice in the region over the kidneys for the 100% count level. The resulting projections were reconstructed using filtered backprojection followed by an order 5 3D Butterworth filter with a cutoff frequency of 0.73 cm−1. This filter was chosen because it subjectively gave the best image quality, but the filter was not otherwise optimized. Sample coronal images for the different body morphometry and count levels are shown in Figure 3.
FIGURE 3
FIGURE 3
Example of set of simulated coronal 99mTc-DMSA SPECT images through kidneys used in CHO study. Images are from 2 phantoms with body weight of 32 kg and with heights of 125 and 147 cm. In this example, images are shown with and without single defect (defect–to–kidney (more ...)
We then extracted 64 × 64 coronal images centered on the defect centroids and used these to perform a channelized Hotelling observer (CHO) study. This CHO methodology has been described previously (24). It is a computer simulation–based procedure for diagnostic task–based evaluation of images to generate receiver-operating-characteristic (ROC) curves. It has been shown to predict human observer performance in the detection of a myocardial perfusion defect (25), including an investigation of the effects of reduced count levels (26). We used difference-of-mesa channels with octave widths in the frequency domain and a starting frequency of 0.5 cycles per pixel, giving a total of 6 frequency domain channels. The CHO for each phantom and noise level was trained and tested using a leave-one-out strategy in which 599 of the images were used to train the observer; the resulting CHO was then tested on the remaining phantom, yielding a single test statistic value. We applied this strategy using each of the 600 images, thus producing 600 CHO test statistic values for each phantom and noise-level combination. The test statistics were then analyzed using LAB-ROC4 (27) and CLABROC (28). The LABROC4 code performs maximum-likelihood estimation of the ROC curve parameters and computes the area under the ROC curve (AUC) and its SD from continuously distributed rating data assuming a binormal model. The method used does not require that rating data themselves be binormally distributed but rather that there exists a monotonic transformation (which may be unknown) of the rating data that produces binormally distributed data. The program CLABROC performs hypothesis testing of the AUC values assuming paired rating data. It takes advantage of the correlations in the data to improve the statistical power, and thus the p-values are not estimated directly from the SD of the AUC produced by LABROC4.
The results of the CHO study and area under the ROC curve (AUC) values for both phantoms at the various count levels are shown in Figure 4. The error bars for the AUC values were obtained from LABROC4. We tested the significance of differences in the AUC for the 2 phantoms at the same noise levels, taking into account the paired nature of the data (using CLABROC). The 1-tailed p-levels for the percentage AAs were 0.017, 0.046, 0.28, 0.008, 0.024, and 0.017 (ordered by increasing AA). Thus, the only AA not significant at a p-level of 0.05 was for the 75% AA. For the remaining cases, we can reject the null hypothesis that the AUC for the tall, thin patient is less than or equal to the AUC for the shorter–stouter patient at the p-level of 0.05.
FIGURE 4
FIGURE 4
Plot of AUC as function of relative AA for 2 different simulated phantoms. AUC provides measure of detectability of regions with reduced kidney function, with higher AUC indicating better performance. AUC curve for taller (147 cm) patient is approximately (more ...)
In interpreting the importance of the difference in AUC values, one should note that the curves would be roughly equivalent if the curve for the 147-cm phantom were shifted to the left by 50% of the standard AA. In other words, for the 147-cm-tall phantom, the same performance in terms of identifying the defect could have been obtained with half of the standard weight-based injected activity normally used for the 125-cm-tall patient. This result indicates that body morphometry, and not just weight, is a significant factor in determining image quality for a given AA. The difference in image quality was due to a combination of 3 factors: the 125-cm patient had about 10% fewer photons escaping the body, required a 10% larger average camera orbit radius resulting in poorer resolution, and a 10% higher scatter-to-primary ratio, resulting in slightly higher noise, poorer resolution, and poorer contrast for the shorter patient.
To evaluate the dose implications, 99mTc S values were generated for the same 2 phantoms. Model-derived residence times for each source organ of a 10-y-old were used to calculate organ-absorbed doses for selected target tissues listed in ICRP 103. The tissue-weighting factors were applied to the absorbed doses, and tissue-factor-weighted equivalent doses were obtained for each phantom. The resulting effective doses per unit AA were 1.42 × 10−2 mSv/MBq and 9.84 × 10−3 mSv/MBq for the short, stout and tall, thin phantoms, respectively. Organ LAR values were also applied to obtain the overall risk for cancer incidence in these 2 hypothetical pediatric patients. The consensus weight-based AA of 1.85 MBq/kg for 99mTc-DMSA (2) was applied to obtain the organ and overall risk values for the short, stout patient; half of this AA was used for the tall, thin patient to take into account the fact that the same image quality could be obtained with 50% of the AA. Individual organ results are summarized in Table 3, and overall cancer incidence risks are compared in Table 4. For this particular agent, and its corresponding organ dose distribution, a difference between age- and sex-averaged risk and patient-specific risk can be noted. The lower AA needed for the tall, thin child translates into overall cancer incidence risk values that are less than half of those for the short, stout child.
TABLE 4
TABLE 4
Summary of Risk Calculations
Consensus guidelines for pediatric AA in nuclear medicine imaging have been recently published (2). The guidelines describe the consensus of expert practitioners in the field. The challenge in identifying the appropriate AA in pediatric nuclear medicine imaging is to balance the long-term concern of radiation risk with the immediate need to obtain images that are diagnostically adequate. The Dosimetry and Pediatrics Committees of the European Association of Nuclear Medicine (8,9,29) developed a method that addresses this challenge. In that approach, each radiopharmaceutical considered was assumed to be uniformly distributed in one of the age-specific, Oak Ridge National Laboratory–stylized phantoms (30), and the complement of the whole-body absorbed fraction (1 minus the whole-body absorbed fraction) was taken as proportional to the photon flux that would be observed by an imaging device. The ratio of this value to the equivalent value obtained for the adult male phantom was used as the normalization factor for photon flux.
The methodology described in the present study is similar in principle to that used by the European Association of Nuclear Medicine Dosimetry Committee but with the following differences. First, in place of calculated photon flux as the measure of image quality, simulated patient images were generated and evaluated by computer observers to obtain a task-based measure of image quality. Second, instead of a uniform whole-body distribution, pharmacokinetic model–based realistic activity distributions were used in the simulated images. As shown in this work, this approach makes it possible to consider anatomic variations for a specific age and weight. Such variations can substantially influence the activity required to obtain good-quality diagnostic images. The approach could ultimately be used to obtain a dataset to help establish pediatric dosing guidelines across a wide spectrum of ages, weights, and body types, bringing us closer to achieving the goal of deriving the maximum diagnostic information in the shortest time and with the lowest patient radiation exposure.
For this rather extreme example (only 10% of the pediatric population at that weight falls into these 2 categories), the effective dose per AA to the short, stout pediatric patient is 44% greater than that to the tall, thin patient. This is reflected in the lower values of the entries in the D(rT)/A0 column (Table 3) for the tall, thin patient than those corresponding to the short, stout patient. The distances between source and target organs in the short phantom are less than those in the tall phantom. There are also organ mass differences that will affect the absorbed dose. Larger target organs will absorb more energy, but the energy will be divided by a larger mass; consequently, it is difficult to predict in advance the impact on each organ.
To better understand the cancer incidence risks for pediatric patients associated with imaging at different AAs, organ-absorbed doses were used to calculate LAR values. Because effective dose calculations use sex- and age-averaged tissue-weighting factors, we examined the effect of using averaged values versus sex- and age-specific values. In the case of 99mTc-DMSA, the dose distribution gave overall risk values that were greater than those obtained by averaging (Table 4). As shown in Table 2, however, the effect of age and sex averaging is strongly dependent on the organ dose profile, and, in particular, for agents that localize to marrow the LAR can differ by almost 2-fold. In most cases, age- and sex-specific risk estimation would allow less AA for equal cancer risk than age- and sex-averaged risk estimation, suggesting that both sex and age, in addition to body morphometry, should be considered when defining optimal pediatric dosing.
We present a methodology that balances image quality with cancer risk in evaluating the appropriate AA level for pediatric patients. Using this method, we show that weight-based adjustment alone can lead to unnecessarily high AAs in a certain portion of the pediatric patient population. We also demonstrate age- and sex-specific evaluations of overall cancer incidence risk and examine the implications of using sex- and age-averaged risk parameters. The techniques developed in this work can lead to pediatric imaging AA guidelines that rigorously and transparently balance image quality with associated excess radiation-induced cancer risks. Extension of this methodology to other radiopharmaceuticals would yield the data required to generate objective and well-founded AA guidelines for pediatric and other patients.
1. Treves ST. Pediatric Nuclear Medicine/PET. 3rd ed. New York, NY: Springer Science+Business Media LLC; 2007.
2. Gelfand MJ, Parisi MT, Treves ST. Pediatric radiopharmaceutical administered doses: 2010 North American consensus guidelines. J Nucl Med. 2011;52:318–322. [PubMed]
3. Fahey FH, Palmer MR, Strauss KJ, Zimmerman RE, Badawi RD, Treves ST. Dosimetry and adequacy of CT-based attenuation correction for pediatric PET: phantom study. Radiology. 2007;243:96–104. [PubMed]
4. Gelfand MJ, Lemen LC. PET/CT and SPECT/CT dosimetry in children: the challenge to the pediatric imager. Semin Nucl Med. 2007;37:391–398. [PubMed]
5. Fahey FH. Dosimetry of pediatric PET/CT. J Nucl Med. 2009;50:1483–1491. [PubMed]
6. Alessio AM, Kinahan PE, Manchanda V, Ghioni V, Aldape L, Parisi MT. Weight-based, low-dose pediatric whole-body PET/CT protocols. J Nucl Med. 2009;50:1570–1577. [PubMed]
7. Accorsi R, Karp JS, Surti S. Improved dose regimen in pediatric PET. J Nucl Med. 2010;51:293–300. [PubMed]
8. Lassmann M, Biassoni L, Monsieurs M, Franzius C. The new EANM paediatric dosage card: additional notes with respect to F-18. Eur J Nucl Med Mol Imaging. 2008;35:1666–1668. [PubMed]
9. Lassmann M, Biassoni L, Monsieurs M, Franzius C, Jacobs F. The new EANM paediatric dosage card. Eur J Nucl Med Mol Imaging. 2007;34:796–798. [PubMed]
10. Smith T, Evans K, Lythgoe MF, Anderson PJ, Gordon I. Radiation dosimetry of technetium-99m-DMSA in children. J Nucl Med. 1996;37:1336–1342. [PubMed]
11. Evans K, Lythgoe MF, Anderson PJ, Smith T, Gordon I. Biokinetic behavior of technetium-99m-DMSA in children. J Nucl Med. 1996;37:1331–1335. [PubMed]
12. Bolch WE, Eckerman KF, Sgouros G, Thomas SR. MIRD pamphlet No. 21: a generalized schema for radiopharmaceutical dosimetry—standardization of nomenclature. J Nucl Med. 2009;50:477–484. [PubMed]
13. Loevinger R, Budinger TF, Watson EE. MIRD Primer for Absorbed Dose Calculations. Revised ed. New York, NY: The Society of Nuclear Medicine, Inc.; 1991.
14. International Commission on Radiological Protection (ICRP) Pergamon, Oxford: ICRP; 1991. 1990 Recommendations of the International Commission on Radiological Protection. Publication 60.
15. International Commission on Radiological Protection (ICRP) ICRP publication 103: recommendations of the International Commission on Radiological Protection. Ann ICRP. 2007;37:1–332. [PubMed]
16. National Research Council. Washington, DC: National Research Council; 2005. Health Risks from Exposure to Low Levels of Ionizing Radiation: BEIR VII—Phase 2.
17. Lee C, Lodwick D, Hurtado J, Pafundi D, Williams JL, Bolch WE. The UF family of reference hybrid phantoms for computational radiation dosimetry. Phys Med Biol. 2010;55:339–363. [PMC free article] [PubMed]
18. Johnson P, Whalen S, Wayson M, Juneja B, Lee C, Bolch W. Hybrid patient-dependent phantoms covering statistical distributions of body morphometry in the US adult and pediatric population. Proc IEEE. 2009;97:2060–2075.
19. Lee C, Lodwick D, Hasenauer D, Williams JL, Lee C, Bolch WE. Hybrid computational phantoms of the male and female newborn patient: NURBS-based whole-body models. Phys Med Biol. 2007;52:3309–3333. [PubMed]
20. Pelowitz DB. MCNPX User’s Manual Version 2.5.0. Los Alamos, NM: Los Alamos National Laboratory; 2005. LA-CP-05-0369.
21. Eckerman KF, Endo A. MIRD: Radionuclide Data and Decay Schemes. 1st ed. Reston, VA: The Society of Nuclear Medicine; 2008.
22. Frey EC, Ju ZW, Tsui BMW. A fast projector-backprojector pair modeling the asymmetric, spatially varying scatter response function for scatter compensation in SPECT imaging. IEEE Trans Nucl Sci. 1993;40:1192–1197.
23. Frey EC, Tsui BMW. A new method for modeling the spatially-variant, object shape dependent scatter response function in SPECT; IEEE Nucl Sci Symp Conf Rec; 1996. pp. 1082–1086.
24. Frey EC, Gilland KL, Tsui BMW. Application of task-based measures of image quality to optimization and evaluation of three-dimensional reconstruction-based compensation methods in myocardial perfusion SPECT. IEEE Trans Med Imaging. 2002;21:1040–1050. [PubMed]
25. Sankaran S, Frey EC, Gilland KL, Tsui BMW. Optimum compensation method and filter cutoff frequency in myocardial SPECT: a human observer study. J Nucl Med. 2002;43:432–438. [PubMed]
26. He X, Links JM, Frey EC. An investigation of the trade-off between the count level and image quality in myocardial perfusion SPECT using simulated images: the effects of statistical noise and object variability on defect detectability. Phys Med Biol. 2010;55:4949–4961. [PMC free article] [PubMed]
27. Metz CE, Herman BA, Shen JH. Maximum likelihood estimation of receiver operating characteristic (ROC) curves from continuously-distributed data. Stat Med. 1998;17:1033–1053. [PubMed]
28. Metz CE, Wang P, Kronman HB. Information Processing in Medical Imaging. The Hague, The Netherlands: Martinus Nijhoff; 1984. A new approach for testing the significance of differences between ROC curves measured from correlated data.
29. Jacobs F, Thierens H, Piepsz A, et al. Optimised tracer-dependent dosage cards to obtain weight-independent effective doses. Eur J Nucl Med Mol Imaging. 2005;32:581–588. [PubMed]
30. Cristy M, Eckerman KF. Specific Absorbed Fractions of Energy at Various Ages for Internal Photon Sources. Oak Ridge, TN: Oak Ridge National Laboratory; 1987. ORNL/TM-8381.