Entering edit mode
23 months ago
Synanth
•
0
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
I'm not totally clear on what it is that you a looking for, but if you want the coefficients for each gene model you can use the
coef
function likecoef(dds)