大家好，又见面了，我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全家桶1年46，售后保障稳定
群体遗传学中测的很多个个体，得到了最终的SNP vcf文件，需要将其分成群体，看那几个物种聚在一起，一般使用的软件就是STRUCTURE，但是STREUTURE运行速度极慢，后面frappe软件提升了速度，但是也不是很快；admixture凭借其运算速度，成为了主流的分析软件。
admixture 软件一共分为5步：
# step 1
/USER/zhusitao/Software/vcftools-master/bin/vcftools –vcf total.final.snp.vcf –plink –out xj
# step 2
/USER/zhusitao/Software/plink/plink –noweb –file xj –geno 0.05 –maf 0.05 –hwe 0.0001 –make-bed –out QC
# step 3
for K in 1 2 3 4 5 6 7;do /USER/zhusitao/Software/admixture_linux-1.3.0/admixture –cv QC.bed $K|tee log${K}.out;done
# step 4 besk K
grep -h CV log*.out
# step 5 画图
tbl = read.table(“/USER/zhusitao/Project/xj/reseq/result/11.admixture/QC.7.Q”);
pdf(“/USER/zhusitao/Project/xj/reseq/result/11.admixture/Q7.pdf”)
barplot(t(as.matrix(tbl)),col= rainbow(7),xlab=”Individual”, ylab=”Ancestry”,border = NA,space = 0);
dev.off()
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/219165.html原文链接：https://javaforall.net
