CMplot绘制曼哈顿图

lfmm=read.table("lfmm.txt",  sep = "\t", header=T)
head(lfmm)
CMplot(lfmm, plot.type="m", multracks=F, threshold=c(18.4342),threshold.lty=c(2),
       
       threshold.lwd=c(1), threshold.col=c("black"), amplify=TRUE,LOG10=F,  ylim = c(0, 50), chr.labels=NULL,
       
       signal.col=c("red","green"),signal.cex=c(1,1), file="pdf",dpi=300,file.output=TRUE,verbose=TRUE)
fst=read.table("FST.txt",  sep = "\t", header=T)
head(fst)
CMplot(fst, plot.type="m", multracks=F, threshold=c(0.428999),threshold.lty=c(2),
       
       threshold.lwd=c(1), threshold.col=c("black"), amplify=TRUE,LOG10=F,  ylim = c(0, 1), chr.labels=NULL,
       
       signal.col=c("red","green"),signal.cex=c(1,1), file="pdf",dpi=300,file.output=TRUE,verbose=TRUE)
XPEHH=read.table("XPEHH.txt",  sep = "\t", header=T)
head(XPEHH)
CMplot(XPEHH, plot.type="m", multracks=F, threshold=c(0.276622),threshold.lty=c(2),
       
       threshold.lwd=c(1), threshold.col=c("black"), amplify=TRUE,LOG10=F,  ylim = c(0, 0.5), chr.labels=NULL,
       
       signal.col=c("red","green"),signal.cex=c(1,1), file="pdf",dpi=300,file.output=TRUE,verbose=TRUE)
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容