复现文章中检验基因的overlap
使用R中的phyper函数实现:
-log10(phyper(58, 555, 19738-555,660, lower.tail = F))
#[1] 14.44417
扩展:
phyper(x, m, N-m, n, lower.tail = TRUE, log.p = FALSE)
文章
[1] Moncada R , Barkley D , Wagner F , et al. Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas[J]. Nature Biotechnology.