as a diagnostic method and model validation, the distribution
of the residuals was evaluated by using the ordinary residuals for the
model with identity link function and the deviance residuals for those
models with log link function, to look for heteroscedasticity (trends in
the degree of scatter of the residuals)