To sample the prior predictive correctly I think you’d need to use gp.Latent. The gp.Marginal implementation marginalizes out the actual GP, which makes it more efficient, but the actual GP f isn’t explicitly included in the model.
1 Like