Entering edit mode
Hello,
I am new to using DEseq2 and multiple factor design. I have gone through forum discussions that deal with DEseq2 and multifactorial designs. However, I would like to understand if my design for the following matrix is correct?
| Sample | Celltype | Condition |
| 1 | 1 | KO |
| 2 | 1 | WT |
| 3 | 2 | KO |
| 4 | 2 | WT |
| 5 | 3 | KO |
| 6 | 3 | WT |
| 7 | 1 | KO |
| 8 | 1 | WT |
| 9 | 2 | KO |
| 10 | 2 | WT |
| 11 | 3 | KO |
| 12 | 3 | WT |
| 13 | 1 | KO |
| 14 | 1 | WT |
| 15 | 2 | KO |
| 16 | 2 | WT |
| 17 | 3 | KO |
| 18 | 3 | WT |
| 19 | 1 | KO |
| 20 | 1 | WT |
| 21 | 2 | KO |
| 22 | 2 | WT |
| 23 | 3 | KO |
| 24 | 3 | WT |
That is we have 3 cell types for each condition (knock out (KO) and wild type(WT)). I would like to perform teh following differential analysis
1. For KO condition, I would like to perform a comparision between each cell type
2. For each cell type, I would like to compare between KO and WT
Should I set up the design as follows? ~CellType+CellType:Condition
Thank you in advance and any pointers or help is appreciated.
Regards
kavi
Hi! Sorry to resurrect this post, but I have a similar doubt.
I have this dataset:
I'll need to test the differential expression in "status" 0 against 1, later "tissue" all against all, "gender" M against F and finally all groups between all groups.
The solution for group is the same of this post right? But, for spliced analysis, should I use a different design formula like: status + tissue + gender?
Best regards,
I'd recommend speaking with a local statistician to help you formulate your design. There are many ways to make these comparisons.
Dear Michael,
I have a similar experiment to this one: Design formula and Design matrix in DESeq
I have 2 groups for genotype (Mut vs. Ctr), and 4 time points, with 8 biological replicates each. As I am interested in the interactions, I ran my model as you suggested in the manual and in similar posts:
dge <- DESeqDataSetFromMatrix(countData = counts, colData = phenomodel1, design = ~group), where group is Ctr_time1, Ctr_time2, Ctr_time3, Ctr_time4, Mut_time1, Mut_time2, Mut_time3, Mut_time4.
dgeWALD <- DESeq(dge)
This works great in order to extract the genotype effect specifically for each age (example of genotype effect for age 1: Waldresults1 <- results(dgeWALD, contrast=c("group", "Mut_time1", "Ctr_time1")) .
In order to extract the overall effect of genotype (an average of the effect in time 1, 2, 3, and 4), I used the following:
Waldresultsgenotype <- results(dgeWALD, contrast=list(c("Mut_time1", "Mut_time2", "Mut_time3", "Mut_time4"), c("Ctr_time1", "Ctr_time2", "Ctr_time3", "Ctr_time4")), listValues=c(1,-1))
Now I am also interested in the effect of genotype across aging... And the aging effect in general really... And so I am trying to use the likelihood ratio. However, after reading the manuals and searching for similar posts, I still can't understand how I can run this. So here are my specific questions:
From my understanding, I should run something like this
1) dgeLRT_genotype_aging <- DESeq(dge, test="LRT", full = ~???, reduced = ~age) # for the effect of genotype accross aging (meaning I take out the effect of age and thus stay with the effect of genotype)
2) dgeLRT_genotype_aging <- DESeq(dge, test="LRT", full = ~???, reduced = ~genotype) # for the effect of age, independently of genotype (meaning I take out the effect of genotype and thus stay with the effect of aging)
But my full model was ~group (interaction!), so I really don't understand very well how I define the full model and the reduced model. In every case I see the LRT being used the full model is something like ~condition+genotype, but in this example one would be looking at the effect of genotype whilst controlling for the effect of condition, right? Which is not my interest.
Therefore, my question is how exactly should I define (for the LRT):
a) DESeqDataSetFromMatrix (colData and design)
c) DESeq
d) DESeq LRT (full and reduced)
If you know of any similar posts to similar studies I would really appreciate if you could help me find them.
Thank you in advance!
Isabel
hi Isabel,
Can you post this as a new question? It will help to organize the thread better. And in the new post, can you post your colData and can you say specifically which variable in the colData you are referring to when you say "aging".
Done. Thank you!