##计算
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())
Structure 计算与绘图
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- pophelper(R 包) Usage:Rscript ../../script/draw_admixture....
- 近年来,为了加快药物发现的进程,人工智能和大数据技术不断地被应用到药物研发的各个领域。其中,在化合物结构解...
- 提到科学计算,可能大家马上就想到matlab、maple等,或者scilab、octave、R、gnuplot等...
- 本节来介绍如何R来自定义构建函数来进行多样性指数计算及绘图 加载R包 导入数据 定义函数计算多样性指数 导出指数计...