作图:
这个代码是qqman画的,具体教程参考:GWAS分析中可视化:qqman和cmplot以及颜值即正义 | 只知道qqman而不知道cmplot是不专业的
但是如何设置上面的红线和蓝线的值呢?下面用代码演示一下:
1,示例代码
library(qqman)
manhattan(gwasResults)
2,设置红线和蓝线
蓝色线设置为10的-2次方,红线设置为10的-4次方:
manhattan(gwasResults,suggestiveline = -log10(1e-2),genomewideline = -log10(1e-4))
3,去掉红线和蓝线
manhattan(gwasResults,suggestiveline = FALSE,genomewideline = FALSE)
其实,这个函数有帮助文档:
R中键入:
?manhattan
如果想要在曼哈顿图中显示snp的信息,可以看这篇博客:GWAS中曼哈顿图如何显示snp的信息,以及多性状曼哈顿图绘制的方法:多性状或者多个模型的QQ和曼哈顿重叠图。