2
vcftools --vcf ALL.chr1.phase3_shapeit2_mvncall_integrated_v5.20130502.genotypes.vcf --weir-fst-pop POP1.txt --weir-fst-pop POP2.txt --out fst.POP1.POP2

上記のスクリプトは、Weir と Cokerham の 1984 年の式を使用して、1000 のゲノム人口データの Fst 距離を計算します。この式は、3 つの分散成分、すなわち a、b、c (集団間、集団内の個体間、集団内の個体内の配偶子間) を使用します。

出力は式の結果を直接提供しますが、最終結果に到達するためにプログラムが計算したコンポーネントは提供しません。Vcftools に a、b、c の値を出力させるにはどうすればよいですか?

4

1 に答える 1