inst/simulations/simu_fwerPowerFdrPower.R
In mshasan/OPWpaper1: Reproduce figures and simulations of OPWeight
# part - I (continuous: for different effects)
#===============================================================================
library(snow)
library(mvnfast) # fast generate multi variate normal
library("IHW") # independent hypotheis weight
library(tibble) # tibble data frame
library(OPWeight)
library(OPWpaper)
cl <- makeCluster(10, type = "MPI") # start zcluster
clusterExport(cl,"rmvn")
clusterExport(cl,"ihw")
clusterExport(cl,"adj_pvalues")
clusterExport(cl,"rejections")
clusterExport(cl,"tibble")
clusterExport(cl, "weightSum_by_c")
clusterExport(cl, "roeder_wasserman_weight")
clusterExport(cl, "runif_by_mean")
clusterExport(cl, "test_by_block")
clusterExport(cl, "fwerPowerFdrPower")
filterEffectVec <- c(seq(0, 1, .2), 2, 3, 5, 8)
datWeight_cont <- read.csv("weight_byEffect_cont_m10000.csv", h = TRUE)
simuVal = 1000 # should use at least 1000
clusterExport(cl, "datWeight_cont")
clusterExport(cl, "filterEffectVec")
clusterExport(cl, "simuVal")
# FwerPowerFdrPower1e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
FwerPowerFdrPower1f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
FwerPowerFdrPower1h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
FwerPowerFdrPower2e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower2f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower2g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower2h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower2i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower3e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
FwerPowerFdrPower3f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
FwerPowerFdrPower3h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
FwerPowerFdrPower4e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
# FwerPowerFdrPower4f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
# FwerPowerFdrPower4g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
# FwerPowerFdrPower4h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
# FwerPowerFdrPower4i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower5e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
# FwerPowerFdrPower5f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
# FwerPowerFdrPower5g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
# FwerPowerFdrPower5h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
# FwerPowerFdrPower5i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
stopCluster(cl)
#-------------------------------------------------------------
# save the workspace to the file .RData in the cwd
save.image("simu_fwerPowerFdrPower_cont.RData")
#===============================================================================
#===============================================================================
# part - II (continuous: for different coefficient of variances)
library(snow)
library(mvnfast) # fast generate multi variate normal
library("IHW") # independent hypotheis weight
library(tibble) # tibble data frame
library(OPWeight)
library(OPWpaper)
cl <- makeCluster(10, type = "MPI") # start zcluster
clusterExport(cl,"rmvn")
clusterExport(cl,"ihw")
clusterExport(cl,"adj_pvalues")
clusterExport(cl,"rejections")
clusterExport(cl,"tibble")
clusterExport(cl, "weightSum_by_c")
clusterExport(cl, "roeder_wasserman_weight")
clusterExport(cl, "runif_by_mean")
clusterExport(cl, "test_by_block")
clusterExport(cl, "fwerPowerFdrPower")
filterEffectVec <- c(seq(0, 1, .2), 2, 3, 5, 8)
datWeight_cont <- read.csv("weight_byEffect_cont_m10000.csv", h = TRUE)
simuVal = 1000 # should use at least 1000
clusterExport(cl, "datWeight")
clusterExport(cl, "filterEffectVec")
clusterExport(cl, "simuVal")
FwerPowerFdrPower2a2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 1, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2a3 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 3, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2a10<- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 10,effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower4a2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 1, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4a3 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 3, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4a10<- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 10,effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
stopCluster(cl)
#-------------------------------------------------------------
# save the workspace to the file .RData in the cwd
save.image("simu_fwerPowerFdrPower_missVar_cont.RData")
#===============================================================================
#------------end: continuous section-----------------------
#===============================================================================
# part - III (Binary: for different effects)
library(snow)
library(mvnfast) # fast generate multi variate normal
library("IHW") # independent hypotheis weight
library(tibble) # tibble data frame
library(OPWeight)
library(OPWpaper)
cl <- makeCluster(10, type = "MPI") # start zcluster
clusterExport(cl,"rmvn")
clusterExport(cl,"ihw")
clusterExport(cl,"adj_pvalues")
clusterExport(cl,"rejections")
clusterExport(cl,"tibble")
clusterExport(cl, "weightSum_by_c")
clusterExport(cl, "roeder_wasserman_weight")
clusterExport(cl, "runif_by_mean")
clusterExport(cl, "test_by_block")
clusterExport(cl, "fwerPowerFdrPower")
filterEffectVec <- c(seq(0, 1, .2), 2, 3, 5, 8)
datWeight_bin <- read.csv("weight_byEffect_bin_m10000.csv", h = TRUE)
simuVal = 1000 # should use at least 1000
clusterExport(cl, "datWeight_bin")
clusterExport(cl, "filterEffectVec")
clusterExport(cl, "simuVal")
# FwerPowerFdrPower1e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
FwerPowerFdrPower1f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
FwerPowerFdrPower1h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
FwerPowerFdrPower2e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower2f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower2g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower2h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower2i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower3e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
FwerPowerFdrPower3f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
FwerPowerFdrPower3h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
FwerPowerFdrPower4e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
# FwerPowerFdrPower4f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
# FwerPowerFdrPower4g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
# FwerPowerFdrPower4h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
# FwerPowerFdrPower4i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower5e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
# FwerPowerFdrPower5f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
# FwerPowerFdrPower5g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
# FwerPowerFdrPower5h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
# FwerPowerFdrPower5i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
stopCluster(cl)
#-------------------------------------------------------------
# save the workspace to the file .RData in the cwd
save.image("simu_fwerPowerFdrPower_bin.RData")
#===============================================================================
#===============================================================================
# part - IV (binary: for different coefficient of variances)
library(snow)
library(mvnfast) # fast generate multi variate normal
library("IHW") # independent hypotheis weight
library(tibble) # tibble data frame
library(OPWeight)
library(OPWpaper)
cl <- makeCluster(10, type = "MPI") # start zcluster
clusterExport(cl,"rmvn")
clusterExport(cl,"ihw")
clusterExport(cl,"adj_pvalues")
clusterExport(cl,"rejections")
clusterExport(cl,"tibble")
clusterExport(cl, "weightSum_by_c")
clusterExport(cl, "roeder_wasserman_weight")
clusterExport(cl, "runif_by_mean")
clusterExport(cl, "test_by_block")
clusterExport(cl, "fwerPowerFdrPower")
filterEffectVec <- c(seq(0, 1, .2), 2, 3, 5, 8)
datWeight <- read.csv("weight_byEffect_bin_m10000.csv", h = TRUE)
simuVal = 1000 # should use at least 1000
clusterExport(cl, "datWeight")
clusterExport(cl, "filterEffectVec")
clusterExport(cl, "simuVal")
FwerPowerFdrPower2a2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 1, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2a3 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 3, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2a10<- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 10,effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower4a2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 1, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4a3 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 3, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4a10<- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 10,effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
stopCluster(cl)
#-------------------------------------------------------------
# save the workspace to the file .RData in the cwd
save.image("simu_fwerPowerFdrPower_missVar_bin.RData")
#===============================================================================
mshasan/OPWpaper1 documentation built on Feb. 22, 2021, 10:22 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
# part - I (continuous: for different effects)
#===============================================================================
library(snow)
library(mvnfast) # fast generate multi variate normal
library("IHW") # independent hypotheis weight
library(tibble) # tibble data frame
library(OPWeight)
library(OPWpaper)
cl <- makeCluster(10, type = "MPI") # start zcluster
clusterExport(cl,"rmvn")
clusterExport(cl,"ihw")
clusterExport(cl,"adj_pvalues")
clusterExport(cl,"rejections")
clusterExport(cl,"tibble")
clusterExport(cl, "weightSum_by_c")
clusterExport(cl, "roeder_wasserman_weight")
clusterExport(cl, "runif_by_mean")
clusterExport(cl, "test_by_block")
clusterExport(cl, "fwerPowerFdrPower")
filterEffectVec <- c(seq(0, 1, .2), 2, 3, 5, 8)
datWeight_cont <- read.csv("weight_byEffect_cont_m10000.csv", h = TRUE)
simuVal = 1000 # should use at least 1000
clusterExport(cl, "datWeight_cont")
clusterExport(cl, "filterEffectVec")
clusterExport(cl, "simuVal")
# FwerPowerFdrPower1e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
FwerPowerFdrPower1f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
FwerPowerFdrPower1h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
# FwerPowerFdrPower1i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,1:10])
FwerPowerFdrPower2e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower2f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower2g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower2h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower2i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
# FwerPowerFdrPower3e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
FwerPowerFdrPower3f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
FwerPowerFdrPower3h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
# FwerPowerFdrPower3i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,21:30])
FwerPowerFdrPower4e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
# FwerPowerFdrPower4f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
# FwerPowerFdrPower4g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
# FwerPowerFdrPower4h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
# FwerPowerFdrPower4i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower5e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr= 0, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr= 0, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.3, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
# FwerPowerFdrPower5f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.3, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.5, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
# FwerPowerFdrPower5g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.5, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.7, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
# FwerPowerFdrPower5h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.7, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
FwerPowerFdrPower5i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.9, cv= 0, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
# FwerPowerFdrPower5i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.9, cv=.5, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,41:50])
stopCluster(cl)
#-------------------------------------------------------------
# save the workspace to the file .RData in the cwd
save.image("simu_fwerPowerFdrPower_cont.RData")
#===============================================================================
#===============================================================================
# part - II (continuous: for different coefficient of variances)
library(snow)
library(mvnfast) # fast generate multi variate normal
library("IHW") # independent hypotheis weight
library(tibble) # tibble data frame
library(OPWeight)
library(OPWpaper)
cl <- makeCluster(10, type = "MPI") # start zcluster
clusterExport(cl,"rmvn")
clusterExport(cl,"ihw")
clusterExport(cl,"adj_pvalues")
clusterExport(cl,"rejections")
clusterExport(cl,"tibble")
clusterExport(cl, "weightSum_by_c")
clusterExport(cl, "roeder_wasserman_weight")
clusterExport(cl, "runif_by_mean")
clusterExport(cl, "test_by_block")
clusterExport(cl, "fwerPowerFdrPower")
filterEffectVec <- c(seq(0, 1, .2), 2, 3, 5, 8)
datWeight_cont <- read.csv("weight_byEffect_cont_m10000.csv", h = TRUE)
simuVal = 1000 # should use at least 1000
clusterExport(cl, "datWeight")
clusterExport(cl, "filterEffectVec")
clusterExport(cl, "simuVal")
FwerPowerFdrPower2a2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 1, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2a3 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 3, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower2a10<- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 10,effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,11:20])
FwerPowerFdrPower4a2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 1, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4a3 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 3, effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
FwerPowerFdrPower4a10<- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 10,effectType = "continuous", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_cont[,31:40])
stopCluster(cl)
#-------------------------------------------------------------
# save the workspace to the file .RData in the cwd
save.image("simu_fwerPowerFdrPower_missVar_cont.RData")
#===============================================================================
#------------end: continuous section-----------------------
#===============================================================================
# part - III (Binary: for different effects)
library(snow)
library(mvnfast) # fast generate multi variate normal
library("IHW") # independent hypotheis weight
library(tibble) # tibble data frame
library(OPWeight)
library(OPWpaper)
cl <- makeCluster(10, type = "MPI") # start zcluster
clusterExport(cl,"rmvn")
clusterExport(cl,"ihw")
clusterExport(cl,"adj_pvalues")
clusterExport(cl,"rejections")
clusterExport(cl,"tibble")
clusterExport(cl, "weightSum_by_c")
clusterExport(cl, "roeder_wasserman_weight")
clusterExport(cl, "runif_by_mean")
clusterExport(cl, "test_by_block")
clusterExport(cl, "fwerPowerFdrPower")
filterEffectVec <- c(seq(0, 1, .2), 2, 3, 5, 8)
datWeight_bin <- read.csv("weight_byEffect_bin_m10000.csv", h = TRUE)
simuVal = 1000 # should use at least 1000
clusterExport(cl, "datWeight_bin")
clusterExport(cl, "filterEffectVec")
clusterExport(cl, "simuVal")
# FwerPowerFdrPower1e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
FwerPowerFdrPower1f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
FwerPowerFdrPower1h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
# FwerPowerFdrPower1i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.2, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,1:10])
FwerPowerFdrPower2e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower2f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower2g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower2h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower2i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
# FwerPowerFdrPower3e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
FwerPowerFdrPower3f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
FwerPowerFdrPower3h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
# FwerPowerFdrPower3i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.75, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,21:30])
FwerPowerFdrPower4e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
# FwerPowerFdrPower4f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
# FwerPowerFdrPower4g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
# FwerPowerFdrPower4h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
# FwerPowerFdrPower4i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower5e1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr= 0, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5e2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr= 0, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5f1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.3, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
# FwerPowerFdrPower5f2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.3, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5g1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.5, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
# FwerPowerFdrPower5g2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.5, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5h1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.7, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
# FwerPowerFdrPower5h2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.7, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
FwerPowerFdrPower5i1 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.9, cv= 0, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
# FwerPowerFdrPower5i2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.99, corr=.9, cv=.5, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,41:50])
stopCluster(cl)
#-------------------------------------------------------------
# save the workspace to the file .RData in the cwd
save.image("simu_fwerPowerFdrPower_bin.RData")
#===============================================================================
#===============================================================================
# part - IV (binary: for different coefficient of variances)
library(snow)
library(mvnfast) # fast generate multi variate normal
library("IHW") # independent hypotheis weight
library(tibble) # tibble data frame
library(OPWeight)
library(OPWpaper)
cl <- makeCluster(10, type = "MPI") # start zcluster
clusterExport(cl,"rmvn")
clusterExport(cl,"ihw")
clusterExport(cl,"adj_pvalues")
clusterExport(cl,"rejections")
clusterExport(cl,"tibble")
clusterExport(cl, "weightSum_by_c")
clusterExport(cl, "roeder_wasserman_weight")
clusterExport(cl, "runif_by_mean")
clusterExport(cl, "test_by_block")
clusterExport(cl, "fwerPowerFdrPower")
filterEffectVec <- c(seq(0, 1, .2), 2, 3, 5, 8)
datWeight <- read.csv("weight_byEffect_bin_m10000.csv", h = TRUE)
simuVal = 1000 # should use at least 1000
clusterExport(cl, "datWeight")
clusterExport(cl, "filterEffectVec")
clusterExport(cl, "simuVal")
FwerPowerFdrPower2a2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 1, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2a3 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 3, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower2a10<- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.5, cv= 10,effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,11:20])
FwerPowerFdrPower4a2 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 1, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4a3 <- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 3, effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
FwerPowerFdrPower4a10<- parSapply(cl, 1:length(filterEffectVec), fwerPowerFdrPower, simu=simuVal, null=.9, cv= 10,effectType = "binary", filterEffectVec=filterEffectVec, datWeightByNull=datWeight_bin[,31:40])
stopCluster(cl)
#-------------------------------------------------------------
# save the workspace to the file .RData in the cwd
save.image("simu_fwerPowerFdrPower_missVar_bin.RData")
#===============================================================================
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.