关于时变copula估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我有个程序，但是总是提醒这个错误。请大家帮我看看到底是什么问题
??? Error using ==> mtimes
Inner matrix dimensions must agree.


Error in ==> clayton_tvp_CL at 19
    CL = log(data(:,1).*data(:,2))*(1+clay);

Error in ==> fmincon at 517
      initVals.f = funfcn{3}(X,varargin{:});

Caused by:
    Failure in initial user-supplied objective function evaluation. FMINCON cannot continue.

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：时变Copula Copula opula Dimensions Evaluation objective function continue initial cannot

1. function [CL, clay] = clayton_tvp_CL(theta,data,clabar)
   
2. 
   
3. u = data(:,1);
   
4. v = data(:,2);
   
5. T  = size(data,1);
   
6. 
   
7. clay = -999.99*ones(T,1);
   
8. clay(1) = clabar;                       
   
9. for jj = 2:T;
   
10.    if jj<=10
    
11.       psi1 = theta(1)+ theta(2)*clay(jj-1) + theta(3)*(mean(abs(u(1:jj-1)-v(1:jj-1))));
    
12.    else
    
13.       psi1 = theta(1)+ theta(2)*clay(jj-1) + theta(3)*(mean(abs(u(jj-10:jj-1)-v(jj-10:jj-1))));
    
14.    end
    
15.    clay(jj) = 0.0001 + exp(psi1);
    
16. end
    
17. if clay~=0;
    
18.     CL = log(data(:,1).*data(:,2))*(1+clay);
    
19.     CL = CL + (2+1/clay).*log( (data(:,1).^(-clay)) + (data(:,2).^(-clay)) -1);
    
20.     CL = log(1+clay) - CL;
    
21.     CL = sum(CL);
    
22.     CL = -CL;
    
23. else
    
24.     CL = 0;  % under independence the copula pdf=1, so the log-likelihood=0
    
25. end

复制代码

1. lower = -5*ones(3,1);  % in theory there are no constraints, but setting loose constraints sometimes helps in the numerical optimisation
   
2. upper =  5*ones(3,1);
   
3. theta0 = [log(kappa2);0;0];
   
4. [ kappa13 LL13] = fmincon('clayton_tvp_CL',theta0,[],[],[],[],lower,upper,[],options,[u,v],kappa2);

复制代码

运行到第4行总是出现那个错误，到底是哪错了，我是根据其他程序模仿写的

您好，调出来了吗？我也做这个能帮帮我吗？
不会时变copula编程

您好，我也做时变copula，但目前没有相关程序代码，能帮帮我吗？
或是相关文章含有相关代码的也可以，谢谢了