Vectorise equicoordinate quantiles for multivariate normal

mitchelloharawild · Sep 15, 2024 · 1410cd1 · 1410cd1
1 parent 4f5b42c
commit 1410cd1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/dist_multivariate_normal.R b/R/dist_multivariate_normal.R
@@ -62,7 +62,7 @@ quantile.dist_mvnorm <- function(x, p, type = c("equicoordinate", "marginal"),
                  sd = rep(diag(sqrt(x[["sigma"]])), each = length(p)), ...)
   } else {
     require_package("mvtnorm")
-    mvtnorm::qmvnorm(p, mean = x[["mu"]], sigma = x[["sigma"]], ...)$quantile
+    vapply(p, function(...) mvtnorm::qmvnorm(...)$quantile, numeric(1L), mean = x[["mu"]], sigma = x[["sigma"]], ...)
   }
 
   matrix(q, nrow = length(p), ncol = dim(x))