Add compute_variance parameter to slow version of stats function.