24小时热门版块排行榜    

查看: 6427  |  回复: 3

逆刃刀客

金虫 (正式写手)

[求助] 关于费舍尔精确检验单侧还是双侧检验选择的请教 已有1人参与

我是做生物信息学的研究生,现在我希望通过费舍尔精确检验去判断两个基因的突变是否是同时发生的。例如:1表示基因突变,0表示基因未突变。假设有基因A=[0,1,1,0,0],基因B=[0,1,1,0,0].也就是A和B同时发生,我用MATLAB测试两个基因是否是同时突变的。
当做右侧检验,p=0.1,
当做左侧检验,p=1,
当做双侧检验,p=0.1。
如果两个基因为A=[0,1,1,0,0],B=[1,0,0,1,1]。也就是A和B不同时发生,使用MATLAB测试,则
当做右侧检验,p=1,
当做左侧检验,p=0.1,
当做双侧检验,p=0.1。
因此我看使用右侧检验更好一些,p值小,更显著,因此我使用了右侧检验。
但是,审稿人说我这样做是不正确的,请问各位,这是为什么?怎样判断使用双侧还是单侧?
回复此楼
为天地立心,为生民立命,为往圣继绝学,为万世开太平!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

退化or进化

木虫 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
逆刃刀客: 金币+3, ★有帮助 2015-05-07 09:06:03
我觉得应该是分别计算两个基因都突变、A突变、B突变、都不突变的个数,你的例子中为2,0,0,3和0,2,3,0,这样直接用Fisher exact双尾检验,都是0.1。
2楼2015-05-06 12:28:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

逆刃刀客

金虫 (正式写手)

引用回帖:
2楼: Originally posted by 退化or进化 at 2015-05-06 12:28:08
我觉得应该是分别计算两个基因都突变、A突变、B突变、都不突变的个数,你的例子中为2,0,0,3和0,2,3,0,这样直接用Fisher exact双尾检验,都是0.1。

谢谢你,还有些疑问,那我希望去看两个基因是否是同时突变的,使用单侧的检验是否正确?是否需要先验条件才能使用?如果需要,那条件是什么啊?
为天地立心,为生民立命,为往圣继绝学,为万世开太平!
3楼2015-05-07 09:05:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

退化or进化

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
逆刃刀客: 金币+12, ★★★★★最佳答案, 谢谢你,我明白了,豁然开朗,非常感谢。 2015-05-07 10:23:31
引用回帖:
3楼: Originally posted by 逆刃刀客 at 2015-05-07 09:05:47
谢谢你,还有些疑问,那我希望去看两个基因是否是同时突变的,使用单侧的检验是否正确?是否需要先验条件才能使用?如果需要,那条件是什么啊?...

这要看你希望test的hypothesis是什么,如果alternative hypothesis是两个基因的突变不独立,则应该用双尾检验;如果alternative hypothesis是两个基因的突变正相关(倾向同时突变),则应该用上单尾检验;同理,负相关(倾向互斥突变)应该用下单尾检验。
4楼2015-05-07 10:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 逆刃刀客 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 X5bHxOULxFqb 2026-10-30 27/1350 2026-11-02 05:28 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 4/200 2026-11-02 03:50 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 03:39 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 8/400 2026-11-02 03:24 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 X5bHxOULxFqb 2026-10-30 24/1200 2026-11-02 02:39 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:05 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 21:52 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 21:45 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
信息提示
请填处理意见