Minor fix to dist_sample methods dimensions for multivariate dist

NAMESPACE

@@ -167,9 +167,9 @@ S3method(format,dist_f)
 S3method(format,dist_gamma)
 S3method(format,dist_geometric)
 S3method(format,dist_gev)
-S3method(format,dist_gpd)
 S3method(format,dist_gh)
 S3method(format,dist_gk)
+S3method(format,dist_gpd)
 S3method(format,dist_gumbel)
 S3method(format,dist_hypergeometric)
 S3method(format,dist_inflated)
@@ -486,9 +486,9 @@ export(dist_f)
 export(dist_gamma)
 export(dist_geometric)
 export(dist_gev)
-export(dist_gpd)
 export(dist_gh)
 export(dist_gk)
+export(dist_gpd)
 export(dist_gumbel)
 export(dist_hypergeometric)
 export(dist_inflated)

R/dist_sample.R

@@ -72,7 +72,7 @@ quantile.dist_sample <- function(x, p, type = "marginal", ..., na.rm = TRUE){
   if(is.matrix(x$x)) {
     # Marginal quantiles
     return(
-      matrix(apply(x$x, 2, quantile, p = p, ..., na.rm = na.rm), nrow = 1)
+      matrix(apply(x$x, 2, quantile, p = p, ..., na.rm = na.rm), nrow = length(p))
     )
   }
   quantile(x$x, probs = p, ..., na.rm = na.rm, names = FALSE)
@@ -106,7 +106,7 @@ mean.dist_sample <- function(x, ...){
 #' @export
 median.dist_sample <- function(x, na.rm = FALSE, ...){
   if(is.matrix(x$x))
-    apply(x$x, 2, median, na.rm = na.rm, ...)
+    matrix(apply(x$x, 2, median, na.rm = na.rm, ...), nrow = 1L)
   else
     median(x$x, na.rm = na.rm, ...)
 }