inst/BNR/averagePower/00.setup.R

In pneuvial/sanssouci: Post Hoc Multiple Testing Inference

library("sansSouci")

sname0 <- "averagePower"
ptag <- sprintf("sansSouci_%s", packageVersion("sansSouci"))
sname <- sprintf("%s,%s", sname0, ptag)

deps <- c(0, 0.2, 0.4) ## correlation coefficient
pi0s <- c(0.8, 0.9, 0.99, 0.999)
alphas <- c(0.01, 0.02, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5)
sf <- c(1, 2, 4)

## calibrate SNR according to pi0 (default scaling factor of 1 places us on the estimation boundary)
getSNR <- function(pi0, sf=1) {
    sqrt(2*sf*log(1/(1-pi0)))
}

m <- 1e3
B <- 1e3
nbSimu <- 1e3

pname <- sprintf("m=%s,B=%s,nbSimu=%s", m, B, nbSimu)
resPath <- "resData"
path <- file.path(resPath, sname, pname)
path <- R.utils::Arguments$getWritablePath(path)

configs <- expand.grid(pi0=pi0s, dep=deps, sf=sf)
configs