如题,u是前面步骤产生的矩阵2列,不知道是什么算法问题吗不太懂
L<-function(pa){
a<-pa
norm.cop<-normalCopula(a,dim=2)
l<-1/nT*sum(log(dCopula(u,norm.cop)))
return(l)
}
est2<-maxLik(L,start=0.5)
报错:Error in maxNRCompute(fn = function (theta, fnOrig, gradOrig = NULL, hessOrig = NULL, :
NA in the initial gradient