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)
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语言入门必备——立刻拥有一个Rstudio账号开启升级模式
-
使用R进行数据提取
-
看看差距有多大,R语言edgeR差异表达分析
-
这个调控因子到底调不调控我的基因啊!——小果带你来学习FIMO数据库
-
小果教你解读WGCNA分析中的cluster dendrogram图谱