面板向量自回归(VAR)模型在应用研究中的应用越来越多。虽然专门用于估计时间序列VAR模型的程序通常作为标准功能包含在大多数统计软件包中,但面板VAR模型的估计和推断通常用通用程序实现,需要一些编程技巧。在本文中,我们简要讨论了广义矩量法(GMM)框架下面板VAR模型的模型选择、估计和推断,并介绍了一套Stata程序来方便地执行它们。
时间序列向量自回归 (VAR) 模型起源于宏观计量经济学文献,作为多元联立方程模型的替代品 (Sims, 1980)。VAR 系统中的所有变量通常都被视为内生变量,尽管可能会根据理论模型或统计程序来确定限制,以解决外生冲击对系统的影响。随着 VAR 在面板数据设置中的引入(Holtz-Eakin、Newey 和 Rosen,1988),面板 VAR 模型已在跨领域的多个应用中使用。
在本文中,我们简要概述了广义矩量法 (GMM) 框架中面板 VAR 模型的选择、估计和推理,并提供了一组 Stata 程序,我们使用国家纵向调查和投资、收入和消费数据。包括实现Granger(1969)因果关系检验的子程序,以及按照Andrews和Lu(2001)进行的最佳时刻和模型选择。
我们考虑具有特定面板固定效应的阶数 -变量面板 VAR,由以下线性方程组表示:
其中,
上面的参数可以与固定效应联合估计,或者在一些转换后独立于固定效应,使用普通最小二乘法 (OLS)。然而,由于方程组右侧存在滞后因变量,即使是大的估计也会有偏差(尼克尔,1981)。尽管偏差随着变大而趋近于零,但 Judson 和 Owen (1999) 的模拟发现即使在 = 30 时也存在显着偏差。
已经提出了基于 GMM 的各种估计器来计算上述方程的一致估计。4在我们假设误差是连续不相关的情况下,一阶差分变换可以通过用较早时期的差异和水平检测滞后差异,如安德森和萧 (1982) 所提出的那样,逐个方程地一致估计。然而,这个估计会带来一些问题。一阶差分变换放大了不平衡面板中的间隙。例如,如果某些不可用,则时间和 − 1 处的一阶差分同样缺失。此外,观察每个面板的必要时间段随着面板 VAR 的滞后顺序而变大。例如,对于二阶面板 VAR,
Arellano 和 Bover (1995) 提出前向正交偏差作为替代变换,它不具有一阶差分变换的缺点。它不使用与过去实现的偏差,而是减去所有可用的未来观察的平均值,从而最大限度地减少数据丢失。可能只有最近的观察不会用于估计。由于过去的实现不包括在这个转换中,它们仍然是有效的工具。例如,在二阶面板 VAR 中,只有
虽然逐个方程的 GMM 估计会产生对面板 VAR 的一致估计,但将模型估计为方程组可能会导致效率增益(Holtz-Eakin、Newey 和 Rosen,1988 年)。考虑以下基于等式 (1) 的变换面板 VAR 模型,但以更紧凑的形式表示:
其中星号表示原始变量的某种变换。如果我们把原始变量表示为
假设我们随着时间的推移将观察叠加在面板上。GMM 估计量由下式给出
其中是一个 ( ) 加权矩阵,假定为非奇异、对称和半正定。假设
方程组的联合估计使交叉方程假设检验变得简单明了。可以基于
面板 VAR 分析的前提是在面板 VAR 规范和矩条件中选择最佳滞后阶数。Andrews和Lu(2001)提出了基于Hansen(1982)统计学的过度识别限制的GMM模型的一致时刻和模型选择标准(MMSC)。他们提出的 MMSC 类似于各种常用的基于最大似然的模型选择标准,即 Akaike 信息标准 (AIC) (Akaike, 1969)、贝叶斯信息标准 (BIC) (Schwarz, 1978; Rissanen, 1978; Akaike, 1977 ),以及 Hannan-Quinn 信息标准 (HQIC)(Hannan 和 Quinn,1979)。
将 Andrews 和 Lu 的 MMSC 应用 GMM 估计,他们提出的标准选择最小化的向量对
其中
通过构造,上述 MMSC 仅在
我们删除外生变量,并专注于方程(1)中面板 VAR 的自回归结构。Lutkepohl (2005) 和 Hamilton (1994) 都表明,如果伴随矩阵的所有模都严格小于 1,则 VAR 模型是稳定的,其中伴随矩阵由
稳定性意味着面板 VAR 是可逆的,并且具有无限阶向量移动平均 (VMA) 表示,为估计的脉冲响应函数和预测误差方差分解提供已知的解释。可以通过将模型重写为无限向量移动平均来计算简单的脉冲响应函数,其中
然而,简单的 IRF 没有因果解释。由于创新是同时相关的,一个变量的冲击很可能伴随着其他变量的冲击。假设我们有一个矩阵
脉冲响应函数置信区间可以基于面板 VAR 参数的渐近分布和交叉方程误差方差-协方差矩阵分析导出。或者,也可以使用蒙特卡罗模拟和自举重采样方法来估计置信区间。
ℎ步提前预测误差可以表示为
其中
在应用中,贡献通常相对于变量的ℎ步超前预测误差方差进行归一化。
与脉冲响应函数类似,置信区间可以通过分析得出或使用各种重采样技术进行估计
我们通过分析年工作时间和小时收入之间的关系来说明 pvar使用,Holtz-Eakin、Newey 和 Rosen(1988)之前在他们关于面板向量自回归的开创性论文中对此进行了分析。为了将我们的新程序与 Stata 的内置 var 命令套件进行比较,我们还将新的 pvar 应用于投资、收入和消费数据时间序列数据。
我们通过分析年工作时间和小时收入之间的关系来说明pvar使用,Holtz-Eakin、Newey和Rosen(1988)曾在他们关于面板向量自回归的开创性论文中分析过这种关系。我们还将pvar应用于Lutkephol(1993)的时间序列数据。
我们使用来自 1968 年至 1978 年国家纵向调查的 1968 年 14-26 岁女性子样本。我们的子样本由 2,039 名女性组成,她们在至少三轮调查中报告了工资(wage)和年工作时数(hours),其中两轮是连续年份。使用相同的调查,但具有不同的时间段和不同的工人子样本,因此结果可能不具有直接可比性。
下面是使用模型选择,用于以工时和工资的前四个滞后期为工具的一到三阶面板VARs。
.ua3 vs(ns(1/4))基于 Andrews 和 Lu (2001) 的三个模型选择标准和整体决定系数,一阶面板 VAR 是首选模型,因为它具有最小的 MBIC、MAIC 和 MQIC。虽然我们也想最小化 Hansen 的 J 统计量,但它并没有像 Andrews 和 Lu 的模型和矩选择标准那样修正模型中的自由度。基于选择标准,我们使用由 pvar 实现的 GMM 估计拟合具有与上述相同的一阶面板 VAR 模型。
面板向量自回归请注意,估计中包括的506名妇女明显少于数据中的全部妇女子样本。默认情况下,pvar会从估计中删除任何缺失数据的观察。由于子样本中的所有妇女的工作时间和工资并不是在所有年份都被观察到的,所以被剔除的观察值的数量会随着作为工具变量的滞后阶数而增加。我们可以通过使用Holtz-Eakin等人提出的 "GMM式 "工具来改善估计。这增加了估计样本,从而使估计更加有效。
面板向量自回归尽管可以从上面的 var 输出推断出一阶面板 VAR 的格兰杰因果关系,但我们仍然使用 granger 作为说明进行检验。下面格兰杰因果检验的结果表明,在通常的置信水平下,工资格兰杰导致工时和工时格兰杰导致工资,这与 Holtz-Eakin 等人的发现相似。
. granger面板向量自回归模型估计很少由其自身解释。在实践中,研究人员通常对面板 VAR 系统中每个内生变量的外生变化对其他变量的影响感兴趣。然而,在估计脉冲响应函数 (IRF) 和预测误差方差分解 (FEVD) 之前,我们首先检查估计面板 VAR 的稳定性条件。生成的特征值表和图证实了估计是稳定的。
.table, rph我们认为工资水平的冲击对同期的工作时间有直接影响,而当前的工作努力只影响未来的工资。使用这种因果顺序,我们使用 fevd 计算了隐含的 IRF,使用 fevd 计算了隐含的 FEVD。IRF 置信区间是根据估计模型使用 200 次蒙特卡罗绘制计算的。FEVD 估计值的标准误差和置信区间同样可用。
.根据FEVD的估计,我们看到,在我们的例子中,多达40%的妇女工作时间的变化可以由她们的工资来解释。另一方面,工作时间只解释了妇女未来工资变化的5%。就水平而言,IRF图显示,实际工资的正冲击导致工作努力的减少,这意味着样本中妇女的劳动供给向后弯曲。同样值得注意的是,当前工作强度的冲击对工作时间和工资都有积极但短暂的影响。另一方面,当前冲击对工资的影响对未来工资有持续的积极影响。
联系客服