View source: R/EFDR_functions.R
fdrpower | R Documentation |
Returns the power of the multiple hypothesis test, by finding the proportion of the correctly rejected null hypotheses.
fdrpower(reject.true, reject)
reject.true |
indices of the true alternative hypotheses |
reject |
indices of the rejected null hypotheses |
Single value (proportion)
Shen, X., Huang, H.-C., and Cressie, N. 'Nonparametric hypothesis testing for a spatial signal.' Journal of the American Statistical Association 97.460 (2002): 1122-1140.
set.seed(1)
wf = "la8"
J = 3
n = 64
h = 0.5
Z <- test_image(h = h, r = 14, n1 = n)$z
sig <- wav_th(Z, wf=wf, J=J, th = h)
Z <- Z + rnorm(n^2)*0.5
m1 <- test.bonferroni(Z, wf="la8",J=3, alpha = 0.05)
m2 <- test.fdr(Z, wf="la8",J=3, alpha = 0.05)
cat(paste0("Bonferroni power: ",fdrpower(sig,m1$reject_coeff)))
cat(paste0("FDR power: ",fdrpower(sig,m2$reject_coeff)))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.