全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
4579 10
2015-08-15
代码如下: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请各位大牛回答,真的超感谢看了曹志广老师的书,案例不太一样,不知如何修改,还望各位牛人帮忙解答下,谢谢






二维码

扫码加我 拉你入群

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

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

全部回复
2015-8-16 10:21:21
有没有高手来帮帮忙啊,不胜感激
二维码

扫码加我 拉你入群

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

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

2015-8-16 10:26:09
我觉得这个像是曹志广的例子!
二维码

扫码加我 拉你入群

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

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

2015-8-16 10:30:22
magicsun 发表于 2015-8-16 10:26
我觉得这个像是曹志广的例子!
您好,可以问下是他哪本书或者哪篇文章中谈到的?非常感谢
二维码

扫码加我 拉你入群

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

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

2015-8-16 10:47:59
他在网上有个PDF文件,你搜一下,可能论坛内就有!
二维码

扫码加我 拉你入群

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

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

2015-8-16 10:54:49
magicsun 发表于 2015-8-16 10:47
他在网上有个PDF文件,你搜一下,可能论坛内就有!
能详细说下搜什么吗?我可以以什么关键字去搜索?我在论坛里以关键字曹志广没搜到,以极大似然估计搜的也没有他,在百度也没找到,可以麻烦您详细说下吗?非常感谢
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群