Hi

I'm working on a model matrix for DESeq2, and have used the trick of grouping factors e.g.,

group <- paste0(coldata$time, ".", coldata$condition)
coldata$time.cond <- factor(group)
dds <- DESeqDataSetFromMatrix(countData = counts, colData=coldata, design=~time.cond)

This works when examining individual contrasts, but I was wondering if there were any way to get the results_names to reflect all permutations (mainly for transferring coefficients to a heatmap.)

Thank you so much,
Synanth