-
Notifications
You must be signed in to change notification settings - Fork 2
/
scenarios_test.R
123 lines (111 loc) · 5.22 KB
/
scenarios_test.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
116
117
118
119
120
121
122
123
sourceDir("datamakers")
fileName = 'data_path.txt'
path = gsub("[[:space:]]", "", readChar(fileName, file.info(fileName)$size))
# dsc-gtex-top10000
addScenario(dsc_gtex,name="null-top10000-2vs2", fn=datamaker,
args=list(tissue="Adipose-Subcutaneous", path=path,
Nsamp=2, Ngene=10000,
voom.normalize=TRUE),
seed=1:1)
# addScenario(dsc_gtex,name="null-top10000-10vs10", fn=datamaker,
# args=list(tissue="Adipose-Subcutaneous", path=path,
# Nsamp=10, Ngene=10000,
# voom.normalize=TRUE),
# seed=1:50)
#
# addScenario(dsc_gtex,name="null-top10000-50vs50", fn=datamaker,
# args=list(tissue="Adipose-Subcutaneous", path=path,
# Nsamp=50, Ngene=10000,
# voom.normalize=TRUE),
# seed=1:50)
# dsc-gtex-break-top10000
addScenario(dsc_gtex,name="null-break-top10000-2vs2", fn=datamaker,
args=list(tissue="Adipose-Subcutaneous",
Nsamp=2, Ngene=10000, breaksample=TRUE, path=path,
voom.normalize=TRUE),
seed=1:1)
# addScenario(dsc_gtex,name="null-break-top10000-10vs10", fn=datamaker,
# args=list(tissue="Adipose-Subcutaneous",
# Nsamp=10, Ngene=10000, breaksample=TRUE, path=path,
# voom.normalize=TRUE),
# seed=1:50)
#
# addScenario(dsc_gtex,name="null-break-top10000-50vs50", fn=datamaker,
# args=list(tissue="Adipose-Subcutaneous",
# Nsamp=50, Ngene=10000, breaksample=TRUE, path=path,
# voom.normalize=TRUE),
# seed=1:50)
# dscr-gtex-diff-top10000
addScenario(dsc_gtex,name="alt-top10000-2vs2", fn=datamaker,
args=list(tissue=c("Adipose-Subcutaneous","Lung"),
Nsamp=2, Ngene=10000, breaksample=FALSE, path=path,
voom.normalize=TRUE),
seed=1:1)
# addScenario(dsc_gtex,name="alt-top10000-10vs10", fn=datamaker,
# args=list(tissue=c("Adipose-Subcutaneous","Lung"),
# Nsamp=10, Ngene=10000, breaksample=FALSE, path=path,
# voom.normalize=TRUE),
# seed=1:50)
#
# addScenario(dsc_gtex,name="alt-top10000-50vs50", fn=datamaker,
# args=list(tissue=c("Adipose-Subcutaneous","Lung"),
# Nsamp=50, Ngene=10000, breaksample=FALSE, path=path,
# voom.normalize=TRUE),
# seed=1:50)
# dscr-gtex-break-diff-top10000
addScenario(dsc_gtex,name="alt-break-top10000-2vs2", fn=datamaker,
args=list(tissue=c("Adipose-Subcutaneous","Lung"),
Nsamp=2, Ngene=10000, breaksample=TRUE, path=path,
voom.normalize=TRUE),
seed=1:1)
# addScenario(dsc_gtex,name="alt-break-top10000-10vs10", fn=datamaker,
# args=list(tissue=c("Adipose-Subcutaneous","Lung"),
# Nsamp=10, Ngene=10000, breaksample=TRUE, path=path,
# voom.normalize=TRUE),
# seed=1:50)
#
# addScenario(dsc_gtex,name="alt-break-top10000-50vs50", fn=datamaker,
# args=list(tissue=c("Adipose-Subcutaneous","Lung"),
# Nsamp=50, Ngene=10000, breaksample=TRUE, path=path,
# voom.normalize=TRUE),
# seed=1:50)
# dscr-gtex-break-mix-top10000
addScenario(dsc_gtex,name="mix-break-top10000-2vs2", fn=datamaker,
args=list(tissue=c("Adipose-Subcutaneous","Lung"),
Nsamp=2, Ngene=10000, path=path,
voom.normalize=TRUE,
breaksample=TRUE, nullpi=0.9),
seed=1:1)
# addScenario(dsc_gtex,name="mix-break-top10000-10vs10", fn=datamaker,
# args=list(tissue=c("Adipose-Subcutaneous","Lung"),
# Nsamp=10, Ngene=10000, path=path,
# voom.normalize=TRUE,
# breaksample=TRUE, nullpi=0.9),
# seed=1:5)
#
# addScenario(dsc_gtex,name="mix-break-top10000-50vs50", fn=datamaker,
# args=list(tissue=c("Adipose-Subcutaneous","Lung"),
# Nsamp=50, Ngene=10000, path=path,
# voom.normalize=TRUE,
# breaksample=TRUE, nullpi=0.9),
# seed=1:5)
# dscr-gtex-poisthin-top10000
addScenario(dsc_gtex,name="poisthin-top10000-2vs2", fn=datamaker,
args=list(tissue=c("Adipose-Subcutaneous"),
Nsamp=2, Ngene=10000, path=path,
voom.normalize=TRUE,
poisthin=TRUE),
seed=1:1)
# addScenario(dsc_gtex,name="poisthin-top10000-10vs10", fn=datamaker,
# args=list(tissue=c("Adipose-Subcutaneous"),
# Nsamp=10, Ngene=10000, path=path,
# voom.normalize=TRUE,
# poisthin=TRUE),
# seed=1:5)
#
# addScenario(dsc_gtex,name="poisthin-top10000-50vs50", fn=datamaker,
# args=list(tissue=c("Adipose-Subcutaneous"),
# Nsamp=50, Ngene=10000, path=path,
# voom.normalize=TRUE,
# poisthin=TRUE),
# seed=1:5)