def report(self, summary_json_paths, barcode_summary_h5_path, recovered_cells, cell_bc_seqs):
assert len(cell_bc_seqs) == len(self.matrices)
barcode_summary_h5 = h5.File(barcode_summary_h5_path, 'r')
d = {}
d.update(self._report_genome_agnostic_metrics(
summary_json_paths, barcode_summary_h5, recovered_cells, cell_bc_seqs))
# Compute genome-specific metrics
for i, (genome, matrix) in enumerate(self.matrices.iteritems()):
for key, value in matrix.report(genome,
barcode_summary_h5,
recovered_cells,
cell_bc_seqs=cell_bc_seqs[i],
).iteritems():
key = '_'.join([genome, key])
d[key] = value
return d
评论列表
文章目录