vlambda博客
学习文章列表

      前面一章，介绍了时间序列中涉及到的基本概念，本章将在此基础上介绍如何对时间序列的资料进行分析，怎么选择合适的模型。时间序列的分析方法，已经给大家整理出一个流程图，如下图。本章介绍时间序列经典分析法。

###书籍集为每年国王的寿命kingsts <- ts(kings) ##转成时间序列格式len <- length(kingsts)t <- seq(1,len)/10 ##产生有序变量kingsts.fit <- lm(kingsts~t) ##配合最小二乘法回归模型summary(kingsts.fit)  ##输出主要计算结果par(mfrow=c(3,1),mar=c(3,3,1.6,1),mgp=c(1.8,0.5,0))plot.ts(kingsts) ##绘制时间序列图plot(resid(kingsts.fit),cex.axis=0.8,cex.lab=0.8) ##画出残差图acf(resid(kingsts.fit),cex.axis=0.8,cex.lab=0.8)  ##做残差的acf图Box.test(kingsts.fit.fit$residuals) ##对残差做检验## p-value = 0.08575

kingsts.fit1 <- lm(kingsts~t+I(t^2)) ##二次拟合summary(kingsts.fit1) ##调整R2=0.2542summary(kingsts.fit)   ##调整R2=0.1419 Box.test(kingsts.fit1$residuals) ##p-value = 0.4331acf(kingsts.fit1$residuals) ##无相关性

      根据不同的历史间隔远近在预报中所起的作用不同而给予不同的权重，历史较近的观察值对当前值影响较大，应给予较大的权重；反之，历史较远的观察值对当前值的影响较小，应赋予较小的权重。该方法特别适合变化比较缓慢的时间序列资料的分析。

      简单指数匀滑法适合不含趋势和季节变动资料，Holt-Winters法则是将简单指数匀滑推广到含有趋势和季节变动的时间序列资料的分析。在R中都采用HoltWinter()函数分析，其具有下面的基本结构：

HoltWinters(x, alpha = NULL, beta = NULL, gamma = NULL, seasonal = c("additive", "multiplicative"), start.periods = 2, l.start = NULL, b.start = NULL, s.start = NULL, optim.start = c(alpha = 0.3, beta = 0.1, gamma = 0.1), optim.control = list())

       函数中的x为时间序列资料的名称，alpha为水平参数α，beta为趋势参数β，gamma为季节参数γ。当beta和gamma都赋值为FALSE时，模型就没有趋势分量和季节分量，也就是简单指数匀滑法。当gamma赋值为FALSE时，模型就没有季节分量。

        简单指数匀滑法

 空气中颗粒物PM10数据集为例 pm10.ts <- ts(pm_10$PM10,start = c(2007,1),frequency = 365) pm10.ts.hw <- HoltWinters(pm10.ts,gamma = F,beta = F) coef <- pm10.ts.hw$coefficients with(pm10.ts.hw,c(alpha,beta,gamma,coef,SSE))

       不含季节趋势指数匀滑

pm10.ts.hw.NS <- HoltWinters(pm10.ts,gamma = F)coef <- pm10.ts.hw.NS$coefficientswith(pm10.ts.hw.NS,c(alpha,beta,gamma,coef,SSE))

       包含长期趋势和季节趋势的Holt-Winters匀滑

pm10.ts.hw.S <- HoltWinters(pm10.ts)coef <- pm10.ts.hw.S$coefficients[1:2]with(pm10.ts.hw.S,c(alpha,beta,gamma,coef,SSE)) ##提取匀滑参数和当前基础水平a和趋势分量bpm10.ts.hw.S$coefficients[c(3:6,365+2)]  ##提取前4个和最后一个季节的匀滑参数s

     上述三个不同模型的结果解释：

     1.该序列没有明显的长期趋势，但可见明显的季节周期，这里为了和后面对比，因此姑且以此数据作为例子进行简单指数平滑。其水平参数α=0.73，说明某一天的pm10浓度受近期的影响比较大。其趋势和季节参数为0，SSE为0.57。

      2. 第二部分，仍以该例子，假使其没有季节分量gamma=F。看到β=0.02，说明其长期趋势不明显。SSE=0.59

      3.第三部分，考虑长期趋势和季节变动，看到α=0.61，表明近期污染水平对当前值影响大于远期污染水平。β=0表明趋势分量不需要随时间更新，γ=0.7，该序列存在明显季节周期。a=0.082是当前的基础水平，b=-1.1699e-06是当前的趋势，s1，s2等则是每天的季节指数。输出的SSE=0.52，误差低于前面两个模型，该模型构建比较合适。

    采用Holt-winters法的加性模型配合时间序列资料的图形见下图:

plot(pm10.ts.hw.S$fitted) ##加性模型拟合图

     图中xhat 为用模型参数估计值回代所得的期望值；level为水平值图，即稳定部分的估计值；trend为时间趋势图，结果未显示趋势，故为一条平行线；season为季节分量图，显示两个周期。2008-2009和2009到2010，预报未来一年的pm10含量图形见下图。

pm10.ts.hw.S.pre <- predict(pm10.ts.hw.S,365)ts.plot(pm10.ts,pm10.ts.hw.S.pre,lty=1:2)

上图中2007-2010年的实线部分为观察值，2010-2011年的虚线部分为预测值。