Hello all,

I have fold changes computed from various pairwise comparisons done using DESeq2. I would like to compare the relative induction of a few genes from each comparison by plotting their Fold-change and associated 95% confidence interval. I believe I can get a 95% confidence interval from the lfcSE values output from DESeq2, but I'm unsure of exactly how to do this. I've read that multiplying the lfcSE value with normal quantiles, such as qnorm(0.025) can do this, but I would like a bit more clarity on this procedure. Is it as simple as saying 95% confidence is [log2FoldChange - qnorm(0.025)lfcSE, log2FoldChange + qnorm(0.025)lfcSE]?

Thanks in advance,

Answered here by the developer: Error bars DESeq or DESeq2 fold change