[求助]请教一下如何用R来看是什么ARIMA...紧急！！谢谢

zihau87

5567

收藏 2009-03-25

请好心人+高手帮忙指教一下。。。。
功课上需要到。。。。
谢谢哦~~~
麻烦请教教如何用 R （或 eview 3.1 ）来看是什么 arima
谢谢^^

附上一个附件。。。。

原版的 acf & pacf

经过了 First Difference 后的 acf & pacf

可是却不会看是 ARIMA (p,d,q)

能教教我吗？？
谢谢 ^^

或者请教我用 eview 3.1 来看也可以。。。谢谢。。。

[此贴子已经被作者于2009-3-25 1:16:29编辑过]

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

ruiqwy

2009-3-25 09:12:00

从你的一阶差分图来看，似乎还不平稳，先做个平稳性检验，等平稳了再考虑ARIMA！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zihau87

2009-3-25 11:41:00

请问 ruiqwy 大大说的平稳性检验是这个？

Box.test(z,lag=24)
Box.test(z,lag=24,type="Ljung")

z 是已经一阶差分了

可是答案的确都小过 "0.10", 所以要二阶差分咯？？

附加小问题：这种做法里，一定是用＝"24" ？？

谢谢ruiqwy 大大
谢谢＾＾

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ruiqwy

2009-3-25 13:12:00

Box.test(z,lag=24)
Box.test(z,lag=24,type="Ljung")

这些是做序列自相关检验，平稳性检验可以用ADF检验、KPSS检验等！！你先找本时间序列的书看看

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

vrooadk

2009-3-25 22:21:00

this does not seem right to me. Doesn't LZ need to fit an arima first and then perform all those tests to the residuals? The ACF and PACF of the 1st differenced series look ok to me. Go ahead fitting an arima model using auto.arima in forecast package like auto.arima(z, d=1,D=0,max.p=5,max.q=3,max.P=0,max.Q=0,ic="aic") and then carry those tests to the residuals of this model.

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zihau87

2009-3-26 03:32:00

回复：（ruiqwy）Box.test(z,lag=24)Box.test(z,lag=2...

哦哦。。。。那么我还真没学到 ADF 耶。。。
毕竟上课只学到基础

那么又好奇。。。我用eviews 做出 adf test
那么是不是要大过那些 % 的才算稳定叻？？


        t-Statistic    Prob.*

Augmented Dickey-Fuller test statistic            -4.106224    0.0014
Test critical values:                      1% level        -3.489659
   5% level        -2.887425
                                                 10% level        -2.580651

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

zihau87

2009-3-26 04:18:00

thank for your suggestion about using
auto.arima in forecast package like auto.arima(z, d=1,D=0,max.p=5,max.q=3,max.P=0,max.Q=0,ic="aic")

and because of using new data , my acf(z) & pacf(z) will be change

and.. can u pls help me hv a look on my new graph
the auto.arima is (2,1,0)
but y isnt (2,1,3) ???

thanks a lot ^^

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

vrooadk

2009-3-27 04:58:00

oh, sorry. did not look very carefully. it should be auto.arima(y, d=1,D=0,max.p=5,max.q=3,max.P=0,max.Q=0,ic="aic"). run this again.

fit<-auto.arima(y, d=1,D=0,max.p=5,max.q=3,max.P=0,max.Q=0,ic="aic")
tsdiag(fit)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ruiqwy

2009-3-27 09:45:00

我觉得上面这样做是有问题，对time series做ARIMA，首先要判断几阶差分平稳，只有平稳后才能做ARMA。否则毫无意义。

auto.arima(y, d=1,D=0,max.p=5,max.q=3,max.P=0,max.Q=0,ic="aic")

这个函数是用来识别最优的ARIMA，但是你把d=1，那就告诉他是一阶差分平稳。你做过了一阶差分平稳检验了吗？

d 的参数如下：

d ：Order of first-differencing. If missing, will choose a value based on KPSS test.

如果d值不设为1，auto.arima会用KPSS去检验平稳性。也就是我说的，先用ADF、KPSS检验平稳后，才可以去做ARMA模型！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zihau87

2009-3-27 10:05:00

谢谢两位......

然后我试了

data <-read.table
y<-ts(data)
z<-ts(diff(y,lag=1))

然后用 KPSS test 时
ndiffs(y,alpha=0.05)   答案是 1
ndiffs(z,alpha=0.05)    答案是 0

所以表示我的 z （一阶差分）是已经平稳了吗？？

这时候,假如我做
auto.arima(y, d=1,D=0,max.p=5,max.q=3,max.P=0,max.Q=0,ic="aic")
答案是出现 (0,1,2)     的确是最小的 AIC

可是又好奇......从第7楼的 ACF & PACF 图来看...
不是应该是 (2,1,0) 或 (2,1,3) 吗？？？

************
笨笨的疑问，为了确实一些事......
请问在算 AIC 时

AIC<-matrix(0,6,6)

for (p in 0:5)

for (q in 0:5)

{

mod.fit<-arima(y,order=c(p,1,q))

AIC[p+1,q+1]<-mod.fit$aic

}

AIC

上面的arima (这个应该用y，还是z呢？？)

[此贴子已经被作者于2009-3-27 10:07:26编辑过]

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

vrooadk

2009-3-31 03:46:00

agreed with fl9. the first order differenced series look fine to me. so I suggest doing auto.arima. maybe visually examining the stationarity is not good enough to ruiqwy.

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

回复：（ruiqwy）Box.test(z,lag=24)Box.test(z,lag=2...

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群