Change zygosity column info from het=1 to hom=1.
authorTom Clegg <tom@curii.com>
Fri, 4 Mar 2022 19:48:35 +0000 (14:48 -0500)
committerTom Clegg <tom@curii.com>
Fri, 4 Mar 2022 19:48:35 +0000 (14:48 -0500)
refs #18581

Arvados-DCO-1.1-Signed-off-by: Tom Clegg <tom@curii.com>

slicenumpy.go

index 8c35e4e425608c592886936227432228dadd6a1c..49b81ca6e02ae278d272d666ee71f094bffd6e0a 100644 (file)
@@ -1236,7 +1236,7 @@ func allele2homhet(colpair [2][]int8) {
 type onehotXref struct {
        tag     tagID
        variant tileVariantID
-       het     bool
+       hom     bool
        pvalue  float64
 }
 
@@ -1306,7 +1306,7 @@ func (cmd *sliceNumpy) tv2homhet(cgs map[string]CompactGenome, maxv tileVariantI
                xref = append(xref, onehotXref{
                        tag:     tag,
                        variant: tileVariantID(col >> 1),
-                       het:     col&1 == 1,
+                       hom:     col&1 == 0,
                        pvalue:  p,
                })
        }
@@ -1335,7 +1335,7 @@ func onehotXref2int32(xrefs []onehotXref) []int32 {
        for i, xref := range xrefs {
                xdata[i] = int32(xref.tag)
                xdata[xcols+i] = int32(xref.variant)
-               if xref.het {
+               if xref.hom {
                        xdata[xcols*2+i] = 1
                }
                xdata[xcols*3+i] = int32(xref.pvalue * 1000000)