我认为这里有些奇怪。例如,以下代码为残差和创新提供了相同的值:
fit <- us_change %>%
model(ARIMA(Consumption ~ Income)) %>%
augment()
似乎该augment()
函数仅提取创新值并将其用于回归的残差。当我们使用 提取残差和创新时可以看到这一点residuals()
:
bind_rows(
`Regression Errors` = as_tibble(residuals(fit, type = "regression")),
`ARIMA Errors` = as_tibble(residuals(fit, type = "innovation")),
.id = "type"
)
那么残差和创新应该是不同的。