fix lv extraction for unequal numbers of level 1 & level 2 variables

ecmerkle · Aug 21, 2023 · ce61b8a · ce61b8a
1 parent 8eed7c8
commit ce61b8a
Show file tree

Hide file tree

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

diff --git a/R/lvgqs.R b/R/lvgqs.R
@@ -300,6 +300,8 @@ samp_lvs_2lev <- function(mcobj, lavmodel, lavsamplestats, lavdata, lavpartable,
         standata$startrow <- c(1, standata$endrow[-length(standata$endrow)] + 1)
         standata$YX <- cbind(clusmns, matrix(0, nrow(clusmns), 2))
         standata$Ntot <- sum(standata$nclus[,2])
+        standata$Nobs <- with(standata, N_between + N_both)
+        standata$Obsvar <- with(standata, matrix(1:standata$Nobs, Np, N_between + N_both, byrow = TRUE))
         tmpmat2[j,,] <- lvgqs(modmat2, standata)
 
         ## now level 1