gravatar for cwwong13

2 hours ago by

I have an experiment that has 3 treatment groups (A, B and C) using 3 different cell types (WT, OE, KD).

I am interested to get the result of:
within WT, whether the treatment effect is significant? I am not trying to know whether B vs A or C vs B etc are significant, but want to get the result similar to ANOVA.

How to get the result for OE cells and KD cells?

Similarly, I wonder within treatment A, B, and C:
whether the 3 cell types are different?

I have been searching for so long but cannot find a good result for my question.
I have tried (and probably made the question more complicated):

dds <- DESeqDataSetFromMatrix(countData = cts,
                              colData = coldata_filtered,
                              design = ~ Treatment + Cell + Treatment:Cell)
dds <- DESeq(dds)

I found that many forums said it can be tested by the LRT of the reduced model: ~ Treatment + Cell.
But how to extract the overall Cell effect within treatment A, B, and C?



Source link