美文网首页
2019-01-06

2019-01-06

作者: 対MS特技兵 | 来源:发表于2019-01-06 15:32 被阅读0次

vcftools --vcf my.vcf --plink --out plink

plink--noweb--fileplink--geno0.1--maf0.05--hwe0.0001--make-bed--outout

forKin1 2 3 4 5;doadmixture --cv XXXXXX.bed $K| teelog${K}.out;done

grep -h CVlog*.out

取最小K值

tbl=read.table("hapmap3.3.Q")

barplot(t(as.matrix(tbl)), col=rainbow(3),xlab="Individual #", ylab="Ancestry", border=NA)

相关文章

网友评论

      本文标题:2019-01-06

      本文链接:https://www.haomeiwen.com/subject/zbesrqtx.html