We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent e3166e8 commit 0dac590Copy full SHA for 0dac590
malariagen_data/anoph/fst.py
@@ -82,10 +82,11 @@ def _fst_gwss(
82
).compute()
83
84
with self._spinner(desc="Compute Fst"):
85
- fst = allel.moving_hudson_fst(ac1, ac2, size=window_size)
86
- # Sometimes Fst can be very slightly below zero, clip for simplicity.
87
- fst = np.clip(fst, a_min=clip_min, a_max=1)
88
- x = allel.moving_statistic(pos, statistic=np.mean, size=window_size)
+ with np.errstate(divide="ignore", invalid="ignore"):
+ fst = allel.moving_hudson_fst(ac1, ac2, size=window_size)
+ # Sometimes Fst can be very slightly below zero, clip for simplicity.
+ fst = np.clip(fst, a_min=clip_min, a_max=1)
89
+ x = allel.moving_statistic(pos, statistic=np.mean, size=window_size)
90
91
results = dict(x=x, fst=fst)
92
0 commit comments