Learn R Programming
mixtools (version 1.0.1)

plotFDR: Plot False Discovery Rate (FDR) estimates from output by EM-like strategies

Description

Plot FDR$(p_i)$ estimates against index of sorted p-values from, e.g., normalmixEM or the semiparametric mixture model posterior probabilities output by spEMsymlocN01, or any EM-algorithm like normalmixEM which returns posterior probabilities. The function can simultaneously plot FDR estimates from two strategies for comparison. Plot of the true FDR can be added if complete data are available (typically in simulation studies).

Usage

plotFDR (post1, post2 = NULL, lg1 = "FDR 1", lg2 = NULL, title = NULL, 
		compH0 = 1, alpha = 0.1, complete.data = NULL, pctfdr = 0.3)

Arguments

post1
The posterior probabilities from objects such as the output from spEMsymlocN01.
post2
A second object like post1 if comparison is desired.
lg1
Text describing the FDR estimate in post1.
lg2
Text describing the FDR estimate in post2 if provided.
title
Plot title, a default is provided if NULL.
compH0
The component indicator associated to the null hypothesis H0, normally 1 since it is defined in this way in spEMsymlocN01, but in case of label switching in other algorithms it can be set to
alpha
The target FDR level; the index at which the FDR estimate crosses the horizontal line for level alpha gives the maximum number of cases to reject.
complete.data
An array with $n$ lines and 2 columns, with the component indicator in column 1 and the p-values in column 2, sorted by p-values.
pctfdr
The level up to which the FDR is plotted, i.e. the scale of the vertical axis.

Value

  • A plot of one or two FDR estimates, with the true FDR if available

References

  • Chauveau, D., Saby, N., Orton, T. G., Lemercier B., Walter, C. and Arrouys, D. Large-scale simultaneous hypothesis testing in soil monitoring: A semi-parametric mixture approach, preprint (2013).

See Also

spEMsymlocN01