简朴指数平滑(SES)

打印 上一主题 下一主题

主题 857|帖子 857|积分 2571

内容泉源
8.1 简朴的指数平滑 | 预测:方法与实践(第三版)

适用于预测没有显着趋势或季节因素的数据
指数平滑法是介于两个极端方法(质朴法和均匀法)之间的方法
质朴法

                                                                y                                  ^                                                      T                                  +                                  h                                  ∣                                  T                                                 =                                       y                               T                                            \hat{y}_{T+h|T}=y_T                     y^​T+h∣T​=yT​
即对将来的所有预测都等于该序列末了一个观测值
这可以看作是一个加权均匀值,其中所有的权重都被赋予了末了一个观测值。
均匀法

                                                                y                                  ^                                                      T                                  +                                  h                                  ∣                                  T                                                 =                                       1                               T                                                 ∑                                           t                                  =                                  1                                          T                                                 y                               t                                            \hat{y}_{T+h|T}=\frac{1}{T}\sum^T_{t=1} y_t                     y^​T+h∣T​=T1​t=1∑T​yt​
即所有汗青观测值具有雷同的权重
简朴指数平滑

                                                                y                                  ^                                                      T                                  +                                  1                                  ∣                                  T                                                 =                            α                                       y                               T                                      +                            α                            (                            1                            −                            α                            )                                       y                                           T                                  −                                  1                                                 +                            α                            (                            1                            −                            α                                       )                               2                                                 y                                           T                                  −                                  2                                                 +                            ⋯                                  \hat{y}_{T+1|T}=\alpha y_T+\alpha(1-\alpha)y_{T-1}+ \alpha(1-\alpha)^2y_{T-2}+\cdots                     y^​T+1∣T​=αyT​+α(1−α)yT−1​+α(1−α)2yT−2​+⋯
其中                                    0                         ≤                         α                         ≤                         1                              0\leq\alpha\leq1                  0≤α≤1 是平滑参数
即将更大的权重赋予到时间较近的观测值,权重随着时间指数下降,下降速度由                                    α                              \alpha                  α 控制
加权均匀形式

                                         {                                                                                                                              y                                                 ^                                                                               2                                                 ∣                                                 1                                                                          =                                           α                                                           y                                              1                                                          +                                           (                                           1                                           −                                           α                                           )                                                           ℓ                                              0                                                                                                                                                                                       y                                                 ^                                                                               t                                                 +                                                 1                                                 ∣                                                 t                                                                          =                                           α                                                           y                                              t                                                          +                                           (                                           1                                           −                                           α                                           )                                                                            y                                                 ^                                                                               t                                                 ∣                                                 t                                                 −                                                 1                                                                                                                             \begin{cases} \hat{y}_{2|1}=\alpha y_1+(1-\alpha)\ell_0\\ \hat{y}_{t+1|t}=\alpha y_t+(1-\alpha)\hat{y}_{t|t-1} \end{cases}                     {y^​2∣1​=αy1​+(1−α)ℓ0​y^​t+1∣t​=αyt​+(1−α)y^​t∣t−1​​

                                                                y                                  ^                                                      T                                  +                                  1                                  ∣                                  T                                                 =                                       ∑                                           i                                  =                                  0                                                      T                                  −                                  1                                                 α                            (                            1                            −                            α                                       )                               i                                                 y                                           T                                  −                                  i                                                 +                            (                            1                            −                            α                                       )                               T                                                 ℓ                               0                                            \hat{y}_{T+1|T}=\sum^{T-1}_{i=0} \alpha(1-\alpha)^iy_{T-i}+ (1-\alpha)^T\ell_0                     y^​T+1∣T​=i=0∑T−1​α(1−α)iyT−i​+(1−α)Tℓ0​
例:阿尔及利亚的出口

绘制时序图

  1. library(fpp3)
  2. algeria_economy <- global_economy |>
  3.   filter(Country == "Algeria")
  4. algeria_economy |>
  5.   autoplot(Exports) +
  6.   labs(y = "GDP百分比", title = "出口:阿尔及利亚", x = "年份")
复制代码

最小化残差平方和估计参数

  1. # 估计参数
  2. fit <- algeria_economy |>
  3.   #error("A")表示误差加在预测值上的
  4.   #trend("N")表示不存在趋势成分
  5.   #season("N")表示不存在季节成分
  6.   model(ETS(Exports ~ error("A") + trend("N") + season("N")))
  7. tidy(fit)
  8. # A tibble: 2 × 4
  9.   Country .model                                       term  estimate
  10.   <fct>   <chr>                                        <chr>    <dbl>
  11. 1 Algeria "ETS(Exports ~ error("A") + trend("N") … alpha    0.840
  12. 2 Algeria "ETS(Exports ~ error("A") + trend("N") … l[0]    39.5  
复制代码
预测

  1. fc <- fit |>
  2.   forecast(h = 5)
  3. fc |>
  4.   autoplot(algeria_economy) +
  5.   geom_line(aes(y = .fitted), col="#D55E00",
  6.             data = augment(fit)) +
  7.   labs(y="占GDP的百分比", x="年份",title="出口:阿尔及利亚") +
  8.   guides(color = "none"))
复制代码

橙色曲线表示拟合值

免责声明:如果侵犯了您的权益,请联系站长,我们会及时删除侵权内容,谢谢合作!更多信息从访问主页:qidao123.com:ToB企服之家,中国第一个企服评测及商务社交产业平台。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?立即注册

x
回复

使用道具 举报

0 个回复

倒序浏览

快速回复

您需要登录后才可以回帖 登录 or 立即注册

本版积分规则

何小豆儿在此

金牌会员
这个人很懒什么都没写!
快速回复 返回顶部 返回列表