updated BW table simuls

R/bw_table.R

@@ -1,15 +1,15 @@
 # Script for generating table 1 in Cattaneo, Chandak, Jansson and Ma (2022)
 
 #parameters
-n = 1000
-num_sims = 1000
+n = 2000
+num_sims = 100
 y = 0.0
 x = 0.0
 #true density function
 true_dens = stats::dnorm(y, mean = 0, sd=1)
 
 #grid for bandwidth range
-h_grid = c(0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.1, 1.2, 1.3, 1.4, 1.5)*0.6
+h_grid = c(0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.1, 1.2, 1.3, 1.4, 1.5)*0.5
 
 bw_testing = function(s){
   print(s)
@@ -57,7 +57,7 @@ table_data = table_data[2:nrow(table_data), ]
 table_data[, 2] = abs(table_data[,2])
 colnames(table_data) = c("BW", "bias", "sd", "rmse", "WBC CE", "RBC CE",
                          "WBC AL", "RBC AL")
-#rownames(table_data) = h_grid/0.6
+rownames(table_data) = h_grid/0.5
 #print table
 table_data
-#print(xtable(table_data), file='bw_table.tex')
+print(xtable(table_data), file='bw_table.tex')

R/bw_table.tex

@@ -1,22 +1,22 @@
 % latex table generated in R 4.1.2 by xtable 1.8-4 package
-% Tue Mar 26 14:26:35 2024
+% Wed Mar 27 13:55:05 2024
 \begin{table}[ht]
 \centering
 \begin{tabular}{rrrrrrrrr}
   \hline
  & BW & bias & sd & rmse & WBC CE & RBC CE & WBC AL & RBC AL \\ 
   \hline
-0.5 & 0.30 & 0.01 & 0.14 & 0.15 & 100.00 & 100.00 & 0.56 & 1.90 \\ 
-  0.6 & 0.36 & 0.01 & 0.08 & 0.10 & 98.00 & 100.00 & 0.32 & 1.08 \\ 
-  0.7 & 0.42 & 0.01 & 0.05 & 0.06 & 96.00 & 100.00 & 0.20 & 0.67 \\ 
-  0.8 & 0.48 & 0.02 & 0.03 & 0.05 & 89.00 & 99.00 & 0.13 & 0.44 \\ 
-  0.9 & 0.54 & 0.02 & 0.02 & 0.04 & 77.00 & 96.00 & 0.10 & 0.32 \\ 
-  1 & 0.60 & 0.02 & 0.02 & 0.03 & 70.00 & 94.00 & 0.07 & 0.23 \\ 
-  1.1 & 0.66 & 0.03 & 0.01 & 0.04 & 38.00 & 86.00 & 0.05 & 0.18 \\ 
-  1.2 & 0.72 & 0.03 & 0.01 & 0.03 & 37.00 & 79.00 & 0.04 & 0.14 \\ 
-  1.3 & 0.78 & 0.04 & 0.01 & 0.04 & 13.00 & 74.00 & 0.03 & 0.11 \\ 
-  1.4 & 0.84 & 0.04 & 0.01 & 0.04 & 10.00 & 69.00 & 0.03 & 0.09 \\ 
-  1.5 & 0.90 & 0.04 & 0.01 & 0.04 & 4.00 & 76.00 & 0.02 & 0.07 \\ 
+0.5 & 0.25 & 0.00 & 0.12 & 0.13 & 100.00 & 100.00 & 0.46 & 1.57 \\ 
+  0.6 & 0.30 & 0.01 & 0.07 & 0.08 & 100.00 & 100.00 & 0.27 & 0.90 \\ 
+  0.7 & 0.35 & 0.01 & 0.04 & 0.06 & 93.00 & 99.00 & 0.17 & 0.56 \\ 
+  0.8 & 0.40 & 0.01 & 0.03 & 0.04 & 89.00 & 99.00 & 0.11 & 0.38 \\ 
+  0.9 & 0.45 & 0.01 & 0.02 & 0.03 & 78.00 & 97.00 & 0.08 & 0.27 \\ 
+  1 & 0.50 & 0.02 & 0.01 & 0.03 & 64.00 & 95.00 & 0.06 & 0.20 \\ 
+  1.1 & 0.55 & 0.02 & 0.01 & 0.03 & 47.00 & 87.00 & 0.04 & 0.15 \\ 
+  1.2 & 0.60 & 0.02 & 0.01 & 0.03 & 35.00 & 81.00 & 0.04 & 0.12 \\ 
+  1.3 & 0.65 & 0.02 & 0.01 & 0.03 & 34.00 & 70.00 & 0.03 & 0.09 \\ 
+  1.4 & 0.70 & 0.03 & 0.01 & 0.03 & 9.00 & 69.00 & 0.02 & 0.07 \\ 
+  1.5 & 0.75 & 0.03 & 0.00 & 0.03 & 9.00 & 64.00 & 0.02 & 0.06 \\ 
    \hline
 \end{tabular}
 \end{table}