這兩天在看一篇論文，文中用到了Fisher's exact p值的計(jì)算。作者給的原始數(shù)據(jù)是（14，8606，21，10197659），計(jì)算出的p值是2.12E-34。由于我學(xué)R沒(méi)有多長(zhǎng)時(shí)間，便試著在R 2.15.2中再次驗(yàn)算，居然算出結(jié)果是2.2E-16。我覺(jué)得這兩個(gè)結(jié)果之間差距太大了，于是請(qǐng)教了一個(gè)師兄，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