remove message + initialize stan matrices with valid input

ecmerkle · ecmerkle · commit 8411de31e93d · 2025-10-07T14:23:32.000-05:00
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: blavaan
 Title: Bayesian Latent Variable Analysis
-Version: 0.5-8.1363
+Version: 0.5-8.1365
 Authors@R: c(person(given = "Edgar", family = "Merkle",
                     role = c("aut", "cre"),
                     email = "merklee@missouri.edu",
diff --git a/R/blavaan.R b/R/blavaan.R
@@ -57,7 +57,6 @@ blavaan <- function(...,  # default lavaan arguments
 
     # multilevel functionality
     if("cluster" %in% dotNames) {
-        cat("blavaan NOTE: two-level models are new, please report bugs!\nhttps://github.com/ecmerkle/blavaan/issues\n\n")
         if(!(target == "stan")) stop("blavaan ERROR: two-level functionality is not available for ", target, ".")
     }
   
diff --git a/R/stanmarg_data.R b/R/stanmarg_data.R
@@ -42,7 +42,8 @@ group_sparse_skeleton <- function(skeleton) {
     tmpv <- c(tmpv, list(parts$v))
   }
 
-  vdat <- wdat <- matrix(0, Ng, max(g_len))
+  vdat <- matrix(1, Ng, max(g_len))
+  wdat <- matrix(0, Ng, max(g_len))
   udat <- matrix(0, Ng, max(u_len))
 
   for (g in 1:Ng) {

Original file line number	Diff line number	Diff line change
`@@ -57,7 +57,6 @@ blavaan <- function(..., # default lavaan arguments`
`57`	`57`
`58`	`58`	`# multilevel functionality`
`59`	`59`	`if("cluster" %in% dotNames) {`
`60`		`- cat("blavaan NOTE: two-level models are new, please report bugs!\nhttps://github.com/ecmerkle/blavaan/issues\n\n")`
`61`	`60`	`if(!(target == "stan")) stop("blavaan ERROR: two-level functionality is not available for ", target, ".")`
`62`	`61`	`}`
`63`	`62`
Original file line number	Diff line number	Diff line change
`@@ -42,7 +42,8 @@ group_sparse_skeleton <- function(skeleton) {`
`42`	`42`	`tmpv <- c(tmpv, list(parts$v))`
`43`	`43`	`}`
`44`	`44`
`45`		`- vdat <- wdat <- matrix(0, Ng, max(g_len))`
	`45`	`+ vdat <- matrix(1, Ng, max(g_len))`
	`46`	`+ wdat <- matrix(0, Ng, max(g_len))`
`46`	`47`	`udat <- matrix(0, Ng, max(u_len))`
`47`	`48`
`48`	`49`	`for (g in 1:Ng) {`