lvclark / polyRAD

Genotype Calling with Uncertainty from Sequencing Data in Polyploids 🍌🍓🥔🍠🥝
24 stars 8 forks source link

Betabinomial #3

Closed lvclark closed 5 years ago

lvclark commented 5 years ago

Updated polyRAD to use the beta-binomial distribution for estimating genotype likelihood rather than the binomial distribution. The depthSamplingPermutations slot also now contains the log of n choose k, rather than simply n choose k, making this version of polyRAD backwards-incompatible with RADdata objects from previous versions. An overdispersion argument has been added to all pipelines, and there is a function to assist users in optimizing the value for this argument. A function for adding even priors has also been added; it is not currently incorporated into any pipeline, but may be useful for future testing purposes.