search

neurorestore / Augur

Cell type prioritization in single-cell data
MIT License
94 stars 10 forks source link

How could I access the AUC per cell? #18

Open sanchezy opened 2 years ago

sanchezy commented 2 years ago

Hi all,

I would like to subset my Seurat object by the AUC score (to keep only the cells with higher AUC). To add the AUC per cell into the metadata of my Seurat object, I need a list of the barcode and the corresponding AUC scores. How could I access the AUC per cell?

Many thanks for your help!

skinnider commented 2 years ago

Hi @sanchezy, you could try something like this (assuming you have a Seurat object sc, cell type in column cell_type of its metadata, and Augur results augur):

sc@meta.data %<>% left_join(augur$AUC, by = 'cell_type')