Question: DESeq interactions: What genes that are different between two experiments affects the conditions
1
gravatar for raquelgarza95
10 months ago by
raquelgarza9550 wrote:

Hello!

I'm super noob on DESeq interactions so sorry in advance if this question is duplicated. I just need to ask directly to see if I understand correctly.

I have two experiments performed in young and old mice. I want to test in both for knock down vs control of a certain gene. I want to make a differential expression analysis between conditions (ko vs ctrl) and between experiments (young vs old).

My coldata looks like

    sample     condition    experiment
    mouse1     ctrl            young
    mouse2     ctrl            young
    mouse3     ctrl            young
    mouse4     ko              young
    mouse5     ko              young
    mouse6     ko              young
    mouse7     ctrl            adult
    mouse8     ctrl            adult
    mouse9     ctrl            adult
    mouse10    ko              adult
    mouse11    ko              adult
    mouse12    ko              adult

But I am not really interested on what are all the differences between young and old mice, but only those that have an impact on the ko vs ctrl test.

So my design for the DESeq object is ~ condition + experiment + condition:experiment

> resultsNames(dds)
[1] "Intercept"                "condition_ko_vs_ctrl"     "experiment_young_vs_adult"   "conditionko.experimentyoung"

I am not sure, but the intuitive thing to do would be to take experiment_young_vs_adult. Is this correct?

Also, could someone explain why do I not have "conditionko.experimentadult", "conditionctrl.experimentadult", "conditionko.experimentyoung" in resultsNames(dds)?

Thank you in advance!

rna-seq deseq2 interactions dea • 440 views
ADD COMMENTlink modified 10 months ago by swbarnes27.2k • written 10 months ago by raquelgarza9550
3
gravatar for swbarnes2
10 months ago by
swbarnes27.2k
United States
swbarnes27.2k wrote:

I found this website helpful in working out the right way to deal with interactions

https://rpubs.com/ge600/deseq2

I think you likely want "conditionko.experimentyoung", or what's likely more useful is to relevel age so that young is the control, and then the contrast will be "conditionko.experimentadult"

ADD COMMENTlink modified 10 months ago • written 10 months ago by swbarnes27.2k

Thank you! This was very helpful.

ADD REPLYlink written 10 months ago by raquelgarza9550
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: 2056 users visited in the last hour