www.cmm.bristol.ac.uk/forum

Many thanks, Chris!

I estimated an unconditional cross-classified model with 3 levels. The intercept is quite a bit off from the sample mean. Any ideas for why that is? Thanks in advance.

I estimated a cross-classified model with 4 chains. I want to get a fit statistic called the Multivariate Effective Sample Size using the multiESS command in R.

I tried the following code after running the model and got this error message:
> multiESS(Model1@chains)
Error in multiESS(Model1@chains ...

Thanks, Chris. This was super helpful.

I want to obtain gelman-rubin statistics after running models with multiple chains. I believe this can be done with the gelman.diag() function in R.

The syntax is:

gelman.diag(x, confidence = 0.95, transform=FALSE, autoburnin=TRUE, multivariate=TRUE)

where x is a mcmc.list object with more ...