## load data and transform with empirical cdf
data(wdbc)
udat <- apply(wdbc[, -1], 2, function(x) rank(x)/(length(x)+1))
## estimation of copula density of variables 5 and 6
dens.est <- kdecop(udat[, 5:6])
## compute fit statistics
logLik(dens.est)
AIC(dens.est)
BIC(dens.est)Run the code above in your browser using DataLab