楼主: pengfeizou
6460 2

[Stata高级班] 请教连老师有关异方差检验的不一致性问题 [推广有奖]

  • 0关注
  • 0粉丝

VIP

硕士生

75%

还不是VIP/贵宾

-

威望
0
论坛币
3459 个
通用积分
29.5469
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
2583 点
帖子
122
精华
0
在线时间
206 小时
注册时间
2006-2-11
最后登录
2023-10-23

相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
我使用李子奈的3.2节中的数据
,比较异方差检验的BP检验和WHITE检验,结果发现BP检验如果使用拟合值(见命令行 estat hettest  //异方差的BP检验),则结果与white检验相反,但使用原始值时结果与WHITE检验相同,在此先谢过?附件是数据。以下是我做的过程。
**例题3.2.2城镇居民2006年的消费支出Y对人均可支配收入X1和
***  其前一年的人均消费X2的回归。X2表示消费惯性。
use "D:\StataSE12\zpf\lizinai\yifangcha.dta", clear
ren var1 y
ren var2 x1
ren var3 x2
regress y x1  x2        //做y对X1和X2的回归
estat hettest  //异方差的BP检验
estat imtest, white  //WHITE TEST
predict e,resid
gen e2=e*e
reg e2 x1  x2
dis "BP=" e(N)*e(r2)  //BP TEST
二维码

扫码加我 拉你入群

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

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

关键词:异方差检验 方差检验 一致性 连老师 异方差 检验 大虾 下载次数 white 李子

附件: 你需要登录才可以下载或查看附件。没有帐号?我要注册
沙发
arlionn 在职认证  发表于 2011-11-22 20:02:38 |只看作者 |坛友微信交流群
只是视频中的命令,你写的好像有些问题:
* B-P test (Breusch-Pagan,1979)
*  sigma_i^2 = sigma^2*(a0 + a1*z1 + a2*z2 + ...)
*  H0: a1=a2=...=0
*  检验思路:用 ei^2/avg(ei^2) 对一系列可能导致异方差的变量作回归
*    LM= 0.5*MSS -- chi2(k)  (k 为解释变量的个数,不包括常数项)
   sysuse auto, clear
   reg price weight mpg turn foreign
   estat hettest, normal   /*B-P检验的原意,同方差假设*/
      * 手动计算
      qui predict yhat
      qui predict e, res
      qui gen e2 = e^2
      qui sum e2
      qui gen e2_nor = e2/r(mean)
      qui reg e2_nor yhat
      dis 0.5*e(mss)
   reg price weight mpg turn foreign
   estat hettest, iid      /*N*R2, 无需同方差假设*/
      qui reg e2 yhat
      dis e(N)*e(r2)
   reg price weight mpg turn foreign
   estat hettest, rhs
   estat hettest length trunk
   estat hettest length trunk, rhs

使用道具

藤椅
pengfeizou 发表于 2011-11-28 22:25:22 |只看作者 |坛友微信交流群
谢谢你辛勤的回复。我想我没有表达清楚我的问题。我在论坛https://bbs.pinggu.org/thread-1256870-1-1.html也发了这个问题。我不知道是因为我的软件是破解版还是其它原因,按理检验结果应该相同。
我使用李子奈的3.2节中的数据  liti030202.rar (584 Bytes, 下载次数: 5)
2011-11-18 09:34:09 上传下载次数: 5
,比较异方差检验的BP检验和WHITE检验,结果发现BP检验如果使用拟合值(回归后用命令estat hettest  检验),则结果与white检验(回归后用命令estat imtest, white  检验,见下方红色部分)相反,但使用原始值时结果与WHITE检验(见最后几行)相同,哪位大虾解释下,在此先谢过?附件是数据。以下是我做的过程。
**例题3.2.2城镇居民2006年的消费支出Y对人均可支配收入X1和
***  其前一年的人均消费X2的回归。X2表示消费惯性。

regress y x1  x2        //做y对X1和X2的回归
estat hettest  //异方差的BP检验
. estat hettest, normal   /*B-P检验的原意,同方差假设*/
Breusch-Pagan / Cook-Weisberg test for heteroskedasticity
         Ho: Constant variance
         Variables: fitted values of y
         chi2(1)      =     0.23
         Prob > chi2  =   0.6348


estat imtest, white  //WHITE TEST
White's test for Ho: homoskedasticity
         against Ha: unrestricted heteroskedasticity
         chi2(5)      =     25.31
         Prob > chi2  =    0.0001
Cameron & Trivedi's decomposition of IM-test
---------------------------------------------------
              Source |       chi2     df      p
---------------------+-----------------------------
  Heteroskedasticity |      25.31      5    0.0001
            Skewness |      21.37      2    0.0000
            Kurtosis |       0.47      1    0.4928
---------------------+-----------------------------
               Total |      47.15      8    0.0000
---------------------------------------------------
上述结论表明两者检验结果相反。
但使用伍德里奇介绍的原理来做bp检验时结果与white检验相同。
predict e,resid
gen e2=e*e
reg e2 x1  x2
dis "BP=" e(N)*e(r2)  //BP TEST

本文来自: 人大经济论坛 Stata专版 版,详细出处参考: https://bbs.pinggu.org/forum.php? ... ge=1&from^^uid=119229

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-6-17 17:02