异方差的最大似然估计 哪儿出了问题？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

根据李子奈老师的《高级应用计量经济学》P39—41，用Mata的option()函数与了一个异方差条件下的最大似然估计，初始值设教材估计结果而设，但stata提示：初始值不可行。不知是哪儿出错了,特向大家请教！
出错提示、数据与我编写的程序如下：


data_do.rar (1.03 KB) 本附件包括：

• data.dta
• mymle_het.do

2012-10-30 16:24:42 上传

set more off
use data,clear
drop if ex==.
gen income2=income^2
gen cons=1
mata
      mata clear
      st_view(y=., ., " expense")
      st_view(X=., ., ("income","income2","cons"))
      void mymax(todo,b,y,X,lnf,g,H)
      {
              real matrix omega,e
  
              omega=b[cols(X)+1]*(diag(X[.,1]):^b[cols(X)+2])
              e=(y-X*b[1..cols(X)]')
              lnf = -0.5*rows(X)*(1+ln(2*pi())-ln(rows(X)))-0.5*ln(det(omega))-0.5*rows(X)*ln(e'*luinv(omega)*e)
      }
       s=optimize_init()
       optimize_init_evaluator(s,&mymax())
       optimize_init_evaluatortype(s,"v0")
       optimize_init_argument(s, 1, y)
       optimize_init_argument(s, 2, X)
       optimize_init_params(s, (-0.2,0,831,50,3)) //初始值设定
       b = optimize(s)      //参数估计结果
       v = optimize_result_V(s)   //系数估计结果

      b
     v
end

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：最大似然估计 最大似然 似然估计 异方差 Evaluator matrix income 李子

xuexi le

看样子又得自己寻找答案了！

已经解决！

楼主，请问初始值设定时时有什么技巧吗，我一直都是0或者1，但就是不知道有什么意义……

请问楼主是哪里出现了问题？我最近在做异方差的极大似然估计，非常迷茫，谢谢你的回答。