X-Git-Url: https://code.communitydata.science/ml_measurement_error_public.git/blobdiff_plain/b52b4f7daaba8a877b041ddb24c8f36b466ddc5b..82fe7b0f482a71c95e8ae99f7e6d37b79357506a:/simulations/summarize_estimator.R?ds=sidebyside

diff --git a/simulations/summarize_estimator.R b/simulations/summarize_estimator.R
index e0e7622..3e4209f 100644
--- a/simulations/summarize_estimator.R
+++ b/simulations/summarize_estimator.R
@@ -31,8 +31,8 @@ summarize.estimator <- function(df, suffix='naive', coefname='x'){
                           var.est = var(.SD[[paste0('B',coefname,'y.est.',suffix)]]),
                           est.upper.95 = quantile(.SD[[paste0('B',coefname,'y.est.',suffix)]],0.975,na.rm=T),
                           est.lower.95 = quantile(.SD[[paste0('B',coefname,'y.est.',suffix)]],0.025,na.rm=T),
-                          mean.ci.upper = mean(.SD[[paste0('B',coefname,'y.ci.upper.',suffix)]]),
-                          mean.ci.lower = mean(.SD[[paste0('B',coefname,'y.ci.lower.',suffix)]]),
+                          mean.ci.upper = mean(.SD[[paste0('B',coefname,'y.ci.upper.',suffix)]],na.rm=T),
+                          mean.ci.lower = mean(.SD[[paste0('B',coefname,'y.ci.lower.',suffix)]],na.rm=T),
                           ci.upper.975 = quantile(.SD[[paste0('B',coefname,'y.ci.upper.',suffix)]],0.975,na.rm=T),
                           ci.upper.025 = quantile(.SD[[paste0('B',coefname,'y.ci.upper.',suffix)]],0.025,na.rm=T),
                           ci.lower.975 = quantile(.SD[[paste0('B',coefname,'y.ci.lower.',suffix)]],0.975,na.rm=T),