How to get info of mcmc behaviour from the result of mcmc

#1

Dear all,
I was just wondering if we could get the info of the MCMC behaviour (especially the % of bad updates in both warning and sampling phases) from the R object that countains the result of the MCMC. That would help.

Sincerely,
Frédéric

#2

No clear reply to date. Seeing other messages I guess the info could be in :
attr(draws, “model_info”)$samplers[[1]]$numerical_rejections
where draws is the output of hmc. Am I right?

Sincerely,
Frédéric

#3

Hi Frédéric,

Thanks for your post - I think that looks about right, but I also think we should make this sort of information easier to get into.

I’ve linked an issue to this for greta, hopefully we can make a nice set of functions to get out this kind of information.

Cheers!