This function calculates CIs for a population correlation coefficient.
For Pearson correlation, "normal" CIs are available (by stats::cor.test()).
Also bootstrap CIs are supported (by default "bca", and the only option for
rank correlations).
ci_cor(iris[1:2])
ci_cor(iris[1:2], type = "bootstrap", R = 999) # Use larger Rci_cor(iris[1:2], method = "spearman", type = "bootstrap", R = 999) # Use larger R