求高手回答matlab求极大似然估计的问题 - MATLAB等数学软件专版

0关注
0粉丝

硕士生

90%

还不是VIP/贵宾

-

0%

威望: 0 级
论坛币: 7 个
通用积分: 3.0000
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 1291 点
帖子: 76
精华: 0
在线时间: 211 小时
注册时间: 2005-8-29
最后登录: 2022-9-27

楼主

我是学士

发表于 2015-8-15 22:41:13 |只看作者 |坛友微信交流群|倒序 |AI写论文

相似文件

换一批

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

代码如下：function [para,fv]=my_mle1(fun,para0,varargin)
para0=para0(:);
[para,fv]=fminsearch(fun,para0,[],2,varargin{:});
fv=-fv;
end

function f=pin(para,num,B,S)
alpha=para(1);
mu=para(2);
eb=para(3);
es=para(4);
delta=para(5);
M=(min(B,S)+max(B,S))/2;
xs=es/(es+mu);
xb=eb/(eb+mu);
yy=(-eb-es + M*(log(xb) +log(xs)) + ...
B*log(mu+eb)+S*log(mu+es) )+ ...
log( alpha*(1-delta)*exp(-mu)*(xs.^(S-M)) .*(xb.^sym(-M)) + ...
alpha* delta*exp(-mu) *(xb.^(B-M)).*(xs.^(-M)) + ...
(1-alpha)*(xs.^(S-M)).*(xb.^(B-M)));
if num==1
f=log(yy);
else
f=-sum(log(yy));
end
end

输入指令：B=x1PIN(:,3)
S=x1PIN(:,4)
[para,fv]=my_mle1('pin',[0.5;400;1000;1000;0.5],B,S)

结果出错，详细如下：
??? Conversion to double from sym is not possible.

Error in ==> fminsearch at 175
fv(:,1) = funfcn(x,varargin{:});

Error in ==> my_mle1 at 3
[para,fv]=fminsearch(fun,para0,[],2,varargin{:});

有没有会的大牛帮我解答下，万分感谢。详细的代码和数据我都放在求助区这个帖子里https://bbs.pinggu.org/thread-3855084-1-1.html请各位大牛回答，真的超感谢看了曹志广老师的书，案例不太一样，不知如何修改，还望各位牛人帮忙解答下，谢谢