Given a SeuratObject created from a 10X file, pbmc, how can I access the expression values across ALL features for a particular cell? I know I can access the cell names with rownames(pbmc@meta.data) and the matrix with GetAssayData(object = pbmc, slot = "counts"), but I'm really confused as to how to index the matrix such that I can extract a column/row from it. (I am not sure if the data I want is a row or column in the matrix - I believe it is a column).
1 answer
The column names of the output from GetAssayData will be the cell UMIs, so one way you can select the counts for a single cell would be:
cell_umi <- "AAACCCACAATCTCGA"
counts <- Seurat::GetAssayData(scrna_object)
cell_counts <- counts[, cell_umi]
You can also use Seurat::subset with the cells parameter to subset your Seurat object to a specific cell or subset of cells via cell UMIs.
Log in to answer this question.
-
What is the difference between the SCT data stored in data and scale.data in Seurat?
written by zhang616123 •
-
Adata.raw.X in LIANA, something wrong with conversion from Seurat to adata in python.
written by beginner123 •
-
Finding differentially expressed genes between Seurat clusters
written by Nona •
-
Cellcyclescoring did not work in Seurat V5
written by synat.keam •
-
How to introduce normalized and scaled seurat data into monocle 3?
written by Sun •
-
How do I project known tSNE coordinates onto SeuratObject
written by samuel.storey •
-
Some cells in meta.data not present in provided counts matrix.
written by Chise •
-
Convert very large sparse matrix to matrix (scRNA-seq analysis)
written by daianehemerich •
-
get counts by gene for single cell rna seq
written by das2000sidd •
-
Manually logging TPM count matrix in Seurat
written by ww22runner •
