[未答复] 用最大似然法求解三维copula函数的参数

[复制链接]
canbaobao 发表于 6 天前
想用matlab最大似然法mle函数求三维copula函数的参数
X=[43370000000 42910000000 44850000000 44160000000 45150000000 40190000000 42380000000
43890000000 45310000000 44190000000 45180000000 45600000000 46570000000;3650800000 3703800000
3922200000 3991700000 4123700000 4172500000 4463900000 4495200000 5245600000
4785600000 5214700000 5733500000 6924300000;180.61 183.6 188.75 190.05 192.5 195.52
215.5 215.55 214.9 218.28 239.87 261.96 264.2];
>> a1=4.4839*10^10;
>> b1=35.42394798;
>> a2=22.2429;
>> b2=0.1868;
>> a3=5.3507;
>> b3=0.1282;
>> pdffun=@(X,thd)(2*thd+1)*(thd+1)*((a1*b1*(X(1,:)*a1)^(b1-1)*exp((-X(1,:)*a1)^b1))*(1/(X(2,:)*b2*(2*pi)^1/2)*exp
(-(log(X(2,:))-a2)^2/(2*b2^2)))*(1/(X(3,:)*b3*(2*pi)^1/2)*exp(-(log(X(3,:))-a3)^2/(2*b3^2))))^(-thd-1)*(((a1*b1*(X
(1,:)*a1)^(b1-1)*exp((-X(1,:)*a1)^b1))^(-thd)+(1/(X(2,:)*b2*(2*pi)^1/2)*exp(-(log(X(2,:))-a2)^2/(2*b2^2)))^(-
thd)+(1/(X(3,:)*b3*(2*pi)^1/2)*exp(-(log(X(3,:))-a3)^2/(2*b3^2)))^(-thd)-2)^(-1/thd-3));
>> [phat,pci]=mle(X,'pdf',pdffun,'start',0)
错误使用 mle (line 226)
DATA must be a vector.
代码如上,一直显示数据必须为向量 想询问一下怎么解决
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
热门教程
站长推荐
快速回复 返回顶部 返回列表