/
supplementary_plots_pesel_homo_hetero.R
76 lines (58 loc) · 2.84 KB
/
supplementary_plots_pesel_homo_hetero.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
#' supplementary materials
#' comparison of pesel n/p homo/hetero on robustness
library(dplyr)
source("plot_templates.R")
date = NULL
scheme.names <- c("Student noise", "Additional variables", "Lognormal noise")
method.names <- c("Minka BIC", "Laplace evidence", "Properly dervied BIC",
"GCV", "Rajan - penalized version", "Passemier", "CSV")
method_labels = expression(paste("PESEL", phantom()["n"]^{"hetero"}), "Laplace evidence",
paste("PESEL", phantom()["p"]^{"hetero"}), "GCV",
paste("PESEL", phantom()["p"]^{"homo"}),"Passemier", "CSV")
y <- NULL
for(batch in c("A", "B", "C", "D", "E")){
f <- paste0("data/robustness_fixed_simulations", batch, "_n_100_", date, ".Rdata")
load(f)
for(sim.scenario in c(1,2,3)){
result <- results[[sim.scenario]]
characteristics <- lapply(strsplit(names(result), "_"), as.numeric)
for(i in 1:length(result)){
y <- rbind(y, data.frame(pc=as.vector(result[[i]]), method=rep(method.names, numb.repetitions),
snr=characteristics[[i]][1], variables=characteristics[[i]][2],
scenario=scheme.names[sim.scenario]))
}
}
}
y$method <- factor(y$method, levels = method.names[c(1:2,4,3,5:7)])
###############################################
############## VARS PLOTS #####################
###############################################
selected.variables.number <- c(150, 800)
selected.snrs <- c(0.25, 0.5, 0.75, 1, 2, 4)
sim.scenario <- 1
selected.method <- method.names[c(1,3,5)]
selected.variables.number <- c(150, 800)
selected.snrs <- c(0.25, 0.5, 0.63, 0.75, 1, 2.00, 8)
sim.scenario = 2
selected.method <- method.names[c(1,3,5)]
selected.variables.number <- c(150, 800)
selected.snrs <- c(0.25, 0.375, 0.5, 0.63, 1, 2.00, 8)
sim.scenario <- 3
selected.method <- method.names[c(1,3,5)]
new.labels <- method_labels[c(1,3,5)]
plot_facet_vars(y, scheme.names, sim.scenario, selected.method, selected.snrs, selected.variables.number, new.labels)
plot_name <- paste0("supplementary_", gsub(" ", "_", scheme.names[sim.scenario]), "_",
paste0(selected.variables.number, collapse = "_"), ".png")
ggsave(filename = plot_name, height = 8, width = 16)
###############################################
############### SNR PLOTS #####################
###############################################
selected.variables.number <- c(50, 75, 150, 400, 800)
selected.snrs <- c(1,4)
sim.scenario = 1
selected.method <- method.names[c(1,3,5)]
new.labels <- method_labels[c(1,3,5)]
plot_facet_snr(y, scheme.names, sim.scenario, selected.method, selected.snrs, selected.variables.number)
plot_name <- paste0("supplementary_", gsub(" ", "_", scheme.names[sim.scenario]), "_",
paste0(selected.snrs, collapse = "_"), ".png")
ggsave(filename = plot_name, height = 8, width = 16)