偽の発見

マス目のひとつひとつが、この分野で検証される仮説です。黒点がついているものが真の仮説を表します。真の仮説の割合や研究の検出力、紛れ込むバイアスを設定したら、研究現場の歩みを順に進めてみてください。研究を行い、バイアスが入り込み、有意な結果が出版され、出版された研究がすべて追試されます。

インタラクティブシミュレーションを読み込んでいます...

レッスン

理論 — 偽の発見

ある研究分野における知見の 陽性的中率(PPV)とは、統計的に有意とされた結果のうち真実であるものの割合を指します。検証される仮説のうち真であるものの割合を π、有意水準 α における各研究の検出力を P、有意差の出なかった分析が発見として報告される割合を u とすると、真の仮説が発見となる確率は P + u(1 − P)、偽の仮説が発見となる確率は α + u(1 − α) となります。したがって PPV = π[P + u(1 − P)] / (π[P + u(1 − P)] + (1 − π)[α + u(1 − α)]) が成り立ちます。これはジョン・イオアニディスが2005年に発表した公式であり、原論文のオッズ R = π/(1 − π) の代わりに π を用いて表したものです。同じ検出力で誠実に行われた追試は、全発見のうち PPV·P + (1 − PPV)·α を再現します。

各記号の意味

π
その分野で検証される仮説のうち、真であるものの割合
P
各研究の検出力(1 − β):真の効果を検出する確率
α
有意水準(0.05):偽の仮説を検証した研究が p < 0.05 に達する確率
u
バイアス:有意差が出なかったにもかかわらず発見として報告される分析の割合
PPV
その分野の発見のうち、真実であるものの割合
前提
同一分野のすべての研究が同じ検出力を持つこと。各仮説は出版前に一度だけ検証されること。バイアスは有意差の出なかった結果にのみ作用すること。すべての発見が出版され、かつ発見のみが出版されること。追試は同一の検出力を持ち、バイアスのない同一の研究であること。なお、ここでの発見には効果量や効果の向きは考慮されていません。有意な研究が効果をどれほど過大評価するかという問題は、ひとつ前の単元である『勝者の呪い』で扱いました。

同じ p < 0.05 でも、検出力の低い研究では価値が薄れる 🖖

検出力80% を押してください。仮説の真偽は10に1つが真で、分析の歪曲もないため、発見とは純粋に p < 0.05 を満たしたものを指します。仮説1,000件あたり、分野からは80件の真の発見と45件の偽の発見が報告され、出版されたものの64.0%が本物となります。次に 検出力20% を押してみましょう。偽の発見は45件のまま動きません。偽である900件の仮説から α = 5% で生じるものであり、検出力の影響を一切受けないからです。今回の試行でも、黒点のない同じ48個のマス目が、どちらの設定でもまったく同じ位置に残っています。一方で真の発見は20件へと激減し、本物の割合も30.8%まで落ち込みます。検出力の不足は、よく『研究者が本物の効果を見落とす損失』として語られます。しかし損をするのは読者も同じです。脆弱な研究から出た有意な結果は、たとえ同じ基準をクリアしていても価値が低いのです。その基準を通過した結果のうち、本物が占める割合自体が小さいからです。ジョン・イオアニディスは2005年、この条件を1行で言い表しました。(1 − β)R > α のとき、発見は偽であるよりも真である公算が高くなります。ここで 1 − β は検出力、R は検証された仮説が真であるオッズです。10に1つが真なら R = 1/9 となり、検出力は45%を超えていなければなりません。

検出力はバイアスを帳消しにできない 🖖

バイアス10% を押してください。研究自体の検出力は先ほどと同じ80%ですが、有意差の出なかった分析の10に1つが、何らかの形で発見として報告されてしまいます。別の評価項目へのすり替え、都合の良いサブグループの抽出、不都合なデータの除外などがそれです。イオアニディスはこの割合を u と呼びました。このバイアスは、もともと大半が検出されていた真の仮説にはほとんど恩恵をもたらしません(80%が82%になる程度です)。深刻なのは偽の仮説への影響です。誠実な検証であれば95%の確率で退けられるはずのところ、その95%の10%にあたる9.5ポイントが漏れ出します。偽の仮説が発見と化す確率は14.5%へと跳ね上がり、実に3倍近い頻度です。本物の発見の割合は64.0%から38.6%へと急落します。下の結果欄を見れば、その深刻さがわかるでしょう。バイアスのない誠実な分野なら、検出力28.3%の研究でも同じ38.6%を達成できるのです。スライダーで検出力を99%まで引き上げても、バイアスのある分野のPPVは43.2%までしか届きません。100%の行を見ても43.4%であり、どんな研究でもこれ以上には達しないのです。ここで 誠実・検出力35% を押してみてください。43.8%となります。検出力35%の誠実な分野は、u = 10% のバイアスを抱えた完全な検出力の分野と同等以上に信頼できる発見を世に送り出せるのです。その境界点である α/(α + u(1 − α)) = 34.5% の式には π が含まれていません。最初のスライダーをどう動かしても、この値はびくともしないのです。

本物の発見であっても追試には失敗する。それも決して稀ではない 🖖

検出力80% を押し、最初から を押して最後のステージまで順に進めてみてください。出版されたすべての発見が、同じ検出力で誠実に追試されます。真の発見が再び有意となる確率は80%、偽の発見は5%ですから、出版された知見のうち再現されるのは53.0%です(今回の試行では129件中69件)。黒点の上に重なった×印を数えてみてください。再現に失敗した60件のうち12件は真の仮説であり、理論的な割合は27.2%にのぼります。誰も不正を働いておらず、研究の手続きに不備があったわけでもありません。検出力80%である以上、本物の効果であっても5回に1回は見逃される運命にあり、それは1回目でも2回目の追試でもまったく同じだからです。したがって、追試の再現率をそのまま『その分野における真実の割合』と読み取ることはできません。真の発見であっても追試の検出力の確率でしか再現されず、偽の発見はそれより低い α の確率でしか再現されないため、再現率は追試の検出力を決して超えられないのです。仮にすべての発見が真実であったとしても、追試の検出力が50%であれば半数の発見が再現に失敗します。バイアス10% を押すと、同じ誠実な追試を行っても再現率は33.9%にまで落ちます。再現を試みた対象の多くが、もともと存在すらしていなかったからです。

全プロセスの詳細解説

  1. p = 0.04 が次の研究に約束するもの 4 ステップ

    ある研究が両側検定で p = 0.04 を報告したとします。真の効果が、その研究で測定された値と完全に一致していると仮定しましょう。同一の実験計画、同一のサンプルサイズで行う全く同じ追試が、同一の方向で p < 0.05 に達する確率はどれほどでしょうか。 得られた検出力を本ツールに入力し、真の仮説が10に1つ、バイアスなしとして、そうした分野の発見のうち本物が占める割合と追試の再現率を読み取ってみてください。

    この問題をツールで開く

    1. p値を、それが算出された元のz値へと逆算します。両側検定なので左右の裾はそれぞれ0.02であり、p < 0.05 の境界となるのは上側確率が0.025となる点です:

    2. 真の効果が観測値と等しいなら、追試におけるz値は平均 2.054、標準偏差 1 の正規分布に従います。追試が成功するのは、この値が 1.960 を上回るときです:

    3. −1.960 を下回った場合は逆向きの有意差となり、再現とは言えません。いずれにしても、その確率は無視できるほどわずかです:

    4. p = 0.04 の研究分野 を押してください。検出力54%(スライダーは1%刻み)、真の仮説は10に1つ、バイアスなしの設定です。

    解答

    53.7%。コイン投げよりわずかにましな程度であり、p = 0.05 であれば全く同じ計算手順でちょうど 50% になります。しかもこれは甘めの見積もりです。ひとつ前の単元で学んだ『勝者の呪い』が示す通り、有意差の出た研究は効果を過大評価しがちであるため、実際の追試における検出力はこれよりもさらに低くなります。検出力54%の研究が行われる分野では、真の発見は全体の54.5%にとどまり、追試で再現されるのは31.7%です。同一の仮定のもとで追試が5回に4回の確率で成功するためには、元の研究で p = 0.005 が得られている必要があり、その場合の検出力は80.2%となります。

学習の道すじ

p < 0.05 が保証してくれないこと

参考文献 (1)
  • The page's model and formula: the positive predictive value of a claimed finding from power, the pre-study odds R and the bias u (Tables 1 and 2), the condition that a finding is more likely true than false when (1 − β)R > α, which the first insight quotes, the definition of u the second insight paraphrases, and Table 4, whose first row is the 'Ioannidis's trial' preset and all nine of whose rows check:compute reproduces: J. P. A. Ioannidis, "Why Most Published Research Findings Are False," PLoS Medicine 2(8): e124, 2005.