In src/likelihoods/likelihoods.jl the predict_obs uses Gaussian quadrature to find the mean and variance of predictions when the likelihood is not normal. There really should be tests for this function (and the Gaussian and Bernoulli cases), which looks quite fiddly to me.
In
src/likelihoods/likelihoods.jlthepredict_obsuses Gaussian quadrature to find the mean and variance of predictions when the likelihood is not normal. There really should be tests for this function (and the Gaussian and Bernoulli cases), which looks quite fiddly to me.