Structure 计算与绘图

##计算
library(LEA)
output = ped2lfmm("51.ped")
obj.snmf = snmf("51.lfmm", K = 1:5, alpha = 100, project = "new")
plot(obj.snmf, col = "blue4", cex = 1.4, pch = 19)
##绘图
sfiles <- list.files("tu3",full.names=TRUE)
slist <- readQ(sfiles)
labset <- read.delim("51lab.txt", header=F,stringsAsFactors=F)
labset
plotQ(slist,imgoutput = "join", 
      outputfilename="plotq9",
      imgtype="PDF", 
      panelspacer =0.05, 
      height = 0.5, 
      splabangle = 0, 
      showindlab=F,
      grplab=labset[,1,drop=FALSE],
      grplabangle = 90,
      grplabsize = 1,
      grplabheight = 5,
      exportpath=getwd())