diff --git a/bilby/core/prior.py b/bilby/core/prior.py
index 9ea0b1934ddf35e4d021329fe6d4a7bd2036f955..5b64546ff315818fb6b551afb1e53e16428e3348 100644
--- a/bilby/core/prior.py
+++ b/bilby/core/prior.py
@@ -2922,7 +2922,7 @@ class MultivariateGaussianDist(object):
                 # sample the multivariate Gaussian keys
                 vals = np.random.uniform(0, 1, len(self))
 
-                samp = self.rescale(vals, mode=mode)
+                samp = np.atleast_1d(self.rescale(vals, mode=mode))
                 samps[i, :] = samp
 
                 # check sample is in bounds (otherwise perform another draw)