斯皮爾曼等級相關係數

斯皮爾曼等級相關係數（簡稱等級相關係數，或稱秩相關係數，Spearman's rank correlation coefficient或Spearman's ρ），在統計學中，常以希臘字母$\rho$（rho）或以$r_s$表示，這一相關係數以查爾斯·斯皮爾曼之名命名。它是衡量兩個變量的相關性的無母數指標。它利用單調函數評價兩個統計變量的相關性。若數據中沒有重複值，且當兩變量完全單調相關時，斯皮爾曼相關係數為+1或−1。

## 定義和計算

斯皮爾曼相關係數的定義為等級變量之間的皮爾遜相關係數。

對於樣本容量為n的樣本，將n個原始數據$X_i, Y_i$轉換成等級數據$\operatorname{R}({X_i}), \operatorname{R}({Y_i})$，則相關係數$r_s$為

$$<pre> r_s = \rho_{\operatorname{R}(X),\operatorname{R}(Y)} = \frac{\operatorname{cov}(\operatorname{R}(X), \operatorname{R}(Y))} {\sigma_{\operatorname{R}(X)} \sigma_{\operatorname{R}(Y)}}, </pre>$$

其中

$\rho$是皮爾遜積矩相關係數，但使用等級變量來計算。

$\operatorname{cov}(\operatorname{R}(X), \operatorname{R}(Y))$為等級變量的協方差。

$\sigma_{\operatorname{R}(X)}$和$\sigma_{\operatorname{R}(Y)}$為等級變量的標準差。

通常，對於數據中相同的值，其等級數等於它們按值升序排列的所處位置的平均值。如下表所示：

變量$X_i$｜升序位置（僅示意，不使用）｜升序位置的平均等級數（使用）

18｜1｜1

2.3｜2｜2

1.2｜3｜$\frac{4+3}{2}=3.5\ $

1.2｜4｜$\frac{4+3}{2}=3.5\ $

0.8｜5｜5

當所有的等級數值都為整數時，可以通過以下簡單的步驟計算等級相關係數：

$$r_s = 1 - \frac{6 \sum d_i^2}{n(n^2 - 1)},$$

其中

$d_i = \operatorname{R}(X_i) - \operatorname{R}(Y_i)$為每組觀測中兩個變量的等級差值。

n為觀測數。

考慮一個雙變量樣本$(x_i, y_i),\, i=1\dots, n$，其相應的位次為$(R(X_i), R(Y_i)) = (R_i, S_i)$。則$x,y$的斯皮爾曼等級相關係數為：

$$r_{s} = <pre> \frac{ \frac{1}{n}\sum_{i=1}^{n}R_{i}S_{i} - \overline{R}\,\overline{S} }{ \sqrt{\sigma_{R}} \sqrt{\sigma_{S}} }, </pre>$$

其中：

$\overline{R} = \textstyle\frac{1}{n}\textstyle\sum_{i=1}^{n}R_{i}$。

$\overline{S} = \textstyle\frac{1}{n}\textstyle\sum_{i=1}^{n}S_{i}$。

$\sigma_{R}^{2} = \textstyle\frac{1}{n}\textstyle\sum_{i=1}^{n}(R_{i} - \overline{R})^{2}$。

$\sigma_{S}^{2} = \textstyle\frac{1}{n}\textstyle\sum_{i=1}^{n}(S_{i} - \overline{S})^{2}$。

若假定樣本中兩變量均沒有重複數值，則$r_{s}$可只用$d_{i} := R_{i} - S_{i}$來給出。

在此假定下，$R,S$可視為隨機變量，其分布類似於均勻分布隨機變量，$U$，其自變量取值為$\{1,2,\ldots,n\}$。

因此$$\overline{R}=\overline{S}=\mathbb{E}[U]$$ 且$\sigma_{R}^{2}=\sigma_{S}^{2}=\mathrm{Var}(U)=\mathbb{E}[U^{2}]-\mathbb{E}[U]^{2}$。

其中$\mathbb{E}[U] = \textstyle\frac{1}{n}\textstyle\sum_{i=1}^{n} i = \textstyle\frac{(n+1)}{2}$。

$\mathbb{E}[U^{2}] = \textstyle\frac{1}{n}\textstyle\sum_{i=1}^{n} i^{2} = \textstyle\frac{(n+1)(2n+1)}{6}$。

故有$\mathrm{Var}(U) = \textstyle\frac{(n+1)(2n+1)}{6} - \left(\textstyle\frac{(n+1)}{2}\right)^{2} = \textstyle\frac{n^{2}-1}{12}$。

（這些求和可以用三角形數和四角錐數的公式來計算，也可以用離散數學的基本求和結果來計算。）

既然

$$\begin{align} <pre> \frac{1}{n}\sum_{i=1}^{n}R_{i}S_{i} - \overline{R}\overline{S} &= \frac{1}{n}\sum_{i=1}^{n}\frac{1}{2}(R_{i}^{2} + S_{i}^{2} - d_{i}^{2}) - \overline{R}^{2}\\ &= \frac{1}{2}\frac{1}{n}\sum_{i=1}^{n}R_{i}^{2} + \frac{1}{2}\frac{1}{n}\sum_{i=1}^{n}S_{i}^{2} - \frac{1}{2n}\sum_{i=1}^{n}d_{i}^{2} - \overline{R}^{2}\\ &= (\frac{1}{n}\sum_{i=1}^{n}R_{i}^{2} - \overline{R}^{2}) - \frac{1}{2n}\sum_{i=1}^{n}d_{i}^{2}\\ &= \sigma_{R}^{2} - \frac{1}{2n}\sum_{i=1}^{n}d_{i}^{2}\\ &= \sigma_{R}\sigma_{S} - \frac{1}{2n}\sum_{i=1}^{n}d_{i}^{2}\\ </pre> \end{align}$$

則綜上可得

$$<pre> r_s = \frac{ \sigma_{R}\sigma_{S} - \frac{1}{2n}\sum_{i=1}^{n}d_{i}^{2} }{ \sigma_{R}\sigma_{S} } = 1 - \frac{ \sum_{i=1}^{n}d_{i}^{2} }{ 2n \cdot \frac{n^{2}-1}{12} } = 1 - \frac{ 6\sum_{i=1}^{n}d_{i}^{2} }{n(n^{2} - 1)}. </pre>$$

當數據中存在相等的數值時，使用該簡化公式會得到錯誤結果：只有在兩組變量中所有數值不重複時，才有$\sigma_{\operatorname{R}(X)} \sigma_{\operatorname{R}(Y)} = \operatorname{Var}{(\operatorname{R}(X))} = \operatorname{Var}{(\operatorname{R}(Y))} = (n^2 - 1)/12$（根據有偏方差計算）。第一個方程（通過標準差進行歸一化）即使在排名標準化為0, 1（「相對排名」）的情況下仍可使用，因為它對平移和線性縮放都不敏感。

對於截取的數據也不應使用簡化公式。即，當希望計算前X條記錄的等級相關係數時，應當使用前述的皮爾遜積矩相關係數公式。

## 相關度量

度量一對觀測數據的統計相關性還有其他的幾種度量指標。其中最常用的是皮爾遜積矩相關係數。

斯皮爾曼相關也可稱為「級別相關」（grade correlation）；也就是說，被觀測數據的「等級」被替換成 「級別」。在連續的分布中，被觀測數據的級別，通常總是小於等級的一半。然而，在這個案例中，級別和等級相關係數是一致的。更一般的，被觀測數據的「級別」與估計的總體樣本的比值小於給定的值，即被觀測值的一半。也就是說，它是相應的等級係數的一種可能的解決方案。雖然不常用，「級別相關」還是仍然有被使用。

## 解釋

正的斯皮爾曼相關係數反映兩個變量X和Y之間單調遞增的趨勢。

斯皮爾曼相關係數表明X（自變量）和Y（因變量）的相關方向。如果當X增加時，Y趨向於增加，則斯皮爾曼相關係數為正。如果當X增加時，Y趨向於減少，則斯皮爾曼相關係數為負。斯皮爾曼相關係數為0表明當X增加時Y沒有任何趨向性。當X和Y越來越接近完全的單調相關時，斯皮爾曼相關係數會在絕對值上增加。當X和Y完全單調相關時，斯皮爾曼相關係數的絕對值為1。完全的單調遞增關係意味著對任意兩對數據Xi, Yi和Xj, Yj，有Xi − Xj和Yi − Yj總是同號。完全的單調遞減關係意味著對任意兩對數據Xi, Yi和Xj, Yj，有Xi − Xj和Yi − Yj總是異號。

斯皮爾曼相關係數經常被稱作「非參數」的，其中有兩層含義。首先，當X和Y的關係由任意單調函數描述時，則它們是完全皮爾遜相關的。與此相應的，皮爾遜相關係數只能給出由線性方程描述的X和Y的相關性。其次，斯皮爾曼不需要先驗知識（也就是說，知道其參數）便可以準確獲取X和Y的採樣機率分布。

## 示例

在此例中，我們要使用下表所給出的原始數據計算一個人的智商和其每周看電視的小時數的相關性（數據為虛構）。

智商, $X_i$｜每周看電視小時數, $Y_i$

106｜7

86｜0

100｜27

101｜50

99｜28

103｜29

97｜20

113｜12

112｜6

110｜17

首先，我們必須根據以下步驟計算出$d^2_i$，如下表所示。

- 排列第一列數據（$X_i$）。創建新列 $x_i$ 並賦以等級值1、2、3......n。
- 然後，排列第二列數據（$Y_i$）。創建第四列 $y_i$ 並相似地賦以等級值1、2、3......n。
- 創建第五列$d_i$，填入兩個等級列（$x_i$和$y_i$）的差值。
- 創建最後一列$d^2_i$填入$d_i$的平方。

智商, $X_i$｜每周看電視小時數, $Y_i$｜$x_i$的排名｜$y_i$的排名｜$d_i$｜$d^2_i$

86｜0｜1｜1｜0｜0

97｜20｜2｜6｜−4｜16

99｜28｜3｜8｜−5｜25

100｜27｜4｜7｜−3｜9

101｜50｜5｜10｜−5｜25

103｜29｜6｜9｜−3｜9

106｜7｜7｜3｜4｜16

110｜17｜8｜5｜3｜9

112｜6｜9｜2｜7｜49

113｜12｜10｜4｜6｜36

<image>

根據$d^2_i$計算$\sum d_i^2 = 194$。樣本容量n為10。將這些值帶入方程

$$\rho = 1- {\frac {6\times194}{10(10^2 - 1)}}$$

得ρ = −0.175757575...，p-value = 0.627188（使用t分布）

該數值接近0，表明儘管看電視時間和智商似乎呈負相關，但兩個變量之間的關係很弱。在原始數據中存在相同數值的情況下，不應使用此公式，而應當用排名計算皮爾遜相關係數（如上文所述）。

## 顯著性的確定

一種確定被觀測數據的ρ值是否顯著不為零（r總是有1 ≥ r ≥ −1）的方法是計算它是否大於r的機率，作為零假設，並使用排列檢驗。這種方法的優勢在於它考慮了樣本中的重複出現的數據個數，以及在計算等級相關性時處理它們的方式。

另一種方法是使用皮爾遜積矩中使用到的費雪變換。也就是，ρ的置信區間和假說檢定可以通過費雪變換獲得

$$F(r) = {1 \over 2}\ln{1+r \over 1-r} = \operatorname{arctanh}(r).$$

如果F(r)是r的費雪變換，則

$$z = \sqrt{\frac{n-3}{1.06}}F(r)$$

是r的z-值，其中，r在統計獨立性（1=ρ = 0）的零假設下近似服從標準正態分布。

顯著性為

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

其在零假設下近似服從自由度為n − 2的t分布。A justification for this result relies on a permutation argument.

一般地，斯皮爾曼相關係數在有三個或更多條件的情況下是有用的。並且，它預測觀測數據有一個特定的順序。例如，在同一任務中，一系列的個體會被嘗試多次，並預測在多次嘗試過程中，性能會得到提升。在這種情況下，對條件間趨勢的顯著性檢驗由E. B. Page發展了，並通常稱為給定序列下的Page趨勢檢驗。

## 基於斯皮爾曼相關係數的一致性分析

經典的一致性分析是一種統計方法，它給兩個標稱變量賦給一個分數。通過這種方法，兩個變量間的皮爾遜相關係數被最大化了。

有一種被稱為級別相關分析的等價方法，它能夠最大化斯皮爾曼相關係數或肯德爾等級相關係數。

## 參考文獻

- G.W. Corder, D.I. Foreman, "Nonparametric Statistics for Non-Statisticians: A Step-by-Step Approach", Wiley (2009)
- C. Spearman, "The proof and measurement of association between two things" Amer. J. Psychol., 15 (1904) pp. 72–101
- M.G. Kendall, "Rank correlation methods", Griffin (1962)
- M. Hollander, D.A. Wolfe, "Nonparametric statistical methods", Wiley (1973)
- J. C. Caruso, N. Cliff, "Empirical Size, Coverage, and Power of Confidence Intervals for Spearman's Rho", Ed. and Psy. Meas., 57 (1997) pp. 637–654

## 外部連結

- "Understanding Correlation vs. Copulas in Excel" by Eric Torkia, Technology Partnerz 2011
- Table of critical values of ρ for significance with small samples
- A calculator that shows the working out for Spearman's correlation
- Spearman's rank online calculator
- Chapter 3 part 1 shows the formula to be used when there are ties
- Spearman's rank correlation : Simple notes for students with an example of usage by biologists and a spreadsheet for Microsoft Excel for calculating it (a part of materials for a Research Methods in Biology course).