时间序列预测专项——SARIMAX介绍
时间序列预测专项——SARIMAX介绍
- 一、延迟算子
- 二、AR model
- 三、MA model
- 四、ARMA model
- 五、ARIMA model
- 六、加法季节ARIMA model
- 七、乘法季节ARIMA model
- 参考文档
本文主要介绍如何用延迟算子(Backshift operator ,记为B)或滞后算子(Lag Operator,记为L)表示 ARIMA model。
一、延迟算子
1、一阶差分
Δ y t = y t − y t − 1 = ( 1 − B ) y t \Delta{y_t}=y_t-y_{t-1}=(1-B)y_t Δyt=yt−yt−1=(1−B)yt
2、二阶差分
Δ 2 y t = ( 1 − B ) 2 y t = ( y t − y t − 1 ) − ( y t − 1 − y t − 2 ) \Delta^{2}{y_t}=(1-B)^2y_t=(y_t-y_{t-1})-(y_{t-1}-y_{t-2}) Δ2yt=(1−B)2yt=(yt−yt−1)−(yt−1−yt−2)
3、n阶差分
Δ n y t = ( 1 − B ) n y t \Delta^{n}{y_t}=(1-B)^ny_t Δnyt=(1−B)nyt
二、AR model
1、AR(1)模型,如下:
y t = ϕ y t − 1 + ε t y_t=\phi y_{t-1}+\varepsilon_t yt=ϕyt−1+εt
用延迟算子表示AR(1)模型,如下:
( 1 − ϕ B ) y t = ε t (1-\phi B)y_t=\varepsilon_t (1−ϕB)yt=εt
2、AR( p )模型,如下:
y t = ϕ 1 y t − 1 + ϕ 2 y t − 2 + . . . + ϕ p y t − p + ε t y_t=\phi_1y_{t-1}+\phi_2y_{t-2}+...+\phi_py_{t-p}+\varepsilon_t yt=ϕ1yt−1+ϕ2yt−2+...+ϕpyt−p+εt
用延迟算子表示AR( p )模型,如下:
( 1 − ϕ 1 B − ϕ 2 B 2 − . . . − ϕ p B p ) y t = ε t (1-\phi_1B-\phi_2B^2-...-\phi_pB^p)y_t=\varepsilon_t (1−ϕ1B−ϕ2B2−...−ϕpBp)yt=εt
三、MA model
MA(q)模型,如下:
y t = ε t + θ 1 ε t − 1 + θ 2 ε t − 2 + . . . + θ q ε t − q y_t=\varepsilon_t+\theta_1\varepsilon_{t-1}+\theta_2\varepsilon_{t-2}+...+\theta_q\varepsilon_{t-q} yt=εt+θ1εt−1+θ2εt−2+...+θqεt−q
用延迟算子表示MA(q)模型,如下:
y t = ( 1 + θ 1 B + θ 2 B 2 + . . . + θ q B q ) ε t y_t=(1+\theta_1B+\theta_2B^2+...+\theta_qB^q)\varepsilon_t yt=(1+θ1B+θ2B2+...+θqBq)εt
四、ARMA model
令 ϕ ( B ) = 1 − ϕ 1 B − ϕ 2 B 2 − . . . − ϕ p B p , θ ( B ) = 1 + θ 1 B + θ 2 B 2 + . . . + θ q B q \phi(B)=1-\phi_1B-\phi_2B^2-...-\phi_pB^p,\theta(B)=1+\theta_1B+\theta_2B^2+...+\theta_qB^q ϕ(B)=1−ϕ1B−ϕ2B2−...−ϕpBp,θ(B)=1+θ1B+θ2B2+...+θqBq,于是ARMA(p,q)模型可以表示为:
ϕ ( B ) y t = θ ( B ) ε t \phi(B)y_t=\theta(B)\varepsilon_t ϕ(B)yt=θ(B)εt
五、ARIMA model
ARIMA(p,d,q)模型可以表示为:
ϕ ( B ) y t ∗ = θ ( B ) ε t \phi(B)y_t^{*}=\theta(B)\varepsilon_t ϕ(B)yt∗=θ(B)εt
其中 y t ∗ = Δ d y t = ( 1 − B ) d y t y_t^{*}=\Delta^{d}y_t=(1-B)^{d}y_t yt∗=Δdyt=(1−B)dyt.
六、加法季节ARIMA model
考虑加法季节效应时,令$p=1,d=1,q=(1,0,0,1),模型可表示为:
Δ y t = ϕ Δ y t − 1 + θ 1 ε t − 1 + θ 2 ε t − 4 \Delta y_t=\phi \Delta y_{t-1}+\theta_1\varepsilon_{t-1}+\theta_2 \varepsilon_{t-4} Δyt=ϕΔyt−1+θ1εt−1+θ2εt−4
七、乘法季节ARIMA model
1、乘法季节模型SARIMA(p,d,q) × \times ×(P,D,Q,s)
ϕ p ( B ) ϕ ~ P ( B s ) y t ∗ = θ q ( B ) θ ~ Q ( B s ) ε t \phi_p(B)\widetilde{\phi}_P(B^s)y_t^{*}=\theta_q(B)\widetilde{\theta}_Q(B^s)\varepsilon_t ϕp(B)ϕ P(Bs)yt∗=θq(B)θ Q(Bs)εt
其中, y t ∗ = A ( t ) + Δ d Δ s D y t = ( 1 − B ) d ( 1 − B s ) D y t y_t^{*}=A(t)+\Delta^d\Delta_s^{D}y_t=(1-B)^d(1-B^s)^Dy_t yt∗=A(t)+ΔdΔsDyt=(1−B)d(1−Bs)Dyt.
- ϕ p ( B ) \phi_p(B) ϕp(B):非季节性自回归滞后多项式
- ϕ ~ P ( B s ) \widetilde{\phi}_P(B^s) ϕ P(Bs):季节性自回归滞后多项式
- θ q ( B ) \theta_q(B) θq(B):非季节性移动平均滞后多项式
- θ ~ Q ( B s ) \widetilde{\theta}_Q(B^s) θ Q(Bs):季节性移动平均滞后多项式
- A ( t ) A(t) A(t):趋势多项式,可以是常数
2、特别地,SARIMA(1,1,1) × \times ×(1,1,0,12)可以表示为:
( 1 − ϕ B ) ( 1 − ϕ ~ B 12 ) y t ∗ = ( 1 − θ B ) ε t (1-\phi B)(1-\widetilde{\phi}B^{12})y_t^{*}=(1-\theta B)\varepsilon_t (1−ϕB)(1−ϕ B12)yt∗=(1−θB)εt
其中, y t ∗ = Δ Δ 12 y t = ( 1 − B ) ( 1 − B 12 ) y t y_t^{*}=\Delta\Delta_{12}y_t=(1-B)(1-B^{12})y_t yt∗=ΔΔ12yt=(1−B)(1−B12)yt.
参考文档
[1] http://devdoc.net/python/statsmodels-0.8.0/examples/notebooks/generated/statespace_sarimax_stata.html
时间序列预测专项——SARIMAX介绍相关推荐
- 时间序列预测专项——基于Prophet的业务预测
时间序列预测专项 基于Prophet的业务预测 1. 趋势模型 1.1 饱和增长模型 1.2 分段线性模型 1.3 变点自动选择 2. 季节模型 3. 节假日模型 4. Prophet模型参数解读 5 ...
- SLS机器学习介绍(05):时间序列预测
00系列文章目录 0.1 算法原理目录 SLS机器学习介绍(01):时序统计建模 SLS机器学习介绍(02):时序聚类建模 SLS机器学习介绍(03):时序异常检测建模 SLS机器学习介绍(04):规 ...
- Python数据分析案例-分别使用时间序列ARIMA、SARIMAX模型与Auto ARIMA预测国内汽车月销量
1. 前言 模型: ARIMA模型(英语:Autoregressive Integrated Moving Average model),差分整合移动平均自回归模型,又称整合移动平均自回归模型(移动也 ...
- 使用Flow forecast进行时间序列预测和分类的迁移学习介绍
ImageNet首次发表于2009年,在接下来的四年里,它成为了大多数计算机视觉模型的基础.到目前为止,无论您是在训练一个模型来检测肺炎还是对汽车模型进行分类,您都可能从在ImageNet或其他大型( ...
- 理论加实践,终于把时间序列预测ARIMA模型讲明白了
上篇我们一起学习了一些关于时间序列预测的知识.而本文将通过一段时间内电力负荷波动的数据集来实战演示完整的ARIMA模型的建模及参数选择过程,其中包括数据准备.随机性.稳定性检验.本文旨在实践中学习,在 ...
- 基于SARIMA、XGBoost和CNN-LSTM的时间序列预测对比
来源:DeepHub IMBA 本文约6800字,建议阅读10+分钟 本文将讨论通过使用假设测试.特征工程.时间序列建模方法等从数据集中获得有形价值的技术. 利用统计测试和机器学习分析和预测太阳能发电 ...
- 时间序列预测,非季节性ARIMA及季节性SARIMA
Python 3中使用ARIMA进行时间序列预测的指南 在本教程中,我们将提供可靠的时间序列预测.我们将首先介绍和讨论自相关,平稳性和季节性的概念,并继续应用最常用的时间序列预测方法之一,称为ARIM ...
- 11种常见的时间序列预测方法
参考内容:4大类11种常见的时间序列预测方法总结和代码示例 代码地址: https://github.com/SeafyLiang/machine_learning_study/blob/master ...
- 时间序列预测任务的模型选择最全总结
在第一部分,将了解多种时间序列的模型,如 经典的时间序列模型 监督学习模型 基于深度学习的模型 在第二部分,将建立几个时间序列模型来预测股市的应用案例,并了解一些时间序列建模技术.这些模型将相互比较, ...
最新文章
- AVL树入门(转载)
- python快速加引号_在python中如何快速地将一串字符串首尾加上双引号?
- zs040蓝牙模块使用方法_三菱CCLINK-IE网络模块使用方法
- Windows上的单个进程所能访问的最大内存量是多少?它与系统的最大虚拟内存一样吗?这对于系统设计有什么影响?...
- 【MySQL】查看 mysql 数据库最后修改时间、创建时间、表大小
- zabbix设置mysql登陆免报警_zabbix3.0.2 监控mysql服务down的时候不触发报警的问题修复...
- Numpy random生成不重复的随机数
- 2022年C语言程序设计题库习集带答案(史上最全版)
- Azkaban Flow 2.0的使用
- 合天网安就业班_【合天网安实验室】SQL注入入门一
- 用C++编写一个走迷宫游戏(1.0版)
- 【Unity3D开发小游戏】《青蛙过河》Unity开发教程
- 应用内更新,直接跳转到对应手机应用市场
- ActiveMQ简介
- js中的四种常用数组排序方法(冒泡、选择、插入、快排)及sort排序
- 特拉华大学计算机专业怎么样,2019上海软科世界一流学科排名计算机科学与工程专业排名特拉华大学排名第201-300...
- NTL库在Win上基于MinGW的安装
- 【论文导读】- Subgraph Federated Learning with Missing Neighbor Generation(FedSage、FedSage+)
- 架构师培训:aop是什么
- python 自动化测试基础