phyper(inter-1, a, 20000-a, b, lower.tail = F)
> phyper(900-1, 2800, 20000-2800, 3300, lower.tail = F)
[1] 4.116406e-110
- 这里的a为A数据集的基因数,b为B数据集的基因数,inter为两者交集的基因数。
- 这里需要解释的是代码phyper(inter-1, a, 20000-a, b, lower.tail = F)中的20000代表的是背景基因总数,如果是mRNA我这里设置为20000。
网友评论