爾云間? 一個專門做科研的團(tuán)隊

生信果

對于Admixture軟件生成的遺傳結(jié)構(gòu)分析結(jié)果進(jìn)行可視化如果只使用boxplot函數(shù)，會比較單調(diào)如下圖

如果想要繪圖更賞心悅目一些可以使用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語言完成序列比對及進(jìn)化樹美化

關(guān)于R包安裝的一些小Tips

你不知道的PCA及在R中的實(shí)現(xiàn)

小果教你三分鐘看懂多條線共存的ROC圖的R語言畫法

和小果學(xué)習(xí)如何一鍵轉(zhuǎn)換gtf和gff3文件