powered by
Compute the pointwise deviance for all the observations in the sample
pointwise.deviance(mu, y, family = gaussian())