全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
3237 9
2013-12-04
在计算某一时间序列的自相关系数时,发现了如下问题,思而不解,来求教各位大侠,程序很简单:> 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)#只显示自相关的值,不显示图形。

图形如下:
Rplot03.png
自相关系数值如下:
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个,什么情况?
想不通,请明白的各位给在下解下惑吧~~

附件列表
QQ截图20131204101617.png

原图尺寸 5.92 KB

QQ截图20131204101617.png

二维码

扫码加我 拉你入群

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

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

全部回复
2013-12-4 11:41:32
那问题出现在pmi_tcdiff1上咯    pmi_tc这个结果在哪?   没结果数据,也看不出你哪步错啊
二维码

扫码加我 拉你入群

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

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

2013-12-4 13:17:24
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)#只显示自相关的值,不显示图形。
二维码

扫码加我 拉你入群

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

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

2013-12-4 16:36:56
目测这个好像是12列的数  表格里有多少数就显示多少个值
横坐标应该是2009到2011吧
二维码

扫码加我 拉你入群

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

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

2013-12-6 22:31:49
厉害呀
二维码

扫码加我 拉你入群

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

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

2013-12-7 09:10:24
acf(...,xaxt="n")
axis(1,at=seq(0,1.9,0.1),1:20)
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群