eviews分布滞后模型和自回归模型.ppt

第六章滞后变量模型 科克分布滞后模型 科克模型 在估计的过程中存在以下问题 1 由于作为解释变量 因此模型中包含随机解释变量 2 即使原模型中的不存在序列相关 然而是序列相关的 3 解释变量和误差项存在序列相关 因此 使用OLS估计将导致估计量不仅是有偏的而且非一致的 可以采用工具变量法来估计 有学者建议用作为的工具变量 例1 table8 1 wf1工作文件中 给出的是1978 2006年北京市城镇家庭平均每人全年消费性支出 PPCE 单位元 和城镇家庭平均每人可支配收入 PPDI 单位元 由于人们消费习惯等原因 使得收入对消费支出的影响存在时间滞后 因此建立消费函数的分布滞后模型 本实验打算建立如下模型 这里以做为滞后解释变量的工具变量 虽然工具变量法可以消除科克模型中解释变量的随机性以及解释变量与误差项之间的序列相关等问题 但由于引入的工具变量是 其与存在高度相关性 因此模型估计存在多重共线性问题 这样 虽然工具变量方法给出了方程的一致性估计 但是这些估计量很可能是低效的 有限分布滞后模型 一般模型为 对于滞后长度的确定 可以根据实际经济问题的需要和经验进行判定 也可以利用一些判定方法和准则 如赤池 Akaike AIC准则与施瓦兹 Schwarz SC准则等 对于滞后长度为已知的分布滞后模型 修正的估计方法有经验加权法 阿尔蒙 Almon 多项式滞后法等 各种方法的基本思想大致相同 都是通过对各滞后变量加权 组成线性组合变量 即滞后变量的线性组合 作为新解释变量引入方程 有目的地减少滞后变量的数目 缓解多重共线性 保证自由度 1 经验加权估计法 所谓经验加权法 是根据实际经济问题的特点及经验判断 对滞后变量赋予一定的权数 利用这些权数构成各滞后变量的线性组合 以形成新的变量 再应用最小二乘法进行估计 由于随机误差项与解释变量不相关 从而也与滞后解释变量的线性组合变量不相关 因此可直接应用最小二乘法对该模型进行估计 经验加权法具有简单易行 不损失自由度 避免多重共线性干扰及参数估计具有一致性等优点 缺陷是设置权数的主观随意性较大 要求分析者对实际问题的特征有比较透彻的了解 通常的做法是 多选几组权数 分别估计多个模型 然后根据样本决定系数 F检验值 t检验值 估计标准误差以及DW值 从中选出最佳估计方程 例 已知某地区制造业部门1955 1974年期间的资本存量Y和销售额X的统计资料如下表 金额单位 百万元 设定有限分布滞后模型为 运用经验加权法 选择下列三组权数 1 1 1 2 1 4 1 8 2 1 4 1 2 2 3 1 4 3 1 4 1 4 1 4 1 4 分别估计上述模型 并从中选择最佳的方程 数据见case25 记新的线性组合变量分别为 分别估计如下经验加权模型 YT 66 52294932 1 071395456 Z1 3 662182 50 96149 R squared 0 994257DW 1 439440F 2597 074YT 133 1722303 1 366668187 Z2 5 029746 37 37033 R squared 0 989373DW 1 042713F 1396 542YT 121 7394467 2 237930494 Z3 4 813143 38 68578 R squared 0 990077DW 1 158530F 1496 590 从上述回归分析结果可以看出 模型一的扰动项无一阶自相关 模型二和模型三扰动项存在一阶正相关 在综合判断可决系数 F 检验值 t检验值 可以认为 最佳的方程式模型一 即权数为1 1 2 1 4 1 8的分布滞后模型 2 阿尔蒙法 主要思想 针对有限滞后期模型 通过阿尔蒙变换 定义新变量 以减少解释变量个数 然后用OLS法估计参数 主要步骤为 第一步 阿尔蒙变换对于分布滞后模型假定其回归系数 i可用一个关于滞后期i的适当阶数的多项式来表示 即 i 0 1 s其中 m s 1 阿尔蒙变换要求先验地确定适当阶数k 例如取m 2 得 将 代入分布滞后模型定义新变量将原模型转换为 第二步 模型的OLS估计对变换后的模型进行OLS估计 得再计算出 求出滞后分布模型参数的估计值 由于m 1 s 可以认为原模型存在的自由度不足和多重共线性问题已得到改善 需注意的是 在实际估计中 阿尔蒙多项式的阶数m一般取2或3 不超过4 否则达不到减少变量个数的目的 例 case26是某水库1998年至2000年各旬的流量 降水量数据 分别建立水库流量与降水量序列 命名为vol和ra 试对其建立多项式分布滞后模型 Eviews操作 在主窗口命令行键入如下命令建立PDL模型 Lsyx1x2pdl seriesname lags order options 其中 lags代表滞后期s order表示多项式次数m options指定约束类型 有下面三个选项 1近端约束 使x对y的一期前导作用为02远端约束 使大于滞后期p后x对y的作用为03同时采用近端和远端两种约束如果模型中没有约束条件 则options缺省 本例 假定降水量对水库流量滞后3月的影响仍然显著 即滞后期p 9 若采用4阶多项式 m 4 且不施加端点限制条件 则输入命令Lsvolcpdl ra 9 4 模型输出窗口的上半部分给出了各参数估计值及检验的t统计量 下半部分模型检验所需的各个统计量 这里用PDL01 PDL02 PDL03等代表式中的w1t w2t等变量 本例m 4 所以除常数项外共有5个参数估计值 该命令还同时绘制出估计值的分布图 表中最后一行的SumofLags是系数估计值的总和 在序列平稳的假设下 它反映了分布滞后变量 本例即降水量ra 对因变量的长期作用大小 表中系数即为 若认为降水量对水库流量的作用在3月之后几乎消失 则可利用远端限制条件 即输入命令Lsvolcpdl ra 9 4 2 比较发现 远端约束模型的调整后的决定系数略高于无约束模型 AIC和SC信息量略低于无约束模型 因此认为加入远端约束条件后的多项式分布滞后模型较优 但二者差异不大 系数估计值差异也不大 说明滞后期为3月时降水量对水库流量的作用本身已衰减接近0 根据需要 可以为模型增加ARMA项 比如对某商品销售额 sale 价格 price 和顾客流量 customer 建立分布滞后模型的同时 加入AR和MA项 在主窗口命令行输入 Lssalecpricepdl customer 10 3 ar 1 ma 1 这里 变量customer的系数取决于无端点约束的次数为3的多项式 PDL模型也可以用两阶段最小二乘法估计参数 命令基本格式为tslsyx1x2pdl seriesname lags orde options zlZ2 格兰杰因果检验 先估计当前的y值被其自身滞后期取值所能解释的程度 然后验证通过引入序列x的滞后值是否可以提高y的被解释程度 如果是 则称序列x是y的格兰杰成因 此时x的滞后期系数具有统计显著性 一般地 还应该考虑问题的另一方面 即y是否是x的格兰杰成因 Eviews计算如下的双变量回归 其中k是最大滞后阶数 通常可以取稍大一些 检验的原假设是序列x y 不是序列y x 的格兰杰成因 即 1999年1月4日至2001年10月15日深圳成指和上海综指序列数据见case9