レッスン
理論 — 勝者の呪い
統計学における勝者の呪いとは、統計的に有意であるという理由で選別された推定値に生じる偏り(バイアス)のことです。真の効果量 δ に対する研究の推定値 d̂ は、z = 1.96 として |d̂| > z·SE のとき α = 0.05 で有意となりますが、有意であるという条件のもとでの期待値 E(|d̂| | 有意) は真の δ を上回ります。Andrew Gelman と John Carlin は、これを誇大比 E(|d̂| | 有意)/δ(M型過誤)と、有意な推定値の符号が逆になる確率(S型過誤)によって定義しました。λ = δ/SE と置いた正規分布の推定値では、検出力は Φ(λ − z) + Φ(−λ − z) で与えられ、S型過誤率は Φ(−λ − z) を検出力で割った比率となります。
各記号の意味
δ- 結果変数の標準偏差で表した真の効果量
d̂- 一つの研究から得られた δ の推定値
n- 二つの群それぞれの人数
SE- 推定値の標準誤差 √(2/n)
λ- 標準誤差単位で測った真の効果量 δ/SE
z- 標準誤差単位の有意基準(両側 p < 0.05 では 1.96)
- 前提
- 各群 n 人の二群において、結果変数が既知の標準偏差 1 を持つ正規分布に従うと仮定します。したがって推定値は厳密に正規分布 N(δ, 2/n) に従い、検定は有意水準 α = 0.05 の両側z検定となります。真の効果量は正の値であらゆる研究で同一、各研究は独立であり、有意となった研究のみが報告・採択されるものとします。
全プロセスの詳細解説
-
0.45を報告した研究 4 ステップ
各群50人の研究が、0.45標準偏差の効果量(p = 0.024)を報告しました。しかし類似の介入研究から、この結果変数の真の変化はおよそ0.10であることが知られています。もし真の効果量が0.10であるならば、この研究計画は効果を何倍過大評価せざるを得なかったでしょうか。また符号を誤る確率はどれほどで、80%の検出力を得るには各群何人が必要になるでしょうか。0.45を報告した研究を押すと、ページ上で δ = 0.10、n = 50 に設定されます。
-
まず報告された数値から、標準誤差と z 統計量を算出します。
-
研究外の知見を取り入れ、真の効果量を δ = 0.10 と仮定します。標準誤差単位では λ = 0.5 となり、検出力は次の通りです。
-
全研究のうち有意となるのは7.91%で、0.69%は逆の符号で有意になります。これは有意となった研究全体の8.8%に当たります。有意な推定値の大きさは、平均して0.479となります。
-
80%の検出力を達成するには、反対側の裾を無視できるため λ が z + Φ⁻¹(0.80) に達する必要があり、λ = δ√(n/2) より次式が成り立ちます。
解答
各群1,570人、実に31.4倍。各群50人の研究計画では検出力は7.9%にとどまり、有意になるには推定値が少なくとも0.392以上でなければなりませんでした。そして有意となった推定値の平均は0.479と、真の値の4.79倍に膨張し、その8.8%で符号まで逆転します。論文に報告された0.45という値は、真の効果量が0.10であるときにこの小規模な研究計画がごく自然に生み出してしまう数値であり、0.10という真実を否定する根拠には全くなりません。ページ上で1570と入力してみてください。検出力は80.0%となり、有意な推定値の平均は0.112、真の値の1.12倍へと収束します。
-
学習の道すじ
p < 0.05 が保証してくれないこと
参考文献 (1)
- Type S and Type M errors and the exaggeration ratio, the formulas for power and the Type S rate under a normal estimate (note 2), the R function retrodesign() whose figures check:compute reproduces, the beauty-and-sex-ratio example the third insight quotes, and the circle the second insight describes (a significant estimate from a small study fed back into a power analysis makes the study look reasonable): A. Gelman and J. Carlin, "Beyond Power Calculations: Assessing Type S (Sign) and Type M (Magnitude) Errors," Perspectives on Psychological Science 9(6), 641–651, 2014.