DESeq2 genes which vary between 4 groups
I am trying to identify which genes are significantly variable in four groups. How I was doing this was through:
Samples_results<- results(Samples_dds, alpha = 0.05, lfcThreshold = 1)
Which then provided me a list of a few hundred genes which were over the threshold. However, I am not sure if this was the correct way of testing this, and I was wondering if anyone had any advice?
Thanks!
• 837 views
•
link
updated
by
Kevin Blighe
•
written
by
j.obrien2013 •
•
written
by
j.obrien2013 •
1 answer
It sounds like you want to do a likelihood ratio test in order to capture genes whose variances differ across your four groups. If this is the case, then you will likely want to perform, for example:
dds <- DESeqDataSetFromMatrix(..., design= ~ group)
dds <- DESeq(dds, test = 'LRT', reduced = ~1)
res <- results(dds)
Please take a look at the 'likelihood ratio test' section of the vignette.
Kevin
• 0 views
•
link
written
by
Kevin Blighe
Log in to answer this question.
More posts like this
-
Interpretaion of tensorQTL output for ieQTL analysis
written by cool.aditya286 •
-
PCA interpretation
written by nuriagaralon •
-
How to select differential expressed genes using DeSeq2?
written by Aynur •
-
Single model matrix_design formula_RNASeq
written by bhanu.chandra1
-
DESeq2 LFC thresholding strategy and its affect on p-value adjustment
written by ivn25 •
-
RNA-Seq analysis problem
written by Lepomis_8 •
-
Differential gene expression analysis by DESeq2
written by umeshtanwar2 •
-
Many Differentially expressed genes but few GO terms
written by firestar
-
DESeq2 with different replicates
written by robertorun •
-
Why is the significance threshold on my cummeRbund volcano plot so high?
written by keisha.carlson •
