diff --git a/R/is_multivariate.R b/R/is_multivariate.R
index a7a9d1717..5f3bc4402 100644
--- a/R/is_multivariate.R
+++ b/R/is_multivariate.R
@@ -54,5 +54,10 @@ is_multivariate <- function(x) {
     return(isTRUE(ncol(x$coefficients) > 1L))
   }
 
-  return(FALSE)
+  vgam_classes <- c("vglm", "vgam")
+  if (inherits(x, vgam_classes)) {
+    return(isTRUE(x@extra$multiple.responses))
+  }
+
+  FALSE
 }
diff --git a/R/model_info.R b/R/model_info.R
index 612b57914..106d0350c 100644
--- a/R/model_info.R
+++ b/R/model_info.R
@@ -39,7 +39,7 @@
 #' * `is_hurdle`: model has zero-inflation component and is a hurdle-model (truncated family distribution)
 #' * `is_dispersion`: model has dispersion component (not only dispersion _parameter_)
 #' * `is_mixed`: model is a mixed effects model (with random effects)
-#' * `is_multivariate`: model is a multivariate response model (currently only works for _brmsfit_ objects)
+#' * `is_multivariate`: model is a multivariate response model (currently only works for _brmsfit_ and _vglm/vgam_ objects)
 #' * `is_trial`: model response contains additional information about the trials
 #' * `is_bayesian`: model is a Bayesian model
 #' * `is_gam`: model is a generalized additive model
@@ -1092,6 +1092,7 @@ model_info.vgam <- function(x, ...) {
     fitfam = faminfo@vfamily[1],
     logit.link = any(.string_contains("logit", faminfo@blurb)),
     link.fun = link.fun,
+    multi.var = is_multivariate(x),
     ...
   )
 }
diff --git a/inst/WORDLIST b/inst/WORDLIST
index e9ca0ae20..0f21e6da1 100644
--- a/inst/WORDLIST
+++ b/inst/WORDLIST
@@ -178,6 +178,7 @@ unstandardizing
 variates
 vectorized
 vgam
+vglm
 visualisation
 warmup
 warmups
diff --git a/man/model_info.Rd b/man/model_info.Rd
index aad1e90c5..fdd9bc615 100644
--- a/man/model_info.Rd
+++ b/man/model_info.Rd
@@ -53,7 +53,7 @@ is returned, where all values starting with \code{is_} are logicals.
 \item \code{is_hurdle}: model has zero-inflation component and is a hurdle-model (truncated family distribution)
 \item \code{is_dispersion}: model has dispersion component (not only dispersion \emph{parameter})
 \item \code{is_mixed}: model is a mixed effects model (with random effects)
-\item \code{is_multivariate}: model is a multivariate response model (currently only works for \emph{brmsfit} objects)
+\item \code{is_multivariate}: model is a multivariate response model (currently only works for \emph{brmsfit} and \emph{vglm/vgam} objects)
 \item \code{is_trial}: model response contains additional information about the trials
 \item \code{is_bayesian}: model is a Bayesian model
 \item \code{is_gam}: model is a generalized additive model
diff --git a/tests/testthat/test-gam.R b/tests/testthat/test-gam.R
index df96440cd..5430cc45f 100644
--- a/tests/testthat/test-gam.R
+++ b/tests/testthat/test-gam.R
@@ -40,7 +40,7 @@ test_that("model_info", {
 
 test_that("n_parameters", {
   expect_identical(n_parameters(m1), 5L)
-  expect_identical(n_parameters(m1, component = "conditional"), 1)
+  expect_identical(n_parameters(m1, component = "conditional"), 1L)
 })
 
 test_that("clean_names", {
diff --git a/tests/testthat/test-glmmTMB.R b/tests/testthat/test-glmmTMB.R
index a4a861bd4..673c1d05f 100644
--- a/tests/testthat/test-glmmTMB.R
+++ b/tests/testthat/test-glmmTMB.R
@@ -968,9 +968,6 @@ test_that("model_info, ordered beta", {
   out <- model_info(m)
   expect_true(out$is_orderedbeta)
   expect_identical(out$family, "ordbeta")
-  skip_on_cran()
-  out <- get_variance(m)
-  expect_equal(out$var.distribution, 1.44250604187634, tolerance = 1e-4)
 })
 
 
@@ -987,3 +984,25 @@ test_that("model_info, recognize ZI even without ziformula", {
   expect_true(out$is_zero_inflated)
   expect_true(out$is_hurdle)
 })
+
+
+skip_if_not_installed("withr")
+
+withr::with_environment(
+  new.env(),
+  test_that("get_variance, ordered beta", {
+    skip_if_not_installed("glmmTMB", minimum_version = "1.1.8")
+    skip_if_not_installed("datawizard")
+    skip_if_not_installed("lme4")
+    skip_on_cran()
+    data(sleepstudy, package = "lme4")
+    sleepstudy$y <- datawizard::normalize(sleepstudy$Reaction)
+    m <- glmmTMB::glmmTMB(
+      y ~ Days + (Days | Subject),
+      data = sleepstudy,
+      family = glmmTMB::ordbeta()
+    )
+    out <- get_variance(m)
+    expect_equal(out$var.distribution, 1.44250604187634, tolerance = 1e-4)
+  })
+)
diff --git a/tests/testthat/test-vgam.R b/tests/testthat/test-vgam.R
index 3d96e0208..bf60ed485 100644
--- a/tests/testthat/test-vgam.R
+++ b/tests/testthat/test-vgam.R
@@ -48,10 +48,11 @@ test_that("find_response", {
 })
 
 test_that("get_response", {
-  expect_equal(get_response(m1), hunua$agaaus)
+  expect_identical(get_response(m1), hunua$agaaus)
   expect_equal(
     get_response(m2),
-    data.frame(agaaus = hunua$agaaus, kniexc = hunua$kniexc)
+    data.frame(agaaus = hunua$agaaus, kniexc = hunua$kniexc),
+    ignore_attr = TRUE
   )
 })
 
@@ -195,7 +196,7 @@ test_that("find_parameters", {
 
 test_that("is_multivariate", {
   expect_false(is_multivariate(m1))
-  expect_false(is_multivariate(m2))
+  expect_true(is_multivariate(m2))
 })
 
 test_that("find_statistic", {