When the data are not phased, we have n9 = n9a + n9b. This then gives∏i∈diseased,notphasedP(gmi,gfi,gci|θ,disi=1,Ψi=0)P(Ψi=0|disi=1)=[(1−r)R2S2ImIpγ22μ1]n1×…×[(1−r)R2S1ImIpγ12μ4]n8×[(1−r)R1S1(Ip+Im)γ11μ4]n9×[(1−r)S1μ4]n10×[(1−r)μ6]n15.Therefore,∏i∈diseasedP(gmi,gfi,gci,Ψi|θ,disi=1)=rn1′(1−r)n1×…×rn8′(1−r)n8×rn10′(1−r)n10×…×rn15′(1−r)n15×[R2S2ImIpγ22μ1]N1×…×[R2S1ImIpγ12μ4]N8×[rR1S1Ipγ11μ4]n9a′×[rR1S1Imγ11μ4]n9b′×[(1−r)R1S1(Ip+Im)γ11μ4]n9×[S1μ4]N10×…×[μ6]N15.