確かに、@ cardinalは、論文がそれと同じくらい明確であることはまったく正しいです。そのため、この記事にアクセスできない場合に備えて、Benjamini–Hochbergが主張する方法を少し詳しく説明します。
FDR は、すべての拒否rに対する偽拒否vの割合の期待値です。ここで、rは明らかに、偽の拒否と正しい拒否の合計です。後者の呼び出し秒。Qevrrs
要約すると、(ランダム変数には大文字を、実現値には小文字を使用して)、
Qe=E(VR)=E(VV+S)=:E(Q).
R = 0の場合、取ります。Q=0R=0
現在、2つの可能性があります。すべての nullが真であるか、それらのm 0 < mのみが真です。前者の場合、正しい拒否はできないため、r = vです。したがって、任意の拒絶(存在する場合、R ≥ 1)、Q = 1、そうでなければQ = 0。したがって、mm0<mr=vr≥1q=1q=0
FDR=E(Q)=1⋅P(Q=1)+0⋅P(Q=0)=P(Q=1)=P(V≥1)=FWER
したがって、この場合はであり、F D Rを簡単に制御するプロシージャはF W E Rも制御し、その逆も同様です。FDR=FWERFDRFWER
ここで第二の場合には場合、V > 0(そう少なくとも一つの本人拒否がある場合)、我々は明らかに(これも有する画分である持っているVこと分母における)V / R ≤ 1。これは、少なくとも一つの本人拒否がある場合は値1をとるインジケータ機能、ということを意味1 V ≥ 1が未満になることはありませんQ、1 V ≥ 1 ≥ Q。ここで、Eの単調性により、不平等の両側に期待を寄せてください。m0<mv>0vv/r≤11V≥1Q1V≥1≥QE 不平等をそのまま残し、
E(1V≥1)≥E(Q)=FDR
インジケータにおける事象の確率であるインジケータ関数の期待値は、我々は再び、F W E Rを。E(1V≥1)=P(V≥1)FWER
我々が制御することを手順を有する場合したがって、の意味でそのF W E R ≤ α、我々はそれを持っている必要がありF D R ≤ α。FWERFWER≤αFDR≤α
逆に、あるαで制御すると、かなり大きなF W E Rが発生する場合があります。直観的には、テストされた仮説の潜在的に大きな合計からゼロ以外の予想される誤棄却率(F D R)を受け入れることは、少なくとも1つの誤棄却(F W E R)の非常に高い確率を意味します。FDRαFWERFDRFWER
したがって、制御のみが必要な場合、手順はそれほど厳しくはなく、電力にも適しています。これは、基本的な仮説検定と同じ考え方です。5%レベルでテストする場合は、1%レベルでテストする場合よりも頻繁に拒否します(正解と偽ヌルの両方)。FDR