|
|
||||||||
Basic Science Investigations |
1 Department of Medical Physics, Memorial Sloan-Kettering Cancer Center, New York, New York
2 CPS Innovations, Knoxville, Tennessee
| ABSTRACT |
|---|
|
|
|---|
Key Words: National Electrical Manufacturers Association acceptance testing PET/CT
| INTRODUCTION |
|---|
|
|
|---|
The National Electrical Manufacturers Association (NEMA) has updated their standard PET performance measurement recommendations. The updated document (NU 2-2001) includes revised measurements for spatial resolution, scatter fraction, sensitivity, counting rate performance, and the accuracy of count loss and randoms corrections (3). This NU 2-2001 document also includes a specification for a new image quality measurement that was not considered in the original NU 2-1994 standard (4). Because of the prevailing interest in clinical whole-body scanning, this measurement is intended to simulate realistic clinical conditions for such scans. A recent publication by Daube-Witherspoon et al. (5) discussed in detail the performance measurements of several PET cameras using the new NU 2-2001 standards.
A brief introduction of each NU 2-2001 test will be followed by results of those tests for an LSO PET/CT scanner (CPS Innovations; marketed as the Reveal RT by CTI, Inc., and as the Biograph LSO by Siemens Medical Solutions). Other sources provide more details, and users wishing to perform the measurements described here are strongly encouraged to consult those publications for a complete description of procedures (3,5). However, it is not possible to follow all of these procedures exactly on LSO systems because of the presence of intrinsic background radiation in the detectors.
LSO contains natural lutetium, of which the radioactive nuclide 176Lu composes 2.6%. The background single-event counting rate resulting from this activity depends on the details of the detectors construction and energy discrimination settings. This rate is typically on the order of 1 million counts per second (cps) in the whole scanner, for a 350- to 650-keV acceptance window, which corresponds to a total random coincidence rate of about 1,600 cps. Singles rates in clinical acquisitions are generally 1020 times higher. Noting that randoms increase approximately as the square of the singles rate (neglecting counting losses), the main effect of the intrinsic events is to increase the rate of random coincidences in the scanner by a few percent in the clinical range (about 8% at peak noise equivalent counting rate [NECR]). There is also a small (approximately 600 cps) background of true coincidences resulting from the intrinsic radiation, which may be compared with typical clinical trues rates of several hundred thousand counts per second. An intrinsic true event may result when a 176Lu nuclide undergoes ß-decay in one detector, depositing energy locally in that detector, and also emits a prompt
-ray that is absorbed in a second detector.
The benefits of the short decay time and high light output of LSO strongly outweigh the deleterious effects of the intrinsic radiation on the counting rate performance of the CPS LSO PET/CT. Nevertheless, one must properly account for the intrinsic contribution when performing the NEMA tests. It is not possible, for example, to achieve a randoms-to-trues ratio of less than 1% in the counting rate test as specified by the NEMA standard. Consequently, one cannot accurately determine the scatter fraction according to its prescription. In another example, the small intrinsic trues contribution must be accounted for when the sensitivity test is performed with a weak source. For these reasons, we have not followed exactly or in all cases the NEMA NU 2-2001 specifications for experimental technique and data processing. Instead, we have implemented the suggestions for adapting the NEMA protocols to LSO systems discussed in some detail by Watson et al. (6). These refined procedures preserve the intent of the NEMA standards but properly account for aspects of the special nature of LSO cameras not addressed in NU 2-2001.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Test Phantom Set and Activity
The NEMA NU 2-2001 tests require 3 sets of phantoms. First is an International Electrotechnical Commission (IEC) body phantom set, which consists of a torso cavity, removable lung insert, and 6 spheres (1037 mm in inner diameter [ID]). All are fillable. Second is a scatter phantom set, including a solid circular cylinder composed of polyethylene with a specific gravity of 0.96, an outside diameter of 203 mm, an overall length of 700 mm, and a plastic tube that is 800 mm long and has an ID of 3.2 mm to hold the source activity. A sensitivity phantom set is the third required and consists of 5 concentric aluminum tubes (each 70 mm long) and a 1.8-mL fillable polyethylene tube inserted into the center sleeve. These NEMA phantoms were manufactured by Data Spectrum Corp. and marketed by Biodex Medical Systems. Detailed information and specifications for the NEMA NU 2-2001 phantoms can be found at www.biodex.com. As specified by NU 2-2001, all measurements were performed with the 18F isotope.
Spatial Resolution
This measurement requires the construction of a point source. NU 2-2001 specifies that the activity be contained in a glass capillary tube with an ID less than 1 mm and that the tube also have an axial extent of less than 1 mm. It is further required that the source (i.e., tube) be placed parallel to the long axis of the tomograph for all measurements. We found this technique problematic for several reasons. The small volume of the source (<0.785 mm3) necessitates a very high specific activity to get a reasonable counting rate. It is difficult to accurately draw up this amount of fluid without wetting the walls of the capillary beyond the desired region. Moreover, positron annihilation is not confined to the location of the 18F in the fluid; positrons may stream up the capillary and annihilate in the glass of the tube. This positron range effect may tend to elongate the source in the direction of the tube. In our attempts at following the prescribed procedure, we measured significantly wider axial resolution when such a source was placed axially than when it was placed transaxially in the scanner, suggesting that the source was indeed asymmetric.
Consequently, we have followed a slightly different procedure for fabricating a point source. A small amount of ion exchange resin was placed in a small plastic tube (about 5 mm in length with an ID of 1 mm and outer diameter [OD] of 2 mm). Light packing material was placed in the tube on either end to hold the resin in place with an axial extent of less than 1 mm, and the tube was placed in a capsule that allowed it to be coupled to a syringe. These sources were manufactured by D and W, Inc. Approximately 300 MBq of 18F ion in 2 mL of water were then flushed 3 times through the resin, followed by a rinse with approximately 50 mL of nonradioactive water. This resulted in an activity of 1.7 MBq residual in the resin. The plastic tube containing the resin was removed from the capsule and placed on a piece of tape, which was then suspended from a needle in the FOV of the scanner (Fig. 1A). Two complete sets of resolution measurements were made with this source, one with the axis of the tube parallel to the axis of the scanner and one with it perpendicular. Although both sets of results are reported, we believe that the axial resolution measurements are more accurate for the transaxially oriented source and that the transverse resolution measurements are more accurate for the axially oriented source.
|
Intrinsic Scatter Fraction
Scatter fraction is the ratio of scattered events to the sum of true (unscattered) and scattered events, excluding randoms. The NU 2-2001 procedure specifies that the scatter fraction be measured as part of the counting rate performance test, which is used to acquire prompt coincidence data only (trues plus scatter and randoms), under the assumption that low-activity frames contain negligible randoms (3). However, for LSO-based detectors such a randoms-free regime does not exist. Consequently, if the NEMA prescription were followed in this case it would result in a significant overestimate of the scatter fraction and a very inaccurate estimate of the randoms. A more accurate approach would be to directly measure the random event rate together with the prompts. From these data the scattered event rate can be estimated very accurately as a function of activity by subtracting the measured randoms and trues event rates from the total event rate (6). The standard NEMA procedure neglects the counting rate dependence of the scatter fraction.
Therefore, for the scatter fraction, count losses, and randoms measurement, separate prompt and delayed sinograms were acquired using a standard delayed coincidence window technique. Activity was injected into the line source (800-mm length, 3.2-mm ID, 4.8-mm OD,
6-mL fillable volume), and the line source was threaded through the 6.4-mm hole in the test phantom (Fig. 1B). The test phantom was placed so that the line source was at the position nearest the patient table and was centered within 5 mm in the transverse and axial FOV.
Data analysis followed the modified procedure described elsewhere (6). Instead of estimating the randoms (Rr) from the prompts sinogram by correcting the totals (RTOT) for trues (Rt) and scatter (Rs), these were computed directly from the delayed window sinogram using the same bins used for the analysis of the prompts sinogram. The scatter event rate for each slice (i) at each activity (j) was then computed from:
![]() | (Eq. 1) |
![]() | (Eq. 2) |
Sensitivity
The sensitivity of a PET scanner is defined as the rate of true coincidence counts per second and per unit radioactivity that are detected for a line source located within the scanner FOV in the absence of attenuating media. To obtain an attenuation-free estimate, successive measurements are made with a uniform line source surrounded by a varying number of known absorbers. The absorbers consist of a set of nested, concentric aluminum tubes. Extrapolation of the response to 0 absorption gives an attenuation-free estimate of sensitivity. A 700-mm portion of the plastic tubing was filled with
2.0 MBq of 18F, a sufficiently low activity to ensure that counting losses were less than 1%. At this activity level, however, the true event rate was on the order of 104 cps, low enough that the intrinsic trues from the LSO itself (about 600 cps) were not negligible and required correction. For this purpose, a separate background scan was acquired with no phantom or activity present, and the total measured intrinsic trues rate was subtracted from the phantom data before computing the sensitivity (6).
After activation, the tube was suspended in the center of the transaxial FOV (Fig. 1C), parallel to the axis of the scanner, with all the absorbers in place. To help equalize counting rates as the source decayed, the measurement was started with all the aluminum tubes on the phantom. After each acquisition, a tube was pulled out and the measurement was repeated. More than 2 million counts were acquired with each configuration to ensure at least 10,000 trues per slice. The acquisitions were performed in the default energy window, which was 350650 keV. The same measurements were repeated at a 10-cm radial offset from the center of the transaxial view. Analysis of the data was performed according to NU 2-2001, using software supplied by the manufacturer.
Counting Rate Performance
To understand the scanners performance under a wide range of imaging conditions, it is necessary to measure its behavior as a function of counting rate. For this purpose the 70-cm polyethylene cylinder with a line source threaded through along its length was used (the same measurement as for the scatter fraction determination previously outlined). NU 2-2001 specifies that the initial amount of radioactivity shall be such as to cause the true event rate to reach 50% dead-time loss value. However, on LSO-based scanners, detector dead time is so low that a 50% loss is not observed until the communications channel of the coincidence processor saturates. With the present generation of electronics, this occurs at a total coincidence rate of about 2.6 million events per second, well beyond the clinical regime and peak NECR of the scanner. Nevertheless, this saturation regime was included in the test reported here.
Following the suggestion of the manufacturer, the line source was filled with 740 MBq of 18F initial activity, enough to achieve counting rates well beyond both the peak trues and the peak NECRs. Thirty-five frames were acquired, each with a 10-min duration, followed by a 15-min gap between the scans. Data were acquired over 15.5 h, and each acquisition had more than 3.3 million coincidence counts. In the analysis, the same process as described for intrinsic scatter fraction (delayed window sinogram) was used to account for randoms resulting from the intrinsic radioactivity of LSO.
Accuracy of Corrections for Count Losses and Randoms
PET scanners must be able to compensate for dead-time losses and randoms to achieve quantitative clinical measurements. Inaccuracy of corrections, especially at high counting rates, can cause bias in clinical scanning results. Thus it is important to measure the accuracy of the corrections for count losses and randoms. For this purpose, the counting rate data described previously were reconstructed with standard corrections and the FORE with ordered-subset expectation maximum (OSEM) iterative algorithm (2 iterations, 8 subsets). The counting rate error (
r) as a function of effective activity concentration was calculated as the deviation of the trues from a linear trend extrapolated from low-activity acquisitions where dead-time losses were negligible. Details of the analysis are available in the NEMA standard (3).
Image Quality Measurement
An image quality measurement has been incorporated into the NEMA testing scheme to better compare the quality of different imaging systems for a standardized situation that simulates clinical conditions. According to NU 2-2001, the large background region of the IEC body phantom is to be filled with an activity of 5.3 kBq/mL. Activated and nonactivated spheres of various sizes are placed in this warm background. The hot spheres (10, 13, 17, and 22 mm in diameter) are to be filled with an activity concentration of n times that of the background, where n = 4 and 8. The 2 largest spheres (28 and 37 mm) are filled with water only. A 5-cm-diameter cylindric insert with a density approximately equal to the average value in lung tissue (0.30 g/mL) is to be placed in the center of the phantom. For the measurements reported here, however, this cylinder, which was packed with foam beads, was inadvertently left dry, and thus had a mean density of only 0.025 g/mL. This provided a somewhat more stringent test of attenuation and scatter correction than was intended, but the CT-based attenuation correction appeared to work quite well for this configuration, with no observable artifacts. Activity is also present outside the FOV of the scanner in the form of a line source in the scatter phantom, which is placed adjacent to the body phantom on one side. Its effective specific activity should be the same as that of the main chamber of the body phantom.
The intention of this measurement is to produce and evaluate images approximating those obtained in a whole-body study with both hot and cold lesions. The background activity level in the torso cavity is supposed to represent the level one would have in a patient weighing 70 kg after an injection of 370 MBq of 18F-FDG (although no allowance is made for uptake time or excretion of part of the dose) (3). In our experience, however, the NECR of this phantom at a given activity concentration was not similar to that found in patients (7,8). The phantoms NECR peaked at a much higher value and at a much higher activity concentration than for average patients, probably as a result of the higher concentration of activity in the scanners FOV and perhaps lower attenuation than in patients. Thus, at the specified 5.3 kBq/mL, the NECR for this study was significantly higher than would typically be achieved in patients but significantly lower than the maximum that could be achieved for this phantom. It is arguable, then, that this test may capture neither typical clinical performance nor the optimal performance of an LSO PET/CT. Nevertheless, we present measurements made at close to the recommended levels.
The actual specific activities we achieved are listed in Table 1. It must be acknowledged that there are relative uncertainties in these assays, which we believe are on the order of ±5%. The hot-sphere contrast ratios depend on the relative accuracy of these assays and thus will exhibit corresponding uncertainties.
|
![]() | (Eq. 3) |
Based on equation 3, we used 7 min for the emission scan and 0.2 min for the CT and positioning scan. Although this conforms to the NU 2-2001 specification, it should be noted that in our clinic patient scan times for the emission acquisitions on an LSO PET/CT scanner range from 1 to 5 min per bed position, depending on patient weight and part of the body scanned, and whole-body acquisitions never require 60 min.
Data acquired were reconstructed with all corrections using the available clinical reconstruction software. NU 2-2001 specifies that results be reported for "standard" reconstruction parameters. However, the CPS LSO PET/CT offers users a wide choice of reconstruction algorithms and parameters, and different clinics prefer different standards. For this reason, we report results for 3 different reconstructions: one FBP (unwindowed) and 2 using OSEM iterative reconstruction, with either 2 iterations and 8 subsets or 4 iterations and 16 subsets. In all cases, FORE was used to reduce the 3D dataset to 2-dimensional, and a 5-mm full-width-at-half-maximum (FWHM) gaussian filter (axial and transaxial) was applied to the image after reconstruction. Image size was 256 x 256, pixel size was 1.77 mm, and the slice thickness was 3.375 mm. To measure the residual error in scatter and attenuation corrections, the relative error (
Clung,i) in percentage units was calculated for each slice (i) by taking the ratio of the average counts in the lung insert ROI to the average counts in the background ROIs. For the analysis of image quality and accuracy of scatter and attenuation corrections, we have followed the guidelines presented in the NU 2-2001 protocol.
NU 2-2001 suggests that the image quality measurement may be repeated and that the results can be averaged to achieve more precise results. We were concerned particularly with the repeatability of the contrast measurement. Rather than repeating scans at approximately the same activity level, we instead performed repeated scans over a range of activity levels for the n = 8 case, to assess the effect of image noise level on contrast measurement. The measured activity levels were 2.3, 1.0, 0.46, and 0.20 times the standard background activity level. For this study, only FBP reconstruction was used. The relative standard deviations of the contrast values across the activity range were computed.
| RESULTS |
|---|
|
|
|---|
|
Counting Rate Performance
Peak true counting rate was measured as 171 kcps at an activity concentration level of 21 kBq/mL. Peak NECR was 44.4 kcps at 14.7 kBq/mL for noiseless randoms correction and 31.4 kcps at 12.6 kBq/mL for online randoms subtraction. A plot of the trues, randoms, and NECRs versus activity is shown in Figure 2. It is important to realize that peak NECR occurs at lower activity concentrations in patients than it does in this phantom (8); thus, it is not accurate to infer optimal patient injected doses directly from these results. Counting rate analyses on clinical data on the LSO PET/CT scanner showed that peak NECR in an average adult patient weighing 70 kg was achieved with a mean activity concentration (assuming no excretion) of 6.7 kBq/mL (7).
|
|
21 kBq/mL or less. Although the intention of this measurement is to determine accuracy, not precision, it seems clear from the fluctuations of the deviations, as well as the symmetry between positive and negative values, that there remains a large element of statistical uncertainty in these results, probably because of the difficulty in determining the precise linear trend, despite the fact that the NU 2-2001 specification was followed here. Therefore, we also showed the average deviation versus aeff in the figure, which may give a better sense of the actual bias in the corrections. In NU 2-2001, bias is defined as the absolute value of the deviations (
ri,j). Over the clinical range of the scanner (the peak NEC activity [aNEC,peak] and below), the maximum bias is less than 8.5% and the average bias is less than 4%.
|
|
|
|
|
The results for the repeatability measurement of the contrast at the 8:1 ratio (Table 4) indicated that, except for the largest sphere, these values were repeatable to within 2%4%. No systematic trend in contrast versus activity was observed for these. The largest sphere, on the other hand, exhibited a monotonic increase in contrast value with activity and a much higher relative deviation than did the other sphere. The reason for this is not yet understood.
Figures 7A and 7B show the value of the lung residual (
Clung, i) for each slice, for the 4:1 and 8:1 contrast ratios, respectively. For each contrast, results for the 3 reconstructions are plotted in percentage units. Because of the construction of this phantom, its active region does not extend beyond plane 42 of the LSO PET/CT when positioned according to NU 2-2001, resulting in the drop in residuals at higher plane numbers. It can be observed that there is a strong dependence of the residual values on the reconstruction algorithm used, although the scatter and attenuation corrections are the same in all cases. The mean residuals over planes 142 are given in Tables 3 and 4. Although the average error in the case of FBP reconstruction was very close to 0, there was significant bias in the iterative reconstructions, probably as a result of incomplete convergence in this large cold region. The plane-to-plane variability of the residuals was likely the result of statistical fluctuations. We observed similar variance across planes in a background ROI equal in size to the lung region. Thus, in our data, the residual deviations shown in Figures 7A and 7B appeared to measure primarily noise and reconstruction convergence errors, not errors in scatter and attenuation correction as was intended by NEMA. We suggest that users exercise care in interpreting the results of this so-called test for attenuation and scatter correction.
|
| DISCUSSION |
|---|
|
|
|---|
The new standards also incorporate the contribution of out-of-field activity in both the counting rate and image quality tests, which is reflected in the increased scatter fraction for a 3D-only system. These new tests attempted to approximate more realistic clinical conditions for whole-body scanning than did the NU 2-1994 standard. Although the new standards are partially successful in this attempt, we find some limitations as well, particularly with regard to the relationship between image quality, NECR, and injected patient dose. Thus, although they certainly provide useful information to the user or the scanner manufacturer and make it possible to compare systems from different vendors, we believe they do not accurately reflect clinical performance in most whole-body studies (8,9).
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Yusuf E. Erdi, DSc, Department of Medical Physics, Memorial Sloan-Kettering Cancer Center, New York, NY 10021.
E-mail: erdiy{at}mskcc.org
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Belakhlef, C. Church, A. Hays, R. Fraser, and S. Lakhanpal Quantitative Assessment of the Influence of Location, Internal Temperature, Idle Time, and Normalization on the Sensitivity of a Mobile PET/CT Scanner J. Nucl. Med. Technol., September 1, 2008; 36(3): 147 - 150. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. C. Ong, H. Schoder, N. Y. Lee, S. G. Patel, D. Carlson, M. Fury, D. G. Pfister, J. P. Shah, S. M. Larson, and D. H. Kraus Clinical Utility of 18F-FDG PET/CT in Assessing the Neck After Concurrent Chemoradiotherapy for Locoregional Advanced Head and Neck Cancer J. Nucl. Med., April 1, 2008; 49(4): 532 - 540. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Belakhlef, C. Church, R. Fraser, and S. Lakhanpal Sensitivity and Daily Quality Control of a Mobile PET/CT Scanner Operating in 3-Dimensional Mode J. Nucl. Med. Technol., December 1, 2007; 35(4): 229 - 231. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. S. Kim, J. S. Lee, K. C. Im, S. J. Kim, S.-Y. Kim, D. S. Lee, and D. H. Moon Performance Measurement of the microPET Focus 120 Scanner J. Nucl. Med., September 1, 2007; 48(9): 1527 - 1535. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. C. Ong, H. Schoder, S. G. Patel, I. M. Tabangay-Lim, I. Doddamane, M. Gonen, A. R. Shaha, R. M. Tuttle, J. P. Shah, and S. M. Larson Diagnostic Accuracy of 18F-FDG PET in Restaging Patients with Medullary Thyroid Carcinoma and Elevated Calcitonin Levels J. Nucl. Med., April 1, 2007; 48(4): 501 - 507. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Are, J. F. Hsu, H. Schoder, J. P. Shah, S. M. Larson, and A. R. Shaha FDG-PET Detected Thyroid Incidentalomas: Need for Further Investigation? Ann. Surg. Oncol., January 1, 2007; 14(1): 239 - 247. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. J. Kemp, C. Kim, J. J. Williams, A. Ganin, and V. J. Lowe NEMA NU 2-2001 Performance Measurements of an LYSO-Based PET/CT System in 2D and 3D Acquisition Modes J. Nucl. Med., December 1, 2006; 47(12): 1960 - 1967. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Schoder, D. L. Carlson, D. H. Kraus, H. E. Stambuk, M. Gonen, Y. E. Erdi, H. W.D. Yeung, A. G. Huvos, J. P. Shah, S. M. Larson, et al. 18F-FDG PET/CT for Detecting Nodal Metastases in Patients with Oral Cancer Staged N0 by Clinical Examination and CT/MRI J. Nucl. Med., May 1, 2006; 47(5): 755 - 762. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Matsumoto, K. Kitamura, T. Mizuta, K. Tanaka, S. Yamamoto, S. Sakamoto, Y. Nakamoto, M. Amano, K. Murase, and M. Senda Performance Characteristics of a New 3-Dimensional Continuous-Emission and Spiral-Transmission High-Sensitivity and High-Resolution PET Camera Evaluated with the NEMA NU 2-2001 Standard J. Nucl. Med., January 1, 2006; 47(1): 83 - 90. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Brambilla, C. Secco, M. Dominietto, R. Matheoud, G. Sacchetti, and E. Inglese Performance Characteristics Obtained for a New 3-Dimensional Lutetium Oxyorthosilicate-Based Whole-Body PET/CT Scanner with the National Electrical Manufacturers Association NU 2-2001 Standard J. Nucl. Med., December 1, 2005; 46(12): 2083 - 2091. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. S. Halpern, M. Dahlbom, M. A. Auerbach, C. Schiepers, B. J. Fueger, W. A. Weber, D. H.S. Silverman, O. Ratib, and J. Czernin Optimizing Imaging Protocols for Overweight and Obese Patients: A Lutetium Orthosilicate PET/CT Study J. Nucl. Med., April 1, 2005; 46(4): 603 - 607. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. C. Watson, M. E. Casey, L. Eriksson, T. Mulnix, D. Adams, and B. Bendriem NEMA NU 2 Performance Tests for Scanners with Intrinsic Radioactivity J. Nucl. Med., May 1, 2004; 45(5): 822 - 826. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |