-
Notifications
You must be signed in to change notification settings - Fork 1
/
plot_simulations_all.R
70 lines (58 loc) · 1.42 KB
/
plot_simulations_all.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
suppressMessages({
source("utils/plot_methods_all_group_corr.R")
library("tools")
library("data.table")
library("ggpubr")
library("scales")
})
file_path <- paste0(getwd(), "/results/results_csv/")
filename <- paste0(file_path, "simulation_results_blocks_100_groups_n_100_p_50_1::1_non_stack_cpi")
N_CPU <- 100
list_func <- c(
"Marg",
"Permfit-DNN",
"CPI-DNN",
"CPI-RF",
"gpfi",
"gopfi",
"dgi",
"goi"
)
run_plot_auc <- TRUE
run_plot_type1error <- TRUE
run_plot_power <- TRUE
run_time <- TRUE
filename_lbl <- strsplit(filename, "_results_")[[1]][2]
if (run_plot_auc) {
plot_method(paste0(filename, ".csv"),
paste0("AUC_", filename_lbl),
compute_auc,
title = "AUC",
list_func = list_func
)
}
if (run_plot_type1error) {
plot_method(paste0(filename, ".csv"),
paste0("type1error_", filename_lbl),
compute_pval,
upper_bound = 0.05,
title = "Type I Error",
list_func = list_func
)
}
if (run_plot_power) {
plot_method(paste0(filename, ".csv"),
paste0("power_", filename_lbl),
compute_power,
upper_bound = 0.05,
title = "Power",
list_func = list_func
)
}
if (run_time) {
plot_time(paste0(filename, ".csv"),
paste0("time_bars_", filename_lbl),
list_func = list_func,
N_CPU = N_CPU
)
}