Я использовал функцию rcorr.cens в пакете R Hmis c. Заметили, что определения для выходов «соответствующие пары» и «согласный» трудно понять. В этом примере используются 4 точки выборки: A (1,0), B (2,1), C (3,0) и D (4,1).
Я ожидаю, что количество пар будет 4 * (4-1) / 2 = 6, но функция дает 8, согласованные пары должны быть AB, AD и C -D, поэтому 3, но функция дает 6. Как эта функция R получается с этими выходами?
library(Hmisc)
x=seq(1,4,1)
s=rep(0,4)
s[c(2,4)]=1
rcorr.cens(x,s)
C Index Dxy S.D. n missing uncensored Relevant Pairs
0.75 0.50 0.50 4.00 0.00 4.00 8.00
Concordant Uncertain
6.00 0.00