plotNull {reconsi} | R Documentation |
Plot the obtained null distribution along with a histogram of observed test statistics
plotNull(fit, lowColor = "yellow", highColor = "blue", dens = TRUE, idDA = NULL, nResampleCurves = length(fit$weights), hSize = 0.5)
fit |
an object returned by the reconsi() (or testDAA()) function |
lowColor, highColor |
The low and high ends of the colour scale |
dens |
a boolean, should fdr and Fdr be plotted? |
idDA |
indices of known null taxa |
nResampleCurves |
The number of resampling null distributions to plot |
hSize |
A double, the size of the line of the collapsed null estimate |
a ggplot2 plot object
p = 175; n = 50; B = 1e2 #Low number of resamples keeps computation time down x = rep(c(0,1), each = n/2) mat = cbind( matrix(rnorm(n*p/10, mean = 5+x),n,p/10), #DA matrix(rnorm(n*p*9/10, mean = 5),n,p*9/10) #Non DA ) fdrRes = reconsi(mat, x, B = B) plotNull(fdrRes)