powered by
Full Gaussian likelihood (includes marginals; only needed for model comparison)
llike_full(data, psi, sigma2)