Hello, I have following as my experimental design in DESeq2 and I am wondering how to retrieve results for all the pairwise comparisons.
Design
Genotype | Days
WT 50
WT 100
WT 200
A 50
A 100
A 200
B 50
B 100
B 200
Default resultsnames are
Intercept
Genotype_A_vs_WT
Genotype_B_vs_WT
Days_100_vs_50
Days_200_vs_50
GenotypeA.Days100
GenotypeB.Days100
GenotypeA.Days200GenotypeB.Days200
How do I retrieve results for
1. WT.50 vs A.50, WT.50 vs B.50, A.50 vs B.50 and similar comparisons for 100 and 200 days?
2. WT.50 vs WT.100, WT.56 vs WT.100, WT.100 vs WT.200 and similar results for Genotypes A and B?
1 answer
I hope you have a few replicates for each genotype and day combination. If you really only have 9 samples you can't do those comparisons. But if have replicates, make a new column that has genotype and days stuck together. Make that column name your design. You can specify whatever contrast you want when calling for the results, even if it's not in resultsNames.
Log in to answer this question.
-
Interaction analysis using DESeq while adjusting for covariates.
written by Sumit Paliwal •
-
P-value Threshold Consideration in Multi-Sample RNA-Seq Experiment
written by Netanel •
-
DESeq2 design for two variables and pairwise comparisons
written by bioinfo8
-
Performing Differential expression analysis using DESeq with 3 conditions
written by bionewbie •
-
DESeq2 interaction padj NA values
written by meshi.barsheshet •
-
DESeq2: design and/or results for 2 groups and 2 treatments
written by asta.tranholm •
-
missing coeeficients in resultsNames(dds)
written by Assa Yeroslaviz
-
multiple comparisons in DESeq2
written by 2405592M
-
DESEQ2 Design formula testing cell type across condition
written by Barry Digby
-
DESeq2 design and interaction terms
written by ATRX
