compute_rsquared function

Compute Bayesian R-squared from matrix of posterior replications