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 345de78 commit 7d10ae4Copy full SHA for 7d10ae4
sgkit/stats/aggregation.py
@@ -1058,7 +1058,7 @@ def individual_heterozygosity(
1058
variables.validate(ds, {call_allele_count: variables.call_allele_count_spec})
1059
1060
AC = da.asarray(ds.call_allele_count)
1061
- K = AC.sum(axis=-1)
+ K = da.sum(AC, axis=-1)
1062
# use nan denominator to avoid divide by zero with K - 1
1063
K2 = da.where(K > 1, K, np.nan)
1064
AF = AC / K2[..., None]
0 commit comments