plotNull {reconsi}R Documentation

Plot the obtained null distribution along with a histogram of observed test statistics

Description

Plot the obtained null distribution along with a histogram of observed test statistics

Usage

plotNull(fit, lowColor = "yellow", highColor = "blue", dens = TRUE,
  idDA = NULL, nResampleCurves = length(fit$weights), hSize = 0.5)

Arguments

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

Value

a ggplot2 plot object

Examples

 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)

[Package reconsi version 1.0.0 Index]