异方差eviews操作.docx

第四章 异方差性一、参数估计进入Eviews软件包，确定时间范围，编辑输入数据；选择估计方程菜单：(1)在Workfile对话框中，由路径：Quick/EstimateEquation，进入EquationSpecification对话框，键入“log(y)clog(x1)log(x2)”，确认ok，得到样本回归估计结果；(2)直接在命令栏里输入“lslog(y)clog(x1)log(x2)”，按Enter,得到样本回归估计结果；(3)在Group的当前窗口，由路径：Procs/MakeEquation，进入EquationSpecification窗口，键入“log(y)clog(x1)log(x2)”，确认ok，得到样本回归估计结果。如表4.1：表4.1二、检验模型的异方差性（一）图形法（1）生成残差平方序列1.在workfile 对话框中，由路径procs/generate series，进入GenerateSeriesbyEquation对话框，键入“e2=resid2”，生成残差平方项序列e2；直接在命令栏里输入“genre2=resid2”，按Enter,得到残差平方项序列e2。（2）绘制散点图直接在命令框里输入“scatlog(x2)e2”，按Enter，可得散点图4.2选择变量名log(x2)与e2（注意选择变量的顺序，先选的变量将在图形中表示横轴，后选的变量表示纵轴），再按路径view/graph/scatter/simplescatter，可得散点图4.2由路径quick/graph进入serieslist窗口，输入“log(x2)e2”，确认并ok，再在弹出的graph窗口把linegraph换成scatterdiagram，再点ok，可得散点图4.2。由图4.2可以看出，残差平方项e2对解释变量log(X2)的散点图主要分布图形中的下三角部分，大致看出残差平方项e2随log(X2)的变动呈增大的趋势，因此，模型很可能存在异方差。但是否确实存在异方差还应通过更进一步的检验。（二）Goldfeld-Quanadt检验（原假设：具有同方差，服从F分布）（1）对变量取值排序（按递增或递减）。在Workfile窗口中，由路径：Procs/SortSeries进入sortworkfileseries对话框，键入“X2”，如果以递增型排序，选Ascending，如果以递减型排序，则应选Descending，点ok。本例选递增型排序，选Ascending。将数据按组打开，鼠标右键，选择sort（2）构造子样本区间，建立回归模型在本例中，样本容量n=31，删除中间1/4的观测值，即大约7个观测值，余下部分平分得两个样本区间：1-12和20-31，它们的样本个数均是12个。在workfile 的Sample菜单里，把sample值改为“112”再用OLS方法进行第一个子样本回归估计，估计结果如表4.2。同样地，在Sample菜单里，把sample值改为“2031”再用OLS方法进行第二个子样本回归估计，估计结果如表4.3。（3）求F统计量值基于表4.2和表4.3中残差平方和RSS的数据，即Sumsquaredresid的值，得到RSS1=0.0702和RSS2=0.1912，根据Goldfeld-Quanadt检验，F统计量为：F=RSS2/RSS1=0.1912/0.0702=2.73。（4）判断在5%与10%的显著性水平下，查F分布表得：自由度为（9,9）（即n-k-1,其中，k为解释变量个数）的F分布的临界值分别为F0.05=3.18与F0.10=2.44。因为F=2.73F0.10(9,9)=2.44，因此10%显著性水平下拒绝两组子样方差相同的假设，即存在异方差。（三）White检验（原假设：具有同方差，服从卡方分布）在equation01中，按路径view/residualtests/heteroskedasticity test(crossterms)/white进入White检验，其中crossterms表示有交叉乘积项。得到表4.4的结果。由表4.4结果得到：怀特统计量nR2=310.6629=20.55，查2分布表得到在5%的显著性水平下，自由度为5（辅助回归式中待估参数个数）的2分布的临界值为20.05=11.07，因为nR2=20.5520.05=11.07，所以拒绝同方差的原假设。（辅助回归式为残差平方和对解释变量及其高次方项进行回归，有时有交叉项，当解释变量较多时可省去交叉项。）三、异方差性的修正（一）加权最小二乘法（WLS）（WLS的关键在于寻找模型中随机干扰项的方差与解释变量间的适当的函数形式）如果不知道异方差的形式，可以用FLGS估计量，也就是可行gls估计量。步骤是这样的：1.将y对x1 x2 x3 x4回归，得到残差e2.将e平方，取自然对数log(e2)3.将log(e2)对x1 x2 x3 x4做回归，得到拟合值gi4.权重就是1/exp(gi)(1/2)（1）生成权数。按路径：Procs/GenerateSeries，进入GenerateSeriesbyEquation对话框，键入“w=1/sqr(exp(93.20-25.981*log(x2)+1.701*(log(x2)2)”或者直接在命令栏输入“genrw=1/sqr(exp(93.20-25.981*log(x2)+1.701*(log(x2)2)”生成权数w。（2）加权最小二乘法估计（WLS）。在表4.1的结果中，由路径：Procs/Specify/Estimate进入EquationSpecification对话框，点击Options按钮，在EstimationOptions对话框的weighted前面打勾并在下面输入栏处输入w，如图4.3，连续两次确认OK后，得到表4.6的估计结果：加权最小二乘法估计（WLS）结果为：可以看出运用加权最小二乘法消除异方差性后，LnX1参数的t检验有了显著的改进，这表明即使在1%显著性水平下，都不能拒绝从事农业生产带来的纯收入对农户人均消费支出有着显著影响的假设。虽然LnX1的参数值有了较大程度的提高，但仍没有LnX2的参数估计值大，说明其他来源的纯收入确实比来自农业经营的纯收入对农户人均消费支出的影响更大一些。（3）检验加权回归模型的异方差性。在命令栏中直接输入“lsw*log(Y)ww*log(X1)w