怎么用wls0命令实现陈强老师书上的异方差命令

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

陈强老师的《高级计量经济学及stata应用》中介绍异方差的WLS操作部分（第二版的96页）提到了一个外部命令wls0，我下载了但是不知道这个命令的调用格式。

试了
wls0 lntc lnq lnpl lnpk lnpf,wvar (lnq) type(abse) noc


wls0 lntc lnq lnpl lnpk lnpf,wvar (lnq) type(loge2) noc


这两条命令都得不到和书上给的操作一样的回归结果。请问谁知道要得到书上96页一样的结果，wls0的命令该怎么写？

扫码加我 拉你入群

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

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

分享0 收藏4 回帖

关键词：Wls 异方差 高级计量经济学及STATA应用 高级计量经济学 stata应用 经济学

可以用type(e2)，得到的结果更好。具体的解释我也在找。

我用这个命令，type（abse）得到的结果跟书上很相近。

这个的确做得跟书上答案不一样

wls0中的命令，到底应该怎么选择呢？type（abse）还是loge2？还是xb2呢？

在plus文档中打开wls0.ado文件，看看do文件的命令就懂了。跟书上写的命令不一样

  quietly regress `varlist' `if' `in'
  quietly predict `p3'
  quietly predict `e', resid
  
  if "`type'"=="abse" {
    generate `ee' = abs(`e')
    local eetype " type: proportional to abs(e)"
  }
  if "`type'"=="e2" {
    generate `ee' = (`e')^2
    local eetype " type: proportional to e^2"
  }
  if "`type'"=="loge2" {
    generate `ee' = log((`e')^2 )
    local eetype " type: proportional to log(e^2)"
  }
  if "`type'"=="xb2" {
    local eetype " type: proportional to xb^2 "
  }
  
  if "`type'"~="xb2" {
    quietly regress `ee' `wvars', `noconst'
    quietly predict `p1'
  }
  else {
    gen `p1' = `p3'
  }

  generate _wls_wgt = 1/`p1'^2

  label variable _wls_wgt "wls weights"

  display
  display in green "WLS regression - `eetype'"
  display
  regress `varlist' `if' `in' [aw = _wls_wgt], `robust' `hc2' `hc3'

应该是：
gen aa=sqrt(e2f)
wls0 intc lnq lnpl lnpk lnpf,wvar (aa) type(abse) noc

我用的wls0 lntc lnq lnpl lnpk lnpf,wvar (lnq) type(loge2) noc
得到书上的结果了哎

wvar是什么意思呢