Classify cells using seurat
1
0
Entering edit mode
3.2 years ago
Apoorva ▴ 290

Hi,

I am working on single cell data and I have a general question. Monocle has a function 'classifyCells' to assign celltype using known marker genes. Here is the bit from the documentation

cth <- newCellTypeHierarchy()

MYF5_id <- row.names(subset(fData(cds), gene_short_name == "MYF5"))
ANPEP_id <- row.names(subset(fData(cds), gene_short_name == "ANPEP"))

cth <- addCellType(cth, "Myoblast", classify_func =
    function(x) { x[MYF5_id,] >= 1 })
cth <- addCellType(cth, "Fibroblast", classify_func =
    function(x) { x[MYF5_id,] < 1 & x[ANPEP_id,] > 1 } )

cds <- classifyCells(cds, cth, 0.1)

Is there a similar function is Seurat. Something that allows to classify cells by cell type and adds it to the meta-data in the seurat object?

scRNA seurat monocle • 3.0k views
ADD COMMENT
1
Entering edit mode

I have used SetIdent() to do something like that.

ADD REPLY
0
Entering edit mode
2.8 years ago
dppb05 ▴ 100

Is there a similar function is Seurat. Something that allows to classify cells by cell type and adds it to the meta-data in the seurat object?

There is ClassifyCells function. But keep in mind that this function works differently from monocle::classifyCells. Monocle's docs, as you have probably read, give a good high-level explanation on their classification approach.

For Seurat::ClassifyCells you need some labelled data that will be used to train Random Forest classifiers. Those trained classifiers will then be used to classify your unlabelled data. You could subset your Seurat object (using SubsetData) based on some marker genes and set the ident (using SetIdent) of this subset according to those markers, then use that as your training set.

ADD COMMENT

Login before adding your answer.

Traffic: 1317 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6