How to fit a copula model in R [heavily revised]. Part 2: fitting the copula. Full article at

copulas_revised_2_4.R

# Estimate copula parameters
cop_model <- claytonCopula(dim = 2)
m <- pobs(as.matrix(mydata))
fit <- fitCopula(cop_model, m, method = 'ml')
coef(fit)
##    param 
## 1.482168 

# Check Kendall's tau value for the Clayton copula with theta = 1.48
tau(claytonCopula(param = 1.48))
## [1] 0.4252874