-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathbenchOne.R
More file actions
125 lines (94 loc) · 8.39 KB
/
Copy pathbenchOne.R
File metadata and controls
125 lines (94 loc) · 8.39 KB
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
124
125
# Written by Micha³ Makowski
# install.packages("MASS")
# install.packages("glasso")
# install.packages("huge")
require(MASS)
require(glasso)
require(huge)
source("auxilaryFunctions.R")
source("admmGLASSO.R")
source("admmGSLOPE.R")
source("benchOneFunctions")
data <- scale(huge.generator(n=100, d=100, graph="scale-free")$data)
str(data)
firstBench_0.1_F_T <- benchOne(data, alpha = .1)
firstBench_0.1_F_F <- benchOne(data, alpha = .1, truncate = F)
firstBench_0.1_T_T <- benchOne(data, alpha = .1, penalizeDiagonal = T)
firstBench_0.1_T_F <- benchOne(data, alpha = .1, truncate = F, penalizeDiagonal = T)
firstBench_0.05_F_T <- benchOne(data)
firstBench_0.05_F_F <- benchOne(data, truncate = F)
firstBench_0.05_T_T <- benchOne(data, penalizeDiagonal = T)
firstBench_0.05_T_F <- benchOne(data, truncate = F, penalizeDiagonal = T)
firstBench_0.01_F_T <- benchOne(data, alpha = .01)
firstBench_0.01_F_F <- benchOne(data, alpha = .01, truncate = F)
firstBench_0.01_T_T <- benchOne(data, alpha = .01, penalizeDiagonal = T)
firstBench_0.01_T_F <- benchOne(data, alpha = .01, truncate = F, penalizeDiagonal = T)
firstBench_0.05_F_T_ric <- benchOne(data, selectCriterion = "ric")
firstBench_0.05_F_T_ebic <- benchOne(data, selectCriterion = "ebic")
ggsave("firstBench_0.05_F_T.png", plotBenchOne(firstBench_0.05_F_T), device = "png", path = "./!01 Plots/01 first/01 criterion/", dpi = 200)
ggsave("firstBench_0.05_F_T_ric.png", plotBenchOne(firstBench_0.05_F_T_ric), device = "png", path = "./!01 Plots/01 first/01 criterion/", dpi = 200)
ggsave("firstBench_0.05_F_T_ebic.png", plotBenchOne(firstBench_0.05_F_T_ebic), device = "png", path = "./!01 Plots/01 first/01 criterion/", dpi = 200)
ggsave("firstBench_0.1_F_T.png", plotBenchOne(firstBench_0.1_F_T), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.1_T_F.png", plotBenchOne(firstBench_0.1_T_F), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.1_F_F.png", plotBenchOne(firstBench_0.1_F_F), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.1_T_T.png", plotBenchOne(firstBench_0.1_T_T), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.05_F_T.png", plotBenchOne(firstBench_0.05_F_T), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.05_T_F.png", plotBenchOne(firstBench_0.05_T_F), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.05_F_F.png", plotBenchOne(firstBench_0.05_F_F), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.05_T_T.png", plotBenchOne(firstBench_0.05_T_T), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.01_F_T.png", plotBenchOne(firstBench_0.01_F_T), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.01_T_F.png", plotBenchOne(firstBench_0.01_T_F), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.01_F_F.png", plotBenchOne(firstBench_0.01_F_F), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
ggsave("firstBench_0.01_T_T.png", plotBenchOne(firstBench_0.01_T_T), device = "png", path = "./!01 Plots/01 first/", dpi = 200)
##################################################################################################
data = scale(huge.generator(n=50, d=75, graph="hub")$data) # Generate data with hub structures
secondHubBench_0.1_F_T <- benchOne(data, alpha = .1)
secondHubBench_0.1_F_F <- benchOne(data, alpha = .1, truncate = F)
secondHubBench_0.1_T_T <- benchOne(data, alpha = .1, penalizeDiagonal = T)
secondHubBench_0.1_T_F <- benchOne(data, alpha = .1, truncate = F, penalizeDiagonal = T)
secondHubBench_0.05_F_T <- benchOne(data)
secondHubBench_0.05_F_F <- benchOne(data, truncate = F)
secondHubBench_0.05_T_T <- benchOne(data, penalizeDiagonal = T)
secondHubBench_0.05_T_F <- benchOne(data, truncate = F, penalizeDiagonal = T)
secondHubBench_0.01_F_T <- benchOne(data, alpha = .01)
secondHubBench_0.01_F_F <- benchOne(data, alpha = .01, truncate = F)
secondHubBench_0.01_T_T <- benchOne(data, alpha = .01, penalizeDiagonal = T)
secondHubBench_0.01_T_F <- benchOne(data, alpha = .01, truncate = F, penalizeDiagonal = T)
ggsave("secondHubBench_0.1_F_T.png", plotBenchOne(secondHubBench_0.1_F_T), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.1_T_F.png", plotBenchOne(secondHubBench_0.1_T_F), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.1_F_F.png", plotBenchOne(secondHubBench_0.1_F_F), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.1_T_T.png", plotBenchOne(secondHubBench_0.1_T_T), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.05_F_T.png", plotBenchOne(secondHubBench_0.05_F_T), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.05_T_F.png", plotBenchOne(secondHubBench_0.05_T_F), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.05_F_F.png", plotBenchOne(secondHubBench_0.05_F_F), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.05_T_T.png", plotBenchOne(secondHubBench_0.05_T_T), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.01_F_T.png", plotBenchOne(secondHubBench_0.01_F_T), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.01_T_F.png", plotBenchOne(secondHubBench_0.01_T_F), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.01_F_F.png", plotBenchOne(secondHubBench_0.01_F_F), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
ggsave("secondHubBench_0.01_T_T.png", plotBenchOne(secondHubBench_0.01_T_T), device = "png", path = "./!01 Plots/02 second/", dpi = 200)
#######################################
data = scale(huge.generator(n=200, d=200, graph="cluster")$data) # Generate data with hub structures
thirdClusterBench_0.1_F_T <- benchOne(data, alpha = .1)
thirdClusterBench_0.1_F_F <- benchOne(data, alpha = .1, truncate = F)
thirdClusterBench_0.1_T_T <- benchOne(data, alpha = .1, penalizeDiagonal = T)
thirdClusterBench_0.1_T_F <- benchOne(data, alpha = .1, truncate = F, penalizeDiagonal = T)
thirdClusterBench_0.05_F_T <- benchOne(data)
thirdClusterBench_0.05_F_F <- benchOne(data, truncate = F)
thirdClusterBench_0.05_T_T <- benchOne(data, penalizeDiagonal = T)
thirdClusterBench_0.05_T_F <- benchOne(data, truncate = F, penalizeDiagonal = T)
thirdClusterBench_0.01_F_T <- benchOne(data, alpha = .01)
thirdClusterBench_0.01_F_F <- benchOne(data, alpha = .01, truncate = F)
thirdClusterBench_0.01_T_T <- benchOne(data, alpha = .01, penalizeDiagonal = T)
thirdClusterBench_0.01_T_F <- benchOne(data, alpha = .01, truncate = F, penalizeDiagonal = T)
ggsave("thirdClusterBench_0.1_F_T.png", plotBenchOne(thirdClusterBench_0.1_F_T), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.1_T_F.png", plotBenchOne(thirdClusterBench_0.1_T_F), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.1_F_F.png", plotBenchOne(thirdClusterBench_0.1_F_F), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.1_T_T.png", plotBenchOne(thirdClusterBench_0.1_T_T), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.05_F_T.png", plotBenchOne(thirdClusterBench_0.05_F_T), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.05_T_F.png", plotBenchOne(thirdClusterBench_0.05_T_F), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.05_F_F.png", plotBenchOne(thirdClusterBench_0.05_F_F), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.05_T_T.png", plotBenchOne(thirdClusterBench_0.05_T_T), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.01_F_T.png", plotBenchOne(thirdClusterBench_0.01_F_T), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.01_T_F.png", plotBenchOne(thirdClusterBench_0.01_T_F), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.01_F_F.png", plotBenchOne(thirdClusterBench_0.01_F_F), device = "png", path = "./!01 Plots/03 third/", dpi = 200)
ggsave("thirdClusterBench_0.01_T_T.png", plotBenchOne(thirdClusterBench_0.01_T_T), device = "png", path = "./!01 Plots/03 third/", dpi = 200)