forward_back_allele function

Forward-backward algorithm for allele HMM

forward_back_allele(hmm)

Arguments

• hmm: HMM object; expect variables x (allele depth), d (total depth), logPi (log transition prob matrix), delta (prior for each state), alpha (alpha for each state), beta (beta for each state), states (states), p_s (phase switch probs)

Returns

numeric matrix; posterior probabilities

Examples

forward_back_allele(pre_likelihood_hmm)