Abstract
Background
Receiver operating characteristic (ROC) curves are useful tools to evaluate classifiers in biomedical and bioinformatics applications. However, conclusions are often reached through inconsistent use or insufficient statistical analysis. To support researchers in their ROC curves analysis we developed pROC, a package for R and S+ that contains a set of tools displaying, analyzing, smoothing and comparing ROC curves in a userfriendly, objectoriented and flexible interface.
Results
With data previously imported into the R or S+ environment, the pROC package builds ROC curves and includes functions for computing confidence intervals, statistical tests for comparing total or partial area under the curve or the operating points of different classifiers, and methods for smoothing ROC curves. Intermediary and final results are visualised in userfriendly interfaces. A case study based on published clinical and biomarker data shows how to perform a typical ROC analysis with pROC.
Conclusions
pROC is a package for R and S+ specifically dedicated to ROC analysis. It proposes multiple statistical tests to compare ROC curves, and in particular partial areas under the curve, allowing proper ROC interpretation. pROC is available in two versions: in the R programming language or with a graphical user interface in the S+ statistical software. It is accessible at http://expasy.org/tools/pROC/ webcite under the GNU General Public License. It is also distributed through the CRAN and CSAN public repositories, facilitating its installation.
Background
A ROC plot displays the performance of a binary classification method with continuous or discrete ordinal output. It shows the sensitivity (the proportion of correctly classified positive observations) and specificity (the proportion of correctly classified negative observations) as the output threshold is moved over the range of all possible values. ROC curves do not depend on class probabilities, facilitating their interpretation and comparison across different data sets. Originally invented for the detection of radar signals, they were soon applied to psychology [1] and medical fields such as radiology [2]. They are now commonly used in medical decision making, bioinformatics[3], data mining and machine learning, evaluating biomarker performances or comparing scoring methods [2,4].
In the ROC context, the area under the curve (AUC) measures the performance of a classifier and is frequently applied for method comparison. A higher AUC means a better classification. However, comparison between AUCs is often performed without a proper statistical analysis partially due to the lack of relevant, accessible and easytouse tools providing such tests. Small differences in AUCs can be significant if ROC curves are strongly correlated, and without statistical testing two AUCs can be incorrectly labelled as similar. In contrast a larger difference can be non significant in small samples, as shown by Hanczar et al. [5], who also provide an analytical expression for the variance of AUC's as a function of the sample size. We recently identified this lack of proper statistical comparison as a potential cause for the poor acceptance of biomarkers as diagnostic tools in medical applications [6]. Evaluating a classifier by means of total AUC is not suitable when the performance assessment only takes place in high specificity or high sensitivity regions [6]. To account for these cases, the partial AUC (pAUC) was introduced as a local comparative approach that focuses only on a portion of the ROC curve [79].
Software for ROC analysis already exists. A previous review [10] compared eight ROC programs and found that there is a need for a tool performing valid and standardized statistical tests with good data import and plot functions.
The R [11] and S+ (TIBCO Spotfire S+ 8.2, 2010, Palo Alto, CA) statistical environments provide an extensible framework upon which software can be built. No ROC tool is implemented in S+ yet while four R packages computing ROC curves are available:
1) ROCR [12] provides tools computing the performance of predictions by means of precision/recall plots, lift charts, cost curves as well as ROC plots and AUCs. Confidence intervals (CI) are supported for ROC analysis but the user must supply the bootstrapped curves.
2) The verification package [13] is not specifically aimed at ROC analysis; nonetheless it can plot ROC curves, compute the AUC and smooth a ROC curve with the binomial model. A Wilcoxon test for a single ROC curve is also implemented, but no test comparing two ROC curves is included.
3) Bioconductor includes the ROC package [14] which can only compute the AUC and plot the ROC curve.
4) Pcvsuite [15] is an advanced package for ROC curves which features advanced functions such as covariate adjustment and ROC regression. It was originally designed for Stata and ported to R. It is not available on the CRAN (comprehensive R archive network), but can be downloaded for Windows and MacOS from http://labs.fhcrc.org/pepe/dabs/rocbasic.html webcite.
Table 1 summarizes the differences between these packages. Only pcvsuite enables the statistical comparison between two ROC curves. Pcvsuite, ROCR and ROC can compute AUC or pAUC, but the pAUC can only be defined as a portion of specificity.
Table 1. Features of the R packages for ROC anaylsis
The pROC package was designed in order to facilitate ROC curve analysis and apply proper statistical tests for their comparison. It provides a consistent and userfriendly set of functions building and plotting a ROC curve, several methods smoothing the curve, computing the full or partial AUC over any range of specificity or sensitivity, as well as computing and visualizing various CIs. It includes tests for the statistical comparison of two ROC curves as well as their AUCs and pAUCs. The software comes with an extensive documentation and relies on the underlying R and S+ systems for data input and plots. Finally, a graphical user interface (GUI) was developed for S+ for users unfamiliar with programming.
Implementation
AUC and pAUC
In pROC, the ROC curves are empirical curves in the sensitivity and specificity space. AUCs are computed with trapezoids [4]. The method is extended for pAUCs by ignoring trapezoids outside the partial range and adding partial trapezoids with linear interpolation when necessary. The pAUC region can be defined either as a portion of specificity, as originally described by McClish [7], or as a portion of sensitivity, as proposed later by Jiang et al. [8]. Any section of the curve pAUC(t_{0}, t_{1}) can be analyzed, and not only portions anchored at 100% specificity or 100% sensitivity. Optionally, pAUC can be standardized with the formula by McClish [7]:
where min is the pAUC over the same region of the diagonal ROC curve, and max is the pAUC over the same region of the perfect ROC curve. The result is a standardized pAUC which is always 1 for a perfect ROC curve and 0.5 for a nondiscriminant ROC curve, whatever the partial region defined.
Comparison
Two ROC curves are "paired" (or sometimes termed "correlated" in the literature) if they derive from multiple measurements on the same sample. Several tests exist to compare paired [1622] or unpaired [23] ROC curves. The comparison can be based on AUC [1619,21], ROC shape [20,22,23], a given specificity [15] or confidence bands [3,24]. Several tests are implemented in pROC. Three of them are implemented without modification from the literature [17,20,23], and the others are based on the bootstrap percentile method.
The bootstrap test to compare AUC or pAUC in pROC implements the method originally described by Hanley and McNeil [16]. They define Z as
where θ_{1 }and θ_{2 }are the two (partial) AUCs. Unlike Hanley and McNeil, we compute sd(θ_{1 } θ_{2}) with N (defaults to 2000) bootstrap replicates. In each replicate r, the original measurements are resampled with replacement; both new ROC curves corresponding to this new sample are built, the resampled AUCs θ_{1,r }and θ_{2,r }and their difference D_{r }= θ_{1,r } θ_{2,r }are computed. Finally, we compute sd(θ_{1 } θ_{2}) = sd(D). As Z approximately follows a normal distribution, one or twotailed pvalues are calculated accordingly. This bootstrap test is very flexible and can be applied to AUC, pAUC and smoothed ROC curves.
Bootstrap is stratified by default; in this case the same number of case and control observations than in the original sample will be selected in each bootstrap replicate. Stratification can be disabled and observations will be resampled regardless of their class labels. Repeats for the bootstrap and progress bars are handled by the plyr package [25].
The second method to compare AUCs implemented in pROC was developed by DeLong et al. [17] based on Ustatistics theory and asymptotic normality. As this test does not require bootstrapping, it runs significantly faster, but it cannot handle pAUC or smoothed ROC curves. For both tests, since the variance depends on the covariance of the ROC curves (Equation 3), strongly correlated ROC curves can have similar AUC values and still be significantly different.
Venkatraman and Begg [20] and Venkatraman [23] introduced tests to compare two actual ROC curves as opposed to their respective AUCs. Their method evaluates the integrated absolute difference between the two ROC curves, and a permutation distribution is generated to compute the statistical significance of this difference. As the measurements leading to the two ROC curves may be performed on different scales, they are not generally exchangeable between two samples. Therefore, the permutations are based on ranks, and ranks are recomputed as described in [20] to break the ties generated by the permutation.
Finally a test based on bootstrap is implemented to compare the ROC curve at a given level of specificity or sensitivity as proposed by Pepe et al. [15]. It works similar to the (p)AUC test, but instead of computing the (p)AUC at each iteration, the sensitivity (or specificity) corresponding to the given specificity (or respectively sensitivity) is computed. This test is equivalent to a pAUC test with a very small pAUC range.
Confidence intervals
CIs are computed with Delong's method [17] for AUCs and with bootstrap for pAUCs [26]. The CIs of the thresholds or the sensitivity and specificity values are computed with bootstrap resampling and the averaging methods described by Fawcett [4]. In all bootstrap CIs, patients are resampled and the modified curve is built before the statistics of interest is computed. As in the bootstrap comparison test, the resampling is done in a stratified manner by default.
Smoothing
Several methods to smooth a ROC curve are also implemented. Binormal smoothing relies on the assumption that there exists a monotone transformation to make both case and control values normally distributed [2]. Under this condition a simple linear relationship (Equation 4) holds between the normal quantile function (φ) values of sensitivities and specificities. In our implementation, a linear regression between all quantile values defines a and b, which then define the smoothed curve.
This is different from the method described by Metz et al. [27] who use maximum likelihood estimation of a and b. Binormal smoothing was previously shown to be robust and to provide good fits in many situations even when the deviation from basic assumptions is quite strong [28]. For continuous data we also include methods for kernel (density) smoothing [29], or to fit various known distributions to the class densities with fitdistr in the MASS package [30]. If a user would like to run a custom smoothing algorithm that is optimized for the analysed data, then pROC also accepts class densities or the customized smoothing function as input. CI and statistical tests of smoothed AUCs are done with bootstrap.
Results and Discussion
We first evaluate the accuracy of the ROC comparison tests. Results in Additional File 1 show that all unpaired tests give uniform pvalues under a null hypothesis (Additional Files 1 and 2) and that there is a very good correlation between DeLong's and bootstrap tests (Additional Files 1 and 3). The relation between Venkatraman's and the other tests is also investigated (Additional Files 1 and 4).
Additional file 1. Assessment of the ROC comparison tests. We evaluate the uniformity of the tests under the null hypothesis (ROC curves are not different), and the correlation between the different tests.
Format: PDF Size: 15KB Download file
This file can be viewed with: Adobe Acrobat Reader
Additional file 2. Histograms of the frequency of 600 test pvalues under the null hypothesis (ROC curves are not different). A: DeLong's paired test, B: DeLong's unpaired test, C: bootstrap paired test (with 10000 replicates), D: bootstrap unpaired test (with 10000 replicates) and E: Venkatraman's test (with 10000 permutations).
Format: PDF Size: 23KB Download file
This file can be viewed with: Adobe Acrobat Reader
Additional file 3. Correlations between DeLong and bootstrap paired tests. X axis: DeLong's test; Yaxis: bootstrap test with number of bootstrap replicates. A: 10, B: 100, C: 1000 and D: 10000.
Format: PDF Size: 101KB Download file
This file can be viewed with: Adobe Acrobat Reader
Additional file 4. Correlation between DeLong and Venkatraman's test. X axis: DeLong's test; Yaxis: Venkatraman's test with 10000 permutations.
Format: PDF Size: 33KB Download file
This file can be viewed with: Adobe Acrobat Reader
We now present how to perform a typical ROC analysis with pROC. In a recent study [31], we analyzed the level of several biomarkers in the blood of patients at hospital admission after aneurysmal subarachnoid haemorrhage (aSAH) to predict the 6month outcome. The 141 patients collected were classified according to their outcome with a standard neurological scale, the Glasgow outcome scale (GOS). The biomarker performances were compared with the well established neurological scale of the World Federation of Neurological Surgeons (WFNS), also obtained at admission.
Case study on clinical aSAH data
The purpose of the case presented here is to identify patients at risk of poor postaSAH outcome, as they require specific healthcare management; therefore the clinical test must be highly specific. Detailed results of the study are reported in [31]. We only outline the features relevant to the ROC analysis.
ROC curves were generated in pROC for five biomarkers (HFABP, S100β, Troponin I, NKDA and UFD1) and three clinical factors (WFNS, Modified Fisher score and age).
AUC and pAUC
Since we are interested in a clinical test with a high specificity, we focused on partial AUC between 90% and 100% specificity.
The best pAUC is obtained by WFNS, with 3.1%, closely followed by S100β with 3.0% (Figure 1). A perfect clinical test within the same region corresponds to a pAUC of 10%, while a ROC curve without any discrimination power would yield only 0.5%. In the case of WFNS, we computed a standardized pAUC of 63.7% with McClish's formula (Equation 1). Of these 63.9%, 50% are due to the small portion (0.5% nonstandardized) of the ROC curve below the identity line, and the remaining 13.9% are made of the larger part (2.6% nonstandardized) above the curve. In the R version of pROC, the standardized pAUC of WFNS can be computed with:
roc(response = aSAH$outcome, predictor = aSAH$wfns, partial.auc = c(100, 90), partial.auc.correct = TRUE, percent = TRUE)
Figure 1. ROC curves of WFNS and S100β. ROC curves of WFNS (blue) and S100β (green). The black bars are the confidence intervals of WFNS for the threshold 4.5 and the light green area is the confidence interval shape of S100β. The vertical light grey shape corresponds to the pAUC region. The pAUC of both empirical curves is printed in the middle of the plot, with the pvalue of the difference computed by a bootstrap test on the right.
In the rest of this paper, we report only not standardized pAUCs.
CI
Given the pAUC of WFNS, it makes sense to compute a 95% CI of the pAUC to assess the variability of the measure. In this case, we performed 10000 bootstrap replicates and obtained the 1.65.0% interval. In our experience, 10000 replicates give a fair estimate of the second significant digit. A lower number of replicates (for example 2000, the default) gives a good estimate of the first significant digit only. Other confidence intervals can be computed. The threshold with the point farthest to the diagonal line in the specified region was determined with pROC to be 4.5 with the coords function. A rectangular confidence interval can be computed and the bounds are 89.098.9 in specificity and 26.054.0 in sensitivity (Figure 1). If the variability of sensitivity at 90% specificity is considered more relevant than at a specific threshold, the interval of sensitivity is computed as 32.868.8. As shown in Figure 1 for S100β, a CI shape can be obtained by simply computing the CI's of the sensitivities over several constantly spaced levels of specificity, and these CI bounds are then joined to generate the shape. The following R code calculates the confidence shape:
plot(x = roc(response = aSAH$outcome, predictor = aSAH$s100, percent = TRUE, ci = TRUE, of = "se", sp = seq(0, 100, 5)), ci.type="shape")
The confidence intervals of a threshold or of a predefined level of sensitivity or specificity answer different questions. For instance, it would be wrong to compute the CI of the threshold 4.5 and report only the CI bound of sensitivity without reporting the CI bound of specificity as well. Similarly, determining the sensitivity and specificity of the cutoff 4.5 and then computing both CIs separately would also be inaccurate.
Statistical comparison
The second best pAUC is that of S100β with 3.0%. The difference to WFNS is very small and the bootstrap test of pROC indicates that it is not significant (p = 0.8, Figure 1). Surprisingly, a Venkatraman's test (over the total ROC curve) indicates a difference in the shape of the ROC curves (p = 0.004), and indeed a test evaluating pAUCs in the high sensitivity region (90100% sensitivity) would highlight a significant difference (p = 0.005, pAUC = 4.3 and 1.4 for WFNS and S100β respectively). However, since we are not interested in the high sensitivity region of the AUC there is no significant difference between WFNS and S100β.
In pROC pairwise comparison of ROC curves is implemented. Multiple testing is not accounted for and in the event of running several tests, the user is reminded that as with any statistical test, multiple tests should be performed with care, and if necessary appropriate corrections should be applied [32].
The bootstrap test can be performed with the following code in R:
roc.test(response = aSAH$outcome, predictor1 = aSAH$wfns, predictor2 = aSAH$s100, partial.auc = c(100, 90), percent = TRUE)
Smoothing
Whether or not to smooth a ROC curve is a difficult choice. It can be useful in ROC curves with only few points, in which the trapezoidal rule consistently underestimates the true AUC [17]. This is the case with most clinical scores, such as the WFNS shown in Figure 2 where three smoothing methods available in pROC are plotted: (i) normal distribution fitting, (ii) density and (iii) binormal. In our case study:
(i) The normal fitting (red) gives a significantly lower AUC estimate (Δ = 5.1, p = 0.0006, Bootstrap test). This difference is due to the nonnormality of WFNS. Distribution fitting can be very powerful when there is a clear knowledge of the underlying distributions, but should be avoided in other contexts.
(ii) The density (green) smoothing also produces a lower (Δ = 1.5, p = 6*10^{7}) AUC. It is interesting to note that even with a smaller difference in AUCs, the pvalue can be more significant due to a higher covariance.
(iii) The binormal smoothing (blue) gives a slightly but not significantly higher AUC than the empirical ROC curve (Δ = +2.4, p = 0.3). It is probably the best of the 3 smoothing estimates in this case (as mentioned earlier we were expecting a higher AUC as the empirical AUC of WFNS was underestimated). For comparison, Additional File 5 displays both our implementation of binormal smoothing with the one implemented in pcvsuite [15].
Figure 2. ROC curve of WFNS and smoothing. Empirical ROC curve of WFNS is shown in grey with three smoothing methods: binormal (blue), density (green) and normal distribution fit (red).
Additional file 5. Binormal smoothing. Binormal smoothing with pcvsuite (green, solid) and pROC (black, dashed).
Format: PDF Size: 66KB Download file
This file can be viewed with: Adobe Acrobat Reader
Figure 3 shows how to create a plot with multiple smoothed curves with pROC in S+. One loads the pROC library within S+, selects the new ROC curve item in the Statistics menu, selects the data on which the analysis is to be performed, and then moves to the Smoothing tab to set parameters for smoothing.
Figure 3. Screenshot of pROC in S+ for smoothing WFNS ROC curve. Top left: the General tab, where data is entered. Top right: the details about smoothing. Bottom left: the details for the plot. Checking the box "Add to existing plot" allows drawing several curves on a plot. Bottom right: the result in the standard S+ plot device.
Conclusion
In this case study we showed how pROC could be run for ROC analysis. The main conclusion drawn from this analysis is that none of the measured biomarkers can predict the patient outcome better than the neurological score (WFNS).
Installation and usage
R
pROC can be installed in R by issuing the following command in the prompt:
install.packages("pROC")
Loading the package:
library(pROC)
Getting help:
?pROC
S+
pROC is available from the File menu, item Find Packages.... It can be loaded from the File menu, item Load Library....
In addition to the command line functions, a GUI is then available in the Statistics menu. It features one window for univariate ROC curves (which contains options for smoothing, pAUC, CIs and plotting) and two windows for paired and unpaired tests of two ROC curves. In addition a specific help file for the GUI is available from the same menu.
Functions and methods
A summary of the functions available to the user in the command line version of pROC is shown in Table 2. Table 3 shows the list of the methods provided for plotting and printing.
Conclusions
The pROC package is a powerful set of tools analyzing and comparing ROC curves in R and S+. Unlike existing packages such as ROCR or verification, it is solely dedicated to ROC analysis, but provides in our knowledge the most complete set of statistical tests and plots for ROC curves. As shown in the case study reported here, pROC features the computation of AUC and pAUC, various kinds of confidence intervals, several smoothing methods, and the comparison of two paired or unpaired ROC curves. We believe that pROC should provide researchers, especially in the biomarker community, with the necessary tools to better interpret their results in biomarker classification studies.
pROC is available in two versions for R and S+. A thorough documentation with numerous examples is provided in the standard R format. For users unfamiliar with programming, a graphical user interface is provided for S+.
Availability and requirements
• Project name: pROC
• Project home page: http://expasy.org/tools/pROC webcite/
• Operating system(s): Platform independent
• Programming language: R and S+
• Other requirements: R ≥ 2.10.0 or S+ ≥ 8.1.1
• License: GNU GPL
• Any restrictions to use by nonacademics: none
List of abbreviations
aSAH: aneurysmal subarachnoid haemorrhage; AUC: area under the curve; CI: confidence interval; CRAN: comprehensive R archive network; CSAN: comprehensive SPLUS archive network; pAUC: partial area under the curve; ROC: receiver operating characteristic.
Authors' contributions
XR carried out the programming and software design and drafted the manuscript. NTu, AH, NTi provided data and biological knowledge, tested and critically reviewed the software and the manuscript. FL helped to draft and to critically improve the manuscript. JCS conceived the biomarker study, participated in its design and coordination, and helped to draft the manuscript. MM participated in the design and coordination of the bioinformatics part of the study, participated in the programming and software design and helped to draft the manuscript. All authors read and approved the final manuscript.
Acknowledgements
The authors would like to thank E. S. Venkatraman and Colin B. Begg for their support in the implementation of their test.
This work was supported by Proteome Science Plc.
References

Swets JA: The Relative Operating Characteristic in Psychology.
Science 1973, 182:9901000. PubMed Abstract  Publisher Full Text

Pepe MS: The statistical evaluation of medical tests for classification and prediction. Oxford: Oxford University Press; 2003.

Sonego P, Kocsor A, Pongor S: ROC analysis: applications to the classification of biological sequences and 3D structures.
Brief Bioinform 2008, 9:198209. PubMed Abstract  Publisher Full Text

Fawcett T: An introduction to ROC analysis.
Pattern Recogn Lett 2006, 27:861874. Publisher Full Text

Hanczar B, Hua J, Sima C, Weinstein J, Bittner M, Dougherty ER: Smallsample precision of ROCrelated estimates.
Bioinformatics 2010, 26:822830. PubMed Abstract  Publisher Full Text

Robin X, Turck N, Hainard A, Lisacek F, Sanchez JC, Müller M: Bioinformatics for protein biomarker panel classification: What is needed to bring biomarker panels into in vitro diagnostics?
Expert Rev Proteomics 2009, 6:675689. PubMed Abstract  Publisher Full Text

McClish DK: Analyzing a Portion of the ROC Curve.
Med Decis Making 1989, 9:190195. PubMed Abstract  Publisher Full Text

Jiang Y, Metz CE, Nishikawa RM: A receiver operating characteristic partial area index for highly sensitive diagnostic tests.
Radiology 1996, 201:745750. PubMed Abstract  Publisher Full Text

Streiner DL, Cairney J: What's under the ROC? An introduction to receiver operating characteristics curves.
Canadian Journal of Psychiatry Revue Canadienne De Psychiatrie 2007, 52:121128. PubMed Abstract

Stephan C, Wesseling S, Schink T, Jung K: Comparison of Eight Computer Programs for ReceiverOperating Characteristic Analysis.
Clin Chem 2003, 49:433439. PubMed Abstract  Publisher Full Text

R Development Core Team: R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2010.

Sing T, Sander O, Beerenwinkel N, Lengauer T: ROCR: visualizing classifier performance in R.
Bioinformatics 2005, 21:39403941. PubMed Abstract  Publisher Full Text

NCAR: verification: Forecast verification utilities v. 1.31. [http://CRAN.Rproject.org/package=verification] webcite

Carey V, Redestig H: ROC: utilities for ROC, with uarray focus, v. 1.24.0. [http://www.bioconductor.org] webcite

Pepe M, Longton G, Janes H: Estimation and Comparison of Receiver Operating Characteristic Curves.
The Stata journal 2009, 9:1. PubMed Abstract  PubMed Central Full Text

Hanley JA, McNeil BJ: A method of comparing the areas under receiver operating characteristic curves derived from the same cases.
Radiology 1983, 148:839843. PubMed Abstract  Publisher Full Text

DeLong ER, DeLong DM, ClarkePearson DL: Comparing the Areas under Two or More Correlated Receiver Operating Characteristic Curves: A Nonparametric Approach.
Biometrics 1988, 44:837845. PubMed Abstract  Publisher Full Text

Bandos AI, Rockette HE, Gur D: A permutation test sensitive to differences in areas for comparing ROC curves from a paired design.
Stat Med 2005, 24:28732893. PubMed Abstract  Publisher Full Text

Braun TM, Alonzo TA: A modified sign test for comparing paired ROC curves.
Biostat 2008, 9:364372. Publisher Full Text

Venkatraman ES, Begg CB: A distributionfree procedure for comparing receiver operating characteristic curves from a paired experiment.
Biometrika 1996, 83:835848. Publisher Full Text

Bandos AI, Rockette HE, Gur D: A Permutation Test for Comparing ROC Curves in Multireader Studies: A Multireader ROC, Permutation Test.
Acad Radiol 2006, 13:414420. PubMed Abstract  Publisher Full Text

Moise A, Clement B, Raissis M: A test for crossing receiver operating characteristic (roc) curves.
Communications in Statistics  Theory and Methods 1988, 17:19852003. Publisher Full Text

Venkatraman ES: A Permutation Test to Compare Receiver Operating Characteristic Curves.
Biometrics 2000, 56:11341138. PubMed Abstract  Publisher Full Text

Campbell G: Advances in statistical methodology for the evaluation of diagnostic and laboratory tests.
Stat Med 1994, 13:499508. PubMed Abstract  Publisher Full Text

Wickham H: plyr: Tools for splitting, applying and combining data v. 1.4. [http://CRAN.Rproject.org/package=plyr] webcite

Carpenter J, Bithell J: Bootstrap confidence intervals: when, which, what? A practical guide for medical statisticians.
Stat Med 2000, 19:11411164. PubMed Abstract  Publisher Full Text

Metz CE, Herman BA, Shen JH: Maximum likelihood estimation of receiver operating characteristic (ROC) curves from continuouslydistributed data.
Stat Med 1998, 17:10331053. PubMed Abstract  Publisher Full Text

Hanley JA: The robustness of the "binormal" assumptions used in fitting ROC curves.
Med Decis Making 1988, 8:197203. PubMed Abstract  Publisher Full Text

Zou KH, Hall WJ, Shapiro DE: Smooth nonparametric receiver operating characteristic (ROC) curves for continuous diagnostic tests.
Stat Med 1997, 16:21432156. PubMed Abstract  Publisher Full Text

Venables WN, Ripley BD: Modern Applied Statistics with S. Fourth edition. New York: Springer; 2002.

Turck N, Vutskits L, SanchezPena P, Robin X, Hainard A, GexFabry M, Fouda C, Bassem H, Mueller M, Lisacek F, et al.: A multiparameter panel method for outcome prediction following aneurysmal subarachnoid hemorrhage.
Intensive Care Med 2010, 36:107115. PubMed Abstract  Publisher Full Text

Ewens WJ, Grant GR: Statistics (i): An Introduction to Statistical Inference. In Statistical methods in bioinformatics. New York: SpringerVerlag; 2005.