Details
type = "response" (default) returns observed Y minus the
response-scale fitted values, so residuals() reconciles with
fitted(). type = "working" returns the Gaussian working-scale
residuals the estimator actually minimised (identical to "response"
for the normal family).