关于序列自相关结果的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在计算某一时间序列的自相关系数时，发现了如下问题，思而不解，来求教各位大侠，程序很简单：> pmi_tcdiff1<-diff(pmi_tc,differences=1)#进行一阶差分
> plot.ts(pmi_tcdiff1)
> acf(pmi_tcdiff1,lag.max=20)#画自相关图，最大滞后阶数20
> acf(pmi_tcdiff1,lag.max=20,plot=FALSE)#只显示自相关的值，不显示图形。

图形如下：

自相关系数值如下：
Autocorrelations of series ‘pmi_tcdiff1’, by lag

0.0000 0.0833 0.1667 0.2500 0.3333 0.4167 0.5000 0.5833 0.6667 0.7500 0.8333 0.9167 1.0000 1.0833 1.1667
1.000  0.052 -0.194 -0.090 -0.073  0.024  0.070  0.002 -0.095 -0.116  0.190  0.164 -0.275 -0.083 -0.057
1.2500 1.3333 1.4167 1.5000 1.5833 1.6667
0.007  0.245  0.254 -0.088 -0.153 -0.090

问题1，图形横坐标不是应该是1,2,3，……20，吗？怎么出现小数的情况，也没有达到20阶。
问题2，自相关系数值应该有20个，可是结果有42个，什么情况？
想不通，请明白的各位给在下解下惑吧~~

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：序列自相关 自相关 correlations correlation relations 小白

那问题出现在pmi_tcdiff1上咯    pmi_tc这个结果在哪？   没结果数据，也看不出你哪步错啊

yywan0913 发表于 2013-12-4 11:41
那问题出现在pmi_tcdiff1上咯    pmi_tc这个结果在哪？   没结果数据，也看不出你哪步错啊

这是全部程序，哪儿有问题呢？
data01<-read.csv("pmi09.csv")
pmits<-ts(data01,frequency=12,start=c(2009,3))
pmitscom<-decompose(pmits)#进行季节调整
pmitscom$seasonal#季节调整因素
pmitsad<-pmits-pmitscom$seasonal#剔除季节变化的影响
pmi_tc<-pmitsad#将剔除季节因素的变量赋给一个新的变量值，相当于重命名
pmi_tcdiff1<-diff(pmi_tc,differences=1)#进行一阶差分
plot.ts(pmi_tcdiff1)
acf(pmi_tcdiff1,lag.max=20)#画自相关图，最大滞后阶数20
acf(pmi_tcdiff1,lag.max=20,plot=FALSE)#只显示自相关的值，不显示图形。

目测这个好像是12列的数  表格里有多少数就显示多少个值
横坐标应该是2009到2011吧

厉害呀

acf(...,xaxt="n")
axis(1,at=seq(0,1.9,0.1),1:20)

nuomin 发表于 2013-12-7 09:10
acf(...,xaxt="n")
axis(1,at=seq(0,1.9,0.1),1:20)

谢谢您的回答。按照您的代码加上以后，图像横坐标正常了，但是自相关值还是老样子，看结果：
Autocorrelations of series ‘pmi_tcdiff1’, by lag

0.0000 0.0833 0.1667 0.2500 0.3333 0.4167 0.5000 0.5833 0.6667 0.7500 0.8333
1.000  0.052 -0.194 -0.090 -0.073  0.024  0.070  0.002 -0.095 -0.116  0.190
0.9167 1.0000 1.0833 1.1667 1.2500 1.3333 1.4167 1.5000 1.5833 1.6667
0.164 -0.275 -0.083 -0.057  0.007  0.245  0.254 -0.088 -0.153 -0.090
不明白哪里不对了，怎么是这样的显示结果？

1. acf_out<-acf(pmi_tcdiff1,lag.max=20,plot=FALSE)
   
2. str(acf_out) #内部结构
   
3. acf_out$acf #这是自相关函数值

复制代码

显示结果是上下为一对

0,1/12,2/12,...,12/12,...,20/12
acf(0),acf(1),acf(2),...,acf(12),...,acf(20)

nuomin 发表于 2013-12-9 19:24
显示结果是上下为一对

0,1/12,2/12,...,12/12,...,20/12

问题解决，多谢！R显示为什么不直观点呢，对初学者来说，不太容易理解。

skytreee 发表于 2013-12-10 08:58
问题解决，多谢！R显示为什么不直观点呢，对初学者来说，不太容易理解。

觉得别扭，才会有动力去改变他