Hi, I try to check on some entries where the IQV entry value is 5. I use filter_entries() to do the filter. It works for -1 but not the other values. what am I doing wrong with the code below?
iqv_values = mt.aggregate_entries(hail.agg.counter(mt.IQV))
print(iqv_values)
{5: 131699, 2: 6, 3: 6, -1: 220296658, 4: 6}
mt_gt_enc_iqv = mt.select_entries(‘GT’, ‘ENC’, ‘IQV’)
mt_gt_enc_iqv.show(5)
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
| locus | alleles | ‘10’.GT | ‘10’.ENC | ‘10’.IQV | ‘100’.GT | ‘100’.ENC | ‘100’.IQV |
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
| locus<deid_reference_genome> | array | call | int32 | int32 | call | int32 | int32 |
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
| D:1 | [“A”,“46548”] | 0|0 | 1 | -1 | 0|0 | 1 | -1 |
| D:2 | [“A”,“21381”] | 1|1 | -1 | -1 | 0|1 | 0 | -1 |
| D:3 | [“A”,“10273”] | 0|1 | 0 | -1 | 1|1 | -1 | -1 |
| D:4 | [“A”,“21708”] | 0|0 | 1 | -1 | 0|1 | 0 | -1 |
| D:5 | [“A”,“10270”] | 0|1 | 0 | -1 | 0|1 | 0 | -1 |
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
mt_iqv = mt_gt_enc_iqv.filter_entries(mt_gt_enc_iqv.IQV == -1)
mt_iqv.show(5)
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
| locus | alleles | ‘10’.GT | ‘10’.ENC | ‘10’.IQV | ‘100’.GT | ‘100’.ENC | ‘100’.IQV |
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
| locus<deid_reference_genome> | array | call | int32 | int32 | call | int32 | int32 |
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
| D:1 | [“A”,“46548”] | 0|0 | 1 | -1 | 0|0 | 1 | -1 |
| D:2 | [“A”,“21381”] | 1|1 | -1 | -1 | 0|1 | 0 | -1 |
| D:3 | [“A”,“10273”] | 0|1 | 0 | -1 | 1|1 | -1 | -1 |
| D:4 | [“A”,“21708”] | 0|0 | 1 | -1 | 0|1 | 0 | -1 |
| D:5 | [“A”,“10270”] | 0|1 | 0 | -1 | 0|1 | 0 | -1 |
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
mt_iqv = mt_gt_enc_iqv.filter_entries(mt_gt_enc_iqv.IQV == 5)
mt_iqv.show(5)
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
| locus | alleles | ‘10’.GT | ‘10’.ENC | ‘10’.IQV | ‘100’.GT | ‘100’.ENC | ‘100’.IQV |
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
| locus<deid_reference_genome> | array | call | int32 | int32 | call | int32 | int32 |
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+
| D:1 | [“A”,“46548”] | NA | NA | NA | NA | NA | NA |
| D:2 | [“A”,“21381”] | NA | NA | NA | NA | NA | NA |
| D:3 | [“A”,“10273”] | NA | NA | NA | NA | NA | NA |
| D:4 | [“A”,“21708”] | NA | NA | NA | NA | NA | NA |
| D:5 | [“A”,“10270”] | NA | NA | NA | NA | NA | NA |
±-----------------------------±--------------±--------±---------±---------±---------±----------±----------+