Hi I am trying to use NUTS to sample poterior distribution of two hyperparameters for Gaussian Process Regression. I have two questions of it:
-
How to get the log probability of samples with pyro? I found an answer for numpyro but that seems do not work for pyro.
-
The results are shown in the figure below, it seems like the posterior samples does not converge in areas with high likelihood. I am wondering is it a nature bahabviour of NUTS or how can I improve it? Note, Empirical Bayes indicate the optimised hyperparameters by SVI.