Fix bug in null sigma_V prior

seananderson · seananderson · commit b73581575aa3 · 2025-03-26T12:21:48.000-07:00
diff --git a/R/fit.R b/R/fit.R
@@ -1091,7 +1091,7 @@ sdmTMB <- function(
 
   # expand sigma_V priors if needed
   if (nrow(priors_sigma_V) == 1L && ncol(X_rw_ik) > 1L) {
-    priors_sigma_V <- replicate(ncol(X_rw_ik), priors_sigma_V, simplify = "matrix")
+    priors_sigma_V <- t(replicate(ncol(X_rw_ik), priors_sigma_V, simplify = "matrix"))
   }
   if (nrow(priors_sigma_V) != ncol(X_rw_ik)) {
     cli_abort("sigma_V (time-varying SD) priors do not match the fitted model.")

Original file line number	Diff line number	Diff line change
`@@ -1091,7 +1091,7 @@ sdmTMB <- function(`
`1091`	`1091`
`1092`	`1092`	`# expand sigma_V priors if needed`
`1093`	`1093`	`if (nrow(priors_sigma_V) == 1L && ncol(X_rw_ik) > 1L) {`
`1094`		`- priors_sigma_V <- replicate(ncol(X_rw_ik), priors_sigma_V, simplify = "matrix")`
	`1094`	`+ priors_sigma_V <- t(replicate(ncol(X_rw_ik), priors_sigma_V, simplify = "matrix"))`
`1095`	`1095`	`}`
`1096`	`1096`	`if (nrow(priors_sigma_V) != ncol(X_rw_ik)) {`
`1097`	`1097`	`cli_abort("sigma_V (time-varying SD) priors do not match the fitted model.")`