MLE for a region switch model

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

The estimation is done with NLMIXED.

Following the link for additional information.
http://support.sas.com/rnd/app/da/new/802ce/ets/chap4/sect13.htm



data tmp;
   do i = 1 to 100;
       x1=rannor(123); x2=rannor(123);x3=rannor(123);
       s=1+1*x1>  rannor(123);
       if s=1 then y= 1+2*x2+ 1.5*rannor(123);
       else  y= 3+1*x3+ 1*rannor(123);
       output;
    end;
run;

   proc nlmixed data=tmp;
      parms a0 a1 b0 b1 c0 c1=0 s1 s2=10  ;
      bounds s1 s2 >0;

      xbeta=a0+a1*x1;
      prob=probnorm(xbeta);

      pdf1=( 1/( s1*sqrt(2*constant('pi')) ) *exp(-0.5*( (y-(b0+b1*x2))/s1) **2) );
      pdf2=( 1/( s2*sqrt(2*constant('pi')) ) *exp(-0.5*( (y-(c0+c1*x3))/s2) **2) );

      ll= min( max(prob *  pdf1 + (1-prob)*pdf2, 1e-20), 1-1e-20);

      loglik=log(ll);

      model y ~ general(loglik);
run;

扫码加我 拉你入群

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

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

关键词：switch Region Witch model mode additional

what's your point here?

这个model X~general() 里面的genera的作用时什么，是不是自定义的loglikelihood就要用这个