我正在使用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 ...)。此列应该是数字还是字符?