In ergm.MCMLE() for missing data, guard against constrained variance …

…being higher than unconstrained on a particular interation by treating the difference as 0 instead.
statnet · Jan 26, 2024 · 3f3ebbf · 3f3ebbf
1 parent c002361
commit 3f3ebbf
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/R/ergm.MCMLE.R b/R/ergm.MCMLE.R
@@ -315,6 +315,8 @@ ergm.MCMLE <- function(init, s, s.obs,
       estdiff <- NVL3(esteq.obs, colMeans(.), 0) - colMeans(esteq)
       pprec <- diag(sqrt(control$MCMLE.MCMC.precision), nrow=length(estdiff))
       Vm <- pprec%*%(cov(esteq) - NVL3(esteq.obs, cov(.), 0))%*%pprec
+      novar <- diag(Vm) <= 0
+      Vm[,novar] <- Vm[novar,] <- 0
       d2 <- xTAx_seigen(estdiff, Vm)
       if(d2<2) last.adequate <- TRUE
     }