DESeq2 genes which vary between 4 groups
1
0
Entering edit mode
12 months ago

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!

RNA-Seq • 216 views
ADD COMMENT
2
Entering edit mode
12 months ago

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

ADD COMMENT

Login before adding your answer.

Traffic: 2124 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