powered by
Compute the overall deviance averaging the contributions of all data
matrix.deviance(mu, y, family = gaussian())