From 292a9459c8dde94b7b9a4245a1d96d10b527f193 Mon Sep 17 00:00:00 2001 From: Tom Clegg Date: Fri, 4 Feb 2022 16:17:23 -0500 Subject: [PATCH] Don't include het just because corresponding hom passed. refs #18581 Arvados-DCO-1.1-Signed-off-by: Tom Clegg --- slicenumpy.go | 13 +++++-------- 1 file changed, 5 insertions(+), 8 deletions(-) diff --git a/slicenumpy.go b/slicenumpy.go index 601cb40eac..77877b3fca 100644 --- a/slicenumpy.go +++ b/slicenumpy.go @@ -1217,20 +1217,17 @@ func (cmd *sliceNumpy) tv2homhet(cgs map[string]CompactGenome, maxv tileVariantI var onehot [][]int8 var xref []onehotXref for homcol := 4; homcol < len(obs); homcol += 2 { - p := [2]float64{ - pvalue(obs[homcol], cmd.chi2Cases), - pvalue(obs[homcol+1], cmd.chi2Cases), - } - if cmd.chi2PValue < 1 && !(p[0] < cmd.chi2PValue || p[1] < cmd.chi2PValue) { - continue - } for het := 0; het < 2; het++ { + p := pvalue(obs[homcol+het], cmd.chi2Cases) + if cmd.chi2PValue < 1 && !(p < cmd.chi2PValue) { + continue + } onehot = append(onehot, bool2int8(obs[homcol+het])) xref = append(xref, onehotXref{ tag: tag, variant: tileVariantID(homcol / 2), het: het == 1, - pvalue: p[het], + pvalue: p, }) } } -- 2.30.2