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
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"))

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


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.

Traffic: 2697 users visited in the last hour
FAQ
API
Stats

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