gumbel.cop <- gumbelCopula(dim=2)
u<-rank(x)/(length(x)+1)
v<-rank(y)/(length(y)+1)
f<-cbind(u,v)
uv<-pobs(f)
fit<- fitCopula(gumbel.cop, uv, method="mpl")
gumbel.cop <- gumbelCopula(2.28458,dim = 2)
2.28458是估计的参数。这样写代码最后得到的copula函数是正确的要估计的函数吗？
然后怎么求u和v的偏导数啊？也就是他们的条件概率？
用deriv求偏导的时候变量名字是什么啊？