R语言pophelper包对 Admixture遗传结构分析可视化

尔云间 一个专门做科研的团队

生信果

对于Admixture软件生成的遗传结构分析结果进行可视化如果只使用boxplot函数，会比较单调如下图

如果想要绘图更赏心悦目一些可以使用R语言pophelper包，效果如下

如果上图成功引起了你的注意，那接下来就和小果一起敲代码吧

代码如下

rm(list=ls());gc();
BiocManager::install('pophelper')
#安装依赖包
install.packages(c("devtools","ggplot2","gridExtra","gtable","label.switching","tidyr"),dependencies=T)


# install pophelper package from GitHub
devtools::install_github('royfrancis/pophelper')


library(pophelper)
library(data.table)
options(stringsAsFactors = F)
alist <- readQ(filetype="auto",list.files(path="D:/test/", full.names=T))




alist<-alignK(alist)
plotQ(alist,exportpath=getwd(),sortind="all" , height = 5,
      width = 5,showyaxis = TRUE,imgtype="pdf",showlegend = TRUE,showdiv = FALSE,outputfilename="plotq7")


plotQMultiline(alist,sortind="all",exportpath=getwd())           

完工[得意]~

推荐阅读

使用R语言完成序列比对及进化树美化

关于R包安装的一些小Tips

你不知道的PCA及在R中的实现

小果教你三分钟看懂多条线共存的ROC图的R语言画法

和小果学习如何一键转换gtf和gff3文件