Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
CellRanger vs Seurat akitrav1 Bioinformatics 10 08-20-2019 07:13 AM
scSeq Analysis - Seurat and OOBE usage joshualpollack Bioinformatics 0 07-24-2017 12:00 PM
What is the meaning of repeat cluster and unique cluster? seungchul Lee Bioinformatics 0 04-13-2014 12:29 AM
Seurat - Analysis tool for tumor/normal DNA and RNA seq alexischr Bioinformatics 5 10-31-2013 11:16 AM
v3: Effect of high cluster densities on cluster PF and %Q30 pmiguel Illumina/Solexa 3 10-05-2011 06:36 AM

Thread Tools
Old 09-19-2019, 11:38 AM   #1
Location: ca

Join Date: Oct 2012
Posts: 12
Default how to calculate DEG in a certain cluster using Seurat

I am a new user of Seurat, and I'd like to know how to correctly calculate the DEG from a certain cluster/sub-population?

Here are some basic description about the data I have and steps I use -
Samples: 3 control samples and 3 case/disease samples
1) In the first several steps, I used FindVariableFeatures() followed by FindIntegrationAnchors(...anchor.features = 1000...) and IntegrateData() to make a Seurat object (1000 gene x 36000 cells from 3+3 samples)
2) In the next steps, I got the 11 clusters in the tsne plot
3) subset() to get the cells in a Cluster1, and Idents(Cluster1)$Group to divide case and control
4) Then shall I used a) or b) or some other method to calculate DEG?
a) deg_1000 <- FindMarkers(hsc, ident.1 = "ctrl", ident.2 = "case", verbose = FALSE)
b) deg_all <- FindMarkers(hsc, ident.1 = "ctrl", ident.2 = "case", verbose = FALSE, assay="RNA")
As in a) the max number of genes I have is 1000, but from the original assay "RNA", I have totally ~20K genes, for which I want to know the pattern of gene expression changes.
miaom is offline   Reply With Quote


Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 05:48 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO