Sensitivity and specificity, (positive and negative) predictive values and (positive and negative) diagnostic likelihood ratios of a single binary diagnostic test.
acc.1test(tab, alpha, testname, method.ci, ...)A list of class acc.1test:
A contingency table (matrix) of test results; the same tab that is supplied as an argument.
| Diseased | Non-diseased | Total | |
| Test pos. | ... | ... | ... |
| Test neg. | ... | ... | ... |
| Total | ... | ... | ... |
A numeric vector containing the estimated sensitivity (est), its standard error (se), lower confidence limit (lcl) and upper confidence limit (ucl).
A numeric vector containing the estimated specificity (est), its standard error (se), lower confidence limit (lcl) and upper confidence limit (ucl).
A numeric vector containing the estimated positive predictive value (est), its standard error (se), lower confidence limit (lcl) and upper confidence limit (ucl).
A numeric vector containing the estimated negative predictive value (est), its standard error (se), lower confidence limit (lcl) and upper confidence limit (ucl).
A numeric vector containing the estimated positive diagnostic likelihood ratio (est), the standard error of the logarithm of the positive diagnostic likelihood ratio (se.ln), the lower confidence limit (lcl) and the upper confidence limit (ucl).
A numeric vector containing the estimated negative diagnostic likelihood ratio (est), the standard error of the logarithm of the negative diagnostic likelihood ratio (se.ln), the lower confidence limit (lcl) and the upper confidence limit (ucl).
The significance level alpha used to compute 100(1-alpha)%-confidence intervals, the default is 0.05.
A character string containing the name of the diagnostic test.
A character string describing the method used to compute the confidence intervals for sensitivity, specificity, and predictive values.
An object of class tab.1test.
Significance level alpha for 100(1-alpha)%-confidence intervals, the default is 0.05.
A character string containing the name of the diagnostic test.
A character string with the name of the function to compute
the confidence intervals for sensitivity, specificity, and predictive values.
The default is ''waldci'' for a Wald asymptotic normal-based confidence
intervals. See 'Details'.
Additional arguments (usually not required).
The calculation of accuracy measures follows standard methodology, e.g. described in Pepe (2003) or Zhou et al. (2011).
The confidence intervals for sensitivity, specificity, and predictive values are
computed using the methodology implemented in the function passed to the
argument method.ci.
Available options are:
waldci - Wald asymptotic normal-based confidence interval; the
default,
logitci - asymptotic normal-based confidence interval on the
logit scale and then back-transformed,
exactci - Clopper-Pearson exact confidence interval,
add4cici - Agresti-Coull add-4 confidence interval,
addz2ci - Agresti-Coull add-z^2/2 confidence interval,
blakerci - Blaker exact confidence interval,
scoreci - Wilson score confidence interval,
midPci - mid-P confidence interval.
Options (3) to (8) are based on the homonymous functions from the {PropCIs}
package. Please see the respective package documentation for more details.
Confidence intervals for diagnostic likelihood ratios are computed according to Simel et al. (1991).
Pepe, M. (2003). The statistical evaluation of medical tests for classification and prediction. Oxford Statistical Science Series. Oxford University Press, 1st edition.
Simel, D.L., Samsa, G.P., Matchar, D.B. (1991). Likelihood ratios with confidence: sample size estimation for diagnostic test studies. J Clin Epidemiol, 44(8):763-70.
Zhou, X., Obuchowski, N., and McClish, D. (2011). Statistical Methods in Diagnostic Medicine. Wiley Series in Probability and Statistics. John Wiley & Sons, Hoboken, New Jersey, 2nd edition.
tab.1test,
print.acc.1test,
acc.paired.
data(Paired1) # Hypothetical study data
a1 <- tab.1test(d=d, y=y1, data=Paired1)
a2 <- acc.1test(a1)
print(a2)
a3 <- acc.1test(a1, method="exactci", conf.level=0.99)
print(a3)
Run the code above in your browser using DataLab