Как работает predic (), если параметр newdata отсутствует?

Question

Недавно я узнал о моделях препятствий. Я нашел блог с этими кодами.

library(AER)
data("NMES1988")
nmes <- NMES1988[, c(1, 6:8, 13, 15, 18)]
plot(table(nmes$visits))
mod1 <- glm(visits ~. , data = nmes, family = "poisson")

mu <- predict(mod1, type = "response")

exp <- sum(dpois(x=0, lambda = mu))

Параметр newdata отсутствует в функции predic (). Будет ли эта функция делать прогноз на основе старых nmes данных?

Answer 1

Да, он использует исходные данные, т.е. возвращает fitted.values (yhat). Проверить:

mu <- predict(mod1, type = "response")
mu2 <- predict(mod1, type = "response", newdata=nmes)
mu3 <- mod1$fitted.values

identical(mu, mu2)
# [1] TRUE
identical(mu, mu3)
# [1] TRUE