Я сейчас занимаюсь моделированием продаж автомобилей в США.Моя модель аддитивная модель.Я пришел к выводу, что лучшая модель аримы - это арима (2,1,3).Тем не менее, я хочу добавить сезонность в модель арима.Итак, я попытался включить сезонные чайники с xreg.Это не работает.Я только выдаю мне сообщение об ошибке.Есть ли другой способ добавить сезонность?
Я пытался увидеть, работает ли модель, когда я не различаю ее, что она делает.
s=12
st=c(1976,1)
styear=1976
data=ts(data,frequency = s,start = st)
y=data[1:504]
T=length(y)
y=ts(y,frequency = s,start=st)
Created my seasonal dummies
S2 = rep(c(0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), T/s)
S3 = rep(c(0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0), T/s)
S4 = rep(c(0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0), T/s)
S5 = rep(c(0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0), T/s)
S6 = rep(c(0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0), T/s)
S7 = rep(c(0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0), T/s)
S8 = rep(c(0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0), T/s)
S9 = rep(c(0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0), T/s)
S10 = rep(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0), T/s)
S11 = rep(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0), T/s)
S12 = rep(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1), T/s)
TrSeas = model.matrix(~ t+S2+S3+S4+S5+S6+S7+S8+S9+S10+S11+S12)
TrSeas
arima213=Arima(y,order = c(2,1,3),xreg = TrSeas,include.mean = FALSE,method = "ML")
Я получаю следующее сообщение: Ошибка в optim (init [маска], armafn, method = optim.method, hessian = TRUE,: неконечное значение, предоставляемое optim