def samples(self, gp, Y_metadata=None):
orig_shape = gp.shape<EOL>gp = gp.flatten()<EOL>Ysim = np.array([sp.stats.fisk.rvs(self.r, loc=<NUM_LIT:0>, scale=self.gp_link.transf(f)) for f in gp])<EOL>return Ysim.reshape(orig_shape)<EOL>
Returns a set of samples of observations based on a given value of the latent variable.

:param gp: latent variable