iamcooler 发表于 2012-6-4 13:10

求助,如何用fminsearch求最大似然估计

如题。L(x;b)是似然函数,求使得L最大化的b,那么用fminsearch时可以理解是求使-L最小的b。
b0=;%b的范围
bt11=fminsearch(@mlfunction,b0,[],x,r,m);%[]求的就是最大似然估计么?
似然函数就是g,但是不知道怎么表示才对,请教各位大侠!
function g=mlfunction(b0,b,x,r,m)
y=x.^b0;
h2=0;h3=0;h4=0;
for i=1:m
    h2=h2+(1+r(i))*y(i)*log(x(i));
    h3=h3+(1+r(i))*y(i);
    h4=h4+log(x(i));
end
g=-m/b0+h4+m*h2/h3;

thfdj007 发表于 2012-6-7 10:49

路过!!!
帮顶……
页: [1]
查看完整版本: 求助,如何用fminsearch求最大似然估计