added script to compute gelman factors