ROC.plot: Estimated ROC curves

Description Usage Arguments Details Value Note Author(s) References See Also Examples

View source: R/ProbYX.R


Plot of ROC curves estimated under parametric model assumptions on the continuous diagnostic marker.


ROC.plot(ydat, xdat, distr = "exp", method = "RPstar", mc = 1)



data vector of the diagnostic marker measurements on the sample of non-diseased individuals (from Y).


data vector of the diagnostic marker measurements on the sample of diseased individuals (from X).


character string specifying the type of distribution assumed for Y and X. Possible choices for distr are "exp" (default) for the one-parameter exponential, "norm_EV" and "norm_DV" for the Gaussian distribution with, respectively, equal or unequal variances assumed for the two random variables.


character string specifying the methodological approach used for estimating the probability R, which is here interpreted as the area under the ROC curve (AUC). The argument method can be set equal to "Wald", "RP" or RPstar" (default), according as inference is based on the Wald statistic, the signed log-likelihood ratio statistic (directed likelihhod, r_p) or the modified signed log-likelihood ratio statistic (modified directed likelihood, r_p^*), respectively. For estimating the ROC curve parametrically, methods "Wald" and "RP" are equivalent and supply maximum likelihood estimation (MLE), whereas, by using method "RPstar", estimate of the ROC curve is based on the modified signed log-likelihood ratio statistic (r_p^*). See rpstar for details on this statistic.


a numeric value indicating single or multiple plots in the same figure. In case mc is equal to 1 (default), only the method specified in method is applied and the corresponding estimated ROC curve is plotted. If mc is different from 1, both MLE and r_p^*-based methods are applied, and two differently estimated ROC curves are plotted.


If mc is different from 1, method does not need to be specified.


Plot of ROC curves


The two independent random variables Y and X with given distribution distr are measurements of the diagnostic marker on the diseased and non-diseased subjects, respectively.
In "Wald" method, or equivalently "RP" method, MLEs for parameters of the Y and X distributions are computed and then used to estimate specificity and sensitivity. These measures are evaluated as P(Y<t) and P(X>t), respectively.
In "RPstar" method, parameters of the Y and X distributions are estimated from the r_p^*-based estimate of the AUC.


Giuliana Cortese


Cortese G., Ventura L. (2013). Accurate higher-order likelihood inference on P(Y<X). Computational Statistics, 28:1035-1059.

See Also



	# data from the non-diseased population
	Y <- rnorm(15, mean=5, sd=1)       
	# data from the diseased population
	X <- rnorm(10, mean=7, sd=1.5)      
 	ROC.plot(Y, X, "norm_DV", method = "RP", mc = 2)

ProbYX documentation built on May 30, 2017, 8:12 a.m.

Search within the ProbYX package
Search all R packages, documentation and source code