lines.thres2 | R Documentation |
The function includes vertical lines for the threshold and confidence interval in a plot created with plot.thres2()
.
## S3 method for class 'thres2'
lines(x, ci = TRUE, which.boot = c("norm", "perc"),
col = 1, lty = c(1, 2), lwd = 1, ...)
x |
an object of class |
ci |
should the confidence interval be plotted? Default, |
which.boot |
in case |
col |
color for the threshold and its corresponding confidence interval. Default, 1. |
lty |
a 2-dimensional vector containing:
Default, |
lwd |
line width for the threshold and its corresponding confidence interval. Default, 1. |
... |
further arguments to be passed to |
With a plot.thres2
open, this function adds lines for the required threshold.
Skaltsa K, Jover L, Carrasco JL. (2010). Estimation of the diagnostic threshold accounting for decision costs and sampling uncertainty. Biometrical Journal 52(5):676-697.
thres2
, plot.thres2
n1 <- 100
n2 <- 100
set.seed(1234)
par1.1 <- 0
par1.2 <- 1
par2.1 <- 2
par2.2 <- 1
rho <- 0.2
k1 <- rnorm(n1, par1.1, par1.2) # non-diseased
k2 <- rnorm(n2, par2.1, par2.2) # diseased
thres <- thres2(k1, k2, rho, method="eq", ci.method="d")
plot(thres, col=c(1, 2, 4), lwd=c(2, 2, 1), leg.pos="topright")
thresunequal <- thres2(k1, k2, rho, method="unequal", ci=FALSE)
lines(thresunequal, col=3) # almost the same; no confidence interval
## Not run:
thresboot <- thres2(k1, k2, rho, method="emp", ci.method="b")
lines(thresboot, col=5, which.boot="norm")
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.