X-Git-Url: https://code.communitydata.science/ml_measurement_error_public.git/blobdiff_plain/d0c5766bdf867a81a2477d2cac1d40812110af90..fa05dbab6bd2c5db6ed4eccf38cff03bb4fd6683:/simulations/summarize_estimator.R?ds=inline 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),