UsageΒΆ

Import COSG as:

import cosg as cosg

amd import Scanpy as:

import scanpy as sc

Next, load the data via:

adata = sc.datasets.pbmc68k_reduced()

then identify marker genes for each cell group by running:

cosg.cosg(adata, key_added='cosg', groupby='bulk_labels')

and the top marker genes can be visualized via:

sc.pl.rank_genes_groups(adata, key='cosg')