Tau Kendalla

Tau Kendalla – statystyka będąca jedną z miar monotonicznej zależności dwóch zmiennych losowych. Służy w praktyce do opisu korelacji między zmiennymi porządkowymi.

Tau Kendalla dwóch zmiennych losowych

Ogólnie tau Kendalla jest definiowane jako:

Szablon:Wzór

Można też wyrazić tau Kendalla w terminach kopuł:

$\tau =4\underset{0}{\overset{1}{\int }}\underset{0}{\overset{1}{\int }}C(u_1,u_2)\mathrm{d}C(u_1,u_2)-1.$

Tau Kendalla z próby

Zwykle estymuje się tau Kendalla dane wzorem Szablon:LinkWzór na podstawie próby statystycznej.

Aby je obliczyć, należy zestawić obserwacje z próby we wszystkie możliwe pary, a następnie podzielić te pary na trzy możliwe kategorie:

• pary zgodne – porównywane zmienne w obrębie tych dwóch obserwacji zmieniają się w tę samą stronę, tzn. albo w pierwszej obserwacji obydwie są większe niż w drugiej, albo obydwie mniejsze. Liczba takich par w próbie będzie dalej oznaczana przez $P.$
• pary niezgodne – zmienne zmieniają się w przeciwną stronę, to znaczy jedna z nich jest większa dla tej obserwacji w parze, dla której druga jest mniejsza. Liczba takich par w próbie będzie oznaczana przez $Q.$
• pary wiązane – jedna ze zmiennych ma równe wartości w obydwu obserwacjach. Liczba takich par w próbie będzie oznaczana przez $T.$

Estymator tau Kendalla obliczany jest wówczas ze wzoru:

$\tau =\frac{P-Q}{P+Q+T}.$

Warto zauważyć, że:

$P+Q+T=(\genfrac{}{}{0ex}{}{N}{2})=\frac{N(N-1)}{2},$

gdzie $N$ to liczność próby, tak więc wzór ten można przedstawić w postaci:

$\tau =2\frac{P-Q}{N(N-1)}.$

Interpretacja

Tau Kendalla stanowi różnicę między prawdopodobieństwem, że porównywane zmienne będą układały się w tym samym porządku dla dwóch obserwacji, a prawdopodobieństwem, że ułożą się w przeciwnym porządku.

Tau Kendalla przyjmuje wartości od –1 do 1 włącznie. +1 oznacza, że każda ze zmiennych rośnie przy wzroście drugiej. –1 oznacza że każda maleje przy wzroście drugiej. Tym samym tau Kendalla, podobnie jak korelacja rangowa i rho Spearmana jest miarą monotonicznej zależności zmiennych losowych.

Ponieważ liczba par zgodnych, niezgodnych i wiązanych, a zatem również wartość tau nie zmienia się pod wpływem rangowania, tau Kendalla można uważać za jedną ze statystycznych metod rangowych.

Istotność statystyczna

Istotność statystyczna wyniku tau Kendalla jest równa jeden minus prawdopodobieństwo, że taki wynik może powstać przypadkowo dla zmiennych niezależnych. Ponieważ tau Kendalla jest metodą rangową, a rozkład rang jest niezależny od rozkładu zmiennych rangowanych (o ile nie ma rang wiązanych), więc określenie istotności nie wymaga żadnych założeń odnośnie do rozkładów w próbie.

Szablon:Wikiźródła Dla małych prób tablice podają Siegel i Castellan (1988). Rozkład tau bardzo szybko dąży do rozkładu normalnego, więc dla prób $N⩾10$ można już stosować rozkład normalny (Hays 1988).

Wersje

Tau Kendalla ma dobrą interpretację i właściwości statystyczne, jednak wartości 1 i –1 nie wystąpią, gdy którakolwiek z porównywanych zmiennych ma tę samą wartość dla dwóch obserwacji. Intuicyjnie wydawałoby się, że miara pokazująca zależność dwóch zmiennych powinna dawać maksymalną wartość (czyli tu 1) jeśli porównamy jakąś zmienną z nią samą. W przypadku tau Kendalla nie będzie to jednak prawda, o ile zmienna przyjmie dla dwóch obserwacji tę samą wartość.

Zaproponowano szereg innych estymatorów z próby tau Kendalla, które nie mają tej właściwości. Dwa z nich, nazywane są tau-b i tau-c. Dla odróżnienia miara opisana wcześniej nazywana bywa tau-a. Miary te różnią się wartością tylko w przypadku występowania par wiązanych.

Współczynnik tau-b (${\tau }_B$) zdefiniowany w następujący sposób:

${\tau }_B=\frac{P-Q}{\sqrt{(n_0-n_1)(n_0-n_2)}}$,

gdzie P to liczba par zgodnych, Q to liczba par niezgodnych, n₀ to liczba wszystkich par ($n_0=N(N-1)/2$), zaś n₁ i n₂ wyznaczone są dla wartości wiązanych za pomocą następujących wzorów:

$n_1=\sum _i{t}_i(t_i-1)/2$, gdzie t_i to liczba powiązanych wartości pierwszej zmiennej dla i-tej grupy wiązań

$n_2=\sum _j{u}_j(u_j-1)/2$, gdzie u_j to liczba powiązanych wartości pierwszej zmiennej dla j-tej grupy wiązań

Współczesne języki programowania stosowane do analizy danych, takie jak R lub Python, stosują domyślnie wersję tau-b.

Tau-c (${\tau }_C$), znane również pod nazwą Tau-c Stewarta-Kendalla, zalecane gdy liczba wiązań jest duża i różna w przypadku obu zmiennych, zdefiniowane jest następująco:

${\tau }_C=\frac{2(P-Q)}{N^2\frac{(m-1)}{m}}={\tau }_A\frac{N-1}{N}\frac{m}{m-1}$,

gdzie m to mniejsza z dwóch liczb: liczby unikalnych wartości pierwszej zmiennej i liczby unikalnych wartości drugiej zmiennej.

Związki z innymi miarami

Współczynnik korelacji

Współczynnik korelacji Pearsona służy do oceny liniowej zależności pomiędzy dwiema zmiennymi. Tau Kendalla mierzy zależność monotoniczną. Współczynniki te mierzą co innego, więc ich wartości nie są bezpośrednio porównywalne.

Tau Kendalla jako metoda rangowa jest odporne na obserwacje odstające, w przeciwieństwie do współczynnika korelacji.

Tau Kendalla jako metoda nieparametryczna nie czyni żadnych założeń dotyczących rozkładu zmiennych w populacji, nadaje się zatem do dowolnych rozkładów. Współczynnik korelacji ma poprawne wyprowadzenie matematyczne tylko dla wielowymiarowego rozkładu normalnego, dla innych rozkładów może nie mieć sensownej interpretacji, niepoprawnie działają również testy istotności statystycznej.

Rho Spearmana

Tau Kendalla i rho Spearmana to dwie miary zależności monotonicznej zmiennych losowych. Inne jest jednak ich wyprowadzenie i interpretacja. Tau Kendalla należy interpretować w kategoriach prawdopodobieństwa, natomiast rho Spearmana należy interpretować w kategoriach procentu wariancji rang jednej zmiennej, wyjaśnianej przez drugą zmienną.

Między tymi wartościami zachodzi nierówność:

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

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

Podawane jest też grubsze oszacowanie: $-1⩾3\tau -2{\rho }_s⩾1.$

Można też pokazać, że jeśli przedstawimy łączny rozkład dwuwymiarowy zmiennych X i Y w postaci unormowanej do 1 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 drugiego (Kowalczyk 2000):

$\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|,$

$\rho =\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|,$

gdzie:

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

a ${\mathrm{Sc}}_{row}$ i ${\mathrm{Sc}}_{col}$ to tzw. skory gradacyjne (grade scores):

${\mathrm{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},$

${\mathrm{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}.$

Analiza odpowiedniości oparta na tau

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. scory), aby przy pewnych założeniach maksymalizować współczynnik korelacji Pearsona między tymi zmiennymi.

Zespół prof. Elżbiety Pleszczyńskiej opracował odpowiedniki klasycznej analizy odpowiedniości, zwane gradacyjną analizą odpowiedniości (GCA), które maksymalizują rho Spearmana oraz tau Kendalla.

Zobacz też

• gamma Kruskala
• korelacja rangowa
• rho Spearmana

Bibliografia

• Szablon:Cytuj książkę
• Szablon:Cytuj książkę
• Szablon:Cytuj książkę
• Szablon:Cytuj książkę
• Szablon:Cytuj książkę Biometrika, 30, 81-89.
• Szablon:Cytuj książkę
• Szablon:Cytuj książkę
• Szablon:Cytuj książkę Journal of the American Statistical Association, 53(284), 814-861.
• Szablon:Cytuj książkę
• Statsoft: Statistica Electronic Manual