tweaked gelman-factors.r