24小时热门版块排行榜    

查看: 6346  |  回复: 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-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 18:22 by 7huP8NPKZCp9
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 6/300 2026-08-15 18:18 by 7huP8NPKZCp9
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] filecode=后面第一个是大写字母 +3 wangze12014 2026-08-14 4/200 2026-08-15 18:05 by babangida
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 7lpolszZVXgi 2026-08-14 7/350 2026-08-15 17:58 by 7huP8NPKZCp9
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 17:58 by 7huP8NPKZCp9
[基金申请] 时间戳又变了8-15 +11 archvillain 2026-08-15 20/1000 2026-08-15 17:47 by 医学老男孩
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 HFw0lei2R37i 2026-08-14 6/300 2026-08-15 17:34 by 7huP8NPKZCp9
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 P4JhWrbJGjSe 2026-08-14 5/250 2026-08-15 15:46 by QXtWNz7PI7wZ
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 7lpolszZVXgi 2026-08-14 3/150 2026-08-15 14:53 by QXtWNz7PI7wZ
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 14:41 by QXtWNz7PI7wZ
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见