rhat function

Compute split Rhat statistic