ability to handle sex chromosomes
jtweir opened this issue · 2 comments
Is your feature request related to a problem? Please describe.
The problem is that pixy seems unable to handle sex chromosomes. For example, the X chromosome is homogametic and thus haploid in males. VCF files in which genotypes have been properly called and coded for the X will have haploid genotype calls for males and diploid calls for females. Such VCF inputs are not currently recognized by pixy. Likewise, the Y chromosome cannot be handled because it is haploid.
Describe the solution you'd like
The ability to calculate Fst, pi, and dxy for sex chromosomes containing a mixture of males and females. Currently the only work around I am aware of is to incorrectly calculate diploid genotypes for both males and females on the X and males on the Y which then results in incorrectly calculated Fst, pi, and dxy (though probably not biased in any way).
Hi @jtweir, the correct handling of sex chromosomes is in progress and will be addressed in the next update (along with general support for non-diploid data). I actually need this for my own work, so it will be a high priority for me. As always, the rate-limiting step will be validating the new feature with simulations. I will reference this issue in the release when it is added.
Hi @jtweir, the correct handling of sex chromosomes is in progress and will be addressed in the next update (along with general support for non-diploid data). I actually need this for my own work, so it will be a high priority for me. As always, the rate-limiting step will be validating the new feature with simulations. I will reference this issue in the release when it is added.
I just wanted to say this will be great for people working in haplodiploid species (like bees) that have haploid and diploid individuals in the same population. There's no other software that cares about them…