Log-Fold Change in RNA-seq
Hello everybody, I am quite new to RNA-seq. I am using DESeq2 for differential expression analysis and I am at a stage where I have to identify for each cell type the 10 highest expressed, specific genes (>5-fold higher). I have 3 conditions (3 cell types) and each condition has ( 3, 5, 2 samples). Can you give me advice on how to do it?
• 1,230 views
•
link
updated
by
padwalmk
•
written
by
thepurplebluesociety •
•
written
by
thepurplebluesociety •
1 answer
Hi, First check the tutorial from the DESEq it has explained everything very nicely.
You can follow the following steps for now.
- First you have to extract all three contrast result using the result <- data.frame(results(dds, contrast = "blah_blah"))
- Then subset according to your criteria subset(result, criteria > "whatever")
Log in to answer this question.
More posts like this
-
Am I accounting for replicates correctly in pseudobulk analysis?
written by Aspire
-
Using Pseudobulk Approach for Identifying Marker Genes Within a Single Condition
written by nhaus •
-
differential gene expression analysis comparing different conditions using scRNAseq
written by Sara
-
Neuronal markers are expressed in all cell types in 10X single-cell RNASeq
written by viktorfeketa •
-
How do I compare gene expression levels between clusters in single-cell RNA-seq dataset?
written by Sergei C •
-
Should I normalize the data before or after filtering to a particular cell type?
written by gt •
-
Comparing two conditions within multiple different groups with unequal sizes in DESeq
written by heather.jackson17 •
-
Quality Check of Samples and Biological Replicates
written by doinelpierrot •
-
Differential Expression and Co-expression analysis from Cell type-specific datasets
written by volvicpellegrino •
-
Identify Upregulated Genes From Ranks
written by Jessica
