3

この ARIMA モデルを考えると、次のようになります。

(1-0.8B)*(1-0.2B^6)*(1-B)Y_t = epsilon_t

乗法モデルは (1,1,0*(1,1,0)_6 (季節成分 = 6) です。このモデルから新しい値 (10 番目または 11 番目の値など) を予測するツールはありますか?次のような値のセット:

y <- c(1,4,5,2,0,8,9,4,-3,-3)

私は試した

arima(y,order=c(1,1,0),seasonal=list(order=c(1,1,0),period=6))

 error: initial value in 'vmmin' is not finite
4

1 に答える 1

2

predict()次の関数で事前に予測できます。

> y=c(1,4,5,2,0,8,9,4,-3,-3)
> mymodel = arima(c(1,4,5,2,0,8,9,4,-3,-3) ,order=c(1,1,0),seasonal=list(order=c(1,1,0), period=2))
> mymodel

Call:
arima(x = c(1, 4, 5, 2, 0, 8, 9, 4, -3, -3), order = c(1, 1, 0), seasonal = list(order = c(1, 
    1, 0), period = 2))

Coefficients:
         ar1     sar1
      0.7368  -0.9169
s.e.  0.3696   0.1089

sigma^2 estimated as 11.25:  log likelihood = -20.23,  aic = 46.46


> predict(mymodel, n.ahead = 5)
$pred
Time Series:
Start = 11 
End = 15 
Frequency = 1 
[1]  -7.763438 -16.104376 -25.686464 -28.419524 -35.086436

$se
Time Series:
Start = 11 
End = 15 
Frequency = 1 
[1]  3.354151  6.722215 10.392430 14.061929 19.640317

モデルのデータベクトルが十分に長くなるように、期間を短縮しました。

于 2013-03-16T19:23:29.140 に答える