紧急！求助！用matlab建立三维copula模型，正用极大似然估计方法求copula函数参数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

紧急！求助！我正在利用matlab建立三维copula函数模型，正利用极大似然估计方法估计copula函数中的参数，但是卡在这里了。这个似然函数非常的长
下面是我的程序：function siranhanshu=myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf)h=(theta^2*(1/theta + 1)*(1/theta + 2))/(U(:).^(theta + 1)*V(:).^(theta +1)*W(:).^(theta + 1)*(1/U(:).^theta + 1/V(:).^theta + 1/W(:).^theta - 2)^(1/theta +3))*Xpdf*Ypdf*Zpdf;
k=prod(h);
siranhanshu=0-log(k);

然后我再命令窗口键入
%在命令窗口输入
[theta,fval]=fminsearch(@(theta) myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf),5.000000001)


其中U,V,W,Xpdf,Ypdf,Zpdf是程序中也求出来的一列多行的矩阵


运行以后报错
Error using  ^
Inputs must be a scalar and a square matrix.
To compute elementwise POWER, use POWER (.^) instead.

Error in myfun (line 3)
h=(theta^2*(1/theta + 1)*(1/theta + 2))/(U^(theta + 1)*V^(theta
+1)*W^(theta + 1)*(1/U^theta + 1/V^theta + 1/W^theta - 2)^(1/theta
+3))*Xpdf*Ypdf*Zpdf;

Error in @(theta)myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf)


Error in fminsearch (line 191)
fv(:,1) = funfcn(x,varargin{:});

Error in P_ET0_IR (line 131)
[theta,fval]=fminsearch(@(theta)
myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf),5.000000001)

而且警告Warning: The value of local variables may have been changed to match the
         globals.  Future versions of MATLAB will require that you declare
         a variable to be global before you use that variable.
> In myfun at 2
  In @(theta)myfun(theta,U,V,W,Xpdf,Ypdf,Zpdf)
  In fminsearch at 191
  In P_ET0_IR at 131



新手上路，能否给予指点！谢谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Copula 极大似然估计 MATLAB opula atlab function matlab Error 模型 程序

看到您还没有邮箱验证 去验证吧
关于新注册用户进行邮箱验证的通告
https://bbs.pinggu.org/thread-1018186-1-1.html

到论坛第一区 您发快速提问  看看  讨论

你能求出三元的就已经很不容易了！

楼主的问题解决了吗？

关于copula的问题可以加我qq535844430，欢迎交流。