compute_Q1 function

Compute the haplotype transition matrices based on the fastPHASE HMM

compute_Q1(r, alpha)