Statistical analysisi in DESeq2 with different condition
Entering edit mode
3 months ago

Hi, I did some analysis in DESeq2 to look at the different gene expression between different condition

condition_breakfast$group2 <-  factor(paste0(condition_breakfast$Braccio, condition_breakfast$Campione)
condition_breakfast$group3 <- factor(paste0(condition_breakfast$Braccio, condition_breakfast$Campione, condition_breakfast$Risposta))

dds <- DESeqDataSetFromMatrix(countData=raw, 
                                             design=~ group2)
dds <- DESe2(dds)
res_AL <- results(dds_tot, contrast=c("group2", "AL", "AL"))

dds_tot <- DESeqDataSetFromMatrix(countData=raw, 
                                             design=~ group3)

dds_tot <- DESeq(dds_tot)
res_totAL <- results(dds_tot, contrast=c("group3", "ALRD", "ALpCR"))

now I want to see if adding "Risposta" in the condition, so in the design, is significant or not to change the gene expression. How I can do it?

I look at the gene and they are different, in particulare in dss_tot they are more but I want to understand if is significant or not.

Thank you

statistics DESeq2 • 265 views
Entering edit mode
3 months ago
Shred ★ 1.4k

What you need to do is a Likelihood Ratio Test, between a full model (including Risposta) and a reduced model (without Risposta).

Quoting the vignette:

DESeq2 offers two kinds of hypothesis tests: the Wald test, where we use the estimated standard error of a log2 fold change to test if it is equal to zero, and the likelihood ratio test (LRT). The LRT examines two models for the counts, a full model with a certain number of terms and a reduced model, in which some of the terms of the full model are removed.


Login before adding your answer.

Traffic: 2530 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6