全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5929 9
2013-03-07
这两天在看一篇论文,文中用到了Fisher's exact p值的计算。作者给的原始数据是(14,8606,21,10197659),计算出的p值是2.12E-34。由于我学R没有多长时间,便试着在R 2.15.2中再次验算,居然算出结果是2.2E-16。我觉得这两个结果之间差距太大了,于是请教了一个师兄,他也说不清楚为什么。还请各位大大给帮忙解释解释啊!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-3-7 17:23:14
怎么没人理我啊?拜托潜水的高手冒个泡给解释一下到底是怎么算的啊!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-3-9 14:35:36
希望大家帮我看看啊
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-3-9 16:05:27
楼主,你应该告诉别人你用的是什么代码?
这样才能更好的验证还有,看作者原文章上写的用哪个软件运行的?
> challenge.df = matrix(c(14,8608,21,10197659), nrow = 2)
> fisher.test(challenge.df)

        Fisher's Exact Test for Count Data

data:  challenge.df
p-value < 2.2e-16
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
  368.2732 1600.0000
sample estimates:
odds ratio
  798.1286
里面还有一些参数可以调整,你可以试试,我用2.14算出来也是。

http://stat.ethz.ch/R-manual/R-patched/library/stats/html/fisher.test.html



二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-3-10 11:51:05
这个只是在输出 summary时的显示问题, p-value < 2.2e-16 跟p值是2.12E-34 不矛盾哈。
直接看真实的p值就可以了。

> fisher.test(matrix(c(14,8608,21,10197659), nrow = 2))

        Fisher's Exact Test for Count Data

data:  matrix(c(14, 8608, 21, 10197659), nrow = 2)
p-value < 2.2e-16
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
  368.2732 1600.0000
sample estimates:
odds ratio
  798.1286

> tmp <- .Last.value
> tmp$p.value
[1] 2.128486e-34
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-3-10 19:26:58
谢谢sein_tao的答复。再请教:你说的那个“真实的p值”是什么意思啊?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群