margloglik: error if cannot fill in unrestricted blocks of covariance…

… matrix
ecmerkle · Apr 15, 2024 · 650c94d · 650c94d
1 parent 6a964a6
commit 650c94d
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 2 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: blavaan
 Title: Bayesian Latent Variable Analysis
-Version: 0.5-3.1261
+Version: 0.5-3.1263
 Authors@R: c(person(given = "Edgar", family = "Merkle",
                     role = c("aut", "cre"),
                     email = "merklee@missouri.edu",

diff --git a/R/margloglik.R b/R/margloglik.R
@@ -88,7 +88,12 @@ margloglik <- function(lavpartable, lavmodel, lavoptions,
       tmpmat <- diag(lavpartable$est[varpars])
 
       if(length(covpars) > 0){
-        tmpmat[lower.tri(tmpmat)] <- lavpartable$est[covpars]
+        if(sum(lower.tri(tmpmat)) != length(covpars)){
+          ## TODO need to evaluate blocks of the matrix at a time
+          stop("blavaan ERROR: margloglik problem evaluating covariance matrix")
+        } else {
+          tmpmat[lower.tri(tmpmat)] <- lavpartable$est[covpars]
+        }
       }
       tmpmat <- tmpmat + t(tmpmat)
       diag(tmpmat) <- diag(tmpmat)/2