原文链接:http://tecdat.cn/?p=20434

本文从实际角度探讨了季节性单位根。咱们思考一些工夫序列 ),例如路线上的交通流量,

 > plot(T,X,type="l")> reg=lm(X~T)> abline(reg,col="red")

如果存在趋势,咱们应该将其删除,而后解决残差 

> Y=residuals(reg)> acf(Y,lag=36,lwd=3)

咱们能够看到这里有一些季节性。第一个策略可能是假如存在季节性单位根,因而咱们思考 Y_t),咱们尝试找到ARMA模型。思考工夫序列的自相干函数,

> Z=diff(Y,12)> acf(Z,lag=36,lwd=3)

或偏自相干函数

第一个图可能倡议MA(1),而第二个图可能倡议AR(1)工夫序列。咱们都尝试。

 arimaCoefficients:          ma1  intercept      -0.2367  -583.7761s.e.   0.0916   254.8805sigma^2 estimated as 8071255:  log likelihood = -684.1,  aic = 1374.2 

能够认为是白噪声(如果您不确定,请尝试 Box-Pierce或Ljung-Box 测试)。

 arimaCoefficients:          ar1  intercept      -0.3214  -583.0943s.e.   0.1112   248.8735sigma^2 estimated as 7842043:  log likelihood = -683.07,  aic = 1372.15 

也能够视为白噪声。到目前为止,咱们有

对于一些白噪声 )。这表明以下的SARIMA构造 )

 arimaCoefficients:          ar1      -0.2715s.e.   0.1130sigma^2 estimated as 8412999:  log likelihood = -685.62,  aic = 1375.25

当初,如果咱们认为咱们没有季节性单位根,而在AR构造中只是一个大的自回归系数。让咱们尝试相似

自然而然的猜想是该系数应该(可能)靠近于1。让咱们尝试一下

 arimaCoefficients:          ar1    sar1  intercept      -0.1629  0.9741  -684.9455s.e.   0.1170  0.0115  3064.4040sigma^2 estimated as 8406080:  log likelihood = -816.11,  aic = 1640.21

这与咱们先前(以某种形式)取得的后果具备可比性,因而咱们能够假如该模型是一个乏味的模型。咱们将进一步探讨:第一个系数可能是不重要的。

这两个模型有什么区别?

从(十分)长期的角度来看,模型是齐全不同的:一个模型是安稳的,因而预测将趋向于平均值,而另一个模型则是按节令的,因而置信区间将减少。咱们失去

> pre(model2,600,b=60000)

对于安稳的

> prev(model3,600,b=60000)

然而,应用这些模型进行的预测仅实用于短期范畴。在这种状况下,这里的预测简直雷同,

> pre(model2,36,b=60000)

> pre(model3,36,b=60000)

当初,如果咱们回到第二个模型,自回归系数可能被认为是不重要的。如果咱们将其删除怎么样?

 Call:seasonal = list(order = c(1, 0, 0)Coefficients:        sar1  intercept      0.9662  -696.5661s.e.  0.0134  3182.3017sigma^2 estimated as 8918630:  log likelihood = -817.03,  aic = 1640.07

如果咱们看一下(短期)预测,咱们失去

> pre(model,36,b=32000)

有什么区别吗?如果咱们看一下预测后果数字,咱们会失去

数字不同,但差别不大(请留神置信区间的大小)。这能够解释为什么在R中,当咱们在自回归过程时 ,失去一个模型要预计的参数,即便其中不重要,咱们通常也会保留它们来预测。


最受欢迎的见解

1.用于NLP的Python:应用Keras的多标签文本LSTM神经网络分类

2.Python中利用长短期记忆模型LSTM进行工夫序列预测剖析 – 预测电力耗费数据

3.python在Keras中应用LSTM解决序列问题

4.Python中用PyTorch机器学习分类预测银行客户散失模型

5.R语言多元Copula GARCH 模型工夫序列预测

6.在r语言中应用GAM(狭义相加模型)进行电力负荷工夫序列剖析

7.R语言中ARMA,ARIMA(Box-Jenkins),SARIMA和ARIMAX模型用于预测工夫序列数

8.R语言预计时变VAR模型工夫序列的实证钻研剖析案例

9.用狭义加性模型GAM进行工夫序列剖析