if (FALSE) { # require("rstanarm", quietly = TRUE)
# \donttest{
library(rstanarm)
data("pbcLong")
model <- suppressWarnings(stan_mvmer(
formula = list(
logBili ~ year + (1 | id),
albumin ~ sex + year + (year | id)
),
data = pbcLong,
chains = 1, cores = 1, seed = 12345, iter = 1000,
show_messages = FALSE, refresh = 0
))
f <- find_formula(model)
is_multivariate(model)
is_multivariate(f)
# }
}
Run the code above in your browser using DataLab