关于用R求EACF的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求各路大牛大神指教！
本人欲用R拟合一个ARMA模型，数据是1978到2010年的我国GDP取对数再做二阶差分。在求EACF时总是报错，如下：
In ar.ols(z, order.max = i, aic = FALSE, demean = FALSE, intercept = FALSE) :
  model order:  15 singularities in the computation of the projection matrix results are only valid up to model order 14

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：eacf ACF Computation Projection Intercept matrix 模型

pleae提供你的数据加代码lz你把AR部分和MA部分的最大阶数改一下就行了，不要用他的缺省值，略大。。。。我随便取了组，具体怎么取lz自己定哈

1. > library(TSA)
   
2. > GDP=ts(read.table("clipboard.txt"),start=1978)
   
3. > eacf(diff(log(GDP),difference=2),ar.max =3, ma.max =5)
   
4. AR/MA
   
5.   0 1 2 3 4 5
   
6. 0 o o o o o o
   
7. 1 o o o o o o
   
8. 2 x x o o x o
   
9. 3 x o o o o o

复制代码

回答完毕，希望对lz有用

求证1加1 发表于 2013-6-2 14:54
pleae提供你的数据加代码

数据：
3645.2
4062.6
4545.6
4891.6
5323.4
5962.7
7208.1
9016
10275.2
12058.6
15042.8
16992.3
18667.8
21781.5
26923.5
35333.9
48197.9
60793.7
71176.6
78973
84402.3
89677.1
99214.6
109655.2
120332.7
135822.8
159878.3
184937.4
216314.4
265810.3
314045.4
340902.8
401512.8
472881.6
代码：
library(TSA)
GDP=ts(read.table("clipboard"),start=1978)
eacf(diff(log(GDP),difference=2))

DongchongT 发表于 2013-6-2 15:06
数据：
3645.2
4062.6

已在沙发回复

求证1加1 发表于 2013-6-2 17:45
已在沙发回复

大神果然超神！
但这个是为什么呢

求证1加1 发表于 2013-6-2 17:45
已在沙发回复

我可不可以理解为 这是取对数后的二阶差分序列是白噪声造成的