R/testScript.R

In blebedenko/thescript2: Boris's thesis scripts

# distribution <- "normal"
# n <- 250
# ME_var_ratio <- 1
# var_eq <- TRUE
# AUC <- 0.8
# theta_squared <- 1
# R <- 1000
# US_ratio <-  c(0.5,0.25)
# 
# 
# AUC.values <- c(0.7,0.8,0.9)
# theta_squared.values <- c(0.5,1,2)
# var_eq.values <- c(TRUE,FALSE)
# ME_var_ratio.values <- c(0.5,1,2)
# distribution.values <- distribution
# 
# MC <- 2
# 
# varList <- simsalapar::varlist(
#   n.sim = list(type = "N", expr = quote(n.sim),  value = MC),
#   theta_squared = list(type = "grid",expr = quote(theta^2), value =theta_squared.values),
#   AUC = list(type = "grid",expr = quote(A), value = AUC.values),
#   var_eq = list(type = "grid",expr = quote("VX=VY"), value = var_eq.values),
#   ME_var_ratio = list(type = "grid",expr = quote(Veps/Veta), value = ME_var_ratio.values),
#   distribution = list(type = "grid", expr = quote(F), value = distribution.values),
#   US_ratio = list(type = "frozen", expr = quote(q),value = c(0.5,0.25)),
#   n = list(type = "frozen", expr = quote(n), value = n)
# 
# )
# cl <- makeCluster(detectCores()-1)
# tictoc::tic()
# RES <- simsalapar::doForeach(vList = varList,
#                              doOne = AtomicSim.US,
#                              monitor = interactive(),
#                              extraPkgs = c("QPdecon","thescript","boot","coxed"),
#                              cluster = cl)
# 
# tictoc::toc()
# 
# array2df(getArray(RES))
# devtools::install_github("blebedenko/thescript")

# library(simsalapar)
# library(QPdecon)