wqt2958 发表于 2012-11-2 22:04 
运行果然很慢……
请问您运行出来了吗?我做的时候,算clayton copula算出来了参数值,算gumble copula 和frank copula时候,程序运行出错了。也不知道哪里出问题了,您能帮我看看吗?
>> [p,l]= M_Copulas('Frank',zdata)
??? Error using ==> reshape
To RESHAPE the number of elements must not change.
Error in ==> sym.mtimes at 17
   X = reshape(maple('scalarmul',B(:),A),size(B));
Error in ==> M_Copulas at 72
        copula_pdf=-(1/a)*log(1+((prod((exp(-a*x)-1)'))')/(exp(-a)-1)^(n-1));