plotNull: Plot the obtained null distribution along with a histogram of...

Description Usage Arguments Value Examples

View source: R/F_plotNull.R

Description

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

Usage

1
2
3
4
5
6
7
8
plotNull(
  fit,
  lowColor = "yellow",
  highColor = "blue",
  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

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

1
2
3
4
5
6
7
8
9
 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)

reconsi documentation built on Nov. 8, 2020, 5:04 p.m.

Related to plotNull in reconsi...