Working Gelman Rubin, better d norm handling

Gelman Rubin diagnostics integrated. d now supports separate d's for each
sample or a dmean for each group (where applicable). Also the MSE has been
upgraded/fixed and bee_moments now gives working first and second moments.