Kod 4.11 - NMSwBE

# Kod 4.11 Funkcja obliczania wartości wielowymiarowego współczynnika korelacji rang

r1<-function(R) {

hd=(d+1)/(2^d-(d+1))

suma=0

for (i in 1:n)  {

ilocz=1

for (j in 1:d)  {ilocz=ilocz*(1-R[i,j])}

suma=suma+ilocz

}

r1=(suma*2^d/n-1)*hd

}