计量经济学Stata软件应用4---【Stata软件之异方差】--1次课概要课件.ppt

《计量经济学Stata软件应用4---【Stata软件之异方差】--1次课概要课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、计量经济软件应用——Stata软件实验之异方差实验目的：能够借助Stata软件诊断异方差的存在(White检验)和修正异方差(加权最小二乘法WLS)，能对软件运行结果进行解释。知识点：异方差检验的最常用方法——White检验出现异方差的最常见原因，是误差项的条件方差与某些解释变量相关，因此检验异方差的基本思路就是看误差项的条件方差是否与解释变量相关。怀特检验考虑到误差项方差与所有解释变量的相关关系，还进一步考虑了误差项方差与所有解释变量的平方及每两个解释变量的交互项的相关关系。怀特检验分为怀特一般检验(White’sgeneraltes

2、t)和怀特特殊检验(White’sspecialtest)。例如对于包含3个解释变量的原模型，怀特一般检验的模型为：原假设实际上是对模型进行回归总体显著性检验(F检验)，如果拒绝原假设，则存在异方差，如果不能拒绝原假设，则不存在异方差。容易看出，用于怀特一般检验的模型会包含很多解释变量；例如如果原模型有3个解释变量，那么怀特一般检验的模型将包含9个解释变量，而如果原模型有6个解释变量，那么怀特一般检验的模型将包含27个解释变量。这样对于样本容量不大的数据，使用怀特一般检验会使得估计的时候自由度偏少。为此，怀特特殊检验使用了一个节省自由度

3、的回归模型，即：其中是原模型的拟合值，是拟合值的平方。由于是所有解释变量的线性函数，而是这些解释变量的平方项和交互项的线性函数，因此用该模型代替怀特一般检验的模型来进行异方差检验是可行的。具体来说，针对模型(1)，同方差原假设为：原假设实际上是对模型(1)进行回归总体显著性检验(F检验);如果拒绝原假设，则存在异方差，如果不能拒绝原假设，则不存在异方差。异方差的修正——加权最小二乘法如果通过White检验发现存在异方差性，可以使用加权最小二乘(WLS)进行估计。1、异方差形式已知时的加权最小二乘估计(以一元为例)假如已经知道异方差的具体

4、形式，如：其中简记为是解释变量的一个已知函数；对于原模型两端乘以权重，得到：(同方差模型)2、异方差形式未知时的加权最小二乘估计(以一元为例)在一般情况下，我们不可能知道的异方差的具体形式，这就需要对异方差的函数形式做出估计，然后再进行加权最小二乘估计。这种方法属于可行的广义最小二乘估计(FGLS)的一种。处理异方差问题的FGLS的步骤是：第1步：对进行OLS回归，得到残差平方及其自然对数；第2步：对以下模型进行OLS回归，并得到拟合值第3步：计算的指数（这里exp是Stata指数函数的命令）第4步：以为异方差函数形式的估计对原模型进行

5、WLS估计，权重为，此时，变换后的模型为：此模型为同方差模型。异方差检验和修正的Stata基本命令whitetst对最近的回归进行怀特一般检验。whitetst,fitted对最近的回归进行怀特特殊检验。wls0yx1x2…xk,wvar(hh)type(abse)nocony对x1,x2,…,xk的WLS回归，wvar(hh)中的hh表示异方差函数形式的开平方，注意这里是异方差函数形式的开平方；选项type(abse)和nocon是使用上述FGLS方法时必须指定的。Stata软件操作实例实验1工资方程中异方差的检验和修正本例使用“工资

6、方程1.dta”数据文件介绍异方差的检验和修正。1、打开数据文件。直接双击“工资方程1.dta”文件；或者点击Stata窗口工具栏最左侧的Open键，然后选择“工资方程1.dta”即可；2、估计工资方程：其中lnwage工资对数，exp工作经验，expsq工作经验的平方；命令及运行结果：reglnwageeduexpexpsqStata软件操作实例3、异方差的检验：white检验怀特检验分为怀特一般检验和怀特特殊检验。怀特一般检验的步骤：第1步：对(2)式进行OLS回归。在这里，我们对这一回归的结果不感兴趣，可以在reg命令前加上qui

7、etly选项，其含义是让Stata进行回归，但不显示结果。(quietly可用于任何Stata命令的前面，表示不在Stata的Result窗口中显示分析结果。)quietlyreglnwageeduexpexpsq第2步：使用predict命令生成残差u(e)，并生成残差的平方usq()。predictu,residualgenusq=u^2Stata软件操作实例第3步：生成所有解释变量的平方项：edusq、expsqsq(原模型的解释变量中已经有exp的平方项expsq，所以不用再生成exp的平方项)；生成每两个解释变量的交互项：ed

8、u_exp、edu_expsq、exp_expsq。genedusq=edu^2genexpsqsq=expsq^2genedu_exp=edu*expgenedu_expsq=edu*expsqgenex