Hi RStudio Community,
I have run a Bayesian GLMM model in SAS (proc MCMC) and have aquired the posterior means and 95% Highest Posterior Density (HPD) intervals. I'm looking to validate this in R but I am struggling to get my results to match (datasets are consistent with each other). I have created models in R using 'brm' and 'jags.model' and then used 'MCMCsummary' to get the HPD intervals, but I just can't get my SAS and R results to match.
Any pointers would be much appreciated. TIA.