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

原创  小果  生信果  关注我们

对于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语言edgeR差异表达分析

• 这个调控因子到底调不调控我的基因啊！——小果带你来学习FIMO数据库

• 小果教你解读WGCNA分析中的cluster dendrogram图谱
  

关注小果，小果将会持续为你带来更多生信干货哦。

“生信果”，生信入门、R语言、生信图解读与绘制、软件操作、代码复现、生信硬核知识技能、服务器、生物信息学的教程，以及基于R的分析和可视化等原创内容，一起见证小白和大佬的成长。