R语言 Holt-Winters法

有增长或者降低趋势并且存在季节性波动的时间序列的预测算法即Holt-Winters，这种序列可以被分解为水平趋势部分、季节波动部分，因此这两个因素应该在算法中有对应的参数来控制。

Holt-Winters算法中提供了alpha、beta和gamma 来分别对应当前点的水平、趋势部分和季节部分，参数的去执法范围都是0-1之间，并且参数接近0时，近期的观测值的影响权重就越小。我们以澳大利亚昆士兰州海滨纪念商品的月度销售日子为分析对象，绘制出该时间序列，捞个整体印象：

[plain] view plaincopy

souvenir <- scan("http://robjhyndman.com/tsdldata/data/fancy.dat")
souvenirtimeseries <- ts(souvenir, frequency=12, start=c(1987,1))
logsouvenirtimeseries <- log(souvenirtimeseries)
plot.ts(logsouvenirtimeseries)

我们采用R中提供的HoltWinters算法进行预测，实现和结果如下：

[plain] view plaincopy

souvenirtimeseriesforecasts <- HoltWinters(logsouvenirtimeseries)

Alpha=0.4,意味着当期预测基于平衡了最近和较远期的观测值。Beta 为0表明趋势部分的斜率在整个时间序列上市不变的，且等于初始值，这个也符合我们的直观感受，水平改变非常多，但是趋势部分斜率是基本不变的，于此相反gamma=0.96表明当期季节部分预测仅仅基于最近的观测值。

我们同时画出预测值和观测值来看下预测的效果如何：

[plain] view plaincopy

plot(souvenirtimeseriesforecasts)

可见Holt-Winters算法非常成功的预测了季节峰值，为了预测未来期数的值，我们同样采用forecast函数包，以预测未来48个月销售为例：

[plain] view plaincopy

library("forecast")
souvenirtimeseriesforecasts2 <- forecast.HoltWinters(souvenirtimeseriesforecasts, h=48)
plot.forecast(souvenirtimeseriesforecasts2)

图中蓝色线条显示的为预测值，深灰色部分为80%的置信区间，浅灰色为95的置信区间。

接下来老套路采用Ljung-Box 和画出直方图来检查预测误差的随机性，就不重复贴代码了，有忘记的童鞋请转前两篇，我这里就直接贴结果图了：

R语言 Holt-Winters法相关推荐

R语言贝叶斯方法在生态环境领域中的高阶技术
贝叶斯统计学即贝叶斯学派是一门基本思想与传统基于频率思想的统计学即频率学派完全不同的统计学方法,它在统计建模中具有灵活性和先进性特点,使其可以轻松应对复杂数据和模型结构. 然而,很多初学者在面对思想. ...
R语言：逆变换法生成随机数
逆变换法生成随机数: 一.概念解释 1.PDF 2.PMF 3.CDF 二.连续型情况举例三.离散型情况举例一.概念解释 1.PDF probability density function 概率 ...
R语言科学计数法详解：digits和scipen设置
控制R语言科学计算法显示有两个option: digitis和scipen.介绍的资料很少,而且有些是错误的.经过翻看R语言的帮助和做例子仔细琢磨,总结如下: 默认的设置是: getOption(&q ...
r语言三倍标准差法去除异常值，再计算平均值标准差
博主自己没能找到好的函数去除异常值,于是自己写好了一个简单实用的包.可以通过三倍标准差法删去每一行的异常值,然后计算出平均值标准差. 函数总共四个参数: file= 要计算的文件路径,在工作目录可以 ...
r语言插补法_R语言用多重插补法估算相对风险
在这里,我将用R中的一个小模拟示例进行说明.首先,我们使用X1和X2双变量法线和Y模拟大型数据集,其中Y遵循给定X1和X2的逻辑模型. 首先,我们模拟一个非常大的完整数据集: #simulate完整数 ...
R语言-随机前沿分析法--SFA
3.1介绍生产函数模型: lnqi=x'i*b+vi-ui (随机生产前沿函数) qi:产出变量向量 x'i:投入变量向量 b:变量参数估计 vi:统计噪声的对 ...
R语言-模糊逻辑控制
4.1模糊逻辑控制一种采用模糊集合理论的控制技术, 将模糊数学应用于控制系统,也是一种非线性的智能控制系统. 采用IF条件判断语句和THEN结果语句的方式呈现,其中条件判断语句是利用人的常识判断对控 ...
R语言简单指数平滑（自编函数）
系列文章目录 R语言利用贝叶斯分类法(klaR程序包),训练数据集,预测数据的分类 R语言载入arules程序包的Epub数据集,使用Apriori算法,进行关联规则分析 R语言进行系统聚类分析并作图 ...
R语言时间序列（time series）分析实战：霍尔特指数Holt‘s平滑法预测
R语言时间序列(time series)分析实战:霍尔特指数Holt's平滑法预测目录
时间序列预测之一：指数平滑法（二）R语言——代码实现
参考:https://www.cnblogs.com/fengzzi/p/10044426.html 指数模型是用来预测时序未来值的最常用模型.这类模型相对比较简单,但是实践证明它们的短期预测能力较好 ...

R语言 Holt-Winters法

R语言 Holt-Winters法相关推荐

最新文章

热门文章