In ergm.estimate(), if the steplen < 1 and the unshrunken Hessian is …

…not SNND, return NA in its place. fixes #551
statnet · Jan 22, 2024 · ff704e4 · ff704e4
1 parent eed8b5a
commit ff704e4
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: ergm
-Version: 4.6-7294
-Date: 2024-01-10
+Version: 4.6-7303
+Date: 2024-01-22
 Title: Fit, Simulate and Diagnose Exponential-Family Models for Networks
 Authors@R: c(
              person(c("Mark", "S."), "Handcock", role=c("aut"), email="[email protected]"),

diff --git a/R/ergm.estimate.R b/R/ergm.estimate.R
@@ -337,9 +337,9 @@ ergm.estimate<-function(init, model, statsmatrices, statsmatrices.obs=NULL,
                         eta0=eta0, etamap=etamap.no
                         )
     }
-    
+
     covar <- matrix(NA, ncol=length(theta), nrow=length(theta))
-    covar[!model$etamap$offsettheta, !model$etamap$offsettheta] <- sginv(-Lout$hessian, tol=.Machine$double.eps^(3/4))
+    covar[!model$etamap$offsettheta, !model$etamap$offsettheta] <- ERRVL(try(sginv(-Lout$hessian, tol=.Machine$double.eps^(3/4)), TRUE), NA)
     dimnames(covar) <- list(names(theta),names(theta))
     He <- matrix(NA, ncol=length(theta), nrow=length(theta))
     He[!model$etamap$offsettheta, !model$etamap$offsettheta] <- Lout$hessian