plotPvals {modnets}

Plot the ECDF of p-values from resampling

Description

Plots the empirical cumulative distribution function of the p-values related to iterated resampling via bootstrapping or multi-sample splitting.

Usage

plotPvals(x, outcome = 1, predictor = 1, title = TRUE, alpha = 0.05)

Arguments

x

Output from resample, given that sampMethod = "bootstrap" or sampMethod = "split".

outcome

Character string or numeric value (in terms of columns in the dataset) to indicate which outcome to plot the p-value distribution for.

predictor

Character string or numeric value (in terms of columns in the dataset) to indicate which predictor to plot the p-value distribution for.

title

If TRUE, then a default title will be given according to the outcome and predictor that are specified. If FALSE, then no title will be plotted. A custom title may also be supplied by the user.

alpha

The false discovery rate. Defaults to .05

Details

See Meinshausen, Meier, & Buhlmann (2009) for details.

Value

Returns a plot based on the relationship between a particular outcome and predictor.

References

Meinshausen, N., Meier, L., & Buhlmann, P. (2009). P-values for high-dimensional regression. Journal of the American Statistical Association. 104, 1671-1681.

See Also

resample

Examples

x <- resample(ggmDat, sampMethod = "bootstrap")
plot(x, what = 'pvals')
plot(x, 'pvals', outcome = 'V2', predictor = 'V1')