- fisher.info <- solve(mod.caroll.lik$hessian)
- coef <- mod.caroll.lik$par
- ci.upper <- coef + sqrt(diag(fisher.info)) * 1.96
- ci.lower <- coef - sqrt(diag(fisher.info)) * 1.96
-
+
+ ## tryCatch({
+ ## mod.calibrated.mle <- mecor(y ~ MeasError(w_pred, reference = x.obs) + z, df, B=400, method='efficient')
+ ## (mod.calibrated.mle)
+ ## (mecor.ci <- summary(mod.calibrated.mle)$c$ci['x.obs',])
+ ## result <- append(result, list(
+ ## Bxy.est.mecor = mecor.ci['Estimate'],
+ ## Bxy.ci.upper.mecor = mecor.ci['UCI'],
+ ## Bxy.ci.lower.mecor = mecor.ci['LCI'])
+ ## )
+
+
+
+ fischer.info <- NA
+ ci.upper <- NA
+ ci.lower <- NA
+
+ tryCatch({fischer.info <- solve(mod.caroll.lik$hessian)
+ ci.upper <- coef + sqrt(diag(fischer.info)) * 1.96
+ ci.lower <- coef - sqrt(diag(fischer.info)) * 1.96
+ },
+
+ error=function(e) {result[['error']] <- as.character(e)
+ })
+
+ coef <- mod.caroll.lik$par