nested archimedean copula 的参数估计

heqingsheng1987

代码如下：

theta0<-copGumbel@tauInv(.2)

theta1<-copGumbel@tauInv(.5)

c3<-onacopula("G",C(theta0,1,C(theta1,c(2,3))))

u<-rnacopula(90,c3)

fit.c3<-enacopula(u,c3,method="mle")

错误如下：

错误: max(cop@comp) == d is not TRUE

求各位大侠帮忙改正啊，完全嵌套和部分嵌套的参数怎么估计和拟合