d

Author: stsmall

project/stat_modules/sequtils.py

@@ -43,7 +43,7 @@ def add_seq_error(pos, haps, length_bp, perfixder):
         ix = bisect_left(pos_ls, mp)
         pos_ls.insert(ix, mp)
         haps = np.insert(haps, ix, mon_err[:, i], axis=1)
-    print(f"{er1}\t{np.sum(seg_err)}\t{np.sum(mon_err)}\t{der_err}")
+    #print(f"{er1}\t{np.sum(seg_err)}\t{np.sum(mon_err)}\t{der_err}")
     return np.array(pos_ls), haps
 
 

project/stat_modules/sumstats.py

@@ -227,7 +227,10 @@ def FST(self):
         for p1, p2 in combinations(self.stats["pop_config"], 2):
             gtpops = gt.take(p1+p2, axis=1)
             flt = pwpopstats.fst(len(p1), pos, gtpops, quants)
-            stats_ls.extend(flt)
+            try:
+                stats_ls.extend(flt)
+            except TypeError:
+                breakpoint()
         return stats_ls
 
     def dXY(self):