Specific DE comparisons with complex design matrix
1
0
Entering edit mode
2.6 years ago

I have RNAseq count data of hybrid yeast at three different temperatures. There are counts from both alleles. I want to find DE genes at different combinations of comparisons.

This is how I model the data in DESeq2:

colData<-data.frame(
            allele=factor(c(rep("SC",9), rep("SU",9), levels(c("SC","SU")))),
            temperature=factor(rep(c(rep("22",3), rep("33",3), rep("37",3)),2)))

dds <- DESeqDataSetFromMatrix(all_temp_data, colData, formula(~allele+temperature+allele*temperature))

Now, for example, if I want to find DE genes between two alleles at 33 degrees, is the following a correct way of doing it?

dds <- DESeq(dds)
results(dds, contrast=list(c("allele_SU_vs_SC","alleleSU.temperature33")), test="Wald")

Thanks

DESeq2 • 469 views
ADD COMMENT
2
Entering edit mode
2.6 years ago

Make a new column for colData that concatenates allele with temperature. Then

design <- ~ group

and then

res <- results(dds, contrast=c("group", "SC_33", "SU_33"))
ADD COMMENT
1
Entering edit mode

I also don't think your design is right. If you are looking for interactions, and not main effects, I think you do

design <- ~ allele*temperature

or

design <- ~ allele + temperature + allele:temperature

But that's not what you are asking about right now.

ADD REPLY
0
Entering edit mode

Hi swbarnes2, thanks for reply. You are right, I have confused * and :. But given the right design that you mention, could you please elaborate a bit why my solution is not correct (or point out to relevant literature)? I feel that I don't completely understand how design formulas work.

ADD REPLY

Login before adding your answer.

Traffic: 2697 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6