Współczynnik korelacji rang Spearmana

Animacja wyjaśnia własności korelacji rangowej. Na osiach odłożone są wartości porównywanych zmiennych. Punkty symbolizują obserwacje w próbie. Współczynnik korelacji rang Spearmana jest obliczany w dwóch etapach: najpierw wykonywane jest rangowanie, czyli zastąpienie każdej zaobserwowanej wartości przez jej numer w zbiorze posortowanym rosnąco. Rangowanie pokazane jest tu w formie animacji, aby umożliwić śledzenie losów każdej z obserwacji; naprawdę żadne stadia pośrednie nie występują. Następnie obliczany jest zwykły współczynnik korelacji liniowej Pearsona. Jest on wprawdzie wrażliwy na obserwacje odstające, jednak rangowanie zbliża je do pozostałych, dzięki czemu niweluje się ich zakłócający wpływ na wynik. Dla pokazanych tu danych, przed rangowaniem współczynnik korelacji jest nieistotny statystycznie $(r=0,01561,\alpha =0,8775).$ Ten sam współczynnik dla rang wykazuje już istotną zależność pomiędzy zmiennymi $(r=0,41702,\alpha <0,0001).$ Monotoniczna nieliniowa zależność przekształca się przy rangowaniu w liniową, w wyniku czego współczynnik korelacji liniowej Pearsona, zastosowany do rang, mierzy siłę zależności monotonicznej.

Korelacja rang Spearmana (lub: korelacja rangowa Spearmana, rho Spearmana) – jedna z nieparametrycznych miar monotonicznej zależności statystycznej między zmiennymi losowymi.

Przyjmuje wartości od -1 do +1 włącznie. Skrajny wynik +1 oznacza, że wzrost jednej ze zmiennych zawsze pociąga za sobą wzrost drugiej, wynik -1 oznacza, że wzrost jednej ze zmiennych pociąga za sobą spadek drugiej, 0 to brak wykrytej zależności. W odróżnieniu od klasycznego współczynnika korelacji Pearsona, zależność ta nie musi być liniowa.

Pierwotny pomysł korelowania rang był już znany wcześniej i pochodził od Bineta i Szablon:Link-interwiki, jednak współczynnik ten został solidnie opisany i rozpropagowany dopiero w 1904 przez angielskiego psychologa Charlesa Spearmana. Zauważył on, że w wielu badaniach nie da się zastosować klasycznego współczynnika korelacji lub daje on fałszywie istotne wyniki ze względu na nadmiar obserwacji odstających.

Spearman zdefiniował swój współczynnik jako współczynnik korelacji Pearsona, liczony dla rang zmiennych (stąd nazwa współczynnik korelacji rang)Szablon:Odn. Obecnie stosowanych jest kilka jego wersji, nieznacznie różniących się od siebie. Ich wartości są identyczne w przypadku, gdy obserwacje każdej zmiennej w próbie nie powtarzają się. Jeśli jednak nie jest to prawdą, to współczynnik korelacji dla rang opisuje jedynie wzór Szablon:LinkWzór i jego odmianySzablon:Odn. Mimo to często używany jest też prostszy wzór Szablon:LinkWzór.

Zastosowanie i interpretacja

Korelacja rangowa przyjmuje wartości z przedziału [-1,+1], ich interpretacja jest podobna do klasycznego współczynnika korelacji Pearsona, z zastrzeżeniem: w odróżnieniu od współczynnika Pearsona, który mierzy liniową zależność między zmiennymi, a wszelkie inne związki traktuje jako zaburzone zależności liniowej, korelacja rangowa pokazuje dowolną monotoniczną zależność (także nieliniową)Szablon:Odn.

Model korelacji rangowej zawiera szerszą klasę zależności niż model klasycznego współczynnika korelacji, nie obejmuje jednak wszystkich możliwych zależności. Na przykład zależność okresowa, spotykana często w analizie szeregów czasowych, gdzie nosi nazwę sezonowości, nie jest wykrywana ani przez korelację Pearsona, ani SpearmanaSzablon:Odn.

Jako metoda rangowa, rho Spearmana jest w niewielkim tylko stopniu wrażliwe na obserwacje odstające, dzięki czemu szczególną użyteczność znajduje w analizie danych niskiej jakościSzablon:Odn.

Współczynnik korelacji Spearmana zależy wyłącznie od uporządkowania zaobserwowanych wartości. Może zatem być stosowany do dowolnych zmiennych, których wartości można uporządkować rosnąco. Klasyczny współczynnik korelacji nie ma sensownej interpretacji dla zmiennych na skali porządkowej, gdyż uzależniony jest od różnic między wartościami zmiennych, które dla cech porządkowych nie są określoneSzablon:Odn.

Współczynnik korelacji Spearmana oraz testy jego istotności mogą być stosowane przy dowolnym rozkładzie porównywanych zmiennychSzablon:Odn.

Korelacja rang Spearmana może być też opisana jako nachylenie (współczynnik kierunkowy) prostej najlepiej dopasowanej (w sensie najmniejszych kwadratów) do zbioru par rangSzablon:Odn. Istnieją inne, bardziej egzotyczne interpretacje, nie mają jednak znaczenia praktycznego.

Zależność między zmiennymi losowymi (niezależnie od tego, jakim wskaźnikiem jest mierzona) nie musi oznaczać związku przyczynowo-skutkowego.

Korelacja rang Spearmana zmiennych losowych

Ta wersja ma znaczenie w statystyce teoretycznej. Wartości dowolnych miar statystycznych wyliczanych z próby wygodnie jest uważać za estymatory (przybliżenia) miar liczonych na podstawie rozkładu zmiennej losowej, z którego próba była losowana. W przypadku miar korelacji, dla zmiennych $X$ i $Y$ będzie to dwuwymiarowy rozkład wektora $(X,Y).$

Korelacja rang Spearmana zmiennych losowych $X$ i $Y$ wyrażona jest wzoremSzablon:OdnSzablon:Odn

Szablon:Wzór

gdzie:

$\mathrm{corr}$ – współczynnik korelacji Pearsona,

${\tilde{\text{F}}}_{\zeta }(u)=\frac{F_{\zeta }(u)+F_{\zeta }(u^{-})}{2},$

${\text{F}}_{\zeta }(u)$ – dystrybuanta zmiennej $\zeta$ w punkcie $u,$

${\text{F}}_{\zeta }(u^{-})$ – lewostronna granica ${\text{F}}_{\zeta }(u)$ w punkcie $u.$

Dla ciągłych zmiennych losowych zachodzi ${\tilde{\text{F}}}_{\zeta }(u)={\text{F}}_{\zeta }(u)$ i wzór ten sprowadza się do

Szablon:Wzór

gdzie

${\text{F}}_X(X),{\text{F}}_Y(Y)$ – dystrybuanty porównywanych zmiennych,

$\mathrm{cov}$ – kowariancja.

Korelacja rang Spearmana z próby

W praktyce współczynnik korelacji rang oblicza się dla próby statystycznej. Używane do tego wzory można uważać za estymatory (przybliżenia) korelacji rang danej wzorem Szablon:LinkWzór lub korelacji rang istniejącej w populacji statystycznejSzablon:Odn. Przybliżenia nie są jednak tym samym, co wartość przybliżana. Ich wyniki będą zatem dla odróżnienia oznaczane przez $r_S,$ podobnie jak w literaturze.

Współczynnik obliczany jest w następujący sposóbSzablon:Odn:

1. Dla każdej porównywanej zmiennej dokonywane jest niezależnie rangowanie, czyli:

1a. Zaobserwowane wartości danej zmiennej porządkowane są rosnąco.

1b. Każdej wartości $x_i$ przypisywana jest ranga $\text{R}{x}_i$ równa pozycji danej wartości w rosnącym porządku (najmniejsza uzyskuje rangę 1, kolejna 2 itd.).

1c. W przypadku gdy dana wartość występuje wielokrotnie, każde z wystąpień ma przypisaną tę samą rangę równą średniej arytmetycznej pozycji w rosnącym porządku (tzw. ranga wiązana lub powiązana, ang. tied rank). Tym samym mogą występować rangi ułamkowe, np. ranga 1,5.

2. Po powrocie do pierwotnego porządku wartości w zmiennych obliczana jest korelacja rangowa z jednego ze wzorów omówionych poniżej (Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór). Istnieją różne wzory, gdyż poszczególne estymatory korelacji rang mają różne własności. Jeśli nie ma rang wiązanych, wszystkie one dają ten sam wynik. Jeśli występują rangi wiązane, to wzór Szablon:LinkWzór i jego odmiany są zgodne z pierwotną koncepcją Spearmana. Pomimo to często używany jest wzór Szablon:LinkWzór, jak twierdzi część autorów, nieprawidłowoSzablon:OdnSzablon:OdnSzablon:OdnSzablon:Odn.

Wzory uwzględniające rangi wiązane

Oryginalna propozycja Spearmana

W oryginalnym ujęciu Spearmana, jego korelacja rang jest współczynnikiem korelacji Pearsona liczonym dla rang zmiennych zamiast ich surowych wartości, co jest bezpośrednim przełożeniem wzoru Szablon:LinkWzór na język rangSzablon:OdnSzablon:Odn.

Szablon:Wzór

gdzie:

$\mathrm{corr}$ – klasyczny współczynnik korelacji,

$\text{R}X$ – rangi zmiennej $X$ w próbie,

$\text{R}Y$ – rangi zmiennej $Y$ w próbie.

Wzór ten można uważać za próbkowy odpowiednik wzoru Szablon:LinkWzór. Rozpisanie wzoru na korelację Pearsona prowadzi doSzablon:Odn

Szablon:Wzór

Wersja oparta na różnicy rang

Ten sam estymator można też zapisać w innej, równoważnej wersji jakoSzablon:OdnSzablon:Odn:

Szablon:Wzór

gdzie:

Szablon:Wzór

Szablon:Wzór

Szablon:Wzór

$t_j$ jest liczbą obserwacji w próbie posiadających tę samą $j$-tą wartość rangi zmiennej $X,$

$u_k$ jest liczbą obserwacji w próbie posiadających tę samą $k$-tą wartość rangi zmiennej $Y,$

sumowanie przebiega po wszystkich wartościach rang – wystarczy zsumować rangi wiązane, bo dla pozostałych $t_j^3-t_j=1^3-1=0$ (analogicznie $u_k$); gdy w danej zmiennej nie ma rang wiązanych, $T_X$ lub $T_Y$ jest równe zeru.

Wersja dla tablicy dwudzielczej

Rozkład porządkowych zmiennych losowych w próbie można przedstawić w formie tablicy dwudzielczej (tablicy kontyngencji), w której kolumny odpowiadają uszeregowanym wartościom jednej zmiennej (oznaczonej przez $X$), wiersze uszeregowanym wartościom drugiej zmiennej (oznaczonej przez $Y$), a w komórkach tablicy znajdują się liczności $n_{ij}.$

Wzór Szablon:LinkWzór przyjmuje wtedy postaćSzablon:Odn:

Szablon:Wzór

gdzie:

$i,j$ – odpowiednio indeksy wierszy oraz kolumn,

$m,k$ – odpowiednio liczba wierszy oraz kolumn,

$F_{k+}=\left(\sum _{i=1}^k{n}_{i+}\right)-\frac{1}{2}{n}_{k+},$

$F_{+k}=\left(\sum _{j=1}^k{n}_{+j}\right)-\frac{1}{2}{n}_{+k},$

$n_{i+}$ – suma $i$-tego wiersza,

$n_{+j}$ – suma $j$-tej kolumny,

$n$ – suma całej tabeli.

Dziś estymator Szablon:LinkWzór jest standardowym wzorem używanym np. przez pakiety statystyczne SAS oraz SPSSSzablon:Odn, a także w uwzględniających rangi wiązane pracach naukowych z dziedziny statystykiSzablon:Odn. W podręcznikach statystyki oraz w pracach naukowych z innych dziedzin nadal jednak popularny jest podany dalej wzór Szablon:LinkWzór, ze względu na stopień komplikacji wzorów Szablon:LinkWzór lub Szablon:LinkWzór, utrudniający ręczne obliczenia, mimo że w obliczeniach wykonywanych na komputerze wzór Szablon:LinkWzór jest nawet prostszy w zastosowaniu.

Niekiedy estymator Szablon:LinkWzór/Szablon:LinkWzór/Szablon:LinkWzór/Szablon:LinkWzór nazywany jest „skorygowaną korelacją rangową”.

Poprawka do wzoru na $r_s$ w przypadku rang wiązanych wynikająca z zastosowania wzorów Szablon:LinkWzór/Szablon:LinkWzór/Szablon:LinkWzór/Szablon:LinkWzór jest obszernie dyskutowana w monografii Kendalla. Lehmann wykazuje asymptotyczną zbieżność do rozkładu normalnego, w przypadku, gdy liczba rang wiązanych jest ograniczona.

Wzór dla rang wiązanych powstałych przez agregację

Powtarzające się wartości zmiennych, a tym samym rangi wiązane, mogą powstawać na dwa sposoby w zależności od natury badanego zjawiska:

1. jako odzwierciedlenie prawdziwych powtarzających się wartości w populacji – przykładem może być zastosowanie korelacji rangowej do obliczania zależności dwóch zmiennych na skali porządkowej, takich jak wykształcenie: dwie osoby mogą mieć to samo wykształcenie, powstanie wówczas ranga wiązana;
2. jako efekt pogrupowania (agregacji) wartości zmiennych w kilka rozłącznych przedziałów, a następnie przypisania jednej wartości (centroid) do każdej z nich – taka sytuacja ma miejsce np. gdy opracowywane są wyniki ankiet, w których badanych poproszono o przedziałowe określenie jakiejś wielkości, takiej jak dochód. Rangi wiązane powstają wówczas w sposób sztuczny.

Podczas agregacji tracona jest informacja o zróżnicowaniu obserwacji wewnątrz każdego przedziału, co sprawia, że zmienne, które przed agregacją nie miały identycznych rang, po agregacji mogą już mieć taki sam porządek. Agregacja jest zwykle zabiegiem wymuszonym warunkami badania, którego wpływ na wyniki powinien być jak najmniejszy. Przydatny byłby więc estymator, szacujący korelację rangową zmiennych przed agregacją na podstawie danych po agregacji. Taki estymator osiągałby wartości $\pm 1$ tylko przy próbie bez rang wiązanych.

Kendall proponuje aby w przypadku rang wiązanych powstałych sztucznie stosować w mianowniku wariancje takie, jak gdyby rang wiązanych nie było (gdyż tak jest w hipotetycznej nieskończonej populacji, dla której korelacja rangowa jest estymowana). Uzyskany w ten sposób estymator jest wartością oczekiwaną współczynnika korelacji rang obliczonego dla tych samych zmiennych przed agregacją (przy założeniu, że każda kombinacja rang prowadząca po agregacji do obserwowanej próby jest jednakowo prawdopodobna)Szablon:Odn.

Postać ogólna

W ogólnym przypadku, po uwzględnieniu rang wiązanych, wariancja rang wynosi:

$\mathrm{var}(\text{R}X)=\frac{n(n+1)}{12}-\frac{T_X}{n-1},$

gdzie współczynnik $T_X$ jest zdefiniowany tak jak wcześniej, wzorem Szablon:LinkWzór.

W przypadku braku rang wiązanych, $T_X=T_Y=0,$ wariancje są zależne tylko od $n,$ w szczególności nie zależą od rozkładu zmiennych przed rangowaniem:

$\mathrm{var}(\text{R}X)=\mathrm{var}(\text{R}Y)=\frac{n(n+1)}{12}.$

Niezależnie od tego, czy pojawiły się rangi wiązane, czy nie, średnia rang jest zależna jedynie od liczności próby:

$\text{E}(\text{R}X)=\text{E}(\text{R}Y)=\frac{n+1}{2}.$

Podstawiając powyższe równania do wzoru Szablon:LinkWzór, uzyskuje się:

Szablon:Wzór

Wersja oparta na różnicy rang

Ten estymator można zapisać w równoważnej postaci jakoSzablon:Odn:

Szablon:Wzór

gdzie $d_i,T_X,T_Y$ są zdefiniowane, jak wyżej, wzorami Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór.

Wersja dla tablicy dwudzielczej

Dla tablic dwudzielczych estymator Szablon:LinkWzór przyjmuje postaćSzablon:Odn

Szablon:Wzór

gdzie:

$U_i=\frac{n_{1+}+\dots +n_{i+}}{n},$

$V_j=\frac{n_{+1}+\dots +n_{+j}}{n},$

$n_{i+}$ – suma $i$-tego wiersza,

$n_{+j}$ – suma $j$-tej kolumny,

$n$ – suma całej tabeli.

Istnieje jeszcze inny estymator dla tablic dwudzielczych, zaproponowany przez StuartaSzablon:OdnSzablon:Odn.

Wzór nieuwzględniający rang wiązanych

Często stosowanym estymatorem jestSzablon:Odn:

Szablon:Wzór

gdzie (tak jak wcześniej)

$d_i=\text{R}{x}_i-\text{R}{y}_i$

to różnica między rangami zmiennych $X$ i $Y$ dla obserwacji $i.$

W przypadku gdy nie ma rang wiązanych (połączonych), czyli wartości nie powtarzają się w obrębie próby dla żadnej ze zmiennych z osobna, wzór Szablon:LinkWzór daje te same wyniki, co każdy z podanych wcześniej estymatorów Szablon:LinkWzór i Szablon:LinkWzór. Jeśli choć jedna ranga jest wiązana, każdy z nich daje inny wynik.

Wzór Szablon:LinkWzór jest używany ze względu na prostotę obliczeń istotną dla kalkulacji wykonywanych bez pomocy komputera i do dziś jest popularny w podręcznikachSzablon:Odn. Estymator ten ma jednak nieoczekiwane własności w przypadku wystąpienia rang wiązanych, np.

• nie jest wówczas prawdą, iż $r_S(X,-Y)=-r_S(X,Y),$
• nie jest wtedy zgodny z pierwotną definicją korelacji rang Spearmana jako zwykłego współczynnika korelacji liczonego dla rangSzablon:Odn,
• dla zmiennych dyskretnych, minimalną wartością jego granicy, przy rozmiarze próby dążącym do nieskończoności, jestSzablon:Odn

Szablon:Wzór

gdzie:

$m_X$ to liczba różnych wartości przyjmowanych przez zmienną $X$

$m_Y$ to liczba różnych wartości zmiennej $Y.$

• Wynika stąd, że estymator ten jest dla zmiennych dyskretnych niezgodny i asymptotycznie obciążony, w odróżnieniu od estymatora Szablon:LinkWzór.

Część autorów uważa, że można ten estymator stosować tylko przy braku rang wiązanych, w przeciwnym wypadku jego stosowanie jest błędemSzablon:OdnSzablon:Odn. Inni autorzy stosują go także wówczasSzablon:OdnSzablon:OdnSzablon:Odn. Niektórzy uważają, że wzór można stosować, jeśli rang wiązanych jest nie więcej niż jedna czwarta ogółu i nie występują rangi wiązane z więcej niż dwóch obserwacjiSzablon:OdnSzablon:Odn. Pakiety statystyczne SAS oraz SPSS używają podanego wcześniej bardziej ogólnego wzoru Szablon:LinkWzórSzablon:Odn.

Niekiedy wzór Szablon:LinkWzór nazywany jest „nieskorygowaną korelacją rangową” w odróżnieniu od „skorygowanej korelacji rangowej” Szablon:LinkWzór. Jest to związane z postacią wzoru Szablon:LinkWzór, który przypomina wzór Szablon:LinkWzór z dodaną „korektą na rangi wiązane”.

Własności

• Współczynnik jest unormowanySzablon:Odn:

$-1⩽r_S(X,Y)⩽+1;$

• Im bardziej wartości oddalone są od zera, tym większa siła związku między zmiennymi;
• Gdy każda zmienna jest ściśle rosnącą funkcją drugiej (np. $x_i=y_i^3$), występuje idealna zgodność rang i ich korelacja przyjmuje wartość +1Szablon:Odn;
  • W szczególności wartość ta jest przyjmowana, gdy zmienna jest korelowana sama ze sobą:

$r_S(X,X)=+1;$

• Gdy każda zmienna jest ściśle malejącą funkcją drugiej zmiennej, występuje maksymalna niezgodność rang i ich korelacja przyjmuje wartość –1Szablon:Odn;
  • W szczególności wartość ta jest przyjmowana, gdy zmienna $X$ korelowana jest z $-X:$

Szablon:Wzór

• Dla niezależnych zmiennych losowych wartością oczekiwaną estymatorów jest 0, a rozkład każdego z nich nie zależy od rozkładu zmiennych przed rangowaniem;
• Zachodzi symetria ze względu na zamianę zmiennych:

$r_S(X,Y)=r_S(Y,X);$

• Zachodzi symetria ze względu na zmianę znaku zmiennej:

Szablon:Wzór

W przypadku wystąpienia rang wiązanych część z tych własności nie jest spełniona dla niektórych estymatorów. Dla estymatora Szablon:LinkWzór nie są prawdziwe własności Szablon:LinkWzór i Szablon:LinkWzór, a estymator Szablon:LinkWzór/Szablon:LinkWzór/Szablon:LinkWzór nie osiąga wartości $\pm 1.$

Matematyczne własności rho Spearmana mają związek z tożsamością Czebyszewa oraz twierdzeniem o ciągach jednomonotonicznych.

Przykład

numery obserwacji $i$	1	2	3	4	5
wartości $x_i$	1,1	1,57	0,51	1,1	1,1
wartości $y_i$	1,2	1	2,3	1	18
rangi $\text{R}{x}_i$	3	5	1	3	3
rangi $\text{R}{y}_i$	3	1,5	4	1,5	5
kwadraty różnic rang $d_i^2$	0	12,25	9	2,25	4

We wzorach Szablon:LinkWzór, Szablon:LinkWzór i Szablon:LinkWzór wykonywane są pośrednie obliczenia:

$\sum _{i=1}^n{d}_i^2=(3-3{)}^2+(5-1,5{)}^2+(1-4{)}^2+(3-1,5{)}^2+(3-5{)}^2=27,5.$

W Szablon:LinkWzór i Szablon:LinkWzór także:

$T_X=\frac{1}{12}(3^3-3)=2$

(jest jedna ranga wiązana, mają ją trzy obserwacje),

$T_Y=\frac{1}{12}(2^3-2)=0,5$

(jest jedna ranga wiązana, mają ją dwie obserwacje).

Po podstawieniu do wzorów otrzymuje się:

estymator	wynik
Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór	$r_S\approx -0,5735$
Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór	$r_S=-0,5$
Szablon:LinkWzór	$r_S=-0,375$

Testowanie istotności statystycznej

Szablon:Wikibooks Aby przetestować istotność statystyczną korelacji rangowej, wykorzystuje się fakt, iż przy założeniu hipotezy zerowej o niezależności zmiennych losowych $X$ i $Y$ oraz niezależności od siebie par $(x_i,y_i)$ rozkład statystyki

$t=r_S\sqrt{\frac{n-2}{1-r_S^2}}$

korelacji rangowej dąży wraz ze wzrostem liczebności próby do rozkładu Studenta o $n-2$ stopniach swobody, gdzie $n$ jest licznością próbySzablon:Odn. Po obliczeniu tej statystyki korzysta się z tablic rozkładu Studenta lub komputera w celu obliczenia poziomu istotności $\alpha .$

Rozkład ten jest wyprowadzany przy założeniu braku rang wiązanych, jednak Kendall twierdzi, że w przypadku istnienia rang wiązanych poprawka do testu nie jest koniecznaSzablon:Odn.

Inni autorzy z kolei zalecają stosowanie poprawki na ciągłość (przez dodanie lub odjęcie $\frac{6}{n^3-n}$)Szablon:Odn.Szablon:Wikibooks Dla liczebności próby dążącej do nieskończoności, rozkład rho Spearmana dąży do rozkładu normalnego o wartości oczekiwanej równej prawdziwej wartości ${\rho }_S$ w populacji i wariancjiSzablon:Odn $\frac{1}{n-1},$ stąd używana jest też inna statystykaSzablon:Odn:

$z=r_S\sqrt{n-1},$

której rozkład przy założeniu hipotezy zerowej dąży wraz ze wzrostem liczności próby do standardowego rozkładu normalnego $\text{N}(0,1).$

Część autorów sugeruje stosowanie poprawki na ciągłość, która jednak nie daje znaczącego wzrostu dokładności przybliżeniaSzablon:Odn.

Plik:Rho Spearmana dla H0, N=5.png

Rozkład rho Spearmana dla 5-elementowej próby z dwóch niezależnych zmiennych losowych

Postulowano również stosowanie średniej arytmetycznej statystyki opartej na rozkładzie $t$ i normalnymSzablon:Odn.

Dla małych prób wszystkie te wzory są niedokładne (statystyka oparta na rozkładzie $t$ jest zazwyczaj nieco lepsza niż statystyka $z$), ale można sprawdzić komputerowo wszystkie permutacje rang lub skorzystać z tablic. Często stosuje się przy tym inną statystykę (nazywaną czasem Hotellinga-Pabst): $\sum _{i=1}^n{d}_i^2.$ Należy przy tym zauważyć, że między tymi statystykami zachodzi związek liniowy, a więc ich stosowanie prowadzi do identycznych wyników. Dokładne tablice rozkładu dla wielkości populacji do 25 można znaleźć w pracySzablon:Odn.

Hipotezą alternatywną może być albo:

$H_1:{\rho }_S\ne 0$

(co prowadzi do dwustronnego obszaru krytycznego) albo

$H_1:{\rho }_S>0$ lub $H_1:{\rho }_S<0$

(co prowadzi do jednostronnego obszaru krytycznego).

Dla omawianego powyżej przykładu, dwustronnego obszaru krytycznego i wyliczeń według trzech estymatorów otrzymuje się następujące wartości:

estymator	$r_S$	$t$	$z$	$\alpha$ (przybliżenie rozkładem Studenta)	$\alpha$ (przybliżenie rozkładem normalnym)	$\alpha$ (wartość dokładna)
Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór	$-0,5735$	$-1,2127$	$-1,2824$	$0,312$	$0,0998$	$0,35$
Szablon:LinkWzór, Szablon:LinkWzór, Szablon:LinkWzór	$-0,5$	$-1,0000$	$-1,118$	$0,391$	$0,1318$	$0,45$
Szablon:LinkWzór	$-0,375$	$-0,7006$	$-0,8385$	$0,534$	$0,2009$	$0,5167$

W tabeli podano wartość $\alpha$ wyliczoną za pomocą przybliżenia rozkładem Studenta, z rozkładu normalnego i wreszcie dokładnie – z tablic. Dla tak małej próby przybliżenie rozkładem Studenta daje różnice rzędu 0,05, co może mieć znaczenie przy określaniu istotności statystycznej. Przybliżenie rozkładem normalnym jest w tym wypadku mniej dokładne. Dla małych prób konieczne jest więc stosowanie tablic lub symulacji komputerowych. Przy liczebności próby dążącej do nieskończoności różnica zmniejsza się i coraz bardziej uzasadnione jest stosowanie rozkładu Studenta, ewentualnie rozkładu normalnego, co jednak z reguły da nieco większy od Studenta błąd wyznaczania istotności.

Znacznie dokładniejsze wyniki przynosi stosowanie przybliżenia opartego na krzywych Pearsona II typu, metodzie opisanej przez OldsaSzablon:Odn. Najdokładniejsze wyniki przynosi przybliżenie oparte na skorygowanym przy pomocy szeregu Edgewortha rozkładzie normalnymSzablon:Odn. Implementację tego przybliżenia jako funkcję w programie Mathematica można znaleźć w pracy. Tam też dokonano porównania różnych metod aproksymacji statystyki SpearmanaSzablon:Odn.

Istnieją też stabelaryzowane rozkłady korelacji rangowej dla innych założeń, np. $\rho =0,4.$ Odpowiednie tabele podaje praca Fritza i Henze’aSzablon:Odn.

Związki z innymi współczynnikami i metodami statystycznymi

Współczynnik korelacji Pearsona

Współczynniki te określają innego rodzaju zależność między zmiennymi (Pearson – zależność liniową, Spearman – dowolną monotoniczną), czasem jednak korelacja rang jest używana jako odporna wersja klasycznego współczynnika korelacji PearsonaSzablon:Odn. W takiej roli widział ją zresztą sam SpearmanSzablon:Odn.

Jest to uzasadnione w przypadku zakładanej liniowej zależności między zmiennymi w warunkach zanieczyszczenia próby obserwacjami odstającymi. Korelacja rangowa jest bowiem znacznie bardziej odporna na obserwacje odstające, które potrafią skrajnie zaburzyć wynik zwykłego współczynnika korelacji Pearsona. Wartości tych dwóch współczynników nie są jednak wtedy równe – korelacja rangowa daje na ogół (nie zawsze) wyniki nieco bliższe zeru.

W szczególności dla dwuwymiarowego rozkładu normalnego zachodziSzablon:Odn:

${\rho }_S=\frac{6\arcsin (\rho /2)}{\pi },$

$\rho =2\sin \frac{\pi {\rho }_S}{6},$

gdzie:

${\rho }_S$ – współczynnik korelacji Spearmana pomiędzy zmiennymi (wzór Szablon:LinkWzór),

$\rho$ – współczynnik korelacji Pearsona pomiędzy zmiennymi.

Zależność ta jest ścisła przy braku rang wiązanych i nieskończonej populacji. Dla skończonej próby zależność między estymatorami Spearmana $r_S$ i Pearsona $r$ różni się od tego wyidealizowanego przypadku. Wartość oczekiwana wynosi wtedy dla każdego z podanych estymatorówSzablon:Odn:

$\text{E}(r_S)=\frac{6}{\pi }(\frac{n-1}{n+1}\arcsin \frac{1}{2}\rho +\frac{1}{n+1}\arcsin \rho ).$

Współczynnik korelacji rang Spearmana jest więc estymatorem obciążonym (także asymptotycznie) i niezgodnym współczynnika korelacji Pearsona. Naturalnie na tej samej zasadzie współczynnik korelacji Pearsona będzie obciążonym, niezgodnym i nieefektywnym estymatorem korelacji rangowej SpearmanaSzablon:Odn.

Rho Spearmana jest też przy założeniu rozkładu dwuwymiarowego normalnego mniej efektywne niż współczynnik korelacji Pearsona liczony klasycznym wzorem, bez rangowania. Dla dwuwymiarowego rozkładu normalnego błąd standardowy korelacji Pearsona liczonej za pomocą wzoru:

Szablon:Wzór

(gdzie $r_S$ to dowolny z estymatorów rho Spearmana), jest ok. 1,88 raza większy od błędu korelacji liczonej za pomocą klasycznego wzoru bez rangowaniaSzablon:Odn:

$r=\frac{\mathrm{cov}(X,Y)}{\sqrt{\mathrm{var}X\mathrm{var}Y}}$

Jednak, gdy obserwacje nie spełniają założenia o normalności rozkładu, szczególnie gdy pojawiają się obserwacje odstające, wzór Szablon:LinkWzór często daje lepsze oszacowanie korelacji liniowej. Jeszcze dokładniejszy jest współczynnik tau KendallaSzablon:Odn.

Współczynnik korelacji Pearsona nie zakłada żadnej postaci rozkładu porównywanych zmiennych, jednak wzory na jego istotność statystyczną zakładają już dwuwymiarowy rozkład normalny. W wielu przypadkach warunek ten nie jest spełniony i nie da się łatwo sprawdzić, czy wyniki korelacji Pearsona są przejawem rzeczywistej zależności. Istotność współczynnika korelacji rangowej daje się zawsze określić, gdyż rozkład rang nie zależy od rozkładu porównywanych zmiennych, o ile nie ma rang wiązanych, a nawet wtedy testy istotności nie są znacząco zaburzoneSzablon:Odn.

Inne miary korelacji rangowej

Korelacja rangowa to szersze pojęcie niż korelacja rang Spearmana. Korelacja to ogólnie w statystyce zależność zmiennych losowych. Miary tej zależności wyliczane na bazie rang zwane są miarami korelacji rangowej.

Wymienione poniżej miary nie są jednak uznawane za estymatory korelacji rang Spearmana – są odrębnymi współczynnikami o odrębnej interpretacji. Istnieją też inne, nie wymienione tutaj, współczynniki korelacji rangowej.

Miara Spearmana

Miara Spearmana (ang. Spearman’s footrule) to współczynnik zaproponowany w tej samej pracy, co rho SpearmanaSzablon:Odn, liczony podobnie jak we wzorze Szablon:LinkWzór, jednak z wartością bezwzględną w miejsce kwadratu i z wynikającą z tego inną normalizacją:

$R=1-\frac{3\sum _{i=1}^n|d_i|}{n^2-1}$

Jak pokazał Pearson, współczynnik ten nie ma dobrych własności statystycznych, w szczególności choć osiąga +1, nie osiąga nigdy wartości -1, z wyjątkiem przypadku $n=2$Szablon:Odn.

Tau Kendalla

Inną miarą korelacji rangowej dwóch zmiennych jest tzw. tau Kendalla. Między tymi wartościami zachodzą nierównościSzablon:Odn:

$\frac{3}{2}\tau -\frac{1}{2}⩽r_S⩽\frac{1}{2}+\tau -\frac{1}{2}{\tau }^2\text{ dla }\tau ⩾0,$

$\frac{1}{2}{\tau }^2+\tau -\frac{1}{2}⩽r_S⩽\frac{3}{2}\tau +\frac{1}{2}\text{ dla }\tau ⩽0.$

Podawane jest też grubsze oszacowanieSzablon:Odn:

$-1⩽3\tau -2r_S⩽1.$

Można też pokazać, że jeśli przedstawić łączny rozkład dwuwymiarowy zmiennych $X$ i $Y$ w postaci unormowanej do jedności macierzy prawdopodobieństwa $m\times k$ o elementach $[p_{ij}],$ wówczas obydwie te wielkości dają się przedstawić za pomocą średnich ważonych z minorów stopnia drugiegoSzablon:Odn. W tym ujęciu rho Spearmana jest ważonym współczynnikiem tau Kendalla.

$r_S=\sum _{i=1}^{m-1}\sum _{q=i+1}^m\sum _{j=1}^{k-1}\sum _{r=j+1}^k{w}_{ijqr}\left|\begin{array}{cc}{p}_{ij}& p_{ir}\\ p_{qj}& p_{qr}\end{array}\right|,$

$\tau =\sum _{i=1}^{m-1}\sum _{q=i+1}^m\sum _{j=1}^{k-1}\sum _{r=j+1}^{k}2\left|\begin{array}{cc}{p}_{ij}& p_{ir}\\ p_{qj}& p_{qr}\end{array}\right|,$

gdzie:

$w_{ijqr}=12({\text{Sc}}_{row}(q)-{\text{Sc}}_{row}(i))({\text{Sc}}_{col}(r)-{\text{Sc}}_{col}(j)),$

${\text{Sc}}_{row}(i)=\left(\sum _{a=1}^{i-1}\sum _{b=1}^k{p}_{ab}\right)+\frac{1}{2}\sum _{b=1}^k{p}_{ib},$

${\text{Sc}}_{col}(j)=\left(\sum _{a=1}^m\sum _{b=1}^{j-1}{p}_{ab}\right)+\frac{1}{2}\sum _{a=1}^m{p}_{aj}.$

Kendall i Stuart pokazali, że dla niezależnych zmiennych korelacja między tau i rho wynosi co najmniej 0,98 i dąży do 1 dla $n\to \mathrm{\infty }.$ (Wspólny rozkład tau Kendalla i rho Spearmana w swojej monografii podaje Kendall.) Nie oznacza to jednak, że ich wyniki są proporcjonalne dla zmiennych zależnych, a dla takich właśnie na ogół liczy się korelację.

Uogólnienia rho Spearmana

Chi kwadrat Friedmana

Rho Spearmana jest znormalizowaną i przeskalowaną do przedziału $[-1,1]$ miarą chi kwadrat Friedmana dla dwóch zmiennych. Jeśli wartość chi kwadrat Friedmana wynosi $X,$ to

$r_S=\frac{X}{n-1}-1.$

L Page’a

Kolejnym uogólnieniem rho Spearmana na przypadek wielu zmiennych jest test L Page'a. Korelację rangową można stosować jako metodę sprawdzania, czy zmienna $X$ ma ten sam porządek rang co zmienna $Y.$ Test L Page’a podaje z jakim prawdopodobieństwem ciąg zmiennych $X_1,X_2,\dots ,X_k$ ma pewne zadane ustawienie. Jego wynik można też podać w formie współczynnika z zakresu $[-1,1],$ który dla $k=1$ sprowadza się do korelacji rang Spearmana.

Analiza odpowiedniości oparta na rho Spearmana

Klasyczna analiza odpowiedniości (inna nazwa: analiza korespondencji) jest metodą statystyczną, która wszystkim możliwym wartościom dwóch zmiennych nominalnych przyporządkowuje takie liczby (tzw. skory), aby przy pewnych założeniach maksymalizować współczynnik korelacji Pearsona między tymi zmiennymi.

Istnieje odpowiednik klasycznej analizy odpowiedniości, zwany gradacyjną analizą odpowiedniości (ang. Grade Correspondence Analysis; GCA), który maksymalizuje rho Spearmana lub tau Kendalla.

Krytyka

Te same własności rho Spearmana, które zwolennicy metod rangowych uważają za zalety, przeciwnicy mają za wady. Sam Spearman, który traktował swój współczynnik wyłącznie jako odporne na obserwacje odstające przybliżenie korelacji Pearsona, uważał za wadę fakt, że mierzy ona także zależność nieliniowąSzablon:Odn.

Twórca klasycznego współczynnika korelacji, Karl Pearson, krytykował niezależność od rozkładu korelacji rang:

Przy okazji tej krytyki pierwszy raz w historii użyto określenia „korelacja rangowa”.

Historia

Pomysł korelowania rang był już znany przed Spearmanem i pochodził od Bineta i Henriego. Redakcja czasopisma Biometrika w przypisie pracy Studenta zaznaczyła, że „ich wywód był bardzo niejasny i chyba nie zauważyli, że korelacja zmiennych różni się od korelacji rang”.

Współczynnik został solidnie opisany, zbadany i rozpropagowany dopiero w 1904 przez angielskiego psychologa Charlesa SpearmanaSzablon:Odn. Praca Spearmana była opisem różnych metod korelacji dla psychologów, m.in. korelacji Pearsona dla rang (choć Spearman nie zapisał swojej metody w postaci wzoru). Autor zauważył też, że w wielu badaniach nie da się zastosować klasycznego współczynnika korelacji Pearsona lub daje on nieistotne albo fałszywie istotne wyniki ze względu na nadmiar obserwacji odstających, natomiast problemy te znikają po rangowaniuSzablon:Odn. Nadal traktował jednak korelację rang jedynie jako poszerzenie możliwości współczynnika korelacji Pearsona, choć znał różnice między nimi.

Koncepcja rang wiązanych nie była jeszcze znana w początkach XX wieku – została ona wprowadzona później przez PearsonaSzablon:Odn. Wówczas znany był już wzór Szablon:LinkWzór, wyprowadzony naturalnie przy założeniu braku rang wiązanych. Student (William Sealy Gosset) w pracy z 1921 zauważył, że wzór Szablon:LinkWzór nie zgadza się z definicją Spearmana w przypadku rang wiązanych (sprowadzającą się wówczas do wzoru Szablon:LinkWzór) i podał wzór Szablon:LinkWzór, wyprowadził też wzór na wariancję korelacji rangowej.

W 1948 Maurice Kendall napisał monografię Rank Correlation Methods, w której szczegółowo zbadał własności rho Spearmana i związki z własnym współczynnikiem tau Kendalla.

Nacisk Spearmana na budowę stabilnych metod statystycznych, niezależnych od konkretnych parametrów rozkładu, został uogólniony w filozofii nauki do tzw. zasady Spearmana (ang. Spearman’s Principle):

Podejście to dało początek całej nowej dziedzinie statystyki, zwanej statystyką odpornościową (ang. robust statistics), zajmującej się budową metod statystycznych odpornych na obserwacje odstające.

Oznaczenia

W literaturze spotyka się różne oznaczenia korelacji rang Spearmana: Szablon:Układ wielokolumnowy

Zobacz też

• metody rangowe
• przegląd zagadnień z zakresu statystyki
• współczynnik korelacji Pearsona

Bibliografia

• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj książkę
• Szablon:Cytuj książkę
• Szablon:Cytuj
• Szablon:Cytuj książkę
• Szablon:Cytuj książkę
• Szablon:Cytuj książkę
• Szablon:Cytuj
• Szablon:Cytuj książkę
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj
• Szablon:Cytuj książkę

Linki zewnętrzne

• Szablon:Cytuj stronę – internetowy kalkulator obliczający korelację rangową według wzorów Szablon:LinkWzór i Szablon:LinkWzór.
• Szablon:Cytuj stronę – przykład zastosowania korelacji rangowej do badania zależności cen w Barcelonie od położenia sklepu w mieście