python时间序列预测不连续怎么办_手把手教你用Python处理非平稳时间序列

1. 平稳简介

“平稳”是处理时间序列数据时遇到的最重要的概念之一：平稳序列是指其特性-均值、方差和协方差不随时间而变化的序列。

让我们用一个直观的例子来理解这一点。考虑以下三个图形：

在第一幅图中，我们可以清楚地看到，均值随时间而变化(增加)，呈现上升的趋势。因此，这是一个非平稳序列。平稳序列不应该呈现出随时间变化的趋势。

第二幅图显然看不到序列的趋势，但序列的变化是一个时间的函数。正如前面提到的，平稳序列的方差必须是一个常数。

再来看第三幅图，随着时间的增加，序列传播后变得更近，这意味着协方差是时间的函数。

上述三个例子均是非平稳时间序列。现在来看第四个图：

在这张图中，均值、方差和协方差都是常数，这就是平稳时间序列。

再想一想，上面的哪一幅图预测未来会更容易呢？第四个图，对吧？大多数统计模型都要求序列是平稳的，这样才能进行有效和精确的预测。

因此，总的来说，平稳时间序列是一个不依赖时间变化 (即均值、方差和协方差不随时间变化)的时间序列。在下一节中，我们将介绍各种检测给定序列是否平稳的方法。

2. 加载数据

在本节和后续几节中，将介绍检测时间序列数据的平稳性的方法，以及如何处理非平稳序列。同时，本文还提供了相应的Python代码。大家可以到：AirPassengers下载文中使用的数据集。

在继续分析数据集之前，首先加载和预处理数据。

好了，看来可以继续了！

3. 检验平稳的方法

下一步是确定给定的序列是否是平稳的，并对它做相应的处理。本节将介绍一些常见的方法，利用这些方法来检测序列是否平稳。

目视检验

看一下我们在上一节中使用的图形，仅需通过目测图形，便能够识别出序列的均值和方差是否随时间变化。同样，通过绘制数据图形，便能确定该序列的属性是否随时间而变化。

显然，通过作图，可以看到序列的趋势(变化的均值)，然而，这种目视检测方法得到的结果可能不准确。最好是用一些统计检验方法来验证观测结果。

统计检验

可以利用统计检验来代替目视检验：比如单位根平稳检验。单位根表名给定序列的统计特性（均值，方差和协方差）不是时间的常数，这是平稳时间序列的先决条件。下面是它的数学解释：

假设我们有一个时间序列：

其中yt是t时刻的数据值，ε t 是误差项。需要利用yt-1的值来计算yt，即：

如果利用所有的观察值，yt 的值将是：

假设在上述方程中a的值为1(单位)，则预测值将等于yt-n 和从t-n到t的所有误差之和，这意味着方差将随着时间的推移而增大，这就是时间序列中的单位根。众所周知，平稳时间序列的方差不能是时间的函数。单元根检验通过检查a=1的值来检查序列中是否存在单位根。以下是两个最常用的单位根平稳检测方法：

ADF（增补迪基-福勒）检验

The Dickey Fuller test is one of the most popular statistical tests. It can be used to determine the presence of unit root in the series, and hence help us understand if the series is stationary or not. The null and alternate hypothesis of this test are:

迪基-福勒（Dickey Fuller）检验是最流行的统计检验方法之一，可以用它来确定序列中单位根的存在，从而帮助判断序列是否是平稳。这一检验的原假设与备择假设如下：

原假设：序列有一个单位根(a=1的值)

备择假设：该序列没有单位根。

如果不能拒绝原假设，则该序列是非平稳的，这意味着序列可以是线性的，也可以是差分平稳的(将在下一节中更多地了解差分平稳)。

Python代码：

ADF检验结果：ADF检验的统计量为1%，p值为5%，临界值为10%，置信区间为10%。我们对本序列的检验结果如下：

平稳性检验：如果检验统计量小于临界值，我们可以拒绝原假设(也就是序列是平稳的)。当检验统计量大于临界值时，不能拒绝原假设(这意味着序列不是平稳的)。

在上面的例子中，检验统计量>临界值，这意味着序列不是平稳的。这证实了我们最初在目视检测中观察的结果。

KPSS（科瓦特科夫斯·基菲利普·斯施密特·辛）检验

KPSS检验是另一种用于检查时间序列的平稳性 (与迪基-福勒检验相比稍逊一筹) 的统计检验方法。KPSS检验的原假设与备择假设与ADF检验的原假设与备择假设相反，常造成混淆。

KPSS检验的作者将原假设定义为趋势平稳，并将备择假设定义为单位根序列。我们将在下一节详细了解趋势平稳。现在，来看一下KPSS检验的实现，并查看KPSS检验的结果。

原假设：序列是趋势平稳的。

备择假设：序列有一个单位根(序列是非平稳的)。

Python代码：

KPSS检验结果：KPSS检验-检验统计量、p-值和临界值和置信区间分别为1%、2.5%、5%和10%。对于航空乘客数据集的检验结果如下：

平稳性检验：如果检验统计量大于临界值，则拒绝原假设(序列不是平稳的)。如果检验统计量小于临界值，则不能拒绝原假设(序列是平稳的)。对于航空乘客数据集来说，在所有置信区间，检验统计量的值都大于临界值，因此可以说该序列是不平稳的。

在为时间序列数据集准备模型之前，通常会同时进行两种检验。有一次，这两种检验显示出相互矛盾的结果：其中一个检验结果表明该序列是平稳的，而另一个则表明该序列是非平稳的！我困惑了好几个小时，想弄清楚这是怎么回事。后来才知道，序列的平稳性有多种类型。

综上所述，ADF检验有线性平稳或差分平稳的备择假设，而KPSS检验则是识别序列的趋势平稳。

4. 平稳的种类

通过了解不同类型的平稳，来解释上述检验的结果。

严格平稳：严格平稳序列满足平稳过程的数学定义。严格平稳序列的均值、方差和协方差均不是时间的函数。我们的目标是将一个非平稳序列转化为一个严格平稳序列，然后对它进行预测。

趋势平稳：没有单位根但显示出趋势的序列被称为趋势平稳序列。一旦去除趋势之后，产生的序列将是严格平稳的。在没有单位根的情况下，KPSS检测将该序列归类为平稳。这意味着序列可以是严格平稳的，也可以是趋势平稳的。

差分平稳：通过差分可以使时间序列成为严格平稳的时间序列。ADF检验也称为差分平稳性检验。

应用两种检验总会更优些，通过两种检验之后，可以确定这个序列是否是平稳的。下面，来看一下应用两种平稳检验后的可能结果：

结果1：两种检验均得出结论：序列是非平稳的->序列是非平稳的

结果2：两种检验均得出结论：序列是平稳的->序列是平稳的

结果3：KPSS =平稳；ADF =非平稳->趋势平稳，去除趋势后序列严格平稳

结果4：KPSS =非平稳；ADF =平稳->差分平稳，利用差分可使序列平稳。5. 时间序列的平稳化

在熟悉了平稳性的概念及其不同的类型之后，接下来可以对序列进行平稳化操作。请记住，为了建立时间序列预测模型，必须首先将任何非平稳序列转换为平稳序列。

差分化

在该方法中，计算序列中连续项的差值。执行差分操作通常是为了消除均值的变化。从数学角度，差分可以写成：

yt‘ = yt – y(t-1)

其中yt 是t时刻的数值。

对序列差分化后，绘制出结果：

季节性差分

在季节性差分中，不计算连续值之间的差异，而是计算观察值与同一季节的先前观察值之间的差异。例如，星期一的观察值将与上星期一的观察值相减。从数学角度，它可以写成：

yt‘ = yt – y(t-n)

变换

变换用于对方差为非常数的序列进行平稳化。常用的变换方法包括幂变换、平方根变换和对数变换。对飞机乘客数据集进行快速对数转换和差分：

可以看出，这个图形比先前的图形有了很大的改善。通过对这个序列进行平方根或幂变换，看看是否得出了更好的结果。欢迎在下面的评论里分享你的发现！

尾注

本文介绍了检验时间序列平稳性的不同方法。但并不止步于此，下一步是对得到的序列应用一个预测模型。可以参考以下文章来构建这样的模型：给初学者的时间序列预测综合指南（Beginner’s Guide to Time Series Forecast）。

如果对本文有任何问题或反馈，可以在下面的评论部分与我联系。

还可以在Analytics Vidhya的Android应用程序上阅读这篇文章

python时间序列预测不连续怎么办_手把手教你用Python处理非平稳时间序列相关推荐

python时间序列预测不连续怎么办_手把手教你用Python处理非平稳时间序列（附代码）...
本文约3600字,建议阅读10分钟. 本文将重点介绍时间序列数据的平稳性检验方法. 简介预测一个家庭未来三个月的用电量,估计特定时期道路上的交通流量,预测一只股票在纽约证券交易所交易的价格--这些问 ...
python正确的输入语句_手把手教你在python中如何使用while True语句
在学习过程中,经常能遇到采用while True的用法.下面以一个例子进行说明: 建立一个用户登录系统,用户输入用户名和密码,如果正确就可以进入系统. 1.我自己最开始的写法:d = {} #数据库字 ...
python爬虫爬取房源_手把手教你用Python网络爬虫爬取新房数据
项目背景大家好,我是J哥. 新房数据,对于房地产置业者来说是买房的重要参考依据,对于房地产开发商来说,也是分析竞争对手项目的绝佳途径,对于房地产代理来说,是踩盘前的重要准备. 今天J哥以「惠民之家」 ...
如何用python抓取qq音乐_手把手教你使用Python抓取QQ音乐数据（第三弹）-阿里云开发者社区...
[一.项目目标] 通过手把手教你使用Python抓取QQ音乐数据(第一弹)我们实现了获取 QQ 音乐指定歌手单曲排行指定页数的歌曲的歌名.专辑名.播放链接. 通过手把手教你使用Python抓取QQ音乐 ...
python网页爬虫循环获取_手把手教你用 Python 搞定网页爬虫
原标题:手把手教你用 Python 搞定网页爬虫编译:欧剃作为数据科学家的第一个任务,就是做网页爬取.那时候,我对使用代码从网站上获取数据这项技术完全一无所知,它偏偏又是最有逻辑性并且最容易获得的 ...
python预测实例教程_手把手教你用Python库Keras做预测（附代码）-阿里云开发者社区...
当你在Keras中选择好最合适的深度学习模型,就可以用它在新的数据实例上做预测了.但是很多初学者不知道该怎样做好这一点,我经常能看到下面这样的问题: "我应该如何用Keras对我的模型作出预 ...
python热搜排行功能_手把手教你用Python+Pyecharts让微博热搜榜动起来
今天教大家如何用pyecharts制作微博热搜榜动态展示视频,先上视频看看效果: 教程主要有2部分: 一是python爬取微博热搜内容二是用pyecharts制作动态视频下面给大家详细介绍一下一 ...
python 函数修饰器父类_手把手教你学python第十四讲（函数装饰器,super用法和时间处理）...
文中有些字在图中是因为每篇文章最多100张图片,我把有的小图片和文字一起截图了,文中所有的引用都会标出原文网址,除此以外都是作者原创. 有时候会在文章最前或者最后补充一些知识或者把前面说的有问题的地方 ...
python图片分类技术介绍_手把手教你用Python构建你的第一个多标签图像分类模型（附案例）！...
介绍你正在处理图像数据吗?我们可以使用计算机视觉算法来做很多事情:对象检测图像分割图像翻译对象跟踪(实时),还有更多-- 这让我思考--如果一个图像中有多个对象类别,我们该怎么办?制作一个图像 ...
python 桑基图地理坐标_手把手教你用Python绘制酷炫的桑基图
最近 ,不止一次收到小伙伴的截图追问: "这个图叫什么 ? ? ?" "这个图真好看 ! ! !怎么画啊 ?" ...... 笔者本没有干货 ,问的人多了 ,也 ...

python时间序列预测不连续怎么办_手把手教你用Python处理非平稳时间序列

python时间序列预测不连续怎么办_手把手教你用Python处理非平稳时间序列相关推荐

最新文章

热门文章