Calculate allele likelihoods
calc_allele_lik(pAD, DP, p_s, theta, gamma = 20)
integer vector Paternal allele counts
integer vector Total alelle counts
numeric vector Phase switch probabilities
numeric Haplotype imbalance
numeric Overdispersion in the allele-specific expression