Missing read depth mean after filter_entries()

[iris-garden]

Note

The following post was exported from discuss.hail.is, a forum for asking questions about Hail which has since been deprecated.

(Dec 18, 2023 at 14:38) barioux said:

Hi all,

Just wondered if you had a way to calculate the distribution of the mean read depth across all loci. I get ‘nan’ for dp_stats.mean at mean and std dev as follows after applying filter_entries.

mt_filtered = mt.filter_entries(
((hl.is_snp(mt.alleles[0], mt.alleles[1]) & (mt.DP >= snv_min_coverage)) |
(hl.is_indel(mt.alleles[0], mt.alleles[1]) & (mt.DP >= indel_min_coverage)))
)

mt_filtered.variant_qc.dp_stats.summarize()

dp_stats.mean (float64 ):

Non-missing	63103 (100.00%)
Missing	0
Minimum	7.58
Maximum	91.90
Mean	nan
Std Dev	nan

Thanks in advance!