pi0: Estimate The Proportion of Nulls
In ntyet/sfdr: False Discovery Rate Estimation in Two Independent Experiments

Description Usage Arguments Value Examples

This function estimates the proportion of true nulls using Storey's FDR method.

1	pi0(p, lambda0)

`p`	a numerical vector of p-value.
`lambda0`	the threshold lambda where all pvalue greater than lambda are null. This is used to calculate estimator of the proportion of true nulls.

estimate of the proportion of true null

n <- 1000
n1 <- 25
mu1 <- rep(0, n)
mu1[1:n1] <- rnorm(n = n1, mean = 6,sd = 1)
t1 <- plyr::laply(1:n, function(i)rt(n = 1, df = 4, ncp = mu1[i]))
p1 <- plyr::laply(1:n, function(i)2*pt(q = abs(t1[i]), df = 4, ncp = 0, lower.tail = FALSE))
lambda0 <- estimate.lambda(p1, B = 20)
pi0(p1, lambda0)