Skip to content

Commit ce733c8

Browse files
committed
fix error when dad not present for segregating dominant
closes #181
1 parent 10177dc commit ce733c8

2 files changed

Lines changed: 14 additions & 1 deletion

File tree

js/js-tests.html

Lines changed: 12 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -221,6 +221,18 @@
221221
kid.mom = mom; kid.dad = dad
222222
assert_not(segregating_denovo, [[kid, dad, mom]], "affected mom is not segregating de novo")
223223

224+
var dom_mom = Sample(HET, {id:"mom_dom", sex:"female", affected: true})
225+
var dom_dad = Sample(HOM_REF, {id:"dad_dom", sex:"male", affected: false})
226+
var dom_kid = Sample(HET, {id:"kid_dom", sex:"male", affected: true})
227+
as_trio(dom_kid, dom_dad, dom_mom)
228+
assert(segregating_dominant, [[dom_kid, dom_dad, dom_mom]], "autosomal dominant trio with affected het mom")
229+
230+
var dom_pair_mom = Sample(HET, {id:"mom_dom_pair", sex:"female", affected: true})
231+
var dom_pair_kid = Sample(HET, {id:"kid_dom_pair", sex:"female", affected: true})
232+
dom_pair_mom.kids = [dom_pair_kid]
233+
dom_pair_kid.mom = dom_pair_mom
234+
assert(segregating_dominant, [[dom_pair_kid, dom_pair_mom]], "autosomal dominant mom+kid affected pair")
235+
224236

225237
var k1 = Sample(HET, {id:"kid", sex:"female", affected: true})
226238
var m1 = Sample(HOM_REF, {id:"mom", sex:"female", affected: false})

js/slivar-functions.js

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -86,7 +86,8 @@ function segregating_dominant_x(s) {
8686
if(s.sex != "female"){return false; }
8787
// this block enforces inherited dominant, but not find de novos
8888
if(("mom" in s) || ("dad" in s)) {
89-
if(!((("mom" in s) && s.mom.affected && s.mom.het) || (s.dad && s.dad.affected))) { return false;}
89+
if(!((("mom" in s) && s.mom.affected && s.mom.het) ||
90+
(("dad" in s) && s.dad.affected && s.dad.het))) { return false;}
9091
if(("dad" in s) && !hq1(s.dad, true)){ return false; }
9192
if(("mom" in s) && !hq1(s.mom, true)){ return false; }
9293
}

0 commit comments

Comments
 (0)