协整检验r语言代码_R语言学习笔记（一）：时间序列分析

ARIMA(p，d，q)模型全称为差分自回归移动平均模型(Autoregressive Integrated Moving Average Model,简记ARIMA)，AR是自回归， p为自回归项； MA为移动平均，q为移动平均项数，d为时间序列成为平稳时所做的差分次数。

所谓ARIMA模型，是指将非平稳时间序列转化为平稳时间序列，然后将因变量仅对它的滞后值以及随机误差项的现值和滞后值进行回归所建立的模型。ARIMA模型根据原序列是否平稳以及回归中所含部分的不同，包括移动平均过程(MA)、自回归过程(AR)、自回归移动平均过程(ARMA)以及ARIMA过程。

*什么时候要取对数：原序列的前后数值相差较大或者数量级相差较大时，一般取对数将指数趋势转化为线性趋势，而且还不会改变变量之间的统计性质，同时得到较平稳的序列；可以消除异方差，让方差恒定，即让波动相对稳定。

http://www.zhihu.com/question/22012482

以数据集robot为例

1.打开数据并转换格式

data(robot)

x

ro

2.模型拟合，参数估计

> air.fit

> m1=arima(x = ro, order = c(1, 0, 0))> m2=arima(x = ro, order = c(0, 1, 1))

3.绘制自相关图和偏自相关图

> acf(ro,80)> pacf(ro,30)

4.导出R的数据

> write.table(robot,"robot.csv",sep=",")

5.绘制时间序列图，一次差分时间序列图，一次差分与季节差分时间序列图

> plot.ts(ro)

> plot.ts(diff(aplog)

> aplog1

> aplog222

6.绘制标准残差图，绘制残差自相关图，