【点】p value、FDR、q value

假设检验:

一般步骤:

1.设定假设后需要验证:假设0一般都是“处理组和对照组无区别”

2.构建检验统计量:根据样本数据计算统计量和检验统计量

3.设定置信区间:根据检验统计量的分布,计算假设发生的概率P

# 分布的意义在于划定误差分布,即划定真实世界的随机误差范围

# 小概率事件的发生概率几乎为零,但是为了满足我们的现实需求,我们一般可以设定P>0.05即认定假设为真,即发生假设情况的概率为5%以上时,即认定假设为真。P<0.05时,即认定假设为假,即处理组和对照组有区别。

# 但在真实世界里,真实事件的概率可能会达到99.9%,而不只是95%。所以我们可能会将原本为真的情况,误判为假,即发生一类错误。

# 反过来说,当假设为假时,我们会认定处理组和对照组有区别,但其实有一定概率是“原本没有区别,但是被我们误判为假”的情况。

为了降低发生一类错误的概率,我们可以降低置信区间,相应的P值范围也会改变。多重检验校正便是通过设定一些规则改变P值范围来降低犯错的概率。


多重检验校正:

随着假设检验的次数增高,犯错误的概率就不断增大,这时候就需要进一步调整P,从而保证最终结果的可靠性。

比较基础的有Bonferroni校正和Benjamini and Hochberg校正

基本定义:

FWER(Family-Wise Error Rate):初现至少一次一类错误(原本为真,判定为假)的概率

FDR(False Discovery Rate):所有判为假的结论中(阳性结果),即发生一类错误的概率

Bonferroni校正:

基本原理是控制FWER

adj.P= α/m

# 其中α 为原定显著性,m 为检验次数;

# 举例:原α=0.05,检验m= 10000次,则adj.P= 0.000005。假设全部结果都是假(阳性结果),并且全部都是一类错误,则发生一类错误的次数 n= adj.P*m= 0.05< 1

Benjamini and Hochberg校正:

基本原理是对FDR设限,如1%,一百个阳性结果,有一个是一类错误。但是也可以根据α来设定FDR,即BH校正。

将所有检验P值进行由低到高排序,找到一个最大正整数k,使得

P(k)<= α* (k/m)

# m为总检验次数

# 1~k均为阳性结果

#该策略假设的是假阳性与真阳性虽然都是阳性,但假阳性与真阳性本质截然不同,故从假阳性到真阳性是斜率陡然上升的。

除了BH外,还可以直接设定FDR。


总结:

adj.P其实就是q value;

各种校正其实只是设定q value的标准不同,核心思想都是一样的;

以上两个校正策略,用得较多的是第二个。因为第一个太严格了,很容易将真阳性结果筛掉。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 205,132评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,802评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,566评论 0 338
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,858评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,867评论 5 368
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,695评论 1 282
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,064评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,705评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 42,915评论 1 300
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,677评论 2 323
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,796评论 1 333
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,432评论 4 322
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,041评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,992评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,223评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,185评论 2 352
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,535评论 2 343

推荐阅读更多精彩内容