-
Notifications
You must be signed in to change notification settings - Fork 0
/
table2.R
57 lines (39 loc) · 1.72 KB
/
table2.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
### table2.R ---
#----------------------------------------------------------------------
## Author: Helene Charlotte Wiesj Rijtgaard, Thomas Alexander Gerds
## Created: Nov 10 2020 (17:58)
## Version:
## Last-Updated: Nov 11 2020 (17:21)
## By: Thomas Alexander Gerds
## Update #: 11
#----------------------------------------------------------------------
##
### Commentary:
##
### Change Log:
#----------------------------------------------------------------------
##
### Code:
source("./examples/load.R")
# be careful setting the number of cores and to watch the memory for large values of K
M <- 1000 # number of simulations
nc <- 40 # number of cores
n <- 1000 # sample size
if (FALSE) { # K=5
table2.K5.conTMLE <- runTMLE(no_cores=nc,K=5,misspecify.init = TRUE,M = M,n = n,max.iter=25,eps=0.0000001,seed=1920,progress.bar=-1)
saveRDS(table2.K5.conTMLE,file="./examples/table2-K5-conTMLE.rds")
}
if (FALSE) { # K=30
table2.K30.conTMLE <- runTMLE(no_cores=nc,K=30,misspecify.init = TRUE,M = M,n = n,max.iter=25,eps=0.0000001,seed=9170,progress.bar=-1)
saveRDS(table2.K30.conTMLE,file="./examples/table2-K30-conTMLE.rds")
}
if (FALSE) { # K=50
table2.K50.conTMLE <- runTMLE(no_cores=nc,K=50,misspecify.init = TRUE,M = M,n = n,max.iter=25,eps=0.0000001,seed=9170,progress.bar=-1)
saveRDS(table2.K50.conTMLE,file="./examples/table2-K50-conTMLE.rds")
}
if (FALSE) { # K=100
table2.K100.conTMLE <- runTMLE(no_cores=nc,K=100,misspecify.init = TRUE,M = M,n = n,max.iter=25,eps=0.0000001,seed=9170,progress.bar=-1)
saveRDS(table2.K100.conTMLE,file="./examples/table2-K100-conTMLE.rds")
}
######################################################################
### table2.R ends here