But, I’m not sure why you’d bother. The “approximate likelihood” is pretty easy to compute (I wouldn’t use logistic regression – it’s a low dimensional integral: use Gauss quadrature) and the resulting posterior satisfies pretty straightforward perturbation bounds w.r.t. this integration error. (This is even true for things like log-Gaussian Cox processes)

]]>