cleaning replication files

nppackages · Mar 27, 2024 · 63bfb76 · 63bfb76
1 parent 9558064
commit 63bfb76
Show file tree

Hide file tree

Showing 5 changed files with 42 additions and 24 deletions.
diff --git a/R/bw_table.R b/R/bw_table.R
@@ -60,4 +60,3 @@ colnames(table_data) = c("BW", "bias", "sd", "rmse", "WBC CE", "RBC CE",
 rownames(table_data) = h_grid/0.5
 #print table
 table_data
-print(xtable(table_data), file='bw_table.tex')
diff --git a/R/lpcde_comparison.R b/R/lpcde_comparison.R
@@ -61,26 +61,3 @@ plot(lpcde_est_q2,xlabel="Sepal length", ylabel="density", title=TeX(r'($f(y|x=4
                                                                                                                                                                                      panel.background = ggplot2::element_blank(), text = ggplot2::element_text(size = 7))
 plot(lpcde_est_q3, xlabel="Sepal length", ylabel="density", title=TeX(r'($f(y|x=5.1)$ with confidence bands)'))+ ggplot2::ylim(-0.8, 1.6) + ggplot2::xlim(4.8, 6.9) + ggplot2::theme(legend.position='none', panel.grid.major = ggplot2::element_blank(), panel.grid.minor = ggplot2::element_blank(),
                                                                                                                                                                                      panel.background = ggplot2::element_blank(), text = ggplot2::element_text(size = 7))
-
-
-
-
-
-#open3d()
-#
-#lines3d(x_evals[2], lpcde_est_q1$Estimate[,1], lpcde_est_q1$Estimate[,3])
-#lines3d(x_evals[3], lpcde_est_q2$Estimate[,1], lpcde_est_q2$Estimate[,3])
-#lines3d(x_evals[4], lpcde_est_q3$Estimate[,1], lpcde_est_q3$Estimate[,3])
-#
-#lines3d(x_evals[2], hdr_est_q1$y, hdr_est_q1$z, color='red')
-#lines3d(x_evals[3], hdr_est_q2$y, hdr_est_q2$z, color='red')
-#lines3d(x_evals[4], hdr_est_q3$y, hdr_est_q3$z, color='red')
-#
-#lines3d(x_evals[2], lpcde_est_q1$Estimate[,1], np_est_q1, color='blue')
-#lines3d(x_evals[3], lpcde_est_q2$Estimate[,1], np_est_q2, color='blue')
-#lines3d(x_evals[4], lpcde_est_q3$Estimate[,1], np_est_q3, color='blue')
-#
-#axes3d()
-#title3d('', '', 'Petal Length', 'Sepal Length', 'density')
-#legend3d("topright", c("lpcde", "hdrcde", "np"), pch = c(1,25, 16), col = c('black', 'red', 'blue'), cex=1)
-
diff --git a/R/pw_table_b.tex b/R/pw_table_b.tex
@@ -0,0 +1,14 @@
+% latex table generated in R 4.1.2 by xtable 1.8-4 package
+% Wed Mar 27 10:34:50 2024
+\begin{table}[ht]
+\centering
+\begin{tabular}{rrrrrrrrr}
+  \hline
+ & BW & bias & sd & rmse & WBC CE & RBC CE & WBC AL & RBC AL \\
+  \hline
+1 & 1.00 & 0.02 & 0.01 & 0.02 & 65 & 93 & 0.04 & 0.12 \\
+  2 & 0.80 & 0.01 & 0.02 & 0.03 & 80 & 95 & 0.07 & 0.23 \\
+  3 & 0.90 & 0.04 & 0.02 & 0.04 & 31 & 94 & 0.06 & 0.17 \\
+   \hline
+\end{tabular}
+\end{table}
diff --git a/R/pw_table_int.tex b/R/pw_table_int.tex
@@ -0,0 +1,14 @@
+% latex table generated in R 4.1.2 by xtable 1.8-4 package
+% Wed Mar 27 10:34:47 2024
+\begin{table}[ht]
+\centering
+\begin{tabular}{rrrrrrrrr}
+  \hline
+ & BW & bias & sd & rmse & WBC CE & RBC CE & WBC AL & RBC AL \\
+  \hline
+1 & 0.51 & 0.02 & 0.01 & 0.03 & 69 & 94 & 0.06 & 0.19 \\
+  2 & 0.58 & 0.00 & 0.01 & 0.02 & 67 & 92 & 0.05 & 0.15 \\
+  3 & 0.80 & 0.01 & 0.01 & 0.02 & 49 & 95 & 0.03 & 0.08 \\
+   \hline
+\end{tabular}
+\end{table}
diff --git a/R/pw_table_nb.tex b/R/pw_table_nb.tex
@@ -0,0 +1,14 @@
+% latex table generated in R 4.1.2 by xtable 1.8-4 package
+% Wed Mar 27 10:34:49 2024
+\begin{table}[ht]
+\centering
+\begin{tabular}{rrrrrrrrr}
+  \hline
+ & BW & bias & sd & rmse & WBC CE & RBC CE & WBC AL & RBC AL \\
+  \hline
+1 & 0.60 & 0.01 & 0.02 & 0.03 & 79 & 96 & 0.06 & 0.20 \\
+  2 & 0.60 & 0.02 & 0.01 & 0.03 & 58 & 93 & 0.06 & 0.18 \\
+  3 & 0.60 & 0.01 & 0.02 & 0.03 & 87 & 98 & 0.07 & 0.22 \\
+   \hline
+\end{tabular}
+\end{table}