0.残差
rvfplot
画残差与拟合值散点图
rvpplot varname
画残差与解释变量散点图,varname 替换为变量名
scatter y x 与与x散点图
predict e,r计算残差
gen e2=e^2
scatter e2 x 计算e^2与x的散点图
1.怀特检验与BP检验
capture log close
日志关闭。
log using food,replace
日志命名为food,并替代原名为food的日志。后缀为smcl
set more off
日志更多记录关闭。
use food.dta,clear
使用 food.dta数据。
twoway scatter food_exp income || lfit food_exp income
绘制food_exp与income的散点图,前者为纵轴,后者为横轴||绘制food_exp与income的折线图,前者为纵轴,后者为横轴。twoway表示两个轴。
reg food_exp income
回归food_exp与income。
estat hettest, iid rhs
e类命令,进行辅助回归以进行BP检验。iid表示假设变量为iid(独立同分布),无需正态假定。rhs表示包含模型中等号右侧所有解释变量。P值大于临界值无显著异方差。hettet表示异方差检验。
estat imtest, white
进行怀特检验。P值大于临界值无显著异方差。
gen food_share = food_exp/income
定义food_share 为food_exp/income。
twoway scatter food_share income || lfit food_share income
对food_share与income做散点图和线性拟合图。
reg food_share income
对food_share与income进行回归。
estat hettest, iid rhs
e类命令,进行辅助回归以进行BP检验。iid表示假设变量为iid(独立同分布),无需正态假定。rhs表示包含模型中等号右侧所有解释变量。P值大于临界值无显著异方差。
estat imtest, white
进行怀特检验。P值大于临界值无显著异方差。
log close
关闭日志
estat hettest [varlist], iid rhs
对指定变量进行BP检验
异方差性检验中P值应尽可能大,小存在异方差
2.WLS(加权最小二乘)
reg y x1 x2 x3 [aw=1/var]
aw为扰动项方差的倒数[aw=1/var] 为权重。
reg y x [aw=1/x]
aw为扰动项方差的倒数[aw=1/x] 为权重。
reg y x [aw=1/x^2]
aw为扰动项方差的倒数[aw=1/x^2] 为权重。
reg y x [aw=1/sqrt(x)]
aw为扰动项方差的倒数[aw=1/sqrt(x)] 为权重。
predict e1,residual
计算残差,记为e1
reg lne2 logq,noc
对lne2与logq进行无常数项的回归
3.OLS+稳健标准误
reg y x,r
普通最小二乘后加”,r”即可得OLS+稳健标准误