% It is assumed that every individual has both father or mother in the pedigree, % or neither (= is a founder). Being a founder therefore is equivalent to not having % a father in the pedigree: @founder(w) = (n-or{sformula(father(u,w))|u: u=u}:0,1); @pchrom_mother_is_A(v)=n-or{pchromA(u)|u: mother(u,v)}; @mchrom_mother_is_A(v)=n-or{mchromA(u)|u: mother(u,v)}; @pchrom_father_is_A(v)=n-or{pchromA(u)|u: father(u,v)}; @mchrom_father_is_A(v)=n-or{mchromA(u)|u: father(u,v)}; @pchrom_mother_is_B(v)=n-or{pchromB(u)|u: mother(u,v)}; @mchrom_mother_is_B(v)=n-or{mchromB(u)|u: mother(u,v)}; @pchrom_father_is_B(v)=n-or{pchromB(u)|u: father(u,v)}; @mchrom_father_is_B(v)=n-or{mchromB(u)|u: father(u,v)}; mchromA(v) = (@founder(v): 0.3, (@pchrom_mother_is_A(v): (@mchrom_mother_is_A(v): 0.9, 0.49 ), (@mchrom_mother_is_A(v): 0.49, (@mchrom_mother_is_B(v): (@pchrom_mother_is_B(v): 0.05, 0.02 ), (@pchrom_mother_is_B(v): 0.02, 0.05 ) ) ) ) ); mchromB(v) = (mchromA(v): 0, (@founder(v): 0.42857, (@pchrom_mother_is_B(v): (@mchrom_mother_is_B(v): 0.94737, (@mchrom_mother_is_A(v): 0.9608, 0.5 ) ), (@pchrom_mother_is_A(v): (@mchrom_mother_is_A(v): 0.5, (@mchrom_mother_is_B(v): 0.9608, 0.0392 ) ), (@mchrom_mother_is_A(v): 0.0392, (@mchrom_mother_is_B(v): 0.5, 0.0526 ) ) ) ) ) ); mchromO(v) = (mchromA(v): 0, (mchromB(v): 0, 1 ) ); pchromA(v) = (@founder(v): 0.3, (@pchrom_father_is_A(v): (@mchrom_father_is_A(v): 0.9, 0.49 ), (@mchrom_father_is_A(v): 0.49, (@mchrom_father_is_B(v): (n-or{pchromB(u)|u: father(u,v)}: 0.05, 0.02 ), (n-or{pchromB(u)|u: father(u,v)}: 0.02, 0.05 ) ) ) ) ); pchromB(v) = (pchromA(v): 0, (@founder(v): 0.42857, (n-or{pchromB(u)|u: father(u,v)}: (@mchrom_father_is_B(v): 0.94737, (@mchrom_father_is_A(v): 0.9608, 0.5 ) ), (@pchrom_father_is_A(v): (@mchrom_father_is_A(v): 0.5, (@mchrom_father_is_B(v): 0.9608, 0.0392 ) ), (@mchrom_father_is_A(v): 0.0392, (@mchrom_father_is_B(v): 0.5, 0.0526 ) ) ) ) ) ); pchromO(v) = (pchromA(v): 0, (pchromB(v): 0, 1 ) ); bloodtypeA(v) = (pchromA(v): (mchromB(v): 0.03, 0.9 ), (pchromB(v): (mchromA(v): 0.03, 0.04 ), (mchromA(v): 0.9, (mchromB(v): 0.04, 0.03 ) ) ) ); bloodtypeB(v) = (bloodtypeA(v): 0, (pchromB(v): (mchromA(v): 0.0309, 0.9375 ), (pchromA(v): (mchromA(v): 0.3, (mchromB(v): 0.0309, 0.4 ) ), (mchromA(v): 0.4, (mchromB(v): 0.9375, 0.0309 ) ) ) ) ); bloodtypeAB(v)= (bloodtypeA(v):0,(bloodtypeB(v): 0, (pchromA(v): (mchromA(v): 0.42857, (mchromB(v): 0.95745, 0.5 ) ), (pchromB(v): (mchromA(v): 0.95745, 0.5 ), (mchromO(v): 0.04255, 0.5 ) ) ) )); bloodtypeO(v) = (bloodtypeA(v):0,(bloodtypeB(v):0,(bloodtypeAB(v):0,1)));