如何在曼哈顿情节中标记点数?

时间:2015-05-19 22:41:37

标签: r plot label

我正在使用qqman包,并对qqman包中的标签有一些疑问。我在R:

中运行了以下命令
manhattan(mydata, chr = "CHR", bp = "MAPINFO", p = "P.Value", snp = "SNP", col = c("blue4", "orange3"), chrlabs = NULL,suggestiveline = -log10(1e-05), genomewideline = -log10(5e-08), highlight =NULL, logp = TRUE )

得到了曼哈顿的情节,但是我想根据“mydata”文件中的“Gene”列,用–log10(p) > 20(在y轴上)标记所有点。该列(Gene)包括基因名称而不是SNP坐标(rs ...)。此列应该是数字还是字符?

0 个答案:

没有答案