FDRCurve {rvalues} | R Documentation |
FDR Curve
Description
Estimates the expected proportion of misclassified units when using a given r-value threshold. If plot=TRUE, the curve is plotted before the estimated function is returned.
Usage
FDRCurve(object, q, threshold = 1, plot = TRUE, xlim, ylim, xlab, ylab, main, ...)
Arguments
object |
An object of class "rvals" |
q |
A value in between 0 and 1; the desired level of FDR control. |
threshold |
The r-value threshold. |
plot |
logical; if TRUE, the estimated FDR curve is plotted. |
xlim , ylim |
x and y - axis limits for the plot |
xlab , ylab |
x and y - axis labels |
main |
the title of the plot |
... |
additional arguments to |
Details
Consider parameters of interest with an effect of size of interest
. That is, a unit is taken to be "null" if
and
taken to be "non-null" if
.
For r-values and a procedure which "rejects" units
satisfying
, the FDR is defined to be
FDRCurve
estimates
for values of
across (0,1) and plots (if
plot=TRUE
)
the resulting curve.
Value
A list with the following two components
fdrcurve |
A function which returns the estimated FDR for each r-value threshold. |
Rval.cut |
The largest r-value cutoff which still gives an estimated FDR less than q. |
Author(s)
Nicholas Henderson and Michael Newton
See Also
Examples
n <- 500
theta <- rnorm(n)
ses <- sqrt(rgamma(n,shape=1,scale=1))
XX <- theta + ses*rnorm(n)
dd <- cbind(XX,ses)
rvs <- rvalues(dd, family = gaussian)
FDRCurve(rvs, q = .1, threshold = .3, cex.main = 1.5)