Merge pull request #7 from DoubleML/p-upd-repl-materials

Update Replication Code

Author: PhilippBach

Code/sim_iivm.R

@@ -12,7 +12,12 @@ library(doRNG)
 # Create a new directory plus subdirectories
 dir.create("simresults")
 dir.create("simresults/iivm")
-dir.create("Figures")
+
+save_plot = FALSE
+
+if (save_plot) {
+  dir.create("Figures")
+}
 
 date = Sys.Date()
 learner_name = "regr.cv_glmnet"
@@ -93,11 +98,13 @@ g_iivm = ggplot(df_iivm, aes(x = coef_resc)) +
 
 g_iivm
 
-ggsave(filename = paste0("Figures/densplot_IIVM_", n_obs, "_", dim_x, "_",
-                         "_", dml_procedure, "_", n_folds,
-                         "_", n_rep_folds, "_", learner_name, "_", R,
-                         "_", alpha, ".pdf"),
-       plot = g_iivm)
+if (save_plot) {
+  ggsave(filename = paste0("Figures/densplot_IIVM_", n_obs, "_", dim_x, "_",
+                           "_", dml_procedure, "_", n_folds,
+                           "_", n_rep_folds, "_", learner_name, "_", R,
+                           "_", alpha, ".pdf"),
+         plot = g_iivm)
+}
 
 results_summary = data.table(
   "learner" = learner_name,

Code/sim_irm.R

@@ -12,7 +12,12 @@ library(doRNG)
 # Create a new directory plus subdirectories
 dir.create("simresults")
 dir.create("simresults/irm")
-dir.create("Figures")
+
+save_plot = FALSE
+
+if (save_plot) {
+  dir.create("Figures")
+}
 
 date = Sys.Date()
 learner_name = "regr.cv_glmnet"
@@ -93,12 +98,13 @@ g_est = ggplot(df, aes(x = coef_resc)) +
 
 g_est
 
-ggsave(filename = paste0("Figures/densplot_IRM_", n_obs, "_", dim_x, "_",
-                         "_", dml_procedure, "_", n_folds,
-                         "_", n_rep_folds, "_", learner_name, "_", R,
-                         "_", alpha, ".pdf"),
-       plot = g_est)
-
+if (save_plot) {
+  ggsave(filename = paste0("Figures/densplot_IRM_", n_obs, "_", dim_x, "_",
+                           "_", dml_procedure, "_", n_folds,
+                           "_", n_rep_folds, "_", learner_name, "_", R,
+                           "_", alpha, ".pdf"),
+         plot = g_est)
+}
 results_summary = data.table(
   "learner" = learner_name,
   "n_folds" = n_folds,

Code/sim_plivX.R

@@ -12,7 +12,12 @@ library(doRNG)
 # Create a new directory plus subdirectories
 dir.create("simresults")
 dir.create("simresults/pliv")
-dir.create("Figures")
+
+save_plot = FALSE
+
+if (save_plot) {
+  dir.create("Figures")
+}
 
 date = Sys.Date()
 learner_name = "regr.cv_glmnet"
@@ -89,13 +94,13 @@ g_est = ggplot(df, aes(x = coef_resc)) +
 
 g_est
 
-
-ggsave(filename = paste0("Figures/densplot_PLIV_", n_obs, "_", dim_x, "_",
-                         "_", dml_procedure, "_", n_folds,
-                         "_", n_rep_folds, "_", learner_name, "_", R,
-                         "_", alpha, ".pdf"),
-       plot = g_est)
-
+if (save_plot) {
+  ggsave(filename = paste0("Figures/densplot_PLIV_", n_obs, "_", dim_x, "_",
+                           "_", dml_procedure, "_", n_folds,
+                           "_", n_rep_folds, "_", learner_name, "_", R,
+                           "_", alpha, ".pdf"),
+         plot = g_est)
+}
 
 results_summary = data.table(
   "learner" = learner_name,

Code/sim_plr.R

@@ -12,7 +12,12 @@ library(doRNG)
 # Create a new directory plus subdirectories
 dir.create("simresults")
 dir.create("simresults/plr")
-dir.create("Figures")
+
+save_plot = FALSE
+
+if (save_plot) {
+  dir.create("Figures")
+}
 
 date = Sys.Date()
 learner_name = "regr.ranger"
@@ -96,11 +101,14 @@ g_est = ggplot(df, aes(x = coef_resc)) +
 
 g_est
 
-ggsave(filename = paste0("Figures/densplot_PLR_", n_obs, "_", dim_x, "_",
-                         "_", dml_procedure, "_", n_folds,
-                         "_", n_rep_folds, "_", learner_name, "_", R,
-                         "_", alpha, ".pdf"),
-       plot = g_est)
+if (save_plot) {
+
+  ggsave(filename = paste0("Figures/densplot_PLR_", n_obs, "_", dim_x, "_",
+                           "_", dml_procedure, "_", n_folds,
+                           "_", n_rep_folds, "_", learner_name, "_", R,
+                           "_", alpha, ".pdf"),
+         plot = g_est)
+}
 
 results_summary = data.table(
   "learner" = learner_name,