library(pophelper)
#匹配后缀名是Q的admixture分析结果
files <- paste0("QC.",2:10,".Q")
#读取样品ID信息
ID <- read.table("all.txt",header = F)
slist <- readQ(files = files)
for(i in 1:length(slist)){
row.names(slist[[i]]) <- ID[,1]
sppos = "left"
splab = paste0("K=", 2:10)
splabangle = 180
dpi = 600
outputfilename = paste0("admixture_barplot.png")
units = "cm"
exportpath = getwd()
}
#绘制图形
plotQ(slist, imgoutput = "join",imgtype = "pdf",width = 20,height = 1.2,
showindlab = T,useindlab = T,indlabhjust = 0.5,indlabangle = 90,
indlabsize = 0.6,sortind = "all",sharedindlab = F,
panelspacer = 0.02, indlabspacer = 0,indlabcol = "black",
showtitle = T,titlelab = "Population Structure",titlecol = "black",
titlehjust = 0.5,showsp = T,exportpath=getwd())
按群体划分
plotQ(slist, imgoutput = "join",imgtype = "pdf",grplab=ID[2],width = 20,height = 1.2,
showindlab = F,useindlab = F,indlabhjust = 0.5,indlabangle = 90,
indlabsize = 0.6,sortind = "all",sharedindlab = F,showdiv = F,
panelspacer = 0.02, indlabspacer = 0,indlabcol = "black",
showtitle = T,titlelab = "Population Structure",titlecol = "black",
titlehjust = 0.5,showsp = F,exportpath=getwd())
