接上文,朋友告诉我,并不是只求一个 p 值,而是要计算许多 p 值,即批量进行 fisher 检验。
我们先看看数据:
library(knitr)
df = read.table('fisher_test.tsv', header = TRUE)
kable(df)
Gene | mMut | mWt | nonMut | nonWt |
---|---|---|---|---|
TP53 | 14 | 24 | 13 | 20 |
IFITM3 | 7 | 31 | 4 | 29 |
MTRNR2L2 | 5 | 33 | 4 | 29 |
USP17L10 | 5 | 33 | 4 | 29 |
CDH1 | 8 | 30 | 6 | 27 |
CLCNKA | 6 | 32 | 6 | 27 |
USP17L18 | 3 |