(Circulation. 2007;115:654-657.)
© 2007 American Heart Association, Inc.
Statistical Primer for Cardiovascular Research |
From Childrens Hospital Boston (K.H.Z.), Harvard Medical School (K.H.Z., A.J.O., L.M.), Brigham and Womens Hospital (L.M.), and Harvard Clinical Research Institute (L.M.), Boston, Mass.
Correspondence to Kelly H. Zou, PhD, Department of Health Care Policy, Harvard Medical School, 180 Longwood Ave, Boston, MA 02115. E-mail kelly.zou{at}childrens.harvard.edu
Key Words: diagnosis ROC curve sensitivity and specificity statistics tests
| Introduction |
|---|
|
|
|---|
ROC analysis is a useful tool for evaluating the performance of diagnostic tests and more generally for evaluating the accuracy of a statistical model (eg, logistic regression, linear discriminant analysis) that classifies subjects into 1 of 2 categories, diseased or nondiseased. Its function as a simple graphical tool for displaying the accuracy of a medical diagnostic test is one of the most well-known applications of ROC curve analysis. In Circulation from January 1, 1995, through December 5, 2005, 309 articles were published with the key phrase "receiver operating characteristic." In cardiology, diagnostic testing plays a fundamental role in clinical practice (eg, serum markers of myocardial necrosis, cardiac imaging tests). Predictive modeling to estimate expected outcomes such as mortality or adverse cardiac events based on patient risk characteristics also is common in cardiovascular research. ROC analysis is a useful tool in both of these situations.
In this article, we begin by reviewing the measures of accuracysensitivity, specificity, and area under the curve (AUC)that use the ROC curve. We also illustrate how these measures can be applied using the evaluation of a hypothetical new diagnostic test as an example.
| Diagnostic Test and Predictive Model |
|---|
|
|
|---|
2+ and >2+), in which case methods for binary outcomes can be used.7 Test data such as serum markers (brain natriuretic peptide11) or physiological markers (coronary lumen diameter,12 peak oxygen consumption13) also may be acquired on a continuous scale. | Gold Standard |
|---|
|
|
|---|
In selection of the gold standard, 2 potential problems arise: verification bias and measurement error. Verification bias results when the accuracy of a test is evaluated only among those with known disease status.1416 Measurement error may result when a true gold standard is absent or an imperfect standard is used for comparison.17,18
| Sensitivity and Specificity |
|---|
|
|
|---|
|
The accuracy of such binary-valued diagnostic tests is assessed in terms of the probability that the test correctly classifies a nondiseased subject as negative, namely the specificity (also known as the true negative rate), and the probability that the test correctly classifies a diseased subject as positive, namely the sensitivity (also known as the true positive rate) (Figure 1).
|
When evaluating a continuous-scale diagnostic test, we need to account for the changes of specificity and sensitivity when the test threshold t varies. One may wish to report the sum of sensitivity and specificity at the optimal threshold (discussed later in greater detail). However, because the optimal value of t may not be relevant to a particular application, it can be helpful to plot sensitivity and specificity over a range of values of interest, as is done with an ROC curve. This inherent tradeoff between sensitivity and specificity also can be demonstrated by varying the choice of threshold.
| ROC Analysis |
|---|
|
|
|---|
|
| Estimation Methods |
|---|
|
|
|---|
|
Parametric Methods
As an alternative to the nonparametric approach, parametric models such as the binormal model may be assumed (Figure 3).2125 The binormal model assumes that both measurements have 2 independent normal distributions with different means and SDs. In our example, the distributions have a mean of 0 and an SD of 1 for the nondiseased population and a mean of 1.87 and an SD of 1.5 for the diseased population. These models have the further advantage of allowing easy incorporation of covariates into the model. By incorporating an optimal transformation, typically a log transformation to normal distributions, the estimated ROC curve may yield a better fit.2628
| Summary Measures |
|---|
|
|
|---|
Area Under the Curve
The AUC is an overall summary of diagnostic accuracy. AUC equals 0.5 when the ROC curve corresponds to random chance and 1.0 for perfect accuracy. On rare occasions, the estimated AUC is <0.5, indicating that the test does worse than chance.31
For continuous diagnostic data, the nonparametric estimate of AUC is the Wilcoxon rank-sum test, namely the proportion of all possible pairs of nondiseased and diseased test subjects for which the diseased result is higher than the nondiseased one plus half the proportion of ties. Under the binormal model, the AUC is a simple function of the mean and variance.21,32
Comparison of AUC Curves
An important problem concerns the comparison of 2 AUCs derived from 2 diagnostic tests administered on the same set of patients. Correlated U statistics may be compared.33 Pearson correlation coefficients were used to estimate the correlation of the 2 AUCs.34 A family of nonparametric comparisons based on a weighted average of sensitivities may be conducted.35
Partial Area
The area under the ROC curve is a simple and convenient overall measure of diagnostic test accuracy. However, it gives equal weight to the full range of threshold values. When the ROC curves intersect, the AUC may obscure the fact that 1 test does better for 1 part of the scale (possibly for certain types of patients) whereas the other test does better over the remainder of the scale.32,36 The partial area may be useful for the range of specificity (or sensitivity) of clinical importance (ie, between 90% and 100% specificity). However, partial area may be more difficult to estimate and compare on the basis of numerical integration methods; thus, full area is used more frequently in practice.37
Optimal Threshold
One criterion for evaluating the optimal threshold of a test is to maximize the sum of sensitivity and specificity. This is equivalent to maximizing the difference between the sensitivity of the test and the sensitivity that the test would have if it did no better than random chance.9 For example, if both sensitivity and specificity are of importance in our example binormal model, the optimal threshold of t would be 0.75, where these 2 accuracy measures equal sensitivity and specificity equal 0.77 (Figure 3).
| Discussion |
|---|
|
|
|---|
| Acknowledgments |
|---|
Sources of Funding
This research was made possible in part by grants R01LM007861, R01GM074068, U41RR019703, and P41RR13218 from the National Institutes of Health (NIH), Bethesda, Md. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the NIH.
Disclosures
None.
| References |
|---|
|
|
|---|
2. Lloyd CJ. Using smooth receiver operating characteristic curves to summarize and compare diagnostic systems. J Am Stat Assoc. 1998; 93: 13561364.[CrossRef]
3. Zhou XH, Obuchowski NA, McClish DK. Statistical Methods in Diagnostic Medicine. New York, NY: Wiley & Sons; 2002.
4. Pepe MS. The Statistical Evaluation of Medical Tests for Classification and Prediction. Oxford, UK: Oxford University Press; 2003.
5. Campbell G. General methodology I: advances in statistical methodology for the evaluation of diagnostic and laboratory tests. Stat Med. 1994; 13: 499508.[Medline] [Order article via Infotrieve]
6. Shapiro DE. The interpretation of diagnostic tests. Stat Methods Med Res. 1999; 8: 113134.
7. Obuchowski NA. Receiver operating characteristic curves and their use in radiology. Radiology. 2003; 229: 38.
8. Eng J. Receiver operating characteristic analysis: a primer. Acad Radiol. 2005; 12: 909916.[CrossRef][Medline] [Order article via Infotrieve]
9. OMalley AJ, Zou KH, Fielding JR, Tempany CMC. Bayesian regression methodology for estimating a receiver operating characteristic curve with two radiologic applications: prostate biopsy and spiral CT of ureteral stone. Acad Radiol. 2001; 8: 713725.[CrossRef][Medline] [Order article via Infotrieve]
10. Lasko TA, Bhagwat JG, Zou KH, Ohno-Machado L. The use of receiver operating characteristic curves in biomedical informatics. J Biomed Inform. 2005; 38: 404415.[CrossRef][Medline] [Order article via Infotrieve]
11. Maisel A, Hollander JE, Guss D, McCollouph P, Nowak R, Green G, Saltzberg M, Ellison SR, Bhalla MA, Bhalla V, Clopton P, Jesse R, for the REDHOT Investigators. A multicenter study of B-type natriuretic peptide levels, emergency department decision making, and outcomes in patients presenting with shortness of breath. J Am Coll Cardiol. 2004; 44: 13281333.
12. Mauri L, Orav J, OMalley AJ, Moses JW, Leon MZB, Holmes DR, Teirstein PS, Schofer J, Breithardt G, Cutlip DE, Kereiakes DJ, Shi C, Firth BG, Donohoe DJ, Kuntz R. Relationship of late loss in lumen diameter to coronary restenosis in sirolimus-eluting stents. Circulation. 2005; 111: 321327.
13. ONeill J, Young JB, Pothier CE, Lauer MS. Peak oxygen consumption as a predictor of death in patient with heart failure receiving ß-blockers. Circulation. 2005; 111: 23132318.
14. Begg CB, Greenes RA. Assessment of diagnostic tests when disease verification is subject to selection bias. Biometrics. 1983; 39: 207215.[CrossRef][Medline] [Order article via Infotrieve]
15. Zhou XH, Higgs RE. Assessing the relative accuracies of two screening tests in the presence of verification bias. Stat Med. 2000; 19: 16971705.[CrossRef][Medline] [Order article via Infotrieve]
16. Toledano AY, Gatsonis C. Generalized estimating equations for ordinal categorical data: arbitrary patterns of missing responses and missingness in a key covariate. Biometrics. 1999; 55: 488496.[CrossRef][Medline] [Order article via Infotrieve]
17. Johnson WO, Gastwirth JL, Pearson LM. Screening without a "gold standard": the Hui-Walter paradigm revisited. Am J Epidemiol. 2001; 153: 921924.
18. Phelps CE, Hutson A. Estimating diagnostic test accuracy using a "fuzzy gold standard." Med Decis Making. 1995; 15: 4457.
19. Hsieh F, Turnbull BW. Nonparametric and semiparametric estimation of the receiver operating characteristic curve. Ann Stat. 1996; 24: 2440.
20. Zou KH, Hall WJ, Shapiro DE. Smooth nonparametric receiver operating characteristic (ROC) curves for continuous diagnostic tests. Stat Med. 1997; 16: 21432156.[CrossRef][Medline] [Order article via Infotrieve]
21. Dorfman DD, Alf E. Maximum likelihood estimation of parameters of signal detection theory: a direct solution. Psychometrika. 1968; 33: 117124.[CrossRef][Medline] [Order article via Infotrieve]
22. Metz CE, Herman BA, Shen J. Maximum-likelihood estimation of receiver operating characteristic (ROC) curves from continuous distributed data. Stat Med. 1998; 17: 10331053.[CrossRef][Medline] [Order article via Infotrieve]
23. Zou KH, Hall WJ. Two transformation models for estimating an ROC curve derived from continuous data. J Appl Stat. 2000; 27: 621631.[CrossRef]
24. Cai T, Moskowitz CS. Semi-parametric estimation of the binormal ROC curve for a continuous diagnostic test. Biostatistics. 2004; 5: 573586.[Abstract]
25. Zou KH, Wells WM 3rd, Kikinis R, Warfield K. Three validation metrics for automated probabilistic image segmentation of brain tumours. Stat Med. 2004; 23: 12591282.[CrossRef][Medline] [Order article via Infotrieve]
26. Hanley JA. The robustness of the "binormal" assumptions used in fitting ROC curves. Med Decis Making. 1988; 8: 197203.
27. Walsh SJ. Goodness-of-fit issues in ROC curve estimation. Med Decis Making. 1999; 19: 193201.
28. Zou KH, Resnic FS, Talos IF, Goldberg-Zimring D, Bhagwat JG, Haker SJ, Kikinis R, Jolesz FA, Ohno-Machado L. A global goodness-of-fit test for receiver operating characteristic curve analysis via the bootstrap method. J Biomed Inform. 2005; 38: 395403.[CrossRef][Medline] [Order article via Infotrieve]
29. Platt RW, Hanley JA, Yang H. Bootstrap confidence intervals for the sensitivity of a quantitative diagnostic test. Stat Med. 2000; 19: 313322.[CrossRef][Medline] [Order article via Infotrieve]
30. Zhou XH, Qin G. Improved confidence intervals for the sensitivity at a fixed level of specificity of a continuous-scale diagnostic test. Stat Med. 2005; 24: 465477.[CrossRef][Medline] [Order article via Infotrieve]
31. Hanley JA, McNeil BJ. The meaning and use of the area under a ROC curve. Radiology. 1982; 143: 2736.
32. McClish DK. Analyzing a portion of the ROC curve. Med Decis Making. 1989; 9: 190195.
33. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988; 44: 837845.[CrossRef][Medline] [Order article via Infotrieve]
34. Hanley JA, McNeil BJ. A method of comparing the areas under receiver operating characteristic curves derived from the same cases. Radiology. 1983; 148: 839843.
35. Weiand S, Gail MH, James BR, James KL. A family of nonparametric statistics for comparing diagnostic makers with paired or unpaired data. Biometrika. 1989; 76: 585592.
36. Dodd LE, Pepe MS. Partial AUC estimation and regression. Biometrics. 2003; 59: 614623.[CrossRef][Medline] [Order article via Infotrieve]
37. Walter SD. The partial area under the summary ROC curve. Stat Med. 2005; 24: 20252040.[CrossRef][Medline] [Order article via Infotrieve]
38. OMalley AJ, Zou KH. Bayesian multivariate hierarchical transformation models for ROC analysis. Stat Med. 2006; 25: 459479.[CrossRef][Medline] [Order article via Infotrieve]
39. Linnett K. Comparison of quantitative diagnostic tests: type I error, power and sample size. Stat Med. 1987; 6: 147158.[Medline] [Order article via Infotrieve]
40. Obuchowski NA, McClish DK. Sample size determination for diagnostic accuracy studies involving binormal ROC curve indices. Stat Med. 1997; 16: 15291542.[CrossRef][Medline] [Order article via Infotrieve]
41. Obuchowski NA. Sample size calculations in studies of test accuracy. Stat Methods Med Res. 1998; 7: 371392.
42. Eng J. Sample size estimation: a glimpse beyond simple formulas. Radiology. 2004; 230: 606612.
43. Moses LE, Shapiro DE, Littenberg B. Combining independent studies of a diagnostic test into a summary ROC curve: data-analytic approaches and some additional considerations. Stat Med. 1993; 12: 12931316.[Medline] [Order article via Infotrieve]
44. Rutter CM, Gatsonis C. A hierarchical regression approach to meta-analysis of diagnostic test accuracy evaluations. Stat Med. 2001; 20: 28652884.[CrossRef][Medline] [Order article via Infotrieve]
This article has been cited by other articles:
![]() |
H Adams, A Tzankov, A Lugli, and I Zlobec New time-dependent approach to analyse the prognostic significance of immunohistochemical biomarkers in colon cancer and diffuse large B-cell lymphoma J. Clin. Pathol., November 1, 2009; 62(11): 986 - 997. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. S. Windle, G. A. Rose, R. Devillers, and M.-J. Fortin Exploring spatial non-stationarity of fisheries survey data using geographically weighted regression (GWR): an example from the Northwest Atlantic ICES J. Mar. Sci., September 4, 2009; (2009) fsp224v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Price Bedside Evaluation of Thienopyridine Antiplatelet Therapy Circulation, May 19, 2009; 119(19): 2625 - 2632. [Full Text] [PDF] |
||||
![]() |
F. Drago, M. S. Russo, M. S. Silvetti, A. De Santis, and M. T. Naso Onofrio 'Time to effect' during cryomapping: a parameter related to the long-term success of accessory pathways cryoablation in children Europace, May 1, 2009; 11(5): 630 - 634. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Nambi, R. C. Hoogeveen, L. Chambless, Y. Hu, H. Bang, J. Coresh, H. Ni, E. Boerwinkle, T. Mosley, R. Sharrett, et al. Lipoprotein-Associated Phospholipase A2 and High-Sensitivity C-Reactive Protein Improve the Stratification of Ischemic Stroke Risk in the Atherosclerosis Risk in Communities (ARIC) Study Stroke, February 1, 2009; 40(2): 376 - 381. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. L. Benowitz, J. T. Bernert, R. S. Caraballo, D. B. Holiday, and J. Wang Optimal Serum Cotinine Levels for Distinguishing Cigarette Smokers and Nonsmokers Within Different Racial/Ethnic Groups in the United States Between 1999 and 2004 Am. J. Epidemiol., January 15, 2009; 169(2): 236 - 248. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Soreide Receiver-operating characteristic curve analysis in diagnostic, prognostic and predictive biomarker research J. Clin. Pathol., January 1, 2009; 62(1): 1 - 5. [Full Text] [PDF] |
||||
![]() |
C. Velik-Salchner, S. Maier, P. Innerhofer, W. Streif, A. Klingler, C. Kolbitsch, and D. Fries Point-of-Care Whole Blood Impedance Aggregometry Versus Classical Light Transmission Aggregometry for Detecting Aspirin and Clopidogrel: The Results of a Pilot Study Anesth. Analg., December 1, 2008; 107(6): 1798 - 1806. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Miller, C. E. Rochitte, M. Dewey, A. Arbab-Zadeh, H. Niinuma, I. Gottlieb, N. Paul, M. E. Clouse, E. P. Shapiro, J. Hoe, et al. Diagnostic Performance of Coronary Angiography by 64-Row CT N. Engl. J. Med., November 27, 2008; 359(22): 2324 - 2336. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. C. Saumarez, M. Pytkowski, M. Sterlinski, J. P. Bourke, J. R. Clague, S. M. Cobbe, D. T. Connelly, M. J. Griffith, P. P. McKeown, K. McLeod, et al. Paced ventricular electrogram fractionation predicts sudden cardiac death in hypertrophic cardiomyopathy Eur. Heart J., July 1, 2008; 29(13): 1653 - 1661. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Rossing, H. Mischak, M. Dakna, P. Zurbig, J. Novak, B. A. Julian, D. M. Good, J. J. Coon, L. Tarnow, P. Rossing, et al. Urinary Proteomics in Diabetes and CKD J. Am. Soc. Nephrol., July 1, 2008; 19(7): 1283 - 1290. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Levrat, A. Gros, L. Rugeri, K. Inaba, B. Floccard, C. Negrier, and J.-S. David Evaluation of rotation thrombelastography for the diagnosis of hyperfibrinolysis in trauma patients Br. J. Anaesth., June 1, 2008; 100(6): 792 - 797. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Hachamovitch and M. F. Di Carli Methods and Limitations of Assessing New Noninvasive Tests: Part I: Anatomy-Based Validation of Noninvasive Testing Circulation, May 20, 2008; 117(20): 2684 - 2690. [Full Text] [PDF] |
||||
![]() |
J. Schwitter, C. M. Wacker, A. C. van Rossum, M. Lombardi, N. Al-Saadi, H. Ahlstrom, T. Dill, H. B.W. Larsson, S. D. Flamm, M. Marquardt, et al. MR-IMPACT: comparison of perfusion-cardiac magnetic resonance with single-photon emission computed tomography for the detection of coronary artery disease in a multicentre, multivendor, randomized trial Eur. Heart J., February 2, 2008; 29(4): 480 - 489. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. L. del Rio, P. I. McConnell, M. Kukielka, R. Dzwonczyk, B. D. Clymer, M. B. Howie, and G. E. Billman Electrotonic remodeling following myocardial infarction in dogs susceptible and resistant to sudden cardiac death J Appl Physiol, February 1, 2008; 104(2): 386 - 393. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. R. Hunsaker, K. H. Zou, A. C. Poh, B. Trotman-Dickenson, F. L. Jacobson, R. R. Gill, and S. Z. Goldhaber Routine Pelvic and Lower Extremity CT Venography in Patients Undergoing Pulmonary CT Angiography Am. J. Roentgenol., February 1, 2008; 190(2): 322 - 326. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. R. Cook Statistical Evaluation of Prognostic versus Diagnostic Models: Beyond the ROC Curve Clin. Chem., January 1, 2008; 54(1): 17 - 23. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2007 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |