大家帮忙看一下这个ks检验的程序哪儿错了，谢谢啦

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

zdata=[a,b,c,d];
a=[0.3;3;4;.3;.4;.3];
x=zdata;
k=length(zdata(:,1));
n=length(zdata(1,:));
copula_cdf=a(4)*(sum((x.^(-a(1)))')'-n+1).^(-1/a(1))+a(5)*exp(-sum(((-log(x)).^a(2))')'.^(1/a(2)))-a(6)/a(3)*log(1+prod((exp(-a(3)*x)-1)')'/(exp(-a(3))-1)^(n-1));
for i=1:k
    empirical_cdf(i)=mean((x(:,1)<=x(i,1)).*(x(:,2)<=x(i,2)).*(x(:,3)<=x(i,3)).*(x(:,4)<=x(i,4)));
end
empirical_cdf=empirical_cdf';
[h,p,ksstat,cv]=kstest(empirical_cdf,copula_cdf)
程序是这样的，最后结果说列向量必须是个矩阵，那位大神帮忙看看怎么改啊

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：ks检验 Empirical length Copula opula 程序