Entering edit mode
4.2 years ago
nitinra
▴
50
Hello all,
I am calculating Fst using VCFtools and here is the code I am running:
vcftools --vcf input.vcf --weir-fst-pop pop1.txt --weir-fst-pop pop2.txt --out pop1_pop2
I am getting the mean fst value in negative: -0.0003 (I did expect it to be low since they are very close), but a lot of the loci has negative values in the pop1_pop2.weir. I know that negative Fst is meaningless, so is there a reason why vcftools has values in the negative? Can I remove those negative value loci and recalculate Fst using awk? Thank you!