data(telomerase)
fit <- fitcopula(data=telomerase,
SID = "ID",
copula="fgm",
iter = 400,
warmup = 100,
seed=1,
cores=1)
printcopula(object=fit)
print(fit$model, pars=c('p_i'), digits=4, prob=c(0.025, 0.975))Run the code above in your browser using DataLab