molecule_counter.py 文件源码

def get_cdna_mol_counts_per_gene(self, gene_index, remove_none_gene=True):
        mol_genes = self.get_column('gene')

        num_genes = len(gene_index.get_genes())
        gene_counts = np.bincount(mol_genes, minlength=num_genes + 1)
        if remove_none_gene:
            gene_counts = gene_counts[:num_genes]

        return gene_counts