Вы можете использовать bootstrap для вычисления доверительного интервала для rho:
1) Выполнить функцию для извлечения оценки cor.test (не забудьте поставить индексы, чтобы загрузчик мог выбирать данные):
rho <- function(x, y, indices){
rho <- cor.test(x[indices], y[indices], method = c("spearman"))
return(rho$estimate)
}
2) Используйте пакет boot
для начальной загрузки вашей оценки:
library(boot)
boot.rho <- boot(x ,y=y, rho, R=1000)
3) Возьмите доверительный интервал:
boot.ci(boot.rho)