I am processing transcriptome-wide expression data using the limma pipeline. I have a two factor experiment (2 treatments crossed). My data are standardized with voom, precision weights are calculated with eBayes, and analysis of differential expression is done though limma modeling.
Does anyone know if an LSMean value and confidence intervals for that expression can be obtained from these tools? Limma offers options for outputting CI's on the log fold change between treatments and average expression across all treatment group combinations, but I have yet to find a way to output CI's or LSMean values for the straight-forward normalized expression values of a particular gene in a particular treatment combination.
The R package lsmeans does not appear to work with MArray model objects generated by limma. Voom and eBayes offer options for outputting precision weight values for each gene, but I am unsure how to transform this information into confidence intervals.
I am primarily interested in this information for plotting purposes, for data on individual genes.
Any help appreciated. Thanks in advance.