foldid
table(foldid)
?cv.glmnet
?`LassoSIR-Package`
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000.R')
foldid
table(foldid)
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000_dat.R')
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000.R')
source('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000.R')
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000.R')
seassir_result
c(seassir_result, seasintra_result, seaspfc_result, lassosir_result, lasso_result)
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000.R')
c(seassir_result, seasintra_result, seaspfc_result, lassosir_result, lasso_result)
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000.R')
c(seassir_result, seasintra_result, seaspfc_result, lassosir_result, lasso_result)
a <- read.table("output/M1p1000-1")
a
a <- read.table("output/M2p1000-1")
a
a <- read.table("output/M3a_p1000-1")
a
a <- read.table("output/M3b_p1000-1")
a
a <- read.table("output/M4p1000-1")
a
a <- read.table("output/M1p1000-1")
a
a <- read.table("output/M2p1000-1")
a
ind1 <- grep("(dist)|(TPR)|(FPR)", names(a))
ind2 <- grep("(d_)", names(a))
ind3 <- grep("(time)", names(a))
m <- apply(a, 2, mean, na.rm=TRUE)
m1 <- sprintf("%.1f", m[ind1]*100)
m2 <- sprintf("%.2f", m[ind2])
m3 <- sprintf("%.3f", m[ind3])
m[ind1] <- m1
m[ind2] <- m2
m[ind3] <- m3
se <- apply(a, 2, function(x){sd(x, na.rm = TRUE)/sqrt(sum(!is.na(x)))})
se1 <- sprintf("%.1f", se[ind1]*100)
se2 <- sprintf("%.2f", se[ind2])
se3 <- sprintf("%.3f", se[ind3])
se[ind1] <- se1
se[ind2] <- se2
se[ind3] <- se3
colnames(a)
cat("SEAS-SIR &", paste0(m[1:4], " (", se[1:4], ") ", collapse = "&"), "\\\\", "\n",
"SEAS-Intra &", paste0(m[5:8], " (", se[5:8], ")", collapse = "&"), "\\\\", "\n",
"SEAS-PFC &", paste0(m[c(9:12)], " (", se[c(9:12)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso-SIR &", paste0(m[c(13:16)], " (", se[c(13:16)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso &", paste0(m[17:20], " (", se[17:20], ") ", collapse = "&"), "\\\\", "\n",
)
cat("SEAS-SIR &", paste0(m[1:4], " (", se[1:4], ") ", collapse = "&"), "\\\\", "\n",
"SEAS-Intra &", paste0(m[5:8], " (", se[5:8], ")", collapse = "&"), "\\\\", "\n",
"SEAS-PFC &", paste0(m[c(9:12)], " (", se[c(9:12)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso-SIR &", paste0(m[c(13:16)], " (", se[c(13:16)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso &", paste0(m[17:20], " (", se[17:20], ") ", collapse = "&"), "\\\\", "\n")
a <- read.table("output/M2p1000-1")
colMeans(a)
a <- read.table("output/M2p1000_test")
colMeans(a)
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000.R')
foldid
lambda
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000.R')
source('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation2p1000_dat.R')
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation2p1000.R')
output
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation2p1000.R')
output
a <- read.table("output/M2p1000_test")
a
colMeans(a)
a <- read.table("output/M2p1000-1")
colMeans(a)
?kmeans
kmeans
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000_dat.R')
count
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000_dat.R')
debugSource('~/Dropbox/shared/SEAS_R1/code/simulation_p1000/simulation1p1000.R')
beta[1:20]
a <- read.table("output/M2p1000_LassoSIR")
a
colMeans(a)
a <- read.table("output/M2p1000-1")
colMeans(a)
a <- read.table("output/M2p1000_test")
colMeans(a)
a <- read.table("output/M2p1000")
colMeans()
colMeans(a)
a <- read.table("output/M2p1000_1")
colMeans(a)
a <- read.table("output/M2p1000")
colMeans(a[1:33,])
a <- read.table("output/M2p1000_2")
colMeans(a[1:33,])
colMeans(a)
a <- read.table("output/M2p1000")
colMeans(a[34:67,])
colMeans(a[68:100,])
a <- read.table("output/M2p1000_3")
colMeans(a[68:100,])
colMeans(a)
a <- read.table("output/M3b_p1000_1")
a_tmp <- a
a <- read.table("output/M3b_p1000_2")
a_tmp <- cbind(a_tmp, a)
a_tmp <- rbind(a_tmp, a)
a <- read.table("output/M3b_p1000_3")
a_tmp <- rbind(a_tmp, a)
write.table(a_tmp, file = "output/M3b_p1000")
a <- read.table("output/M3b_p1000_1")
colMeans(a)
a <- read.table("output/M3b_p1000_2")
colMeans(a)
a <- read.table("output/M3b_p1000_3")
colMeans(a)
a <- read.table("output/M3b_p1000")
colMeans(a)
a <- read.table("output/M1p1000")
ind1 <- grep("(dist)|(TPR)|(FPR)", names(a))
ind2 <- grep("(d_)", names(a))
ind3 <- grep("(time)", names(a))
m <- apply(a, 2, mean, na.rm=TRUE)
m1 <- sprintf("%.1f", m[ind1]*100)
m2 <- sprintf("%.2f", m[ind2])
m3 <- sprintf("%.3f", m[ind3])
m[ind1] <- m1
m[ind2] <- m2
m[ind3] <- m3
se <- apply(a, 2, function(x){sd(x, na.rm = TRUE)/sqrt(sum(!is.na(x)))})
se1 <- sprintf("%.1f", se[ind1]*100)
se2 <- sprintf("%.2f", se[ind2])
se3 <- sprintf("%.3f", se[ind3])
se[ind1] <- se1
se[ind2] <- se2
se[ind3] <- se3
cat("SEAS-SIR &", paste0(m[1:4], " (", se[1:4], ") ", collapse = "&"), "\\\\", "\n",
"SEAS-Intra &", paste0(m[5:8], " (", se[5:8], ")", collapse = "&"), "\\\\", "\n",
"SEAS-PFC &", paste0(m[c(9:12)], " (", se[c(9:12)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso &", paste0(m[c(13:16)], " (", se[c(13:16)], ") ", collapse = "&"), "\\\\", "\n")
a <- read.table("output/M2p1000")
cat("SEAS-SIR &", paste0(m[1:4], " (", se[1:4], ") ", collapse = "&"), "\\\\", "\n",
"SEAS-Intra &", paste0(m[5:8], " (", se[5:8], ")", collapse = "&"), "\\\\", "\n",
"SEAS-PFC &", paste0(m[c(9:12)], " (", se[c(9:12)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso &", paste0(m[c(13:16)], " (", se[c(13:16)], ") ", collapse = "&"), "\\\\", "\n")
ind1 <- grep("(dist)|(TPR)|(FPR)", names(a))
ind2 <- grep("(d_)", names(a))
ind3 <- grep("(time)", names(a))
m <- apply(a, 2, mean, na.rm=TRUE)
m1 <- sprintf("%.1f", m[ind1]*100)
m2 <- sprintf("%.2f", m[ind2])
m3 <- sprintf("%.3f", m[ind3])
m[ind1] <- m1
m[ind2] <- m2
m[ind3] <- m3
se <- apply(a, 2, function(x){sd(x, na.rm = TRUE)/sqrt(sum(!is.na(x)))})
se1 <- sprintf("%.1f", se[ind1]*100)
se2 <- sprintf("%.2f", se[ind2])
se3 <- sprintf("%.3f", se[ind3])
se[ind1] <- se1
se[ind2] <- se2
se[ind3] <- se3
cat("SEAS-SIR &", paste0(m[1:4], " (", se[1:4], ") ", collapse = "&"), "\\\\", "\n",
"SEAS-Intra &", paste0(m[5:8], " (", se[5:8], ")", collapse = "&"), "\\\\", "\n",
"SEAS-PFC &", paste0(m[c(9:12)], " (", se[c(9:12)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso &", paste0(m[c(13:16)], " (", se[c(13:16)], ") ", collapse = "&"), "\\\\", "\n")
a <- read.table("output/M3a_p1000")
ind1 <- grep("(dist)|(TPR)|(FPR)", names(a))
ind2 <- grep("(d_)", names(a))
ind3 <- grep("(time)", names(a))
m <- apply(a, 2, mean, na.rm=TRUE)
m1 <- sprintf("%.1f", m[ind1]*100)
m2 <- sprintf("%.2f", m[ind2])
m3 <- sprintf("%.3f", m[ind3])
m[ind1] <- m1
m[ind2] <- m2
m[ind3] <- m3
se <- apply(a, 2, function(x){sd(x, na.rm = TRUE)/sqrt(sum(!is.na(x)))})
se1 <- sprintf("%.1f", se[ind1]*100)
se2 <- sprintf("%.2f", se[ind2])
se3 <- sprintf("%.3f", se[ind3])
se[ind1] <- se1
se[ind2] <- se2
se[ind3] <- se3
cat("SEAS-SIR &", paste0(m[1:4], " (", se[1:4], ") ", collapse = "&"), "\\\\", "\n",
"SEAS-Intra &", paste0(m[5:8], " (", se[5:8], ")", collapse = "&"), "\\\\", "\n",
"SEAS-PFC &", paste0(m[c(9:12)], " (", se[c(9:12)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso &", paste0(m[c(13:16)], " (", se[c(13:16)], ") ", collapse = "&"), "\\\\", "\n")
a <- read.table("output/M3b_p1000")
ind1 <- grep("(dist)|(TPR)|(FPR)", names(a))
ind2 <- grep("(d_)", names(a))
ind3 <- grep("(time)", names(a))
m <- apply(a, 2, mean, na.rm=TRUE)
m1 <- sprintf("%.1f", m[ind1]*100)
m2 <- sprintf("%.2f", m[ind2])
m3 <- sprintf("%.3f", m[ind3])
m[ind1] <- m1
m[ind2] <- m2
m[ind3] <- m3
se <- apply(a, 2, function(x){sd(x, na.rm = TRUE)/sqrt(sum(!is.na(x)))})
se1 <- sprintf("%.1f", se[ind1]*100)
se2 <- sprintf("%.2f", se[ind2])
se3 <- sprintf("%.3f", se[ind3])
se[ind1] <- se1
se[ind2] <- se2
se[ind3] <- se3
cat("SEAS-SIR &", paste0(m[1:4], " (", se[1:4], ") ", collapse = "&"), "\\\\", "\n",
"SEAS-Intra &", paste0(m[5:8], " (", se[5:8], ")", collapse = "&"), "\\\\", "\n",
"SEAS-PFC &", paste0(m[c(9:12)], " (", se[c(9:12)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso &", paste0(m[c(13:16)], " (", se[c(13:16)], ") ", collapse = "&"), "\\\\", "\n")
a <- read.table("output/M4_p1000")
a <- read.table("output/M4p1000")
ind1 <- grep("(dist)|(TPR)|(FPR)", names(a))
ind2 <- grep("(d_)", names(a))
ind3 <- grep("(time)", names(a))
m <- apply(a, 2, mean, na.rm=TRUE)
m1 <- sprintf("%.1f", m[ind1]*100)
m2 <- sprintf("%.2f", m[ind2])
m3 <- sprintf("%.3f", m[ind3])
m[ind1] <- m1
m[ind2] <- m2
m[ind3] <- m3
se <- apply(a, 2, function(x){sd(x, na.rm = TRUE)/sqrt(sum(!is.na(x)))})
se1 <- sprintf("%.1f", se[ind1]*100)
se2 <- sprintf("%.2f", se[ind2])
se3 <- sprintf("%.3f", se[ind3])
se[ind1] <- se1
se[ind2] <- se2
se[ind3] <- se3
cat("SEAS-SIR &", paste0(m[1:4], " (", se[1:4], ") ", collapse = "&"), "\\\\", "\n",
"SEAS-Intra &", paste0(m[5:8], " (", se[5:8], ")", collapse = "&"), "\\\\", "\n",
"SEAS-PFC &", paste0(m[c(9:12)], " (", se[c(9:12)], ") ", collapse = "&"), "\\\\", "\n",
"Lasso &", paste0(m[c(13:16)], " (", se[c(13:16)], ") ", collapse = "&"), "\\\\", "\n")
a <- read.table("output/M1p1000")
colMeans(a)
a <- read.table("output/M1p1000-1")
colMeans(a)
a <- read.table("output/M1p1000_LassoSIR")
colMeans(a)
a <- read.table("output/M1p1000_LassoSIR-1")
colMeans(a)
a <- read.table("output/M2p1000_test")
colMeans(a)
a <- read.table("output/M2p1000")
colMeans(a)
a <- read.table("output/M3a_p1000")
colMeans(a)
a <- read.table("output/M3a_p1000-1")
colMeans(a)
a <- read.table("output/M3b_p1000")
colMeans(a)
a <- read.table("output/M3b_p1000_1")
a_tem <- a
a <- read.table("output/M3b_p1000_2")
a_tem <- rbind(a_tem, a)
a <- read.table("output/M3b_p1000_3")
a_tem <- rbind(a_tem, a)
colMeans(a_tem)
a <- read.table("output/M4")
a <- read.table("output/M4p1000")
colMeans(a_tem)
colMeans(a)
a <- read.table("output/M4p1000-1")
colMeans(a)
a <- read.table("output/M4p1000_LassoSIR")
colMeans(a)
a <- read.table("output/M4p1000_LassoSIR-1")
colMeans(a)
write.table(a_tem, file = "output/M3b_p1000")
colMeans(a)
colMeans(a_tem)
a <- read.table("output/M1_LassoSIR")
colMeans(a)
a <- read.table("output/M1p1000_LassoSIR")
colMeans(a)
debugSource('~/OneDrive - Florida State Students/ongoing_projects/SEAS_code/code_R1/simulation_p1000/simulation1_dat.R')
save(foldid, file = paste0("dat/sim1/foldid",i,'.rda'))
debugSource('~/OneDrive - Florida State Students/ongoing_projects/SEAS_code/code_R1/simulation_p1000/simulation1_dat.R')
writeMat(paste0("dat/sim1/foldid",i,".mat"), foldid = foldid)
a <- read.table("output/M1p1000")
colMeans(a)
a <- read.table("output/M1")
colMeans(a)
a <- read.table("output/M2p1000")
colMeans(a)
a <- read.table("output/M2")
colMeans(a)
a <- read.table("output/M3a_p1000")
colMeans(a)
a <- read.table("output/M3a")
colMeans(a)
a <- read.table("output/M3b_1")
a_tmp <- a
a <- read.table("output/M3b_2")
a_tmp <- rbind(a_tmp, a)
a <- read.table("output/M3b_3")
a_tmp <- rbind(a_tmp, a)
write.table(a_tmp, file = "output/M3b")
a <- read.table("output/M3b")
colMeans(a)
a <- read.table("output/M3b_p1000")
colMeans(a)
a <- read.table("output/M4")
colMeans(a)
a <- read.table("output/M4p1000")
colMeans(a)
debugSource('~/OneDrive - Florida State Students/ongoing_projects/SEAS_code/code_R1/simulation_p1000/simulation1_dat.R')
foldid
debugSource('~/OneDrive - Florida State Students/ongoing_projects/SEAS_code/code_R1/sliceH/simulation3a_error.R')
beta[ind,2] <- c(runif(3, 0.3, 0.6), runif(3,-0.6,-0.3))
beta[1:11]
writeMat("beta/beta3a.mat", beta = beta)
debugSource('~/OneDrive - Florida State Students/ongoing_projects/SEAS_code/code_R1/sliceH/simulation3b_error.R')
writeMat("beta/beta3b.mat", beta = beta)
debugSource('~/OneDrive - Florida State Students/ongoing_projects/SEAS_code/code_R1/sliceH/simulation4_error.R')
writeMat("beta/beta4.mat", beta = beta)
a <- read.table("output/M1_LassoSIR")
colMeans(a)
a <- read.table("M1_LassoSIR")
colMeans(a)
a <- read.table("M2_LassoSIR")
colMeans(a)
a <- read.table("output/M2_LassoSIR")
colMeans(a)
a <- read.table("M3a_LassoSIR")
colMeans(a)
a <- read.table("output/M3a_LassoSIR")
colMeans(a)
a <- read.table("M3b_LassoSIR")
colMeans(a)
a <- read.table("output/M3b_LassoSIR")
colMeans(a)
a <- read.table("M4_LassoSIR")
colMeans(a)
a <- read.table("output/M4_LassoSIR")
colMeans(a)
a <- read.table("M3b_1")
a_tmp <- a
a <- read.table("M3b_2")
a_tmp <- rbind(a_tmp, a)
a <- read.table("M3b_3")
a_tmp <- rbind(a_tmp, a)
a <- read.table("M3b_4")
a_tmp <- rbind(a_tmp, a)
colMeans(a_tmp)
a <- read.table("output/M3b")
colMeans(a)
a <- read.table("output/M3a")
colMeans(a)
a <- read.table("M3b_test_1")
a_tmp <- a
a <- read.table("M3b_test_2")
a_tmp <- rbind(a_tmp, a)
a <- read.table("M3b_test_3")
a_tmp <- rbind(a_tmp, a)
a <- read.table("M3b_test_4")
a_tmp <- rbind(a_tmp, a)
colMeans(a_tmp)
a <- read.table("M1")
colMeans(a)
a <- read.table("output/M1")
colMeans(a)
a <- read.table("output/M2")
colMeans(a)
a <- read.table("M2")
colMeans(a)
a <- read.table("output/M3a")
colMeans(a)
a <- read.table("M3a")
colMeans(a)
a <- read.table("output/M4")
colMeans(a)
a <- read.table("M4")
colMeans(a)
debugSource('~/OneDrive - Florida State Students/ongoing_projects/SEAS/SEAS_code/code_R1/simulation_p1000/simulation1_dat.R')
y_train
count
foldid
table(foldid)
n
debugSource('~/OneDrive - Florida State Students/ongoing_projects/SEAS/SEAS_code/code_R1/simulation_p1000/simulation1.R')
yclass
a <- read.table("M3b_1")
colMeans(a)
a <- read.table("M3b_2")
colMeans(a)
a <- read.table("M3b_3")
colMeans(a)
a <- read.table("M3b_4")
colMeans(a)
a <- read.table("M3b_1")
a_tmp <- a
a <- read.table("M3b_2")
a_tmp <- rbind(a_tmp,a)
a <- read.table("M3b_3")
a_tmp <- rbind(a_tmp,a)
a <- read.table("M3b_4")
a_tmp <- rbind(a_tmp,a)
colMeans(a_tmp)
a <- read.table("M5p1000_cv_test")
colMeans(a)
a <- read.table("output/M3b_1")
a_tmp <- a
a <- read.table("output/M3b_2")
a_tmp <- rbind(a_tmp,a)
a <- read.table("output/M3b_3")
a_tmp <- rbind(a_tmp,a)
colMeans(a_tmp)
a <- read.table("output/M3b")
colMeans(a)
getwd()
a <- read.table("output/M1")
colMeans(a)
a <- read.table("output/M1-1")
colMeans(a)
a <- read.table("output/M2")
colMeans(a)
a <- read.table("output/M2-1")
colMeans(a)
a <- read.table("output/M3a")
colMeans(a)
a <- read.table("output/M3a-1")
colMeans(a)
a <- read.table("output/M4-1")
colMeans(a)
a <- read.table("output/M4")
colMeans(a)
a <- read.table("output/M1_LassoSIR")
colMeans(a)
a <- read.table("output/M1_LassoSIR-1")
colMeans(a)
a <- read.table("output/M2_LassoSIR-1")
colMeans(a)
a <- read.table("output/M2_LassoSIR")
colMeans(a)
a <- read.table("output/M3a_LassoSIR")
colMeans(a)
a <- read.table("output/M3a_LassoSIR-1")
colMeans(a)
a <- read.table("output/M3b_LassoSIR")
colMeans(a)
a <- read.table("output/M3b_LassoSIR-1")
colMeans(a)
a <- read.table("output/M4_LassoSIR")
colMeans(a)
a <- read.table("output/M4_LassoSIR-1")
colMeans(a)
a <- read.table("../simulation_p3000/output/M2_1")
a_tmp <- a
a <- read.table("../simulation_p3000/output/M2_2")
a_tmp <- rbind(a_tmp,a)
a <- read.table("../simulation_p3000/output/M2_3")
a_tmp <- rbind(a_tmp,a)
colMeans(a_tmp)
a <- read.table("../simulation_p3000/output/M2")
colMeans(a)
getwd()
a <- read.table("../simulation_p3000/output/M3a_1")
a_tmp <- a
a <- read.table("../simulation_p3000/output/M3a_2")
a_tmp <- rbind(a_tmp,a)
a <- read.table("../simulation_p3000/output/M3a_3")
a_tmp <- rbind(a_tmp,a)
colMeans(a_tmp)
a <- read.table("../simulation_p3000/output/M3a")
colMeans(a)
a <- read.table("../simulation_p3000/output/M4_1")
a_tmp <- a
a <- read.table("../simulation_p3000/output/M4_2")
a_tmp <- rbind(a_tmp,a)
a <- read.table("../simulation_p3000/output/M4_3")
a_tmp <- rbind(a_tmp,a)
colMeans(a_tmp)
a <- read.table("../simulation_p3000/output/M4")
colMeans(a)
a <- read.table("../simulation_p3000/output/M1_LassoSIR")
colMeans(a)
a <- read.table("../simulation_p3000/output/M1_LassoSIR-1")
colMeans(a)
a <- read.table("../simulation_p3000/output/M2_LassoSIR")
colMeans(a)
a <- read.table("../simulation_p3000/output/M2_LassoSIR-1")
colMeans(a)
a <- read.table("../simulation_p3000/output/M3a_LassoSIR")
colMeans(a)
a <- read.table("../simulation_p3000/output/M3a_LassoSIR-1")
colMeans(a)
a <- read.table("../simulation_p3000/output/M3b_LassoSIR")
colMeans(a)
a <- read.table("../simulation_p3000/output/M3b_LassoSIR-1")
colMeans(a)
a <- read.table("../simulation_p3000/output/M4_LassoSIR")
colMeans(a)
a <- read.table("../simulation_p3000/output/M4_LassoSIR-1")
colMeans(a)
a <- read.table("../simulation_p3000/output/M3b_1")
a_tmp <- a
a <- read.table("../simulation_p3000/output/M3b_2")
a_tmp <- rbind(a_tmp,a)
a <- read.table("../simulation_p3000/output/M3b_3")
a_tmp <- rbind(a_tmp,a)
a <- read.table("../simulation_p3000/output/M3b_1")
a <- read.table("../simulation_p3000/output/M3b_2")
a <- read.table("../simulation_p3000/output/M3b_3")
library(LassoSIR)
?LassoSIR
