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 回答