时间序列模型c语言,R语言进行简单的时间序列分析

原标题：R语言进行简单的时间序列分析

经典的统计分析都假定数据序列具有独立性，而时间序列分析则侧重研究数据序列的互相依赖关系。后者实际上是对离散指标的随机过程的统计分析，所以又可看作是随机过程统计的一个组成部分。

今天我们就使用R语言自带的co2数据集进行简单的时间序列分析演示。

> View(co2)

> class(co2)

[1] "ts" ##说明为时间序列格式

> plot(co2)

> de_data

> plot(de_data)

建立模型之前，需要进行平稳性检验

> library(forecast)

> library(tseries)

> adf.test(co2)

Augmented Dickey-Fuller

Testdata: co2Dickey-Fuller = -2.8299, Lag order = 7, p-value = 0.2269alternative hypothesis: stationary ##可以看到P值是0.23，说明数据时稳定的。如果不稳定则需要进行差分。

ARIMA模型：

ARIMA：自回归积分滑动平均模型(Autoregressive Integrated Moving Average Model,简记ARIMA)

ARIMA(P,D,Q)

D：代表差分的阶数；差分就是上一节课讲到的diff()函数。一阶差分就是diff()一次，二阶差分就是对一阶差分的结果再做一次差分。

Q：代表滑动平均的阶数。

> acf(co2)

> pacf(co2)

> fit

> co2fit

> plot(co2,col="red")

> par(new=t)

> plot(co2fit)

责任编辑：