全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
10986 4
2015-04-16
多元线性 回归之后 检测到 异方差 , 一个参数 和 残差 存在相关性 。
如何通过WLS驱除 影响?

reg DISTANCE AGE KIDS INCOME

1.png
Holiday.dta
大小:(8.67 KB)

 马上下载


通过生成图片 发现 残差和 income 可能存在 异方差。(下图 左上的图)

twoway (scatter res INCOME) ,saving (1)

(file 1.gph saved)

. twoway (scatter res AGE),saving (2)

(file 2.gph saved)

. twoway  (scatter resKIDS),saving (3)

(file 3.gph saved)

graph combine 1.gph 2.gph 3.gph


2.png
进一步 测定:

Golddeld-Quandt Test

sort INCOME

g index = _n

reg DISTANCE AGE KIDS INCOME if index  <81

scalar RSS1 = e(rss)

scalar list RSS1

RSS1 =  7434484.1


reg DISTANCE AGE KIDS INCOME if index  >112

scalar RSS2 = e(rss)

scalar list RSS2

RSS2 =   24239434

scalar GQ = RSS2/RSS1

scalar DOF = 81-3-1

scalar PV = Ftail(DOF,DOF,GQ)

scalar list GQ PV

        GQ =  3.2604056


        PV =  2.460e-07

对比确定存在异方差 GQ 大于 PV

为了修正影响 已经 使用了. reg DISTANCE AGE KIDS INCOME,robust

但我想 知道 如果我想要用 WLS 进行修正 应该如何做

比如:

用ressidual 的标准差 作为 权重参数 调整

用ressidual 的方差 作为 权重参数 调整


我最糊涂的是 我只知道当 一元回归的时候可以用

g yi = y / sqrt (res)

g ic = 1/sqrt(res)

g ix = x/sqrt(res)

reg iy ic ix, noco

来做

那么 如果是 多元 那么是否应该是:

g yi = y / sqrt (res)

g ic = 1/sqrt(res)

g ix1 = x1/sqrt(res)

g ix2 = x2/sqrt(res)

g ix3 = x3/sqrt(res)

reg iy ic ix1 ix2 ix3, noco

这样对吗?


二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2015-4-16 13:02:17
lililll 发表于 2015-4-16 00:49
多元线性 回归之后 检测到 异方差 , 一个参数 和 残差 存在相关性 。
如何通过WLS驱除 影响?
reg DIST ...
直接用命令reg y x1 x2 xk[aw=1/eles] 就可以了。具体的步骤是:1.先做OLS,导出残差e;2.生成残差的平方得到es;3.将残差的平方取对数得到les;4.将les对x1 x2 xk做OLS回归,得到les的拟合值;5.对拟合值取指数得到eles。最后输入上面的命令就可以了,不需要对每个解释变量进行处理。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2015-4-17 00:59:02
hs4601 发表于 2015-4-16 13:02
直接用命令reg y x1 x2 xk[aw=1/eles] 就可以了。具体的步骤是:1.先做OLS,导出残差e;2.生成残差的平方 ...
为什么我要先做log? 直接用e^2 和自变量 作回归 会有什么影响呢?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2015-4-17 08:06:17
lililll 发表于 2015-4-17 00:59
为什么我要先做log? 直接用e^2 和自变量 作回归 会有什么影响呢?
直接用e2对自变量做回归不能保证e2的拟合值为正值。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-6-10 17:44:06
hs4601 发表于 2015-4-16 13:02
直接用命令reg y x1 x2 xk[aw=1/eles] 就可以了。具体的步骤是:1.先做OLS,导出残差e;2.生成残差的平方 ...
请问拟合值是不是系数?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群