[求助]Maxlik怎麼改成Optmum

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我是Gauss的新手，

最近在研究怎麼把Maxlik改成Optmum

library maxlik;          (改成library optmum;)
#include maxlik.ext (改成#include optmum.ext)
maxset;                 (改成optset;)

z0 = {
     8.3 1,
    10.3 2,
    19.0 3,
    16.0 4,
    15.6 5,
    19.8 7 };

proc lnlk(b,z);
    local dev,s2;
    dev = z[.,1] - b[1]*(1 - exp(-b[2]*z[.,2]));
    s2 = dev'dev/rows(dev);
    retp(lnpdfmvn(dev,s2));
endp;

proc grd(b,z);
    local dev,s2,m,r;
    r = exp(-b[2]*z[.,2]);
    m = 1 - r;
    dev = z[.,1] - b[1]*m;
    s2 = dev'dev/rows(dev);
    dev = dev / s2;
    retp((dev.*m)~(b[1]*dev.*z0[.,2].* r));
endp;

start = { 20, .24 };

output file = maxlik10.out reset;

_max_GradProc = &grd;
_max_Options = { newton };

call maxprt(maxlik(z0,0,&lnlk,start));   (改成call optprt(optmum(&lnlk,start)); 

output off;

跑出來卻出現

c:\gauss6.0\src\optutil.src(136):error G0159: Wrong number of parameters

麻煩Gauss高手幫我看看為什麼參數會出錯

非常感謝大家!

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MAXLIK OPTMUM max OPT PTM 求助 MAXLIK OPTMUM