stata基础–回归,画散点图,异质性分析

利用stata的内部数据来进行回归

代码:

sysuse auto

sysuse dir       /*可以看到所有的数据*/

su price mpg foreign

reg price mpg

predict u,residual       /*新变量u=每一个观测的残差*/

                                   /*生成残差u需要紧接着回归*/

stata基础--回归,画散点图,异质性分析

stata基础--回归,画散点图,异质性分析

mpg和price在0.01显著性水平上负相关。

利用回归得出的结果表计算系数,写出方程

经验回归方程:

stata基础--回归,画散点图,异质性分析

Y:price=6165.257

X:mpg=21.2973

B1=-238.8943

所以B0=11253.06057539

Y=11253.06057539-238.8943x

画出x-y图

graph twoway lfit price mpg || scatter price mpg

stata基础--回归,画散点图,异质性分析

画出残差和y的散点图

graph twoway lfit u mpg || scatter u mpg

stata基础--回归,画散点图,异质性分析

【如果误差满足零均值、同方差且互不相关,则回归系数的最佳线性无偏估计。】

由图可知,误差满足零均值,但为异方差。

异质性分析

异质性分析代码:

su price mpg if foreign ==0

reg price mpg if foreign ==0

回归

stata基础--回归,画散点图,异质性分析

stata基础--回归,画散点图,异质性分析

price对mpg的OLS回归,只使用国产车为样本:

mpg和price在0.01显著性水平上负相关。

经验回归方程:

Y:price=6072.423

X:mpg=19.82692

B1=-329.2551

所以B0=12600.537527292

Y=12600.537527292-329.2551x

画图:

graph twoway lfit price mpg || scatter price mpg

graph twoway lfit u mpg || scatter u mpg

stata基础--回归,画散点图,异质性分析

stata基础--回归,画散点图,异质性分析

由图可知,国产车为样本时,误差均值不为0,且为异方差。

异质性分析报告:

一、经验回归方程对比

1、研究对象为总体时:

Y:price=6165.257

X:mpg=21.2973

B1=-238.8943

所以B0=11253.06057539

Y=11253.06057539-238.8943x

2、研究对象为国产车时:

Y:price=6072.423

X:mpg=19.82692

B1=-329.2551

所以B0=12600.537527292

Y=12600.537527292-329.2551x

二、图形对比

1、研究对象为总体时                                                   2、研究对象为国产车时

stata基础--回归,画散点图,异质性分析                         stata基础--回归,画散点图,异质性分析

三、结论

总体和国产车并无异质性(没有意义,论文中不用出现,论文中要出现的是异质性明显的)

总代码:

sysuse auto
sysuse dir
su price mpg foreign

reg price mpg        /*回归*/
predict u,residual        /*残差*/

graph twoway lfit price mpg || scatter price mpg      

 /*画出price(y)和mpg(x)的散点图和样本回归线*/
graph twoway lfit u mpg || scatter u mpg          

 /*画出残差和mpg(x)的散点图和样本回归线*/

/*异质性分析*/
su price mpg if foreign ==0
reg price mpg if foreign ==0
predict u,residual
graph twoway lfit price mpg || scatter price mpg
graph twoway lfit u mpg || scatter u mpg

文章出处登录后可见!

已经登录?立即刷新

共计人评分,平均

到目前为止还没有投票!成为第一位评论此文章。

(0)
青葱年少的头像青葱年少普通用户
上一篇 2023年2月25日 下午8:45
下一篇 2023年2月25日 下午8:46

相关推荐