Determining up-regulation vs down-regulation DESeq2 likelihood ratio test

Determining up-regulation vs down-regulation DESeq2 likelihood ratio test

0

Hello,

I am hoping to receive some guidance on interpreting an RNA-sequencing project I am conducting. For this project, I have 9 time points, 3 biological replicates per time point, and 2 genotypes: wildtype vs mutant. I am using a likelihood ratio test in DESeq2 using the following code:

dds <- DESeqDataSetFromMatrix(ctdata, coldata, design = ~ genotype + age
dds <- DESeq(dds, test="LRT", reduced = ~ age)

From a DESeq2 training page, the log2 fold change is printed in the results table for consistency with other results table outputs, but is not associated with the actual test. However, I am interested in detecting upregulated and downregulated genes. Seeing as I cannot use the log2 fold change from the DESeq results to do so, does anyone have any advice on how to go about this?


differential


RNA-sequencing


DESeq2


expression

• 15 views

Read more here: Source link