将VCF转换为基因流动分析(Treemix)所需文件出现问题

Question

我用以下流程进行分析：

vcftools --gzvcf  DK.vcf.gz  --plink-tped --out DK
plink --tfile DK --aec --chr-set 31 --freq --noweb --missing --within pop.cov
gzip plink.frq
plink2treemix.py plink.frq.gz treemix.frq.gz

vcftools和plink的流程都没有报错，但plink的结果是plink.frq.strat，与plink.frq不一样，我按照接下来的流程进行分析，但plink2treemix.py的运行结果treemix.frq.gz的**文件大小却是0**，这是什么原因呢，求指导，万分感谢。

Answer 1

http://blog.sina.com.cn/s/blog_bd549ac30102x03j.html