-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path3_run_sims_over_grid.R
132 lines (96 loc) · 4.11 KB
/
3_run_sims_over_grid.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
source("2_run_sims.R")
sizes_irgtt <- list(
c(cluster_size_t = 5, clusters_per_arm_t = 12,
cluster_size_c = 5, clusters_per_arm_c = 12,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 5, clusters_per_arm_t = 30,
cluster_size_c = 5, clusters_per_arm_c = 30,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 5, clusters_per_arm_t = 50,
cluster_size_c = 5, clusters_per_arm_c = 50,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 10, clusters_per_arm_t = 12,
cluster_size_c = 10, clusters_per_arm_c = 12,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 10, clusters_per_arm_t = 30,
cluster_size_c = 10, clusters_per_arm_c = 30,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 10, clusters_per_arm_t = 50,
cluster_size_c = 10, clusters_per_arm_c = 50,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 20, clusters_per_arm_t = 12,
cluster_size_c = 20, clusters_per_arm_c = 12,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 20, clusters_per_arm_t = 30,
cluster_size_c = 20, clusters_per_arm_c = 30,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 20, clusters_per_arm_t = 50,
cluster_size_c = 20, clusters_per_arm_c = 50,
control_arm_irgtt = T, control_arm_no_ranef = T))
nsims <- 2000
core_n <- 9
Sys.time()
results_irgtt <- run_full_set(cores = core_n, seed = 11,
sizes = sizes_irgtt,
link_function = "identity",
nsims = nsims)
saveRDS(results_irgtt, "results/results_irgtt.rds")
Sys.time()
rm(results_irgtt)
gc()
Sys.time()
results_irgtt_binary <- run_full_set(cores = core_n, seed = 11,
sizes = sizes_irgtt,
t_eff = c(.35, 0), # make larger for binary
link_function = "logit",
nsims = nsims)
saveRDS(results_irgtt_binary, "results/results_irgtt_binary.rds")
Sys.time()
rm(results_irgtt_binary)
gc()
Sys.time()
results_irgtt_misspec <- run_full_set(cores = core_n, seed = 11,
tmle_misspec_test = T,
sizes = sizes_irgtt,
link_function = "identity",
nsims = nsims)
saveRDS(results_irgtt_misspec, "results/results_irgtt_misspec.rds")
Sys.time()
rm(results_irgtt_misspec)
gc()
sizes_unbalanced <- list(
# .5x smaller N in control arm
c(cluster_size_t = 5, clusters_per_arm_t = 30,
cluster_size_c = 5, clusters_per_arm_c = 15,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 5, clusters_per_arm_t = 50,
cluster_size_c = 5, clusters_per_arm_c = 25,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 10, clusters_per_arm_t = 30,
cluster_size_c = 10, clusters_per_arm_c = 15,
control_arm_irgtt = T, control_arm_no_ranef = T),
c(cluster_size_t = 10, clusters_per_arm_t = 50,
cluster_size_c = 10, clusters_per_arm_c = 25,
control_arm_irgtt = T, control_arm_no_ranef = T))
nsims <- 2000
core_n <- 9
Sys.time()
results_irgtt_unbalanced <- run_full_set(cores = core_n, seed = 11,
sizes = sizes_unbalanced,
link_function = "identity",
nsims = nsims)
saveRDS(results_irgtt_unbalanced, "results/results_irgtt_unbalanced.rds")
Sys.time()
rm(results_irgtt_unbalanced)
gc()
Sys.time()
results_irgtt_binary_unbalanced <- run_full_set(cores = core_n, seed = 11,
sizes = sizes_unbalanced,
t_eff = c(.35, 0), # make larger for binary
link_function = "logit",
nsims = nsims)
saveRDS(results_irgtt_binary_unbalanced, "results/results_irgtt_binary_unbalanced.rds")
Sys.time()
rm(results_irgtt_binary_unbalanced)
gc()
plan(sequential)