astheeggeggs / lshmm

code to run Li and Stephens
MIT License
4 stars 3 forks source link

Remove all-NONCOPY ancestors #86

Closed szhan closed 3 months ago

szhan commented 3 months ago

Sometimes, get_ancestral_haplotypes returns ancestral haplotypes that are non-ancestral at variable sites (e.g. in simulations where the recombination rate is high enough and the mutation rate is low enough), so that the haplotypes contain only NONCOPY values. They should be removed before returning the genotype matrix.