Вызов остатка (модель) по умолчанию будет соответствовать остаткам отклонения, тогда как модель $ остаток даст вам рабочие остатки. Из-за функции связи нет единого определения того, что является остатком модели. Существуют отклонения, рабочие, частичные, остатки Пирсона и отклика. Поскольку они зависят только от средней структуры (а не от дисперсии), остатки для квазипуассона и пуассона имеют одинаковую форму. Вы можете взглянуть на функцию residuals.glm
для деталей, но вот пример:
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
treatment <- gl(3,3)
glm.D93 <- glm(counts ~ outcome + treatment, family=quasipoisson())
glm.D93$resid
#working
resid(glm.D93,type="working")
(counts - glm.D93$fitted.values)/exp(glm.D93$linear)
#deviance
resid(glm.D93,type="dev")
fit <- exp(glm.D93$linear)
poisson.dev <- function (y, mu)
sqrt(2 * (y * log(ifelse(y == 0, 1, y/mu)) - (y - mu)))
poisson.dev(counts,fit) * ifelse(counts > fit,1,-1)
#response
resid(glm.D93,type="resp")
counts - fit
#pearson
resid(glm.D93,type="pear")
(counts - fit)/sqrt(fit)