Entering edit mode
Hi,
Has anyone any experince of extracting the beta coefficients form a limma object (i.e. fit_2) and calculating the 5% and 95% confidence intervals for the coefficients? Any help would be hugely appreciated.
fit <- lmFit(methylation_matrix, mod, method = "ls") # Use this line for full run, hash out for testing
fit_2 <- eBayes(fit)
table <- limma::topTable(fit_2, coef = 2, number = Inf)
Thanks,
Matt
Unfortunately that returns the confidence intervals for logFC, not the B statistic.
You said beta in your original post, which is the logFC. You want the CI for the log odds?
Apologies, ive been a plank and obviously not read closely enough - in the toptable output logFC is essentially the effect size / B statistic for each individual regression per gene / CpG and B the log-odds?
It's not the effect size, which is a different thing. But it is the beta for the model you fit, estimating the difference between groups. The B statistic in that table is the log odds of being differentially expressed.
Ok much appreciated, thanks!