lfdrPlot.ash | R Documentation |
Make an lFDR plot for an ash object
lfdrPlot.ash(fit, n = 99, ...)
fit |
Output from 'ash()“ |
n |
Number of breaks for histogram |
... |
Further arguments to |
library(ashr)
z <- c(rnorm(900), rnorm(200, sd=2))
fit <- ash(z, rep(1, length(z)), optmethod='cxxMixSquarem')
lfdrPlot.ash(fit)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.