f_hat_1: Estimate of F^1, The Distribution of p-values Under The...
In ntyet/sfdr: False Discovery Rate Estimation in Two Independent Experiments

Description Usage Arguments Value Author(s) Examples

This function estimates the cdf of p-values under the alternative H_1.

1	f_hat_1(p, t, pi_0)

`p`	the vector of p-values.
`t`	real number between 0 and 1.
`pi_0`	the proportion of true nulls.

CDF of p-values under the alternative at the threshold t.

Yet Nguyen ntyet@iastate.edu.

n <- 10000
n1 <- 250
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)
pi_0 <- pi0(p1, lambda0)