astheeggeggs / lshmm

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

Support multiple alleles #11

Open jeromekelleher opened 2 years ago

jeromekelleher commented 2 years ago

In the simplest case, the reference panel can contain values > 1. Initial version can assume that mutation rates among the k existing alleles are equal. We should keep the possibility open of providing arbitrary transition matrix between states.

May be as easy to do in the first pass as a kxk uniform transition matrix that we create, and then worry about making a parameter later.

astheeggeggs commented 11 months ago

we currently consider the simplest case, whereby multiple alleles are allowed at a given site but where mutation probabilities to each are equal.

szhan commented 1 week ago

It should be now easy to switch the emission probability matrix to support multiple alleles.