网站地区分站系统,青岛百度推广多少钱,wordpress+Apache升级,做室内意向图的网站理解平稳性
一般来说#xff0c;平稳时间序列是指随着时间的推移具有相当稳定的统计特性的时间序列#xff0c;特别是在均值和方差方面。平稳性可能是一个比较模糊的概念#xff0c;将序列排除为不平稳可能比说序列是平稳的更容易。通常不平稳序列有几个特征#xff1a;
…理解平稳性
一般来说平稳时间序列是指随着时间的推移具有相当稳定的统计特性的时间序列特别是在均值和方差方面。平稳性可能是一个比较模糊的概念将序列排除为不平稳可能比说序列是平稳的更容易。通常不平稳序列有几个特征
平均值随时间推移发生变化不保持稳定每个周期的波峰和波谷之间的距离在增长缩短也就是序列的方差随时间推移而增加减小整个序列表现出强烈的季节性行为与平稳性相反。
如何判断序列是否平稳
时间序列平稳性检验方法一般可分为两类
图形分析方法假设检验方法 1.图形分析方法
这是最简单的方法即可视化时间序列数据或可视化时间序列的统计特征进行肉眼判断。
1.1 可视化数据
import numpy as np
import pandas as pd
import akshare as ak
import matplotlib.pyplot as plt
import seaborn as sns
import warningssns.set(styledarkgrid)
warnings.filterwarnings(ignore)np.random.seed(123)# -------------- 准备数据 --------------
# 白噪声
white_noise np.random.standard_normal(size1000)# 随机游走
x np.random.standard_normal(size1000)
random_walk np.cumsum(x)# GDP
df ak.macro_china_gdp()
df df.set_index(季度)
# df.index pd.to_datetime(df.index)
gdp df[国内生产总值-绝对值][::-1].astype(float)# GDP DIFF
gdp_diff gdp.diff(4).dropna()# -------------- 绘制图形 --------------
fig, ax plt.subplots(2, 2, figsize(22,16), dpi300)ax[0][0].plot(white_noise)
ax[0][0].set_title(white_noise)
ax[0][1].plot(random_walk)
ax[0][1].set_title(random_walk)ax[1][0].plot(gdp)
ax[1][0].set_title(gdp)
ax[1][1].plot(gdp_diff)
ax[1][1].set_title(gdp_diff)plt.show()白噪声曲线围绕0值上下波动波动幅度前后、上下一致为平稳序列随机游走曲线无确定趋势均值、方差波动较大非平稳序列GDP数据趋势上升均值随时间增加非平稳序列GDP季节差分后数据曲线大致在一条水平线上下波动波动幅度前后变化较小可认为是平稳的。 1.2 可视化统计特征
可视化统计特征是指绘制时间序列自相关图和偏自相关图根据自相关图的表现来判断序列是否平稳。
自相关也叫序列相关是一个信号与自身不同时间点的相关度或者说与自身的延迟拷贝或滞后性的相关性是延迟的函数。不同滞后期得到的自相关系数叫自相关图。这里有一个默认假设是序列是平稳序列平稳序列的自相关性只和时间间隔k有关不随时间t的变化而变化因而可以称自相关函数是延迟k的函数
平稳序列通常具有短期自相关性对于平稳的时间序列自相关系数往往会迅速退化到零滞后期约短相关性越高滞后期为0相关性为1而对于非平稳的数据退化会发生得更慢或存在先减后增或者周期性波动等变动。
自相关的计算公式为根据滞后期 k k k将序列拆成等长的两个序列计算两个序列的相关性得到滞后期为 k k k时的自相关性。举例 X [ 2 , 3 , 4 , 3 , 8 , 7 ] , A [ 2 , 3 , 4 , 3 , 8 ] , B [ 3 , 4 , 3 , 8 , 7 ] X[2,3,4,3,8,7],\quad A[2,3,4,3,8], \quad B[3,4,3,8,7] X[2,3,4,3,8,7],A[2,3,4,3,8],B[3,4,3,8,7] X ˉ 1 6 ∑ i 1 6 X i 4.5 \bar{X} \frac{1}{6} \sum_{i1}^{6}X_i 4.5 Xˉ61i1∑6Xi4.5 s 2 ( X ) 1 6 ∑ i 1 6 ( X i − X ˉ ) ( X i − X ˉ ) 4.916 s^2(X)\frac{1}{6} \sum_{i1}^{6}(X_i- \bar{X})(X_i- \bar{X})4.916 s2(X)61i1∑6(Xi−Xˉ)(Xi−Xˉ)4.916 r ( 1 ) 1 5 ( A i − X ˉ ) ( B i − X ˉ ) 1.75 r(1)\frac{1}{5}(A_i-\bar{X})(B_i-\bar{X})1.75 r(1)51(Ai−Xˉ)(Bi−Xˉ)1.75 A C F ( 1 ) r ( 1 ) s 2 ( X ) 0.356 ACF(1)\frac{r(1)}{s^2(X)}0.356 ACF(1)s2(X)r(1)0.356
import statsmodels.api as sm
X [2,3,4,3,8,7]
sm.tsa.stattools.acf(X, nlags1, adjustedTrue) # 返回值[滞后期为0时的自相关性滞后期为1的自相关系数]array([1. , 0.3559322])根据 A C F ACF ACF求出滞后 k k k自相关系数时实际上得到并不是 X ( t ) X(t) X(t)与 X ( t − k ) X(t-k) X(t−k)之间单纯的相关关系。因为 X ( t ) X(t) X(t)还会受到中间 k − 1 k-1 k−1个随机变量的影响而这 k − 1 k-1 k−1个随机变量又都和 X ( t − k ) X(t-k) X(t−k)具有相关关系如果提出中间 k − 1 k-1 k−1个随机变量的影响所计算得到的是偏自相关系数计算较为复杂。
# 数据生成过程在第一个代码块中
from statsmodels.graphics.tsaplots import plot_acf, plot_pacffig, ax plt.subplots(4, 2, figsize(22, 16), dpi300)
fig.subplots_adjust(hspace1.2)plot_acf(white_noise, axax[0][0])
ax[0][0].set_title(ACF(white_noise))
plot_pacf(white_noise, axax[0][1])
ax[0][1].set_title(PACF(white_noise))plot_acf(random_walk, axax[1][0])
ax[1][0].set_title(ACF(random_walk))
plot_pacf(random_walk, axax[1][1])
ax[1][1].set_title(PACF(random_walk))plot_acf(gdp, axax[2][0])
ax[2][0].set_title(ACF(gdp))
plot_pacf(gdp, axax[2][1])
ax[2][1].set_title(PACF(gdp))plot_acf(gdp_diff, axax[3][0])
ax[3][0].set_title(ACF(gdp_diff))
plot_pacf(gdp_diff, axax[3][1])
ax[3][1].set_title(PACF(gdp_diff))plt.show()白噪音的自相关系数很快就衰减到0附近是明显的平稳序列。滞后期为0时的自相关系数和偏自相关系数其实就是序列自己和自己的相关性故为1滞后期为1时自相关系数为0表示白噪声无自相关性随机游走自相关系数下降非常缓慢故为非平稳序列另从偏自相关系数中可以看到随机游走只和前一项有关GDP数据的自相关图中也可以看到存在一定的周期性滞后4、8、12等自相关系数较大下降较慢差分后下降多一些起到一定效果认为差分后序列是平稳的。同可视化数据一样直观判断带有较强主观性但能让我们对数据有更直观的认识。 2.假设检验
平稳性的假设检验方法当前主流为单位根检验检验序列中是否存在单位根若存在则为非平稳序列不存在则为平稳序列。
2.1 DF检验 A D F ADF ADF检验(Augmented Dickey-Fuller Testing)是最常用的单位根检验方法之一通过检验序列是否存在单位根来判断序列是否是平稳的。 A D F ADF ADF是 D F DF DF检验的增强版。
迪基Dickey和弗勒Fuller1979年基于非平稳序列的基本特征将其大致归为三类并提出DF检验 当序列基本走势呈现无规则上升或下降并反复时将其归为无漂移项自回归过程当序列基本走势呈现明显的随时间递增或递减且趋势并不太陡峭时将其归为带漂移项自回归过程当序列基本走势随时间快速递增时则将其归为带趋势项回归过程。 对应检验回归式为 无漂移项自回归过程 Y t ρ Y t − 1 ε t , Y 0 0 Y_t\rho Y_{t-1} \varepsilon_t, \quad Y_00 YtρYt−1εt,Y00带漂移项自回归过程 Y t μ ρ Y t − 1 ε t , Y 0 0 Y_t\mu \rho Y_{t-1} \varepsilon_t, \quad Y_00 YtμρYt−1εt,Y00带漂移项和趋势项自回归过程 Y t μ β t ρ Y t − 1 ε t , Y 0 0 Y_t\mu \beta_t \rho Y_{t-1} \varepsilon_t, \quad Y_00 YtμβtρYt−1εt,Y00 其中 μ \mu μ是常数项 β t \beta_t βt是趋势项 ε t \varepsilon_t εt为白噪声无自相关性。 原假设 H 0 : ρ 1 H_0: \rho 1 H0:ρ1存在单位根为非平稳时间序列备择假设 H 1 : ρ 1 H_1: \rho 1 H1:ρ1不存在单位根时间序列是平稳的不含截距项和趋势项的平稳/含截距项的平稳/含截距项和趋势平稳 若检验统计量大于临界值 P P P值大于显著性水平 α \alpha α不能拒绝原假设序列是非平稳的 若检验统计量小于临界值 P P P值小于显著性水平 α \alpha α拒绝原假设序列是平稳的。 2.2 ADF检验 D F DF DF检验公式为一阶自回归过程为了能适用于高阶自回归过程的平稳性检验迪基等1984年对 D F DF DF检验进行了一定的修正引入了更高阶的滞后项 A D F ADF ADF的检验回归式修正为 无漂移项自回归过程 Y t ρ Y t − 1 ∑ i 1 k C i △ Y t − i ε t , Y 0 0 Y_t\rho Y_{t-1} \sum_{i1}^{k}C_i\bigtriangleup Y_{t-i} \varepsilon_t, \quad Y_00 YtρYt−1∑i1kCi△Yt−iεt,Y00带漂移项自回归过程 Y t μ ρ Y t − 1 ∑ i 1 k C i △ Y t − i ε t , Y 0 0 Y_t\mu \rho Y_{t-1} \sum_{i1}^{k}C_i\bigtriangleup Y_{t-i} \varepsilon_t, \quad Y_00 YtμρYt−1∑i1kCi△Yt−iεt,Y00带漂移项和趋势项自回归过程 Y t μ β t ρ Y t − 1 ∑ i 1 k C i △ Y t − i ε t , Y 0 0 Y_t\mu \beta_t \rho Y_{t-1} \sum_{i1}^{k}C_i\bigtriangleup Y_{t-i} \varepsilon_t, \quad Y_00 YtμβtρYt−1∑i1kCi△Yt−iεt,Y00 其中 μ \mu μ是常数项 β t \beta_t βt是趋势项 ε t \varepsilon_t εt为白噪声无自相关性。 假设条件不变
原假设 H 0 : ρ 1 H_0: \rho 1 H0:ρ1存在单位根为非平稳时间序列备择假设 H 1 : ρ 1 H_1: \rho 1 H1:ρ1不存在单位根时间序列是平稳的不含截距项和趋势项的平稳/含截距项的平稳/含截距项和趋势平稳
检验流程同 D F DF DF检验一致。若要严格判断序列是否是宽平稳的可以直接检验是否不含截距项和趋势项平稳若不能拒绝原假设如 p 0.05 p0.05 p0.05序列非平稳此时仍有必要检验序列是否是趋势平稳的。非平稳且非趋势平稳可以使用一阶差分等平稳化方法处理后再做检验弱势趋势平稳困于过度差分则不宜使用差分方式平稳化。
import numpy as np
from matplotlib import pyplot as pltnp.random.seed(123)y np.random.standard_normal(size100)
for i in range(1, len(y)):y[i] 1 0.1*i y[i]plt.figure(figsize(12, 6), dpi300)
plt.plot(y)
plt.show()from arch.unitroot import ADF
adf ADF(y)
print(adf.summary().as_text())adf ADF(y, trendct)
print(adf.summary().as_text())Augmented Dickey-Fuller Results Test Statistic -0.739
P-value 0.836
Lags 5
-------------------------------------Trend: Constant
Critical Values: -3.50 (1%), -2.89 (5%), -2.58 (10%)
Null Hypothesis: The process contains a unit root.
Alternative Hypothesis: The process is weakly stationary.Augmented Dickey-Fuller Results Test Statistic -9.963
P-value 0.000
Lags 0
-------------------------------------Trend: Constant and Linear Time Trend
Critical Values: -4.05 (1%), -3.46 (5%), -3.15 (10%)
Null Hypothesis: The process contains a unit root.
Alternative Hypothesis: The process is weakly stationary.arch包中 A D F ADF ADF检验可指定trend n不含截距项和时间趋势项、trend c含截距项、trend ct含截距项和时间趋势项、trend ctt含截距项、时间趋势项和二次型时间趋势项分别对应不同平稳类型的检验。
再来看看GDP季节差分前后数据是否为平稳
# 数据在第一个代码块中
from arch.unitroot import ADF
adf ADF(gdp)
print(adf.summary().as_text())adf ADF(gdp_diff) # 差分后平稳
print(adf.summary().as_text())Augmented Dickey-Fuller Results Test Statistic 2.606
P-value 0.999
Lags 12
-------------------------------------Trend: Constant
Critical Values: -3.55 (1%), -2.91 (5%), -2.59 (10%)
Null Hypothesis: The process contains a unit root.
Alternative Hypothesis: The process is weakly stationary.Augmented Dickey-Fuller Results Test Statistic -3.651
P-value 0.005
Lags 8
-------------------------------------Trend: Constant
Critical Values: -3.55 (1%), -2.91 (5%), -2.59 (10%)
Null Hypothesis: The process contains a unit root.
Alternative Hypothesis: The process is weakly stationary.# 数据在第一个代码块中
from arch.unitroot import ADF
adf ADF(gdp, trendct) # 差分前非趋势平稳
print(adf.summary().as_text())Augmented Dickey-Fuller Results Test Statistic -2.358
P-value 0.402
Lags 4
-------------------------------------Trend: Constant and Linear Time Trend
Critical Values: -4.10 (1%), -3.48 (5%), -3.17 (10%)
Null Hypothesis: The process contains a unit root.
Alternative Hypothesis: The process is weakly stationary.2.3 PP检验
Phillips和Perron(1988)提出了一种非参数检验方法主要是为了解决残差项中潜在的序列相关和异方差问题其检验统计量的渐进分布和临界值与ADF检验相同。同样出现较早假设条件一样用法相似可作为ADF检验的补充。
假设条件不变
原假设 H 0 : ρ 1 H_0: \rho 1 H0:ρ1存在单位根为非平稳时间序列备择假设 H 1 : ρ 1 H_1: \rho 1 H1:ρ1不存在单位根时间序列是平稳的不含截距项和趋势项的平稳/含截距项的平稳/含截距项和趋势平稳
同样构造一个趋势平稳序列看下PP检验结果
import numpy as np
from arch.unitroot import PhillipsPerronnp.random.seed(123)y np.random.standard_normal(size100)
for i in range(1, len(y)):y[i] 1 0.1*i y[i]pp PhillipsPerron(y) # 10%显著水平下平稳
print(pp.summary().as_text())pp PhillipsPerron(y, trendct) # 1%显著水平下趋势平稳
print(pp.summary().as_text())Phillips-Perron Test (Z-tau) Test Statistic -2.825
P-value 0.055
Lags 12
-------------------------------------Trend: Constant
Critical Values: -3.50 (1%), -2.89 (5%), -2.58 (10%)
Null Hypothesis: The process contains a unit root.
Alternative Hypothesis: The process is weakly stationary.Phillips-Perron Test (Z-tau) Test Statistic -10.009
P-value 0.000
Lags 12
-------------------------------------Trend: Constant and Linear Time Trend
Critical Values: -4.05 (1%), -3.46 (5%), -3.15 (10%)
Null Hypothesis: The process contains a unit root.
Alternative Hypothesis: The process is weakly stationary.2.4 DF-GLS检验
DF-GLS检验是Elliot,Rothenberg,and Stock 1996年提出的一种单位根检验方法全称Dickey-Fuller Test with GLS Detredding即“使用广义最小二乘法去除趋势的检验”是目前最有功效的单位根检验。
DF-GLS检验利用广义最小二乘法首先对要检验的数据进行一次“准差分”然后利用准差分的数据对原序列进行去除趋势处理再利用ADF检验的模型形式对去除趋势后的数据进行单位根检验但此时ADF检验模型中不再包含常数项或者时间趋势变量。
原假设序列存在单位根时间序列是非平稳的备择假设序列不存在单位根时间序列是平稳的或趋势平稳的
同样构造一个趋势平稳序列来看下检验效果DF-GLS检验trend只能指定为c或ct
import numpy as np
from arch.unitroot import DFGLSnp.random.seed(123)y np.random.standard_normal(size100)
for i in range(1, len(y)):y[i] 1 0.1*i y[i]dfgls DFGLS(y)
print(dfgls.summary().as_text())dfgls DFGLS(y, trendct) # 趋势平稳
print(dfgls.summary().as_text())Dickey-Fuller GLS Results Test Statistic 1.186
P-value 0.945
Lags 4
-------------------------------------Trend: Constant
Critical Values: -2.77 (1%), -2.15 (5%), -1.83 (10%)
Null Hypothesis: The process contains a unit root.
Alternative Hypothesis: The process is weakly stationary.Dickey-Fuller GLS Results Test Statistic -7.565
P-value 0.000
Lags 0
-------------------------------------Trend: Constant and Linear Time Trend
Critical Values: -3.62 (1%), -3.04 (5%), -2.74 (10%)
Null Hypothesis: The process contains a unit root.
Alternative Hypothesis: The process is weakly stationary.2.5 KPSS检验
KPSS检验是由Kwiatkowski, Phillips, and Shin 1992年提出的与之前的检验方法相比最大的不同点就是它的原假设是平稳序列或趋势平稳序列而备择假设是存在单位根。
原假设序列不存在单位根时间序列是平稳的或趋势平稳的备择假设序列存在单位根时间序列是非平稳的
import numpy as np
from arch.unitroot import KPSSnp.random.seed(123)y np.random.standard_normal(size100)
for i in range(1, len(y)):y[i] 0.1 y[i-1] y[i]kpss KPSS(y)
print(kpss.summary().as_text())kpss KPSS(y, trendct)
print(kpss.summary().as_text())KPSS Stationarity Test Results Test Statistic 1.393
P-value 0.000
Lags 5
-------------------------------------Trend: Constant
Critical Values: 0.74 (1%), 0.46 (5%), 0.35 (10%)
Null Hypothesis: The process is weakly stationary.
Alternative Hypothesis: The process contains a unit root.KPSS Stationarity Test Results Test Statistic 0.114
P-value 0.115
Lags 5
-------------------------------------Trend: Constant and Linear Time Trend
Critical Values: 0.22 (1%), 0.15 (5%), 0.12 (10%)
Null Hypothesis: The process is weakly stationary.
Alternative Hypothesis: The process contains a unit root.KPSS检验中原假设为不存在单位根。默认检验趋势类型下P值为0.000拒绝原假设存在单位根序列非平稳。指定trendct后P值不拒绝原假设认为序列趋势平稳检验错误。以上几种检验均不能100%保证检验正确PP检验可认为是ADF检验的补充KPSS检验同样也可和其他检验一同使用当均认为是平稳或趋势平稳时方判定为平稳。
