When drawing approximate samples from the marginal posterior \pi(\phi, \eta \mid y), we don’t need the normalizing constant in \log \pi(y \mid \theta^*, \eta). Our current unit tests benchmark against GPStuff, which as far as I can tell computes the normalizing constant. Is there a reason for this?
Edit: making corrections to the equations, following Aki’s suggestion.
I think \log \pi(y \mid \theta, \eta) should be \log \pi(y \mid \theta^*, \eta)
I think \pi(\phi \mid y) should be \pi(\phi, \eta \mid y)
I think \log \pi(y \mid \theta, \eta) should be \log \pi(y \mid \theta^*, \eta). As \theta^* is fixed, the normalization term with respect to \theta^* is constant and not needed for MCMC, but in most cases the computational cost is negligible compared to matrix and Laplace operations, and it was easier to make comparisons to full MCMC when not dropping those terms. For models with \eta as parameter, the normalization term depends on \eta and needs to be included.