Make an lFDR plot for an ash object

Make an lFDR plot for an ash object

Usage

lfdrPlot.ash(fit, n = 99, ...)

Arguments

fit

Output from `ash()“

n

Number of breaks for histogram

...

Further arguments to Hist()

Examples

library(ashr)
z <- c(rnorm(900), rnorm(200, sd=2))
fit <- ash(z, rep(1, length(z)), optmethod='cxxMixSquarem')
#> Warning: subscript out of bounds (index 0 >= vector size 0)
lfdrPlot.ash(fit)