diff --git a/GPy/core/gp.py b/GPy/core/gp.py
index 1d2c9e2d..dc1519e1 100644
--- a/GPy/core/gp.py
+++ b/GPy/core/gp.py
@@ -325,7 +325,7 @@ class GP(Model):
         :returns: Ysim: set of simulations, a Numpy array (N x samples).
         """
         fsim = self.posterior_samples_f(X, size, full_cov=full_cov)
-        Ysim = self.likelihood.samples(Ysim, Y_metadata)
+        Ysim = self.likelihood.samples(fsim, Y_metadata)
 
         return Ysim