Under the description of the 'type' argument, ?residuals.gls says
'Defaults to ?"pearson"?.'
But residuals.gls starts
function(object, type = c("response", "pearson", "normalized"), ...)
type <- match.arg(type)
which sure looks to me like it defaults to "response", not "pearson"
(and it behaves that way in tests).
It would seem to make more sense to change the documentation rather than
since anyone who looked at the docs would have been confused already,
whereas someone who had
been happily using the code without looking at the docs would see a
sudden change in the results ...
This is in nlme 3.1-96, from a fresh tools/rsync-recommended. Sending it
to r-devel for comment because r-core is listed as the maintainer.