Determining up-regulation vs down-regulation DESeq2 likelihood ratio test
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?
• 15 views
Read more here: Source link