仅作为个人笔记使用
bgzip -c input.vcf > input.vcf.gz
tabix -p vcf input.vcf.gz
./smc++ vcf2smc input.vcf.gz ctg1.pop1.smc ctg1 pop1:X1,X2,X3 && ./smc++ vcf2smc input.vcf.gz ctg1.pop1.smc ctg2 pop1:X1,X2,X3 #如果有多个contig比如ctg1和ctg2
./smc++ vcf2smc input.vcf.gz ctg1.pop2.smc ctg1 pop2:Y1,Y2,Y3
......
./smc++ estimate -o ctg.pop1 1e-10 ctg1.pop1.smc #ctg2.pop1.smc #如果有多个contig
./smc++ estimate -o ctg.pop2 1e-10 ctg1.pop2.smc
......
./smc++ plot out.pdf ctg.pop1/model.final.json ctg.pop2/model.final.json
