Spearman 等級相關性的檢驗統計量

Question

我的問題是我已經計算了兩個變量之間的 Spearman 等級相關性。 一位評論者問我是否可以為 p < 0.001 的所有系數添加檢驗統計量。

這是一個結果：

> cor.test(pl$data, pl$rang, method= "spearman")


#       Spearman's rank correlation rho

data:  pl$data and pl$rang
S = 911164.6, p-value = 1.513e-05
alternative hypothesis: true rho is not equal to 0 
sample estimates:
rho 
-0.3347658 

測試統計是否等於S = 911164.6 ？ 這么大的數字可以嗎？ 如果問題不是很專業，請提前抱歉，但我花了很多時間在書籍和互聯網上尋找答案。 :( 提前感謝您的幫助。

Answer 1

是的。 ?cor.test幫助頁面（在 Value 部分）將cor.test的返回值描述為：

 A list with class '"htest"' containing the following components:

統計量：檢驗統計量的值。

修改該頁面上的示例，我們看到

x <- c(44.4, 45.9, 41.9, 53.3, 44.7, 44.1, 50.7, 45.2, 60.1)
y <- c( 2.6,  3.1,  2.5,  5.0,  3.6,  4.0,  5.2,  2.8,  3.8)
(result <- cor.test(x, y, method = "spearman"))
#         Spearman's rank correlation rho

# data:  x and y 
# S = 48, p-value = 0.0968
# alternative hypothesis: true rho is not equal to 0 
# sample estimates:
# rho 
# 0.6 
result$statistic
#  S 
# 48 

統計量由(n ^ 3 - n) * (1 - r) / 6 ，其中n是x和r <- cor(rank(x), rank(y))的長度。

Answer 2

我找到的最佳答案在頁面上： Rpubs Spearman Rank Correlation

這個問題也在Cross Validated “Interpreting the Spearman's Rank Correlation Coefficient output in R.What is 'S'?”

雖然論壇(n ^ 3 - n) * (1 - r) / 6 ，如果n等於變量的樣本大小並且r等於用於計算 S 統計量的相關系數 (rho)，則很清楚，有沒有關於如何解釋結果的明確解釋。 我還沒有找到明確的答案:(