powered by
Find optimal theta of a chromosome using forward-backward; uses a 3-state allele HMM
approx_theta_post_s3( pAD, DP, R, p_s, t, upper = 0.45, start = 0.05, gamma = 20, r = 0.015 )