单位根检验

单位根检验
单位根检验

Exogenous: Constant, Linear Trend

Lag Length: 2 (Automatic - based on SIC, maxlag=13)

t-Statistic

Prob.*Augmented Dickey-Fuller test statistic -1.605700 0.7863

Test critical values:1% level

-4.0216915% level -3.44068110% level -3.144830

*MacKinnon (1996) one-sided p-values.

Augmented Dickey-Fuller Test Equation Dependent Variable: D(Y)Method: Least Squares

Date: 11/12/14 Time: 23:32Sample (adjusted): 4 150

Included observations: 147 after adjustments

Variable Coefficient Std. Error t-Statistic Prob. Y(-1)-0.0196270.012223-1.6057000.1106D(Y(-1))0.2535480.081530 3.1098870.0023D(Y(-2))0.2146390.081798 2.6240080.0096C

4.074679 2.403364 1.6954070.0922@TREND("1")0.008884

0.006079

1.461269

0.1462R-squared

0.148873 Mean dependent var 0.430612Adjusted R-squared 0.124898 S.D. dependent var 1.450725S.E. of regression 1.357109 Akaike info criterion 3.482012Sum squared resid 261.5277 Schwarz criterion 3.583727Log likelihood -250.9279 Hannan-Quinn criter. 3.523340F-statistic

6.209410 Durbin-Watson stat 2.054851

Prob(F-statistic)

0.000124

Lag length:2 即滞后阶数为2,则初始估计模型为

0111122t t t t Y c c t Y Y Y λββ---?=+++?+?

因为ADF 的t=-1.6057> 5%level 的t=-3.440681,所以接受H 0,即存在单位根。 (或因为p=0.7863>α=0.05,所以接受H 0,即存在单位根。)

又因为@TREND(“1”)的p=0.1462>α=0.05,所以接受H 0,即c 1显著为0。 则模型改为011122t t t t Y c Y Y Y λββ---?=++?+?

Exogenous: Constant

Lag Length: 2 (Automatic - based on SIC, maxlag=13)

t-Statistic

Prob.*Augmented Dickey-Fuller test statistic -0.663786 0.8513

Test critical values:1% level

-3.4751845% level -2.88112310% level -2.577291

*MacKinnon (1996) one-sided p-values.

Augmented Dickey-Fuller Test Equation Dependent Variable: D(Y)Method: Least Squares

Date: 11/12/14 Time: 23:38Sample (adjusted): 4 150

Included observations: 147 after adjustments

Variable Coefficient Std. Error t-Statistic Prob. Y(-1)-0.0035350.005326-0.6637860.5079D(Y(-1))0.2496220.081808 3.0513070.0027D(Y(-2))0.2042410.081811 2.4964990.0137C

1.052158

1.228690

0.856325

0.3933R-squared

0.136074 Mean dependent var 0.430612Adjusted R-squared 0.117950 S.D. dependent var 1.450725S.E. of regression 1.362485 Akaike info criterion 3.483332Sum squared resid 265.4604 Schwarz criterion 3.564704Log likelihood -252.0249 Hannan-Quinn criter. 3.516394F-statistic

7.507838 Durbin-Watson stat 2.046101

Prob(F-statistic)

0.000106

因为ADF 的t=-0.663786> 5%level 的t=-2.881123,所以接受H 0,即存在单位根。 (或因为p=0.8513>α=0.05,所以接受H 0,即存在单位根。) 又因为C 的p=0.3933>α=0.05,所以接受H 0,即c 0显著为0。 则模型改为11122t t t t Y Y Y Y λββ---?=+?+?

Exogenous: None

Lag Length: 2 (Automatic - based on SIC, maxlag=13)

t-Statistic

Prob.*Augmented Dickey-Fuller test statistic 1.940255 0.9875

Test critical values:1% level

-2.5807885% level -1.94301210% level -1.615270

*MacKinnon (1996) one-sided p-values.

Augmented Dickey-Fuller Test Equation Dependent Variable: D(Y)Method: Least Squares

Date: 11/13/14 Time: 18:53Sample (adjusted): 4 150

Included observations: 147 after adjustments

Variable Coefficient Std. Error t-Statistic Prob. Y(-1)0.0010040.000517 1.9402550.0543D(Y(-1))0.2492900.081732 3.0501040.0027D(Y(-2))

0.201130

0.081655

2.463183

0.0149R-squared

0.131644 Mean dependent var 0.430612Adjusted R-squared 0.119584 S.D. dependent var 1.450725S.E. of regression 1.361223 Akaike info criterion 3.474841Sum squared resid 266.8216 Schwarz criterion 3.535870Log likelihood

-252.4008 Hannan-Quinn criter. 3.499638

Durbin-Watson stat

2.043355

因为ADF 的t=1.940255> 5%level 的t=-1.943012,所以接受H 0,即存在单位根。 (或因为p=0.9875>α=0.05,所以接受H 0,即存在单位根。) 所以模型:11122t t t t Y Y Y Y λββ---?=+?+?

Exogenous: Constant, Linear Trend

Lag Length: 1 (Automatic - based on SIC, maxlag=13)

t-Statistic Prob.*

Augmented Dickey-Fuller test statistic-5.739450 0.0000 Test critical values:1% level-4.021691

5% level-3.440681

10% level-3.144830

*MacKinnon (1996) one-sided p-values.

Augmented Dickey-Fuller Test Equation

Dependent Variable: D(Y,2)

Method: Least Squares

Date: 11/12/14 Time: 23:41

Sample (adjusted): 4 150

Included observations: 147 after adjustments

Variable Coefficient Std. Error t-Statistic Prob.

D(Y(-1))-0.5510030.096003-5.7394500.0000

D(Y(-1),2)-0.2007290.081786-2.4543200.0153 C0.2336790.233654 1.0001070.3189 @TREND("1")8.93E-050.0026530.0336640.9732

R-squared0.371186 Mean dependent var0.004082 Adjusted R-squared0.357994 S.D. dependent var 1.703054 S.E. of regression 1.364577 Akaike info criterion 3.486400 Sum squared resid266.2762 Schwarz criterion 3.567773 Log likelihood-252.2504 Hannan-Quinn criter. 3.519463 F-statistic28.13737 Durbin-Watson stat 2.043400 Prob(F-statistic)0.000000

对初始模型进行1阶差分后进行分析:

因为ADF的p=0.0000<α=0.05,所以拒绝原假设H0,即不存在单位根。

单位根过程和单位根检验

第二章 单位根过程和单位根检验 第一节 单位根过程 从本章开始我们进入时间序列的非平稳分析和建模研究。前面的章节的内容主要考虑的是平稳时间序列的建模和预测问题,但对于非平稳的时间序列,只有先进行差分处理,将其转换为平稳的时间序列模型。这样会损失部分信息。本章从理论上介绍非平稳时间序列的性质,讨论非平稳时间序列数据建模的伪回归问题。 非平稳序列的分析建立在维纳过程(布朗运动)和泛函中心极限定理之上。 一. 若干定义 定义1: (1)白噪声过程(white noise ,如图1)。属于平稳过程。 εε2 t t,t y =~iid(0,σ) 图3是日元兑美元差分序列(收益序列),近似于白噪声序列。 (2)随机游走过程(random walk ,如图2)。属于非平稳过程。 εε+2 t t-1t,t y =y ~iid(0,σ) 随机游走的差分过程是平稳过程(白噪声过程)。?yt =t ε。 -3 -2 -1 012 3 100120140160180200220240260280300 white noise -10 -50 510 20 40 60 80 140160y=y(-1)+u 图 1 白噪声序列(σ2=1) 图2 随机游走序列(σ2=1) 随机游走过程是非平稳的,这是因为: +t 012t y =y +u +u +u +t 012t 0E(y )=E(y +u +u +u )=y →∞22t 012t 12t D(y )=D(y +u +u ++u )=E(u +u ++u )=t σ 定义2:单位根过程

随机过程t,{y t =1,2,} 是一单位根过程,若t t-1t y =y +u t =1,2 t u 为一平稳过程,且t t t-s s E(u )=0,cov(u ,u )=μs =0,1,2 定义3:维纳过程 维纳过程(Wiener Process)也称为布朗运动过程(Brownian Motion Process)。 设W(t)是定义在闭区间[0,1]上一连续变化的随机过程,若该过程满足: (a) W(0)=0; (b) 对闭区间[0,1]上任意一组分割 12k 0≤t

单位根检验

Exogenous: Constant, Linear Trend Lag Length: 2 (Automatic - based on SIC, maxlag=13) t-Statistic Prob.*Augmented Dickey-Fuller test statistic -1.605700 0.7863 Test critical values:1% level -4.0216915% level -3.44068110% level -3.144830 *MacKinnon (1996) one-sided p-values. Augmented Dickey-Fuller Test Equation Dependent Variable: D(Y)Method: Least Squares Date: 11/12/14 Time: 23:32Sample (adjusted): 4 150 Included observations: 147 after adjustments Variable Coefficient Std. Error t-Statistic Prob. Y(-1)-0.0196270.012223-1.6057000.1106D(Y(-1))0.2535480.081530 3.1098870.0023D(Y(-2))0.2146390.081798 2.6240080.0096C 4.074679 2.403364 1.6954070.0922@TREND("1")0.008884 0.006079 1.461269 0.1462R-squared 0.148873 Mean dependent var 0.430612Adjusted R-squared 0.124898 S.D. dependent var 1.450725S.E. of regression 1.357109 Akaike info criterion 3.482012Sum squared resid 261.5277 Schwarz criterion 3.583727Log likelihood -250.9279 Hannan-Quinn criter. 3.523340F-statistic 6.209410 Durbin-Watson stat 2.054851 Prob(F-statistic) 0.000124 Lag length:2 即滞后阶数为2,则初始估计模型为 0111122t t t t Y c c t Y Y Y λββ---?=+++?+? 因为ADF 的t=-1.6057> 5%level 的t=-3.440681,所以接受H 0,即存在单位根。 (或因为p=0.7863>α=0.05,所以接受H 0,即存在单位根。) 又因为@TREND(“1”)的p=0.1462>α=0.05,所以接受H 0,即c 1显著为0。 则模型改为011122t t t t Y c Y Y Y λββ---?=++?+?

单位根检验内容及标准规定样式分析

第八章 单位根检验 由于非平稳过程可能存在严重的伪回归问题,所以在对序列进行估计之前,需要检验序列的平稳性。本章介绍了严格的平稳性的统计检验方法--单位根检验。在简要介绍四种主要的非平稳随机过程以产输出单位根检验原理之后,文章主要介绍ADF 检验及PP 检验法,以及介结构突变和单位根检验。 8.1 四种典型非平稳过程简介 前面我们知道,若一个时间序列含有某种变动趋势,即该序列的均值或自协方差函数随时间而改变,则称该序列为非平稳序列。下面介绍四种典型的非平稳过程。 8.1.1随机游走过程 t t t y y ξ+=-1,t=1,2,... (8.11) 若}{t ξ为独立随机分布,即()0=t E ξ,()∞<=2σξt D 。则称}{t y 为随机游走过程(Random Walk Process )。随机游动过程是单位根过程的特例。在现实经济社会中,如股票价格的走势便是随机游走序列。下图是t t t y y ξ+=-1, ()1,0∈t ξ生成的序列。

图8.11 随机游走过程t t t y y ξ+=-1,()1,0∈t ξ生成的序列图 8.1.2随机趋势过程 t t t y y ξα++=-1,),0(2 σξIID t ∈, (8.12) 其中α称为漂移项,由于序列一阶差分后便趋于平稳,又称随机趋势过程为差分平稳过程。 图8.12 t t t y y ξ++=-11.0,()1,0∈t ξ生成的序列 8.1.3趋势平稳过程 t t t y ξβα++= ,其中t t t νρξξ+=-1,1<ρ,),0(2σν∈t (8.13) 由于t t t y ξαβ+=-,即当减去退势后为平稳过程,故趋势平稳过程又称为退势平稳过程。 由t t t y ξβα++=,t t t νρξξ+=-1知: 11)1(--+-+=t t t y ξβα (8.14) 将(4)两边同时乘以ρ,与(3)两边同时相减,整理可得: t t t y t y νρβα+++=-1'' , ),0(2σν∈t (8.15) 其中,ρβρααα+-=',ρβρβ-=' 这样便得出趋势平稳过程的另一种形式。

面板数据分析简要步骤与注意事项 面板单位根—面板协整—回归分析

面板数据分析简要步骤与注意事项 (面板单位根—面板协整—回归分析)步骤一:分析数据的平稳性(单位根检验) 按照正规程序,面板数据模型在回归前需检验数据的平稳性。李子奈曾指出,一些非平稳的经济时间序列往往表现出共同的变化趋势,而这些序列间本身不一定有直接的关联,此时,对这些数据进行回归,尽管有较高的R平方,但其结果是没有任何实际意义的。这种情况称为称为虚假回归或伪回归(spurious regression)。他认为平稳的真正含义是:一个时间序列剔除了不变的均值(可视为截距)和时间趋势以后,剩余的序列为零均值,同方差,即白噪声。因此单位根检验时有三种检验模式:既有趋势又有截距、只有截距、以上都无。 因此为了避免伪回归,确保估计结果的有效性,我们必须对各面板序列的平稳性进行检验。而检验数据平稳性最常用的办法就是单位根检验。首先,我们可以先对面板序列绘制时序图,以粗略观测时序图中由各个观测值描出代表变量的折线是否含有趋势项和(或)截距项,从而为进一步的单位根检验的检验模式做准备。 单位根检验方法的文献综述:在非平稳的面板数据渐进过程中,Levin andLin(1993) 很早就发现这些估计量的极限分布是高斯分布,这些结果也被应用在有异方差的面板数据中,并建立了对面板单位根进行检验的早期版本。后来经过Levin et al. (2002)的改进,提出了检验面板单位根的LLC 法。Levin et al. (2002) 指出,该方法允许不同截距和时间趋势,异方差和高阶序列相关,适合于中等维度(时间序列介于25~250 之间,截面数介于10~250 之间) 的面板单位根检验。Im et al. (1997) 还提出了检验面板单位根的IPS 法,但Breitung(2000) 发现IPS 法对限定性趋势的设定极为敏感,并提出了面板单位根检验的Breitung 法。Maddala and Wu(1999)又提出了ADF-Fisher和PP-Fisher面板单位根检验方法。 由上述综述可知,可以使用LLC、IPS、Breintung、ADF-Fisher 和PP-Fisher5种方法进行面板单位根检验。其中LLC-T 、BR-T、IPS-W 、ADF-FCS、PP-FCS 、H-Z 分别指Levin, Lin & Chu t* 统计量、Breitung t 统计量、lm Pesaran & Shin W 统计量、ADF- Fisher Chi-square统计量、PP-Fisher Chi-square统计量、Hadri Z统计量,并且Levin, Lin & Chu t* 统计量、Breitung t统计量的原假设为存在普通的单位根过程,lm Pesaran & Shin W 统计量、ADF- Fisher Chi-square统计量、PP-Fisher Chi-square统计量的原假设为存在有效的单位根过程,Hadri Z 统计量的检验原假设为不存在普通的单位根过程。 有时,为了方便,只采用两种面板数据单位根检验方法,即相同根单位根检验LLC(Levin-Lin-Chu)检验和不同根单位根检验Fisher-ADF检验(注:对普通序列(非面板序列)的单位根检验方法则常用ADF检验),如果在两种检验中均拒绝存在单位根的原假设则我们说此序列是平稳的,反之则不平稳。如果我们以T(trend)代表序列含趋势项,以I(intercept)代表序列含截距项,T&I代表两项都含,N(none)代表两项都不含,那么我们可以基于前面时序图得出的结论,在单位根检验中选择相应检验模式。 但基于时序图得出的结论毕竟是粗略的,严格来说,那些检验结构均

面板数据的单位根检验

;. 面板数据的单位根检验 1 LLC (Levin-Lin-Chu ,2002)检验(适用于相同根(common root )情形) LLC 检验原理是仍采用ADF 检验式形式。但使用的却是it y ?和it y 的剔出自相关和确定项影响的、标准的代理变量。具体做法是(1)先从? y it 和y it 中剔出自相关和确定项的影响,并使 其标准化,成为代理变量。(2)用代理变量做ADF 回归,*?ij ε=ρ*ij ε% + v it 。LLC 修正的?()t ρ 渐近服从N(0,1)分布。 详细步骤如下: H 0: ρ = 0(有单位根); H 1: ρ < 0。LLC 检验为左单端检验。 LLC 检验以如下ADF 检验式为基础: ? y it = ρ y i t -1 +∑=i k j j i 1γ? y i t -j + Z it 'φ + εit , i = 1, 2, …, N ; t = 1, 2, …, T (38) 其中Z it 表示外生变量(确定性变量)列向量,φ 表示回归系数列向量。 (1)估计代理变量。首先确定附加项个数k i ,然后作如下两个回归式, ? y it = ∑=i k j j i ? 1 γ? y i t -j + Z it '?φ +t i ε?

;. y i t -1 = ∑=i k j j i ~1 γ ? y i t -j + Z it 'φ%+1 ~-it ε 移项得 t i ε ?= ? y it -∑=i k j j i ?1 γ? y i t -j - Z it '?φ 1 ~-it ε= y it -∑=i k j j i ~1 γ? y i t -j - Z it 'φ% 把t i ε?和1 ~-it ε标准化, * ?ij ε= t i ε?/s i *ij ε%= 1~-it ε/s i 其中s i , i = 1, 2, …, N 是用(38)式对每个个体回归时得到的残差的标准差,从而得到? y it 和y it -1 的代理变量*?ij ε和* ij ε%。

面板数据分析简要步骤与注意事项面板单位根面板协整回归分析

面板数据分析简要步骤与注意事项 面板单位根—面板协整—回归分析) 步骤一:分析数据的平稳性(单位根检验) 按照正规程序,面板数据模型在回归前需检验数据的平稳性。李子奈曾指出,一些非平稳的经济时间序列往往表现出共同的变化趋势,而这些序列间本身不一定有直接的关联,此时,对这些数据进行回归,尽管有较高的R平方,但其结果是没有任何实 际意义的。这种情况称为称为虚假回归或伪回归( spurious regression )。他认为平稳的真正含义是:一个时间序列剔除了不变的均值(可视为截距)和时间趋势以后,剩余的序列为零均值,同方差,即白噪声。因此单位根检验时有三种检验模式:既有趋势又有截距、只有截距、以上都无。 因此为了避免伪回归,确保估计结果的有效性,我们必须对各面板序列的平稳性进行检验。而检验数据平稳性最常用的办法就是单位根检验。首先,我们可以先对面板序列绘制时序图,以粗略观测时序图中由各个观测值描出代表变量的折线是否含有趋势项和(或)截距项,从而为进一步的单位根检验的检验模式做准备。单位根检验方法的文献综述:在非平稳的面板数据渐进过程中 ,Levin andLin(1993) 很早就发现这些估计量的极限分布是高斯分布 , 这些结果也被应用在有异方差的面板数据中,并建立了对面板单位根进行检验的早期版本。后来经过Levin et al. (2002) 的改进, 提出了检验面板单位根的LLC法。Levin et al. (2002)指出,该方法允许不同截距和时间趋势,异方差和高阶序列相关,适合于中等维度(时间序列介于25?250之间,截面数介于10?250之间)的面板单位根检验。Im et al. (1997) 还提出了检验面板单位根的 IPS 法, 但 Breitung(2000) 发现 IPS 法对限定性趋势的设定极为敏感 , 并提出了面板单位根检验的 Breitung 法。Maddala and Wu(1999)又提出了 ADF-Fisher 和 PP-Fisher 面板单位根检验方法。 由上述综述可知,可以使用 LLC、IPS、Breintung 、ADF-Fisher 和 PP-Fisher5 种方法进行面板单位根检验。其中LLC-T 、BR-T、IPS-W 、ADF-FCS、PP-FCS、H-Z 分 别指 Levin, Lin & Chu t* 统计量、 Breitung t 统计量、 lm Pesaran & Shin W 统 量、计 ADF- Fisher Chi-square 统计量、PP-Fisher Chi-square 统计量、Hadri Z 统计 量,并且 Levin, Lin & Chu t* 统计量、 Breitung t 统计量的原假设为存在普通的单位根过程, lm Pesaran & Shin W 统计量、 ADF- Fisher Chi-square 统计量、 PP-Fisher Chi-square 统计量的原假设为存在有效的单位根过程, Hadri Z 统计量的检验原假设为不存在普通的单位根过程。 有时,为了方便,只采用两种面板数据单位根检验方法,即相同根单位根检验 LLC(Levin-Lin-Chu )检验和不同根单位根检验 Fisher-ADF 检验(注:对普通序列(非面板序列)的单位根检验方法则常用 ADF检验),如果在两种检验中均拒绝存在单位根的原假设则我 们说此序列是平稳的,反之则不平稳。 如果我们以 T(trend )代表序列含趋势项,以 I (intercept )代表序列含截距项, T&I 代表两项都含,N (none)代表两项都不含,那么我们可以基于前面时序图得出的结论,在单位根检验中选择相应检验模式。 但基于时序图得出的结论毕竟是粗略的,严格来说,那些检验结构均需一一检验。具体操作可以参照李子奈的说法:ADF检验是通过三个模型来完成,首先从含有截距和趋势项的模型开始,再检验只含截距项的模型,最后检验二者都不含的模型。并且认

单位根检验

一、单位根检验 面板数据增强了稳定性,但是也需要进行单位根检验。 面板数据单位根检验有四种方法: 1、LLC检验需要安装命令search levinlin, net ,要求各截面单元具有同质性, H0:具有单位根 命令:levinlin varname ,lags(n) 2、IPS检验安装命令search ipshin, net,各截面存在异质单位根 H0:具有单位根 命令:ipshinvarname ,lags(n) 3、fisher ADF检验 命令:xtfishervarname ,lags(n) 对统计量样本容量和滞后期较为稳健,并且适用于非平衡面板数据 4、fisher PP检验 命令:xtfishervarname ,lags(n) pp N较大时必须对P进行修正,即为fisher PP test 以上各种,还可以加入trend,时间趋势项。加入存在单位根需要差分后再检验。差分即D.varname 注意:以上各种在使用前均需要xtset设置好面板数据。 help xtunitroot 默认带有截距项 二、协整检验 1、在Stata中对面板数据进行协整检验的命令是xtwest, 命令安装ssc install xtwest 命令:xtwestdepvarvarlist [if exp] [in range] , lags(# [#]) leads(# [#]) 具体使用时可以help 通过了协整检验,说明变量之间存在着长期稳定的均衡关系,其方程回归残差是平稳的。因此可以在此基础上直接对原方程进行回归,此时的回归结果是较精确的。 三、长面板的处理 长面板N相对较小,T相对较大,扰动项不一定服从iid分布,需要估计扰动项的具体形式,然后使用广义最小二乘法(FGLS)进行估计。长面板数据关注的焦点在于设定扰动项相关的具体形式,用于提高估计的效率。在对长面板估计时需要确定是否存在异方差或者自相关,因此需要进行检验。 1、组间异方差的检验 quietlyxtglsladdindu L.lofdi huil other ,igls panel(het)est store hetero quietlyxtgls laddindL.lofdihuil other ,iglsest store homo localdf=e(N_g)-1lrtest hetero homo,df(`df') 2)xttest3也用于组间异方差的检验。 2、组内自相关的检验

单位根与协整检验

一、单位根检验的回顾 1、在实际应用中,何种情况下需要对单位根进行检验? 答:理论上,你在实际应用过程中,如果你遇到的样本是时间序列形式的,都要进行单位根检验。原因是,如果你的时间序列数据是单位根的话,类似于你的数据的变化是很不规则的,好像一个“醉汉”。从计量角度看,它影响了我们假设检验当中的“仪器”的准确性。 2、单位根检验的数学形式,或说你应当用数学方式会表述单位根检验的原假设。 3、学会在eviews上对一个时间序列变量进行单位根检验。 (1)如果一个变量具有单位根的特征,那么表示这个变量经过一次差分,就会变成平稳的。 (2)在eviews中,单位根检验的对象是series object。也就是,你要先打开一个series object,然后,在打开的窗口中点击view来观察这个序列是否具有单位根的特征。(3)要特别注意的是,eviews上如果你不

能拒绝你所检验的变量对象是一个单位根,那么此时并不一定表明你所检验的变量一定是I(1),也可能是I(2)或I(3)等更高阶的单整。要注意的是,只要你检验的变量是非平稳的,都会接受原假设。 (4)在eveiws单位根检验要遵循如下的步骤:第一,先对变量(比如Y)进行水平数据的单位根检验(level);第二,如果水平数据拒绝原假设(即不存在单位根),那么检验停止,说明水平数据是一个平稳的时间序列变量;第三,如果水平数据的检验接受原假设,仅能说明你检验的变量是非平稳的,此时需要继续对这个变量的一阶差分进行单位根检验(1S difference)。如果此时拒绝原假设,那么,检验停止,表明这个变量要经过两次差分才会平稳,否则,继续对二阶差分进行单位根检验(1S difference)。总之,检验的目的是判断,到底你所检验的变量经过几次差分后才会平稳?所以,检验一定要到差分平稳后为止。 (5)对你而言,由于有不同的单位根检验方法,所以一个不错的选择是,你同时用不

单位根过程和单位根检验

第二章单位根过程和单位根检验 第一节单位根过程 从本章开始我们进入时间序列的非平稳分析和建模研究。前面的章 节的内容主要考虑的是平稳时间序列的建模和预测问题,但对于非平 稳的时间序列,只有先进行差分处理,将其转换为平稳的时间序列模 型。这样会损失部分信息。本章从理论上介绍非平稳时间序列的性质, 讨论非平稳时间序列数据建模的伪回归问题。 非平稳序列的分析建立在维纳过程(布朗运动)和泛函中心极限定 理之上。 若干定义 定义1: (1) 白噪声过程(white noise ,如图1 )。属于平稳过程。 2 Y t =也 t ?iid (0,(T ) 图3是日元兑美元差分序列(收益序列),近似于白噪声序列。 (2) 随机游走过程(random walk ,如图2)。属于非平稳过程 2 Y t =Y t-i ;t, i ?iid (0,(T ) 随机游走过程是非平稳的,这是因为: y t =y o + U i + U 2 + W u t E(y t ) = E(y 0 + U 1+ U 2+丨1( u 」= y o 2 2 — D(y t ) = D(y o + U i + U 2 + IH+U t ) = E(u i + U 2 + 1卄+U t ) = t ^一 : 定义2 :单位根过程 随机过程{y t,t = 1,2,|||}是一单位根过程,若y t =y t_i + u t = 1,2||| U t 为一平稳过程,且 E(U t )= 0,cov(U t ,U t-s )= Ms S= 0,1,2||| CT 2 =1 ) 随机游走的差分过程是平稳过程(白噪声过程)。心yt = §

时间序列单位根检验

《计量经济学》 6.采用表5.1.1中列出的1980-2013年中国居民实际可支配收入(t X )时间序列数据,分别对t X 、t X ln 、1/-t t X X 3个序列进行单位根检验。 解:对t X 、t X ln 、1/-t t X X 序列分别进行单位根检验,R 代码为:

setwd("D://计量经济学/madongfe/") w <- read.csv("22.csv",header=T) attach(w);library(lmtest);library(tseries) X <- ts(X, start = 1980) X2 <- log(X) X3 <- X[2:34]/X[1:33] par(mfrow=c(3,1)) plot(X, xlab = "时间", type="o",col=2,lwd=2,main = "Xt 序列的波动图") plot(X2, xlab = "时间", ylab = "lnXt",type = "o", col=1,lwd=2, mian = "lnXt 序 列的波动图") plot(X3, ylab = "Xt/Xt-1",type="o",col=2,lwd=2,main = "Xt/Xt-1序列的波动图") adf.test(X); adf,test(X2); adf.test(X3) 查看三个序列的波动图,看序列图是否有明显的变化趋势,若有明显趋势,则说明该序列非平稳,结合单位根检验,单位根检验的原假设为该序列非平稳。具体结果如下: 图一:3个序列的波动图 图一可见序列t X 与序列t X ln 有明显的的增长趋势,故而两序列非平稳;序列1/ t t X X 没有明显趋势,但依然无法说明该序列平稳。借助单位根检验,结果如下表所示: 表一:3个序列的单位根检验结果一览表 Xt 序列的波动图 时间 X 1980 198519901995200020052010 5000 时间 ln X t 1980 198519901995200020052010 9.010.512. 51015 202530 0.95 1.1 Xt/Xt-1序列的波动图 I ndex X t /X t - 1

单位根过程

单位根过程 1、为什么进行单位根检验 单位根检验是检验序列中是否存在单位根,因为存在单位根就是非平稳时间序列了。单位根就是指单位根过程,可以证明,序列中存在单位根过程就不平稳,会使回归分析中存在伪回归。但是进行单位根检验的序列需服从I(d)过程。当然从变量的自相关图和偏相关图也可以判断序列是否平稳,但准确度不高。而单位根检验平稳性是比较准确的,主要方法是DF检验以及ADF检验。 2、什么是单位根检验 单位根检验是针对宏观经济数据序列、货币金融数据序列中是否具有某种统计特性而提出的一种平稳性检验的特殊方法,单位根检验的方法有很多种,包括ADF检验、PP检验、NP检验等。单位根检验时间序列的单位跟研究是时间序列分析的一个热点问题。 时间序列矩特性的时变行为实际上反映了时间序列的非平稳性质。对非平稳时间序列的处理方法一般是将其转变为平稳序列,这样就可以应用有关平稳时间序列的方法来进行相应得研究。对时间序列单位根的检验就是对时间序列平稳性的检验,非平稳时间序列如果存在单位根,则一般可以通过差分的方法来消除单位根,得到平稳序列。 对于存在单位根的时间序列,一般都显示出明显的记忆性和波动的持续性,因此单位根检验是本书中有关协整关系存在性检验和序列波动持续性讨论的基础。 3、单位根过程 定义2-1 随机序列{x_t },t=1,2,…是一单位根过程,若x_t=ρx_t-1+ε,t=1,2…(1)其中ρ=1,{ε}为一平稳序列,且 E[ε]=0, V(ε)=σ<∞, Cov(ε,ε)=μ<∞这里τ=1,2…。特别地,若{ε}是独立同分布的,且E[ε]=0,V(ε)=σ<∞,则式(1)就变成一个随机游走序列,因此随机游走序列是一种最简单的单位根过程。将式(1)改写为下列形式:( 1-ρL)x_t=ε,t=1,2,…其中L 为滞后算子,1-ρL为滞后算子多项式,其特征方程为1-ρz=0,有根z= 。当ρ=1时,时间序列存在一个单位根,此时{ }是一个单位根过程。当ρ<1时,{x_t}为平稳序列。而当ρ〉1时,{x_t}为一类具有所谓爆炸根的非平稳过程,它经过差分后仍然为非平稳过程,因此不为单整过程。一般情况下,单整过程可以称作单位根过程。在经济、金融时间序列中,常会遇到ρ非常接近1的情况,成为近似单位根现象。近似单位根是介于平稳序列I(0)和单正序列I(1)之间。一般情况下,单整过程可以称作单位根过程。 4、单位根检验的基础 单位根检验是建立ARMA模型、ARIMA模型、变量间的协整分析、因果关系检验等的基础。自Nelson和Plosser利用ADF检验研究了美国名义GNP等14 个历史经济和金融时间序列的平稳性以后,单位根检验业已成为分析经济和金融时间序列变化规律和预测的重要组成部分。因此,单位根检验作为一种特殊的假设检验,其可靠性的研究以及如何寻求可靠性较高的检验方法或统计量多年来一直是时间序列分析中的重要课题。

单位根检验

单位根检验、协整检验和格兰杰因果关系检验三者 之间的关系 实证检验步骤:先做单位根检验,看变量序列是否平稳序列,若平稳,可构造回归模型等经典计量经济学模型;若非平稳,进行差分,当进行到第i次差分时序列平稳,则服从i阶单整(注意趋势、截距不同情况选择,根据P值和原假设判定)。若所有检验序列均服从同阶单整,可构造VAR模型,做协整检验(注意滞后期的选择),判断模型内部变量间是否存在协整关系,即是否存在长期均衡关系。如果有,则可以构造VEC模型或者进行Granger因果检验,检验变量之间“谁引起谁变化”,即因果关系。 一、讨论一 1、单位根检验是序列的平稳性检验,如果不检验序列的平稳性直接OLS容易导致伪回归。 2、当检验的数据是平稳的(即不存在单位根),要想进一步考察变量的因果联系,可以采用格兰杰因果检验,但要做格兰杰检验的前提是数据必须是平稳的,否则不能做。 3、当检验的数据是非平稳(即存在单位根),并且各个序列是同阶单整(协整检验的前提),想进一步确定变量之间是否存在协整关系,可以进行协整检验,协整检验主要有EG两步法和JJ检验 A、EG两步法是基于回归残差的检验,可以通过建立OLS模型检验其残差平稳性 B、JJ检验是基于回归系数的检验,前提是建立VAR模型(即模型符合ADL模式)

4、当变量之间存在协整关系时,可以建立ECM进一步考察短期关系,Eviews 这里还提供了一个Wald-Granger检验,但此时的格兰杰已经不是因果关系检验,而是变量外生性检验,请注意识别 二、讨论二 1、格兰杰检验只能用于平稳序列!这是格兰杰检验的前提,而其因果关系并非我们通常理解的因与果的关系,而是说x的前期变化能有效地解释y的变化,所以称其为“格兰杰原因”。 2、非平稳序列很可能出现伪回归,协整的意义就是检验它们的回归方程所描述的因果关系是否是伪回归,即检验变量之间是否存在稳定的关系。所以,非平稳序列的因果关系检验就是协整检验。 3、平稳性检验有3个作用:1)检验平稳性,若平稳,做格兰杰检验,非平稳,作协正检验。2)协整检验中要用到每个序列的单整阶数。3)判断时间学列的数据生成过程。 三、讨论三 其实很多人存在误解。有如下几点,需要澄清: 第一,格兰杰因果检验是检验统计上的时间先后顺序,并不表示而这真正存在因果关系,是否呈因果关系需要根据理论、经验和模型来判定。 第二,格兰杰因果检验的变量应是平稳的,如果单位根检验发现两个变量是不稳定的,那么,不能直接进行格兰杰因果检验,所以,很多人对不平稳的变量进行格兰杰因果检验,这是错误的。 第三,协整结果仅表示变量间存在长期均衡关系,那么,到底是先做格兰杰还是先做协整呢?因为变量不平稳才需要协整,所以,首先因对变量进行差分,平稳后,可

单位根检验

平稳性的单位根检验:DF检验、ADF检验、DFGLS检验、PP检验、KPSS检验、ERS检验和NP检验 (2011-12-21 12:13:27) ADF检验 作用 检查序列平稳性的标准方法是单位根检验。有6种单位根检验方法:ADF检验、DFGLS检验、PP检验、KPSS检验、ERS检验和NP检验,本节将介绍DF检验、ADF 检验。 比较 ADF检验和PP检验方法出现的比较早,在实际应用中较为常见,但是,由于这2种方法均需要对被检验序列作可能包含常数项和趋势变量项的假设,因此,应用起来带有一定的不便;其它几种方法克服了前2种方法带来的不便,在剔除原序列趋势的基础上,构造统计量检验序列是否存在单位根,应用起来较为方便。 来源 ADF检验是在Dickey-Fuller检验(DF检验)基础上发展而来的。因为DF检验只有当序列为AR(1)时才有效。如果序列存在高阶滞后相关,这就违背了扰动项是独立同分布的假设。在这种情况下,可以使用增广的DF检验方法(augmented Dickey-Fuller test )来检验含有高阶序列相关的序列的单位根。 步骤 一般进行ADF检验要分3步: 1 对原始时间序列进行检验,此时第二项选level,第三项选None.如果没通过检验,说明原始时间序列不平稳; 2 对原始时间序列进行一阶差分后再检验,即第二项选1st difference,第三项选intercept,若仍然未通过检验,则需要进行二次差分变换; 3 二次差分序列的检验,即第二项选择2nd difference ,第四项选择Trend and intercept.一般到此时间序列就平稳了! 在进行ADF检验时,必须注意以下两个实际问题: (1)必须为回归定义合理的滞后阶数,通常采用AIC准则来确定给定时间序列模型的滞后阶数。在实际应用中,还需要兼顾其他的因素,如系统的稳定性、模型的拟合优度等。 (2)可以选择常数和线性时间趋势,选择哪种形式很重要,因为检验显著性水平的 t 统计量在原假设下的渐近分布依赖于关于这些项的定义。 ①若原序列中不存在单位根,则检验回归形式选择含有常数,意味着所检验的序列的均值不为0;若原序列中存在单位根,则检验回归形式选择含有常数,意味着所检验的序列具有线性趋势,一个简单易行的办法是画出检验序列的曲线

单位根检验

单位根检验、协整检验和格兰杰因果关系检验三者之间的关系 实证检验步骤:先做单位根检验,看变量序列是否平稳序列,若平稳,可构造回归模型等经典计量经济学模型;若非平稳,进行差分,当进行到第i次差分时序列平稳,则服从i阶单整(注意趋势、截距不同情况选择,根据P值和原假设判定)。若所有检验序列均服从同阶单整,可构造VAR模型,做协整检验(注意滞后期的选择),判断模型内部变量间是否存在协整关系,即是否存在长期均衡关系。如果有,则可以构造VEC模型或者进行Granger因果检验,检验变量之间“谁引起谁变化”,即因果关系。 一、讨论一 1、单位根检验是序列的平稳性检验,如果不检验序列的平稳性直接OLS容易导致伪回归。 2、当检验的数据是平稳的(即不存在单位根),要想进一步考察变量的因果联系,可以采用格兰杰因果检验,但要做格兰杰检验的前提是数据必须是平稳的,否则不能做。 3、当检验的数据是非平稳(即存在单位根),并且各个序列是同阶单整(协整检验的前提),想进一步确定变量之间是否存在协整关系,可以进行协整检验,协整检验主要有EG两步法和JJ检验 A、EG两步法是基于回归残差的检验,可以通过建立OLS模型检验其残差平稳性 B、JJ检验是基于回归系数的检验,前提是建立VAR模型(即模型符合ADL模式) 4、当变量之间存在协整关系时,可以建立ECM进一步考察短期关系,Eviews这里还提供了一个Wald-Granger检验,但此时的格兰杰已经不是因果关系检验,而是变量外生性检验,请注意识别 二、讨论二 1、格兰杰检验只能用于平稳序列!这是格兰杰检验的前提,而其因果关系并非我们通常理解的因与果的关系,而是说x的前期变化能有效地解释y的变化,所以称其为“格兰杰原因”。 2、非平稳序列很可能出现伪回归,协整的意义就是检验它们的回归方程所描述的因果关系是否是伪回归,即检验变量之间是否存在稳定的关系。所以,非平稳序列的因果关系检验就是协整检验。 3、平稳性检验有3个作用:1)检验平稳性,若平稳,做格兰杰检验,非平稳,作协正检验。2)协整检验中要用到每个序列的单整阶数。3)判断时间学列的数据生成过程。

单位根检验详解

第2节 单位根检验 由于虚假回归问题的存在,因此检验变量的平稳性是一个必须解决的问题。在第十二章中介绍用相关图判断时间序列的平稳性。这一章则给出序列平稳性的严格的统计检验方法,即单位根检验。单位根检验有很多方法,这里主要介绍DF 和ADF 检验。 序列均值为0则无C ,序列无时间趋势则无trend 在介绍单位根检验之前,先认识四种典型的非平稳随机过程。 1、四种典型的非平稳随机过程 (1)随机游走过程。 y t = y t -1 + u t , y 0 = 0, u t ~ IID(0, σ 2) 其均值为零,方差无限大(?),但不含有确定性时间趋势。(见图1a )。 -10 -5 5 10 20 40 60 140160y=y(-1)+u 1200 1400 1600 1800 2000 2200 图1a 由y t = y t -1+ u t 生成的序列 图1b 深证成指 (2)随机趋势过程。 y t = α + y t -1 + u t , y 0 = 0, u t ~ IID(0, σ 2) 其中α称作位移项(漂移项)。由上式知,E(y 1)= α(过程初始值的期望)。将上式作如下迭代变换, y t = α + y t -1 + u t = α+ (α+ y t -2 + u t -1) + u t = … = αt +y 0 +∑-t i i u 1

y t 由确定性时间趋势项αt 和y 0 +∑-t i i u 1 组成。可以把y 0 +∑-t i i u 1 看作随机 的截距项。在不存在任何冲击u t 的情况下,截距项为y 0。而每个冲击u t 都表现为截距的移动。每个冲击u t 对截距项的影响都是持久的,导致序列的条件均值发生变化,所以称这样的过程为随机趋势过程(stochastic trend process ),或有漂移项的非平稳过程(non-stationary process with drift ),见图2,虽然总趋势不变,但随机游走过程围绕趋势项上下游动。由上式还可以看出,α是确定性时间趋势项的系数(原序列y t 的增长速度)。α为正时,趋势向上;α为负时,趋势向下。 20 40 60 80 stochastic trend process -100 -80-60-40-20 020 1002003004005006007008009001000 y=-0.1+y(-1)+u 图2a 由y t =0.1+ y t -1+ u t 生成的序列 图2b 由y t =- 0.1+ y t -1+ u t 生成的序列 因为对y t 作一次差分后,序列就平稳了, ? y t = y t - y t -1 = α + u t (平稳过程) 所以也称y t 为差分平稳过程(difference- stationary process )。α是? y t 序列的均值,原序列y t 的增长速度。 (3)趋势平稳过程 y t = β0 + β1 t + u t , u t = ρu t -1 + v t , (ρ <1, v t ~ IID(0, σ2)) y t 与趋势值 β0+β1t 不同,差值为u t 。因为u t 是平稳的,y t 只会

单位根检验、协整检验和格兰杰因果

实证检验步骤:先做单位根检验,看变量序列是否平稳序列,若平稳,可构造回归模型等经典计量经济学模型;若非平稳,进行差分,当进行到第i次差分时序列平稳,则服从i阶单整(注意趋势、截距不同情况选择,根据P值和原假设判定)。若所有检验序列均服从同阶单整,可构造V AR模型,做协整检验(注意滞后期的选择),判断模型内部变量间是否存在协整关系,即是否存在长期均衡关系。如果有,则可以构造VEC模型或者进行Granger因果检验,检验变量之间“谁引起谁变化”,即因果关系。 一、讨论一 1、单位根检验是序列的平稳性检验,如果不检验序列的平稳性直接OLS容易导致伪回归。 2、当检验的数据是平稳的(即不存在单位根),要想进一步考察变量的因果联系,可以采用格兰杰因果检验,但要做格兰杰检验的前提是数据必须是平稳的,否则不能做。 3、当检验的数据是非平稳(即存在单位根),并且各个序列是同阶单整(协整检验的前提),想进一步确定变量之间是否存在协整关系,可以进行协整检验,协整检验主要有EG两步法和JJ检验 A、EG两步法是基于回归残差的检验,可以通过建立OLS模型检验其残差平稳性 B、JJ检验是基于回归系数的检验,前提是建立V AR模型(即模型符合ADL模式) 4、当变量之间存在协整关系时,可以建立ECM进一步考察短期关系,Eviews这里还提供了一个Wald-Granger检验,但此时的格兰杰已经不是因果关系检验,而是变量外生性检验,请注意识别 二、讨论二 1、格兰杰检验只能用于平稳序列!这是格兰杰检验的前提,而其因果关系并非我们通常理解的因与果的关系,而是说x的前期变化能有效地解释y的变化,所以称其为“格兰杰原因”。

面板数据的单位根检验

面板数据的单位根检验 1 LLC (Levin-Lin-Chu ,2002)检验(适用于相同根(common root )情形) LLC 检验原理是仍采用ADF 检验式形式。但使用的却是it y ?和it y 的剔出自相关和确定项影响的、标准的代理变量。具体做法是(1)先从? y it 和y it 中剔出自相关和确定项的影响,并使 其标准化,成为代理变量。(2)用代理变量做ADF 回归,*?ij ε=ρ*ij ε% + v it 。LLC 修正的?()t ρ 渐近服从N(0,1)分布。 详细步骤如下: H 0: ρ = 0(有单位根); H 1: ρ < 0。LLC 检验为左单端检验。 LLC 检验以如下ADF 检验式为基础: ? y it = ρ y i t -1 +∑=i k j j i 1γ? y i t -j + Z it 'φ + εit , i = 1, 2, …, N ; t = 1, 2, …, T (38) 其中Z it 表示外生变量(确定性变量)列向量,φ 表示回归系数列向量。 (1)估计代理变量。首先确定附加项个数k i ,然后作如下两个回归式, ? y it = ∑=i k j j i ? 1 γ? y i t -j + Z it '?φ +t i ε?

y i t -1 = ∑=i k j j i ~1 γ ? y i t -j + Z it 'φ%+1 ~-it ε 移项得 t i ε ?= ? y it -∑=i k j j i ?1 γ? y i t -j - Z it '?φ 1 ~-it ε= y it -∑=i k j j i ~1 γ? y i t -j - Z it 'φ% 把t i ε?和1 ~-it ε标准化, * ?ij ε= t i ε?/s i *ij ε%= 1~-it ε/s i 其中s i , i = 1, 2, …, N 是用(38)式对每个个体回归时得到的残差的标准差,从而得到? y it 和y it -1 的代理变量*?ij ε和* ij ε%。 (2)用代理变量*?ij ε和* ij ε%作如下回归,

相关文档
最新文档