I got 2-sample matrix.
samples = hl.literal(set([‘sample1’, ‘sample2’]))
mt_test = mt.filter_cols(samples.contains(mt.s), keep = True)
Now, I want to take all reference genotypes to calculate how similar are they. Finally, I hope to get one number (for example ratio of common reference genotypes to all). This is will do to check my PI_HAT metrics. But I got some problems in next …
mt_test.aggregate_rows(hl.agg.count_where(mt_test.GT == hl.literal(set([‘0/0’, ‘0/0’]))))
TypeError: Invalid ‘==’ comparison, cannot compare expressions of type ‘call’ and ‘set’
What am I wrong about?