Likelihood of latent variables in hierarchical model for model comparison

I don’t think so, because the log_likelihood group is added as sampling occurs. What do you need that isn’t captured in that group? I can’t quite parse your post-processing snippet.