prototype of SAS/IML for a probit model

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

You may use it as a prototype for another ML estimation in SAS/IML.
********************************************************************;
  data probit;
      do i = 1 to 50;
         x=rannor(123);
         y=1+1*x> rannor(123);
         output;
        end;
    run;

    proc logistic data=probit;
    model y(event='1')=x/link=probit;
    run;

    proc iml;

       use probit;
       read all var{y x};
       nobs= nrow(y);
       yint=j(nobs,1) -y;

       *define a likelihood fun;

   start loglik(pr) global(y, yint, x, nobs);
        c1 = pr[1]; c2 = pr[2];
        xbeta=c1+c2#x;
        p=probnorm(xbeta);
        pint= j(nobs,1) - p;
        logp=log(p#y + pint#yint);
        ll=-2*logp[+];
        return (ll);
  finish loglik;

    pr = j(1,2,1e-5 );
    optn = {0 2 };
    call nlpdd(rc, xr,"loglik", pr, optn);
    print "Start" pr, "RC=" rc, "Opt Par" xr;

quit;

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：prototype Probit model type mode SAS Probit model iml prototype

不错。学习了。