Question: DESeq : base mean value are same for all genes in all the three combinations in compare all levels design
0
gravatar for aiswaryabioinfo
7 weeks ago by
aiswaryabioinfo20 wrote:

My RNA-seq dataset consist of Resis_treated (1st day, 3rd day and 7th day) samples taken under 3 time point with 2 replicated for each time point. I have a condition with 3 levels (day1, day3, day7). I used HTseq-count to get the read counts.

I am trying to get the DEGs between Day3VsDay1, Day7VsDay1 and Day7Vs_Day3. After running the code, to get all of the pairwise comparisons I used lfcShrink() with the contrast= option 3 times.

    sampleCondition<-c('Day1','Day1','Day3','Day3', 'Day7','Day7')
    sampleTable<-data.frame(sampleName=sampleFiles, fileName=sampleFiles, condition=sampleCondition)
    ddsHTSeq<-DESeqDataSetFromHTSeqCount(sampleTable=sampleTable, directory=directory, design=~condition)
    colData(ddsHTSeq)$condition<-factor(colData(ddsHTSeq)$condition, levels=c('Day1','Day3','Day7'))
    ddsHTSeq <- ddsHTSeq[ rowSums(counts(ddsHTSeq)) > 1, ] 
    dds <- DESeq(ddsHTSeq)

After getting all the three combinations

`day73 <- lfcShrink(dds, contrast=c("condition", "day7", "day3"))
 day71 <- lfcShrink(dds, contrast=c("condition", "day7", "day1"))
 day31 <- lfcShrink(dds, contrast=c("condition", "day3", "day1"))`

In the result table of all the three combinations saved separately, It was observed that the base mean value was same for all of three combinations

And also the log2fc is very less in each combination when compared to the results of all the wald test done for each combination (Day3VsDay1, Day7VsDay1 and Day7Vs_Day3) individually as separate three DESeq run. Why is it so?

base mean log2fc deseq2 R • 118 views
ADD COMMENTlink modified 7 weeks ago by h.mon27k • written 7 weeks ago by aiswaryabioinfo20

Cross-posted

https://support.bioconductor.org/p/122326/

ADD REPLYlink written 7 weeks ago by swbarnes26.2k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1499 users visited in the last hour