Merge pull request #418 from imbs-hl/catch_rInfJack_calibration_error

Catch rinfJack calibration errors

R/infinitesimalJackknife.R

@@ -98,8 +98,18 @@ rInfJack = function(pred, inbag, calibrate = TRUE, used.trees = NULL) {
     sigma2 = (delta^2 + (1 - delta)^2) / (2 * (1 - delta)^2) * sigma2.ss
 
     # Use Monte Carlo noise scale estimate for empirical Bayes calibration
-    vars.calibrated = calibrateEB(vars, sigma2)
-    results$var.hat = vars.calibrated
+    results = tryCatch(
+      expr = {
+        vars.calibrated = calibrateEB(vars, sigma2)
+        results$var.hat = vars.calibrated
+        results
+      }, 
+      error = function(e) {
+        warning(sprintf("Calibration failed with error:\n%sFalling back to non-calibrated variance estimates.", e))
+        results = rInfJack(pred, inbag, calibrate = FALSE, used.trees = used.trees)
+        return(results)
+      }
+    )
   }
 
   return(results)