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

原創(chuàng)? 小果? 生信果? 關(guā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)行數(shù)據(jù)提取

• 看看差距有多大，R語言edgeR差異表達(dá)分析

• 這個調(diào)控因子到底調(diào)不調(diào)控我的基因?。　」麕銇韺W(xué)習(xí)FIMO數(shù)據(jù)庫

• 小果教你解讀WGCNA分析中的cluster dendrogram圖譜
  

關(guān)注小果，小果將會持續(xù)為你帶來更多生信干貨哦。

“生信果”，生信入門、R語言、生信圖解讀與繪制、軟件操作、代碼復(fù)現(xiàn)、生信硬核知識技能、服務(wù)器、生物信息學(xué)的教程，以及基于R的分析和可視化等原創(chuàng)內(nèi)容，一起見證小白和大佬的成長。