动态面板GMM回归结果解读
一、方法概述
动态面板数据模型是在静态面板模型的基础上引入被解释变量的滞后项作为解释变量,从而刻画被解释变量的动态调整过程。在传统的固定效应或随机效应面板模型中,假设解释变量与被解释变量之间不存在动态反馈机制,但在实际经济和社会现象中,被解释变量往往受到其自身历史值的影响,即存在状态依赖性。例如,企业当前的绩效往往与其前期的绩效水平密切相关,经济增长具有惯性特征等。此时若仍使用传统面板模型进行估计,将因遗漏滞后被解释变量而导致模型设定偏误。
动态面板模型的估计方法主要包括差分GMM(Difference GMM)和系统GMM(System GMM)两种。差分GMM由Arellano和Bond(1991)提出,通过对模型进行一阶差分消除个体固定效应,然后使用被解释变量和解释变量的滞后水平值作为差分方程的工具变量,采用广义矩估计(GMM)方法进行参数估计。系统GMM由Arellano和Bover(1995)以及Blundell和Bond(1998)进一步发展,将水平方程与差分方程联立构成方程组进行估计,在提高估计效率的同时能够估计不随时间变化的变量的系数。本研究使用SPSSAU软件,采用系统GMM(TwoStep)方法对动态面板模型进行估计,并使用Robust标准误以提高推断的可靠性。
在SPSSAU【计量经济研究】模块选择【GMM估计】,将变量拖拽至右侧对应分析框,操作如下图:
二、模型基本信息
本研究所构建的动态面板模型基本信息如表1所示。
从表1可以看出,本研究采用系统GMM方法,差分类型为两步系统GMM(TwoStep),相较于一步GMM,两步GMM在有限样本下具有更好的小样本性质,能够提供更有效的估计。模型共包含595个组别(个体),每组平均样本量为5个时间期,时间项数为7。模型共使用49个工具变量,其中被解释变量Y使用其2期及更早期的滞后值作为GMM式工具变量,各解释变量X1至X6使用其1期及更早期的滞后值作为GMM式工具变量,同时对工具变量进行了压缩处理(collapse)以避免工具变量过多导致的过度拟合问题。模型的信息准则值为AIC=27.057、BIC=-224.859、HQIC=-72.327,可作为模型比较的参考依据。
三、动态面板模型估计结果
动态面板系统GMM模型的回归估计结果如表2所示。
从表2的估计结果可以看出,各变量的回归系数及其显著性情况如下:
(1)被解释变量滞后项L1.Y:回归系数为0.720,在1%的水平上高度显著(z=9.243,p=0.000),95%置信区间为[0.567, 0.873]。这一结果表明,被解释变量Y存在显著的正向动态效应,即前一期的Y值每增加1个单位,当期的Y值平均增加0.720个单位。该系数高达0.720,说明Y具有较强的惯性特征和持续性,前期水平对当期水平具有很强的解释力。这一结果也验证了使用动态面板模型的必要性——如果忽略这一滞后效应,模型将存在严重的设定偏误。
(2)X1:回归系数为0.022,在1%的水平上显著(z=3.071,p=0.002),95%置信区间为[0.008, 0.035]。这表明在控制其他变量和动态效应后,X1每增加1个单位,Y平均增加0.022个单位,X1对Y具有显著的正向影响。
(3)X2:回归系数为0.003,在5%的水平上显著(z=2.454,p=0.014),95%置信区间为[0.001, 0.005]。X2对Y具有显著的正向影响,但影响幅度相对较小,X2每增加1个单位,Y仅增加0.003个单位。
(4)X3:回归系数为-0.098,在5%的水平上显著(z=-2.269,p=0.023),95%置信区间为[-0.182, -0.013]。X3对Y具有显著的负向影响,即X3每增加1个单位,Y平均减少0.098个单位。该变量在所有显著变量中影响幅度最大。
(5)X4:回归系数为0.021,z=0.491,p=0.624>0.05,未达到统计显著性水平。这表明在控制其他变量后,X4对Y的影响不具有统计显著性,不能拒绝X4对Y无影响的零假设。
(6)X5:回归系数为-0.356,z=-1.772,p=0.076>0.05,在5%的水平上不显著,但在10%的水平上具有边际显著性。X5对Y呈现负向影响趋势,但证据尚不充分,可能需要更大样本量来确认其效应。
(7)X6:回归系数为-0.075,z=-0.780,p=0.435>0.05,未达到统计显著性水平。X6对Y的影响不显著。
(8)常数项:回归系数为1.574,在1%的水平上显著(z=3.500,p=0.000),95%置信区间为[0.693, 2.455]。
四、Hansen过度识别检验
在GMM估计中,工具变量的有效性是保证估计结果一致性的关键前提。Hansen过度识别检验用于检验所有工具变量是否满足正交性条件,即工具变量是否与误差项不相关。该检验的原假设为[工具变量与误差项不相关],若p值大于0.05则接受原假设,认为工具变量有效;若p值小于0.05则拒绝原假设,意味着工具变量与误差项存在相关性,模型设定可能存在问题。
从表4可知,Hansen过度识别检验的Wald χ2值为111.057,自由度为41,对应的p值为0.000,小于0.05的显著性水平,因此拒绝原假设。这意味着工具变量与误差项之间存在显著相关性,部分工具变量可能不满足外生性条件。该结果提示当前模型的工具变量选择可能需要进一步优化,例如调整工具变量的滞后阶数、减少工具变量数量或重新审视模型设定。需要注意的是,Hansen检验在工具变量数量较多时可能出现过度拒绝的情况,本研究使用了49个工具变量,这在一定程度上可能影响检验的可靠性。
五、Arellano-Bond自相关检验
Arellano-Bond AR检验用于检验差分后模型的残差是否存在序列自相关。对于动态面板GMM模型,差分后的残差允许存在一阶自相关AR(1),但不应存在二阶自相关AR(2)。如果差分后残差存在二阶自相关,则说明原始水平方程的误差项存在序列相关,这将导致GMM估计量不一致。AR检验的原假设为[模型不存在自相关]。
从表5可知,AR(1)检验的z值为-4.666,p值为0.000<0.05,拒绝一阶无自相关的原假设,表明差分后残差存在一阶自相关。这一结果是预期之内的,因为动态面板模型中由于包含了被解释变量的滞后项,差分后残差必然存在一阶自相关,这并不影响GMM估计的一致性。
AR(2)检验的z值为2.199,p值为0.028<0.05,拒绝二阶无自相关的原假设,表明差分后残差存在二阶自相关。这一结果意味着原始水平方程的误差项存在序列自相关,可能导致GMM估计量不一致,模型构建存在一定问题。结合Hansen检验的结果(p=0.000),两项诊断检验均提示当前模型设定有待改进。可能的改进方向包括:调整工具变量的滞后阶数范围、减少工具变量个数以避免过度识别问题、考虑使用差分GMM替代系统GMM、或引入更多的控制变量以改善模型设定。
六、结论
本研究采用两步系统GMM方法对动态面板模型进行估计,主要结论如下:第一,被解释变量的滞后一期项L1.Y的回归系数为0.720且在1%水平上高度显著,表明Y具有强烈的动态持续效应,前期水平对当期水平具有重要的解释作用,验证了使用动态面板模型的必要性。第二,在解释变量中,X1(β=0.022, p<0.01)、X2(β=0.003, p<0.05)对Y具有显著的正向影响,X3(β=-0.098, p<0.05)对Y具有显著的负向影响,而X4、X5、X6的影响未达到5%的显著性水平。第三,模型诊断检验方面,Hansen过度识别检验(p=0.000)和AR(2)检验(p=0.028)的结果均提示模型设定存在改进空间,工具变量的有效性和残差的序列独立性有待进一步优化。后续研究可通过调整工具变量滞后阶数、精简工具变量数量等方式改进模型质量。