Added function Pcontam to create matrix of likelihood of coming from a pair of populations if the individual is contamintated.  Fixed minor issue in P_likelihood.