I hope you can answer my previous question.
Further to my previous question, Can we do pair-wise comparison on a DEseq2 object from likelihood ratio, considering that it is similar to ANOVA?.
dds_LRT <- DESeq(ddsobj, test="LRT", reduced = ~ 1)
use this object to do pairwise comparisons for the 8 levels of the factor
res_LRT <- results(dds_LRT, lfcThreshold = 0, contrasts = c('B-control', 'A-control'))