def test_xrec(): efam = EvalFamily(fam) fam.subjects[0].sex = 'female' fam.subjects[1].sex = 'male' fam.subjects[2].sex = 'female' efam.gt_types = [Family.HET, Family.HOM_REF, Family.HOM_ALT] assert efam.x_rec() fam.subjects[2].sex = 'male' efam.gt_types = [Family.HET, Family.HOM_REF, Family.HOM_ALT] assert efam.x_rec() fam.subjects[2].sex = 'female' # mom is hom_alt, but not affected efam.gt_types = [Family.HOM_ALT, Family.HOM_REF, Family.HOM_ALT] assert not efam.x_rec()
def test_x_rec(): mom = Sample('mom_1239NIH', affected=False, sex='female') dad = Sample('dad_1240NIH', affected=False, sex='male') kid_aff = Sample('kidaff_1238NIH', affected=True, sex='female') kid_aff.mom = mom kid_aff.dad = dad efam = EvalFamily(Family([dad, mom, kid_aff], 'oler-trio')) # mom should be a carrier efam.gt_types = [Family.HOM_REF, Family.HOM_REF, Family.HOM_ALT] assert efam.x_rec()