-
Notifications
You must be signed in to change notification settings - Fork 31
FAQ Frequently asked questions
This is not a bug!
Most statistical tools cannot plot specificity from 1 to 0, so instead they plot 1-specificity from 0 to 1. This results in exactly the same plot, just with a different labelling of the axis.
R doesn't have this limitation, and pROC plots specificity from 1 to 0, saving you from doing a few subtractions when you look at the plot. If you are using pROC in S+, you only have the choice to plot 1-specificity.
If you are concerned that your reviewers will be confused and really wish to have the "old-style" 1-Specificity from 0 to 1, or if you plan to change the axis label to "False Positive Rate" (which is really 1-specificity), you can use the legacy.axis argument to the plot function. See ?plot.roc for more details.
You need to use the predict() function to get predictions that can be passed to pROC (as the predictor argument).
Specifically you need to get a single vector of numeric predictions. In some cases you may have to add pass type="probabilities" or something similar to the predict function (to get probabilities instead of the binary response), sometimes you may to have select one of the column of the returned matrix (if predict returns a matrix of class probabilities). Read the documentation of the package you're using and look at your results to find out what to do exactly.
Not directly, but you can! You need to use a modeling function (such as glm()) or package (such as lme4); fit and check the model, then use the predict() function to get predictions that can be passed to pROC (as the predictor argument).
This is because pROC tries to auto-detect the direction of the decision due to the default value of direction="auto". The decision is based on the median of the groups instead of the AUC directly, so you the effect may not be 1:1, but you should expect to obtain AUC > 0.5 on average.
Whenever you are resampling with pROC you need to explicitly pass an explicit value to direction, either direction="<" or direction=">". I also recommend taking a minute to set the levels (value of the response for controls and cases, respectively) argument explicitly too.
A future version of pROC will be more verbose and print a message when auto-detecting the direction.
This is equivalent to asking whether the median values of mean and controls are significantly different. Therefore you can easily perform a Wilcoxon Rank Sum Test.
A large effort has been made to render pROC more and more efficient with large ROC curves. An algorithm with complexity of O(n) has been rolled out with pROC 1.6 and improved over time. Improvements to the DeLong test and the coords functions have been implemented too.
If you feel that a function is abnormally slow, especially in comparison to other packages, feel free to open a bug report.
Yes you can since pROC 1.6 and its ci.coords function. For instance:
> ci.coords(aSAH$outcome, aSAH$s100b, x="best", input = "threshold", ret="threshold")
95% CI (2000 stratified bootstrap replicates):
2.5% 50% 97.5%
threshold 0.115 0.205 0.4854
There is 95% chance that the optimal threshold lies between 0.115 and 0.4854.
Several packages on CRAN provide alternative roc or auc functions. These packages can interfere with pROC, especially if they are loaded later in the session (and hence appear earlier in the search path).
For instance, here are a few messages you may see if you have the AUC package loaded:
Not enough distinct predictions to compute area under the ROC curve.
Error in roc(outcome ~ ndka, data = aSAH) : unused argument (data = aSAH)
If that happens, you should unload the package by typing detach("package:AUC"). Alternatively you can refer to the pROC version of the function specifically through their namespace:
pROC::auc(pROC::roc(aSAH$outcome, aSAH$ndka))
Here is a list of packages providing roc, auc and ci functions which will hide pROC's functions if the package is loaded after pROC:
| Package | roc | auc | ci |
|---|---|---|---|
| analogue | Generic | ||
| asbio | Not Generic | ||
| AUC | Not Generic | Not Generic | |
| aucm | NA | NA | |
| cutpointr | Not Generic | ||
| distillery | Generic | ||
| enrichvs | Not Generic | ||
| epiDisplay | NA | Generic | |
| esvis | Not Generic | ||
| FFTrees | Not Generic | ||
| flux | Not Generic | ||
| fmsb | Not Generic | ||
| gmodels | Generic | ||
| grpreg | NA | ||
| hutils | Not Generic | ||
| kulife | Not Generic | ||
| longROC | Not Generic | Not Generic | |
| MAMSE | Not Generic | ||
| MESS | Not Generic | ||
| Metrics | Not Generic | ||
| MXM | Not Generic | ||
| ncvreg | NA | ||
| NEArender | Not Generic | ||
| PK | Not Generic | ||
| precrec | Generic | ||
| PresenceAbsence | Not Generic | ||
| radiant.model | Not Generic | ||
| RMediation | Not Generic | ||
| sdm | Generic | ||
| SDMTools | Not Generic | ||
| spatstat | Generic | ||
| StatMeasures | Not Generic | ||
| survMisc | Generic |
To get an updated list, run inst/extra/sos_clashes.R. It may be needed to start R with export R_MAX_NUM_DLLS=1000 to be able to run the script.