本文介绍了如何利用R包`qqman`绘制GWAS（全基因组关联研究）的曼哈顿图和QQ图，以及如何计算膨胀系数（Genomic Inflation Factor）。通过示例代码展示了画图过程，并解释了膨胀系数的意义和其在评估假阳性结果中的作用。
画曼哈顿图和QQ plot 首推R包“qqman”，简约方便。下面具体介绍以下。
一、画曼哈顿图
install.packages("qqman")
library(qqman)
1、准备包含SNP, CHR, BP, P的文件gwasResults（如果没有zscore可以不用管），如下所示：
2、上代码，如下所示：
manhattan(gwasResults)
如果觉得不够美观，考虑添加一下参数：
manhattan(gwasResults, main = "Manhattan Plot", ylim = c(0, 10), cex = 0.6, cex.axis =
最低0.47元/天 解锁文章
1497
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
