avoid complete failure if vcov is NaN

stemangiola · Sep 16, 2023 · 0b395c4 · 0b395c4
1 parent f65d233
commit 0b395c4
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/R/glmmSeq.R b/R/glmmSeq.R
@@ -260,7 +260,14 @@ glmmTMBcore = function (geneList, fullFormula, reduced, data, family, control,
   if (!inherits(fit, "try-error")) {
     singular <- conv <- NA
     stdErr <- suppressWarnings(coef(summary(fit))$cond[, 2])
+
     vcov. <- vcov(fit)$cond
+
+    if(vcov. |> as.numeric() |> is.nan() |> all())
+      return(list(stats = NA, coef = NA, stdErr = NA, chisq = NA,
+                  df = NA, predict = NA, optinfo = NA, ci_random_effect_df = NA,
+                  message = "vcov. failed to calculate", tryErrors = fit[1]))
+
     fixedEffects <- glmmTMB::fixef(fit)$cond
     disp <- glmmTMB::sigma(fit)
     msg <- fit$fit$message