Residuals of models fitted with functions betabin
and negbin
(formal class “glimML”), or
quasibin
and quasipois
(formal class “glimQL”).
1 2 3 4 5 |
object |
Fitted model of formal class “glimML” or “glimQL”. |
type |
Character string for the type of residual: “pearson” (default) or “response”. |
... |
Further arguments to be passed to the function, such as |
For models fitted with betabin
or quasibin
, Pearson's residuals are computed as:
(y - n * p.fit) / (n * p.fit * (1 - p.fit) * (1 + (n - 1) * φ))^{0.5}
where y and n are respectively the numerator and the denominator of the response, p.fit
is the fitted probability and φ is the fitted overdispersion parameter. When n = 0, the
residual is set to 0. Response residuals are computed as y/n - p.fit.
For models fitted with negbin
or quasipois
, Pearson's residuals are computed as:
(y - y.fit) / (y.fit + φ * y.fit^2)^{0.5}
where y and y.fit are the observed and fitted counts, respectively. Response residuals are computed as y - y.fit.
A numeric vector of residuals.
Matthieu Lesnoff matthieu.lesnoff@cirad.fr, Renaud Lancelot renaud.lancelot@cirad.fr
residuals.glm
1 2 3 4 5 6 |
Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.