我的似然函数是用mata编的一个极其复杂的函数，想问mata可以做极大似然估计吗？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请问一下，mata可以做极大似然估计吗？或者如何把我写的这个mata环境下的似然函数引用在stata里的极大似然估计中？？

似然函数如下：

function MLE(bbt1,bbt2,bbt3,bbt4,aalpha,pp1,pp2,aa,bb,cc,dd)

{

         n=abs(y)

         ut=y-x*(bbt1\bbt2\bbt3\bbt4)

         zzt_1=randaepd(aalpha,pp1,pp2,n)

         ww=aalpha*kep(pp1)+(1-aalpha)*kep(pp2)

         zzt_2=zzt_1/ww

         ee_aepd=E_aepd(aalpha,pp1,pp2)

         vvar_aepd=Var_aepd(aalpha,pp1,pp2)

         zzt=(zzt_2-ee_aepd)/sqrt(vvar_aepd)

         

         Ezt=funEzt(zzt,n)

         

         sigmat= funegarch_e( ut, Ezt, n, aa, bb, cc, dd)

         

         alphass=aalpha*kep(pp1)/(aalpha*kep(pp1)+(1-aalpha)*kep(pp2))

         e_aepd=E_aepd(aalpha,pp1,pp2)

         var_aepd=Var_aepd(aalpha,pp1,pp2)

         std_aepd=sqrt(var_aepd)

         A1=0

         A2=0

         for(i=1; i<=n; i++)

         {

                   if(zt(i,1)<=-e_aepd/std_aepd)

                   {

                            A1=A1-log(std_aepd/sigmat(i,1))-log(aalpha/alphass)-log(kep(pp1))+((abs((e_aepd+zt(i,1)*std_aepd)/(2*alphass)))^pp1)/pp1

                   }

                   else

                   {

                      A2=A2-log(std_aepd/sigmat(i,1))-log((1-aalpha)/(1-alphass))-log(kep(pp2))+((abs((e_aepd+zt(i,1)*std_aepd)/(2*(1-alphass))))^pp2)/pp2

                   }

         }

}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：极大似然估计 MATA 似然估计 似然函数 极大似然 function 如何