生統(tǒng)筆記5-多重檢驗(yàn)矯正

1. Bonferroni 矯正:

原理:將顯著性水平(α)除以檢驗(yàn)的數(shù)量 m阳距。每個(gè)檢驗(yàn)的閾值為α/m酬蹋。這種方法簡單糖声,但在檢驗(yàn)數(shù)量較多時(shí)可能過于保守。

2. Holm-Bonferroni 方法:

原理:對(duì) p 值進(jìn)行排序澈歉,并逐個(gè)檢驗(yàn)展鸡。第 i 個(gè) p 值的顯著性水平為 α/(m?i+1)。與 Bonferroni 方法相比埃难,它更加靈活莹弊,減少了對(duì)顯著性水平的過度保守。

3. False Discovery Rate (FDR) 控制方法

  • FDR(錯(cuò)誤發(fā)現(xiàn)率涡尘,F(xiàn)alse Discovery Rate)指的是在多個(gè)假設(shè)檢驗(yàn)中忍弛,被錯(cuò)誤拒絕的零假設(shè)所占的比例。具體來說考抄,它衡量的是在所有被認(rèn)為顯著的結(jié)果中细疚,實(shí)際上是錯(cuò)誤的結(jié)果所占的比例。

  • 舉例:比如有1000個(gè)統(tǒng)計(jì)檢驗(yàn)川梅,其中800個(gè)檢驗(yàn)是顯著的疯兼,但實(shí)際上真正顯著的只有700個(gè),有100個(gè)是錯(cuò)誤的發(fā)現(xiàn)贫途,則FDR = 1/8

一般用Benjamini-Hochberg (BH) 方法:

  1. 進(jìn)行m個(gè)假設(shè)檢驗(yàn)吧彪,得到m個(gè) p 值

  2. 將 p 值按從小到大排列

  3. 設(shè)定期望的假發(fā)現(xiàn)率(FDR)水平 α,通常是 0.05

  4. 對(duì)于排序后的p值丢早,計(jì)算其對(duì)應(yīng)的閾值姨裸。例如對(duì)于第i個(gè)p值,閾值為 i/m×α

  5. 對(duì)排序后的每個(gè)p值及其對(duì)應(yīng)閾值進(jìn)行比較怨酝,小于閾值的認(rèn)為是顯著

FDR與p值的相互轉(zhuǎn)換

為了在同一閾值下控制FDR傀缩,可以對(duì)每個(gè)p值計(jì)算其對(duì)應(yīng)q值:qi = pi * m/i,再將q值與α進(jìn)行比較即可

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末农猬,一起剝皮案震驚了整個(gè)濱河市赡艰,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌斤葱,老刑警劉巖瞄摊,帶你破解...
    沈念sama閱讀 206,013評(píng)論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件勋又,死亡現(xiàn)場離奇詭異,居然都是意外死亡换帜,警方通過查閱死者的電腦和手機(jī)楔壤,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,205評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來惯驼,“玉大人蹲嚣,你說我怎么就攤上這事∷钌” “怎么了隙畜?”我有些...
    開封第一講書人閱讀 152,370評(píng)論 0 342
  • 文/不壞的土叔 我叫張陵,是天一觀的道長说贝。 經(jīng)常有香客問我议惰,道長,這世上最難降的妖魔是什么乡恕? 我笑而不...
    開封第一講書人閱讀 55,168評(píng)論 1 278
  • 正文 為了忘掉前任言询,我火速辦了婚禮,結(jié)果婚禮上傲宜,老公的妹妹穿的比我還像新娘运杭。我一直安慰自己,他們只是感情好函卒,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,153評(píng)論 5 371
  • 文/花漫 我一把揭開白布辆憔。 她就那樣靜靜地躺著,像睡著了一般报嵌。 火紅的嫁衣襯著肌膚如雪虱咧。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 48,954評(píng)論 1 283
  • 那天锚国,我揣著相機(jī)與錄音腕巡,去河邊找鬼。 笑死跷叉,一個(gè)胖子當(dāng)著我的面吹牛逸雹,可吹牛的內(nèi)容都是我干的营搅。 我是一名探鬼主播云挟,決...
    沈念sama閱讀 38,271評(píng)論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢(mèng)啊……” “哼转质!你這毒婦竟也來了园欣?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 36,916評(píng)論 0 259
  • 序言:老撾萬榮一對(duì)情侶失蹤休蟹,失蹤者是張志新(化名)和其女友劉穎沸枯,沒想到半個(gè)月后日矫,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,382評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡绑榴,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 35,877評(píng)論 2 323
  • 正文 我和宋清朗相戀三年哪轿,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片翔怎。...
    茶點(diǎn)故事閱讀 37,989評(píng)論 1 333
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡窃诉,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出赤套,到底是詐尸還是另有隱情飘痛,我是刑警寧澤,帶...
    沈念sama閱讀 33,624評(píng)論 4 322
  • 正文 年R本政府宣布容握,位于F島的核電站宣脉,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏剔氏。R本人自食惡果不足惜塑猖,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,209評(píng)論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望介蛉。 院中可真熱鬧萌庆,春花似錦、人聲如沸币旧。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,199評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽吹菱。三九已至巍虫,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間鳍刷,已是汗流浹背占遥。 一陣腳步聲響...
    開封第一講書人閱讀 31,418評(píng)論 1 260
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留输瓜,地道東北人瓦胎。 一個(gè)月前我還...
    沈念sama閱讀 45,401評(píng)論 2 352
  • 正文 我出身青樓,卻偏偏與公主長得像尤揣,于是被迫代替她去往敵國和親搔啊。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,700評(píng)論 2 345

推薦閱讀更多精彩內(nèi)容