时间序列分析R语言程序
#例2.1 绘制1964——1999年中国年纱产量序列时序图（数据见附录1.2） Data1.2=read.csv(\s\\\\Administrator\\\\Desktop\\\\附录1.2.csv\如果有标题，用T；没有标题用F
plot(Data1.2,type='o') #例2.1续 tdat1.2=Data1.2[,2] a1.2=acf(tdat1.2)
#例2.2绘制1962年1月至1975年12月平均每头奶牛产奶量序列时序图（数据见附录1.3）
Data1.3=read.csv(\s\\\\Administrator\\\\Desktop\\\\附录1.3.csv\tdat1.3=as.vector(t(as.matrix(Data1.3)))[1:168]#矩阵转置转向量
plot(tdat1.3,type='l') #例2.2续
acf(tdat1.3) #把字去掉 pacf(tdat1.3)
#例2.3绘制1949——1998年北京市每年最高气温序列时序图
Data1.4=read.csv(\s\\\\Administrator\\\\Desktop\\\\附录1.4.csv\plot(Data1.4,type='o')
##不会定义坐标轴 #例2.3续 tdat1.4=Data1.4[,2] a1.4=acf(tdat1.4) #例2.3续
Box.test(tdat1.4,type=\g-Box\
Box.test(tdat1.4,type=\g-Box\
#例2.4随机产生1000个服从标准正态分布的白噪声序
列观察值，并绘制时序图 Data2.4=rnorm(1000,0,1) Data2.4
plot(Data2.4,type='l') #例2.4续 a2.4=acf(Data2.4) #例2.4续
Box.test(Data2.4,type=\ng-Box\
Box.test(Data2.4,type=\ng-Box\
#例2.5对1950——1998年北京市城乡居民定期储蓄所占比例序列的平稳性与纯随机性进行检验
Data1.5=read.csv(\s\\\\Administrator\\\\Desktop\\\\附录1.5.csv\plot(Data1.5,type='o',xlim=c(1950,2010),ylim=c(60,100))
tdat1.5=Data1.5[,2] a1.5=acf(tdat1.5)
#白噪声检验
Box.test(tdat1.5,type=\g-Box\
Box.test(tdat1.5,type=\g-Box\
#例2.5续选择合适的ARMA模型拟合序列 acf(tdat1.5) pacf(tdat1.5)
#根据自相关系数图和偏自相关系数图可以判断为AR（1）模型
#例2.5续 P81 口径的求法在文档上 #P83
arima(tdat1.5,order=c(1,0,0),method=\极大似然估计
ar1=arima(tdat1.5,order=c(1,0,0),method=\summary(ar1) ev=ar1$residuals acf(ev) pacf(ev)
时间序列分析R语言程序
