Я использую auto.arima для прогнозирования финансовых данных с использованием годовой сезонности и внешних регрессоров. Продолжительность сезонного периода составляет 251, потому что мы имеем дело с бизнесом, а не с календарными днями. Приведенный ниже код показывает, как я создаю свои объекты MSTS и ARIMA.
msts.f <- msts(daily.data$Total[1:895], seasonal.periods = c(251))
msts.f <- msts.f / 100000000
arima.f <- auto.arima(msts.f, xreg = as.matrix(daily.data[1:895,c(3,6)]), D = 1)
Проблема, с которой я столкнулся, связана с моими остатками. Значения намного меньше для первых 251 наблюдений (совпадающих с продолжительностью моего сезонного периода), чем остальные наблюдения.
autoplot(arima.f$residuals)
Любые идеи чего мне не хватает?