sgkit icon indicating copy to clipboard operation
sgkit copied to clipboard

Generalize `convert_probability_to_call` to multiallelic and polyploid cases

Open timothymillar opened this issue 1 year ago • 0 comments

This is currently only implemented for the diploid biallelic case. Generalizing it to the multiallelic case would require that the number of possible alleles at each variant is known. See relevant discussion in #1006 about counting alleles at each variant. We can preserve backwards compatibility by assuming the diploid biallelic case when the relevant dimensions are absent (validating that the size of the probabilities array is correct).

timothymillar avatar Aug 09 '23 07:08 timothymillar