Spearmans rangcorrelatiecoëfficiënt

Spearmans rangcorrelatiecoëfficiënt, of kortweg Spearmans ρ (rho), is in de statistiek de correlatiecoëfficiënt gebaseerd op de rangnummers van de data in plaats van op de data zelf. Het is daarmee een verdelingsvrije maat voor correlatie, ook geschikt voor data die slechts op ordinale schaal gemeten zijn. De coëfficiënt is genoemd naar z'n bedenker, de psychometricus Charles Spearman.

Spearmans rangcorrelatiecoëfficiënt is een eenvoudig geval van de product-momentcorrelatiecoëfficiënt van Karl Pearson, maar dan berekend voor de rangnummers $(r_{i})$ en $(s_{i})$ van de $n$ dataparen.

$\rho _{S}={\frac {\sum _{i}(r_{i}-{\bar {r}})(s_{i}-{\bar {s}})}{{\sqrt {\sum _{i}(r_{i}-{\bar {r}})^{2}}}{\sqrt {\sum _{i}(s_{i}-{\bar {s}})^{2}}}}}=$

$={\frac {\sum _{i}(r_{i}s_{i}-{\bar {r}}{\bar {s}})}{\sum _{i}(r_{i}-{\bar {r}})^{2}}}={\frac {-{\tfrac {1}{2}}\sum _{i}(r_{i}-s_{i})^{2}+\sum _{i}r_{i}^{2}-n{\bar {r}}^{2}}{\sum _{i}r_{i}^{2}-n{\bar {r}}^{2}}}=$

$=1-{\frac {\sum _{i}D_{i}^{2}}{2(\sum _{i}r_{i}^{2}-n{\bar {r}}^{2})}}$

Omdat

${\bar {r}}={\tfrac {1}{2}}(n+1)$

en

$\sum _{i}r_{i}^{2}={\tfrac {1}{6}}n(n+1)(2n+1)$,

volgt, in het geval dat er geen knopen zijn, de eenvoudiger formule:

$\rho _{S}=1-{\frac {6\sum D^{2}}{n(n^{2}-1)}}$

waarin:

$D_{i}=r_{i}-s_{i}$

Knopen

Als er onder de waarden van de variabelen gelijke voorkomen, zogenaamde knopen, kan de bovenstaande formule niet gebruikt worden. Elke waarde in een knoop krijgt een aangepast rangnummer als het gemiddelde rangnummer in de knoop. Spearmans rho is in dat geval de product-momentcorrelatiecoëfficiënt berekend voor de aangepaste rangnummers.

$\rho _{S}={\frac {1-{\frac {6}{n(n^{2}-1)}}\sum _{i}D_{i}^{2}-{\tfrac {1}{2}}(T_{X}+T_{Y})}{\sqrt {\left(1-T_{X}\right)\left(1-T_{Y}\right)}}}$

Daarin is, voor de X-waarden, en ook voor de Y-waarden:

$T={\tfrac {1}{n(n^{2}-1)}}\sum _{k}t_{k}(t_{k}^{2}-1)$,

met $t_{k}$ het aantal waarnemingen in de betrokken steekproef met hetzelfde rangnummer.

Voorbeeld

waarden X	waarden Y	rangnr. X	rangnr. Y	verschil D	D2	$t_{X}$	$\scriptstyle t_{X}^{\,}(t_{X}^{2}-1)$	$t_{Y}$	$\scriptstyle t_{Y}^{\,}(t_{Y}^{2}-1)$
2,0	1,5	1	2½	−1½	2,25	1	0	2	6
3,0	1,5	2½	2½	0	0	2	6
3,0	4,0	2½	5	−2½	6,25			1	0
5,0	3,0	4	4	0	0	1	0	1	0
5,5	1,0	5	1	4	16	1	0	1	0
8,0	5,0	6	6½	-½	0,25	1	0	2	6
10,0	5,0	7½	6½	1	1	2	6
10,0	9,5	7½	8	-½	0,25			1	0
totalen					26		12		12

Het aantal waarnemingsparen is $n=8$, dus

$n(n^{2}-1)=504$

en

$T_{X}=T_{Y}=12/504=0{,}0238$.

Voor Spearmans rho vinden we:

$\rho _{S}={\frac {1-6\times 26/504-0{,}0238}{0{,}976}}=0{,}68$