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 202293f commit 7fdce78Copy full SHA for 7fdce78
sgkit/stats/pca.py
@@ -325,7 +325,7 @@ def count_call_alternate_alleles(ds: Dataset, merge: bool = True) -> Dataset:
325
# TODO: Add to public API (https://github.com/sgkit-dev/sgkit/issues/282)
326
AC = count_call_alleles(ds)["call_allele_count"]
327
AC = AC[..., 1:].sum(dim="alleles").astype("int16")
328
- AC = AC.where(~ds.call_genotype_mask.any(dim="ploidy"), AC.dtype.type(-1))
+ AC = AC.where(~ds.call_genotype_mask.any(dim="ploidy"), -1)
329
new_ds = Dataset({"call_alternate_allele_count": AC})
330
return conditional_merge_datasets(ds, new_ds, merge)
331
0 commit comments