时间序列分析2 ARIMA(python)

1.数据的平稳性

1.1.平稳性:

1.平稳性是要求经由样本时间序列所得到的拟合曲线在未来一段时间内仍能顺着现有的形态“惯性”地延续下去。
2.平稳性要求序列的均值和方差不发生明显的变化

1.2.严平稳和弱平稳

1 严平稳:严平稳表示的分布不随时间的改变而改变。
如:白噪声(正态),无论怎么取,都是期望为0,方差为1。

  1. 弱平稳:期望与相关系数(依赖性)不变。未来某时刻的t的值Xt就要依赖于它的过去信息,所以需要依赖性。

1.3. 差分法:时间序列在t与t-1时刻的差值。

在一阶差分的结果上再做一阶差分那就是二阶差分。
数据波动特别大的时候,使用差分法,使得数据形态相对稳定。例如股票日浮动波动很大,但是五日均线就平缓一些,十日比五日更平缓,如果看年线就几乎没有什么波动。


2. ARIMA

2.1.自回归模型(AR)

2.1.1描述当前值与历史值之间的关系,用变量自身的历史时间数据对自身进行预测。
2.1.2自回归模型必须满足平稳性的要求。

2.1.3P阶自回归过程的公式定义:

P阶的意思是当前值和前面P个值相关。
2.1.4 yt是当前值, μ是常数项(类似线性回归中的b),P是阶数,γi是自相关系数,εt是误差。γi就是我们需要求的参数(类似线性回归中的w)。

2.2.自回归模型(AR)的限制

2.2.1 自回归模型是用自身的数据来进行预测。
2.2.2 必须具有平稳性。
2.2.3 必须基友自相关性,如果自相关系数小于(φi)0.5,则不宜采用。
2.2.4 自回归只适用于预测与自杀前期相关的现象。

2.3. 移动平均模型(MA)
2.3.1. 移动平均模型关注的是自回归模型中的误差项的累加

2.3.2. q阶自回归过程的公式定义:

2.3.3. 移动平均法能有效地消除预测中的随机波动。

2.4. 自回归移动平均模型(ARMA)
2.4.1 自回归与移动平均的结合
2.4.2 公式:


2.4.3. ARIMA(p,d,q)模型全称为差分自回归移动平均模型(Autoregressive Integrated Moving Average Model),在python统计学包statsmodels中。参数p,d,q,p和d为自回归项和移动平均项数,d代表差分数,例如指定d=1 就是一阶差分。
2.4.4. AR是自回归,p为自回归项;MA为移动平均,q为移动平均项数,d为时间序列成为平稳是所做的差分次数。
2.4.5. 原理:将非平稳时间序列转化为平稳时间序列然将因变量仅对它的滞后值以及随机误差项的现值和滞后值进行回归所建立的模型。例如p = 2 就代表自回归滞后两阶。

3. 建模流程和调参

3.1.自相关系数ACF(autocorrelation function)

3.1.1 有序的随机变量序列与其自身相比较。自相关函数反映了同一序列再不同时序的取值之间的相关性。

3.1.2 公式:

3.1.3 Pk的取值范围为[-1,1]

3.2. 偏自相关函数(PACF)(partial autocorrelation function)

3.2.1 对于一个平稳AR(p)模型,求出滞后k自相关系数p(k)时
实际上得到并不是x(t)与x(t-k)之间单纯的相关关系。
3.2.2 x(t)同时还会受到中间k-1个随机变量x(t-1)、x(t-2)、……、x(t-k+1)的影响,而这k-1个随机变量又都和x(t-k)具有相关关系,所以自相关系数p(k)里实际掺杂了其他变量对x(t)与x(t-k)的影响。
3.2.3 剔除了中间k-1个随机变量x(t-1)、x(t-2)、……、x(t-k+1)的干扰之后x(t-k)对x(t)影响的相关程度。
3.2.4 ACF还包含了其他变量的影响,而偏自相关系数PACF是严格这两个变量之间的相关性。

3.3 ARIMA(p,d,q)参数确定

3.3.1 图表



3.3.2 截尾:落在置信区间内(95%的点都符合该规则)

3.3.3 ARIMA(p,d,q)阶数确定
AR(p)看PACF,p阶后截尾,意思是PACF图中找到一阶,后面的部分的点都落在置信区间(蓝色部分)上,少部分点不在可以忽略。此图中p可以取1。
MA(q)看ACF,q阶后截尾,意思是ACF图中找到一阶,后面的部分的点都落在置信区间(蓝色部分)上,少部分点不在可以忽略。此图中q可以取7。

3.4 ARIMA 建模流程:

  1. 将序列平稳(差分法确定d)
  2. p和q阶数确定:ACF与PACF
  3. ARIMA(p,d,q)

3.5. 模型选择指标AIC与BIC: 选择更简单的模型(AIC与BIC都是越低越好)

  1. AIC:赤池信息准则(AkaikeInformation Criterion,AIC)
    𝐴𝐼𝐶=2𝑘−2ln(𝐿)(在保证模型精度的同时,尽量让k更小,这样模型更简单)
  2. BIC:贝叶斯信息准则(Bayesian Information Criterion,BIC)
    𝐵𝐼𝐶=𝑘𝑙𝑛𝑛−2ln(𝐿)
  3. k为模型参数个数,n为样本数量,L为似然函数

3.6. 模型残差检验:

  1. ARIMA模型的残差是否是平均值为0且方差为常数的正态分布
  2. QQ图:线性即正态分布


    用plot_diagnostics()函数得到诊断图

4. 例子


Reference:
https://www.analyticsvidhya.com/blog/2016/02/time-series-forecasting-codes-python/
https://blog.csdn.net/goldxwang/article/details/73351155
https://github.com/aarshayj/Analytics_Vidhya/tree/master/Articles

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,098评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,213评论 2 380
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 149,960评论 0 336
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,519评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,512评论 5 364
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,533评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,914评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,574评论 0 256
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,804评论 1 296
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,563评论 2 319
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,644评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,350评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,933评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,908评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,146评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,847评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,361评论 2 342

推荐阅读更多精彩内容

  • 不同图像灰度不同,边界处一般会有明显的边缘,利用此特征可以分割图像。需要说明的是:边缘和物体间的边界并不等同,边缘...
    大川无敌阅读 13,815评论 0 29
  • "use strict";function _classCallCheck(e,t){if(!(e instanc...
    久些阅读 2,027评论 0 2
  • 过得还是蛮开心的。但我一直有种错觉,我身边的人,把我当猴耍。尽管我不是很介意,能给人带来快乐就好。也许这就是...
    zsm泡泡阅读 197评论 0 0
  • 今日老友相见! 千言万语,不用多说,亲人一般的老盆友。 一起努力,变成更好的自己吧! 元气少女,继续加油哦~~ 【...
    霁青阅读 145评论 0 0
  • 为大鱼海棠鸣不平 到现在还没有来得及去电影院看大鱼海棠,却已经在简书,豆瓣,朋友圈看到了很多吐槽大鱼海棠的文章,当...
    萧牧寒阅读 699评论 3 5