matlab的最大似然估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我现在已知一组数据和它的概率密度函数(函数中有参数,且这个函数matlab中没有),如何用最大似然估计,估计出参数呢?
请各位大虾帮帮忙,急用啊!谢谢了!!!!!

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：MATLAB 最大似然估计 matla atlab 似然估计 MATLAB

[cry][dizzy]

我也遇到相同问题。如果你的参数不超过7个  也许我可以帮你

我的是参数有12个  不能用fminsearch()函数

我用fminunc老出错  忘达人帮忙啊。

附我出错的程序如下：看起来很复杂，其实密度函数本来是混合正态分布，只是对参数又进行了假设所以变复杂了。

%----likelihood function  M函数文件
function fun3=myfun(x,s,f,b)
r01=x(1);r11=x(2);r02=x(3);r12=x(4);f01=x(5);
nd01=x(6);nd11=x(7);nd02=x(8);nd12=x(9);f11=x(10);f02=x(11);f12=x(12);

pi1=(1+exp(f01+f11*b))/(2+exp(f01+f11*b)+exp(f02+f12*b));

mu1=(s-r01-r11*f).^2;mu2=(s-r02-r12*f).^2;

sd1=sqrt(2*pi*exp(nd01+nd11*b));
sd2=sqrt(2*pi*exp(nd02+nd12*b));

fun3=-sum(log(pi1*exp(-mu1/exp(nd01+nd11*b)/2)/sd1+(1-pi1)*exp(-mu2/exp(nd02+nd12*b)/2)/sd2));

%窗口主程序
clc
mydata1=xlsread('spot.xls');
mydata2=xlsread('future.xls');
mydata3=xlsread('basic.xls');
s=mydata1(:,1);
f=mydata2(:,1);
b=mydata3(:,1);
options=optimset('fminunc');
[x,fval,exitflag,output]=fminunc(@fun3,[0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5],options,s,f,b)

忘了把出错信息写上了。

??? Error using ==> mrdivide
Matrix dimensions must agree.

Error in ==> fun3 at 13
fun3=-sum(log(pi1*exp(-mu1/exp(nd01+nd11*b)/2)/sd1+(1-pi1)*exp(-mu2/exp(nd02+nd12*b)/2)/sd2));

Error in ==> fminunc at 182
   f = feval(funfcn{3},x,varargin{:});

我已经解决了,是我自己的初值有问题,还有程序有占小问题!

lz,能分享下如何解决的么，我也碰到类似问题，可否粘出相关程序供参考？谢谢

pzpenny 发表于 2010-4-6 16:24
我已经解决了,是我自己的初值有问题,还有程序有占小问题!

你好，我最经正在为matlab基于fminsearch的最大似然估计的初值问题困扰，请问这个初值设置有什么要求吗？为什么我设置不同的初值得到的结果都是不一样的呢？

seapig83 发表于 2010-4-6 15:09
我也遇到相同问题。如果你的参数不超过7个  也许我可以帮你

我的是参数有12个  不能用fminsearch()函数
...

为什么说参数不超过7个才行呢?