Laura Pyle

I have been using rjags for a multivariate normal model with unknown mean and variance. If I apply the model multiple times to the same dataset, I get very different results. I thought that I wasn't using enough iterations to get convergence, so I tried using runjags and the autorun.jags function to extend the simulations, and after 800K+ iterations, the Gelman-Rubin statistics was still above 1.05 for many of the parameters. I also get a note "unable to calculate the multi-variate psrf." I have...