Fix for r-hat calculation on multidimensional nodes
Added comments to describe error handling in calc_rhat
Merge branch 'rhat_fix'