次の形式のARIMA(1,1,0)モデルを使用して近似したい時系列あるとします。
これは次のように書き直すことができます。
最後の方程式は、係数がおよび AR(2)モデルを表しています。によっては、このAR(2)モデルが非定常である可能性があることを認識しています。ただし、最初に差分を取っていた場合、モデリングしているシリーズは静止していてはいけません。
モデルが定常的でない場合は、差分を使用する必要があることを知っています。しかし、AR(2)モデルとARIMA(1,1,0)モデルを使用した場合、結果はどのように異なりますか?私は(Rが示唆するように)収束に問題があると思います。ただし、Rに近似を実行するように依頼すると、Rは両方を実行し、係数は(ほとんど)上記の私の観察と一致します。ただし、予測は明らかに異なります。
誰かがこれに光を当てたり、私に良い参照を指摘したりできれば、私はそれを感謝します。
これは、両方のモデルを生成するために使用したRコードです。
> set.seed(2)
> x <- arima.sim(n = 1000, model=list(order=c(1,1,0), ar=c(0.3)))
> plot(x)
> arima(x, order=c(1,1,0))
Call:
arima(x = x, order = c(1, 1, 0))
Coefficients:
ar1
0.3291
s.e. 0.0298
sigma^2 estimated as 1.03: log likelihood = -1433.91, aic = 2871.81
> arima(x, order=c(2,0,0))
Call:
arima(x = x, order = c(2, 0, 0))
Coefficients:
ar1 ar2 intercept
1.3290 -0.3294 50.9803
s.e. 0.0298 0.0299 35.9741
sigma^2 estimated as 1.03: log likelihood = -1438.93, aic = 2885.86
Warning messages:
1: In log(s2) : NaNs produced
2: In log(s2) : NaNs produced
3: In log(s2) : NaNs produced
4: In arima(x, order = c(2, 0, 0)) :
possible convergence problem: optim gave code = 1