Question: Classify cells using seurat
0
gravatar for AB
3 months ago by
AB90
United States
AB90 wrote:

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?

seurat scrna monocle • 246 views
ADD COMMENTlink modified 3 days ago by dppb0550 • written 3 months ago by AB90
1

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

ADD REPLYlink written 3 months ago by Madelaine Gogol5.0k
0
gravatar for dppb05
3 days ago by
dppb0550
dppb0550 wrote:

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 COMMENTlink modified 2 days ago • written 3 days ago by dppb0550
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1234 users visited in the last hour