Я вычисляю ро Спирмена на небольших наборах парных рейтингов. Спирмен хорошо известен тем, что неправильно обращается с галстуками. Например, взяв 2 набора из 8 рейтингов, даже если 6 равны в одном из двух наборов, корреляция все еще очень высока:
> cor.test(c(1,2,3,4,5,6,7,8), c(0,0,0,0,0,0,7,8), method="spearman")
Spearman's rank correlation rho
S = 19.8439, p-value = 0.0274
sample estimates:
rho
0.7637626
Warning message:
Cannot compute exact p-values with ties
Значение p <0,05 кажется довольно высокой статистической значимостью для этих данных. Есть ли исправленная версия Спирмена в R? Какова на сегодняшний день лучшая формула для ее вычисления с большим количеством связей?