スピアマンの順位相関の信頼区間を計算する方法は?


13

ウィキペディアには、スピアマンの順位相関を近似Zスコアにフィッシャー変換しています。おそらく、Zスコアは帰無仮説(ランク相関0)との違いですか?

このページには次の例があります。

4, 10, 3, 1, 9, 2, 6, 7, 8, 5
5, 8, 6, 2, 10, 3, 9, 4, 7, 1
rank correlation 0.684848
"95% CI for rho (Fisher's z transformed)= 0.097085 to 0.918443"

Fisher変換を使用して95%信頼区間を取得する方法

回答:


20

一言で言えば、95%の信頼区間は、与えられます ここで、は相関の推定値で、はサンプルサイズです。

tanh(arctanhr±1.96/n3),
rn

説明:フィッシャー変換はアークタンです。変換されたスケールでは、推定値のサンプリング分布はほぼ正規であるため、変換された推定値を取得し、その標準誤差の1.96倍を加算および減算することで95%CIが見つかります。標準誤差は(およそ)です。1/n3

編集:上記のPythonの例:

import math
r = 0.684848
num = 10
stderr = 1.0 / math.sqrt(num - 3)
delta = 1.96 * stderr
lower = math.tanh(math.atanh(r) - delta)
upper = math.tanh(math.atanh(r) + delta)
print "lower %.6f upper %.6f" % (lower, upper)

与える

lower 0.097071 upper 0.918445

これは、小数点以下4桁までの例に一致します。


1つの質問:en.wikipedia.org/wiki / ...の1.06は、あなたの答えにどのように関連していますか?
dfrankow

あなたは私をそこに持っています!正直であることはわかりません。私はそれを試してみて、試してみましたが、試してみた結果と一致しました。
ワンストップ

1
@dfrankow私はその編集を受け入れましたが、これはこの機能の完全な使用ではありません-より良いアイデアは、質問にそのようなテキストを追加することです。

6
ζ^=タン1θ^θ^σζ^21.06/n3
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.