This is a test version of Biostars. For the public version, visit https://www.biostars.org.
scanpy - subset AnnData to remove samples with low numbers of cells

Hi,

I am using scanpy for scRNA-seq analysis. Please could someone help with how to subset an AnnData object to remove samples with a number of cells below a threshold? Could I use an extension of the value_counts() function?

e.g. adata.obs['Sample'].value_counts() < 500

Thank you!

scanpy python scrnaseq
• 6,853 views • link
written by ES • ES

You can filter your cells and genes count with sc.pp.filter_cells(adata, min_genes=200) and sc.pp.filter_genes(adata, min_cells=3), respectively.

• 0 views • link
written by zorbax zorbax

Thank you for your reply, I do use the filters you mentioned. But what I'd like to do is remove a whole sample based on the total number of cells in that sample, rather than just filtering genes and cells (not from a specific sample).

• 0 views • link
written by ES • ES

1 answer

You need to keep only the value counts higher than your threshold value:

cluster_counts = adata.obs['Sample'].value_counts() 
keep = cluster_counts.index[cluster_counts >= 500] 
filtered_adata = adata[adata.obs['Sample'].isin(keep)].copy()
• 0 views • link
written by zorbax zorbax

Great, thank you very much!

• 0 views • link
written by ES • ES

Log in to answer this question.

More posts like this