這兩天在看一篇論文，文中用到了Fisher's exact p值的計算。作者給的原始數(shù)據(jù)是（14，8606，21，10197659），計算出的p值是2.12E-34。由于我學(xué)R沒有多長時間，便試著在R 2.15.2中再次驗算，居然算出結(jié)果是2.2E-16。我覺得這兩個結(jié)果之間差距太大了，于是請教了一個師兄，nice的幫我解答了一下，分享給大家參考。

> 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