Run tmod enrichment tests directly on a limma object
Order the genes according to each of the coefficient found in a limma object and run an enrichment test on the ordered list.
result of linear model fit produced by limma functions lmFit and eBayes
Either the name of the column in f$genes which contains the gene symbols corresponding to the gene set collection used, or a character vector with gene symbols
How the gene names should be ordered: "msd" (default), "pval" or "lfc"
The function to run the enrichment tests. Either tmodCERNOtest or tmodUtest
If not NULL, only run tmod on these coefficients
Further parameters passed to the tmod test function
A list with length equal to the number of coeffients. Each element is the value returned by tmod test function. The list can be directly passed to the functions tmodSummary and tmodPanelPlot.
For each coefficient in the fit returned by the eBayes / lmFit functions from the limma package, tmodLimmaTest will order the genes run an enrichment test and return the results.
The ordering of the genes according to a certain metric is the fundament for gene enrichment analysis. tmodLimmaTest allows three orderings: p-values, "MSD" and log fold changes. The default MSD ("minimal significant difference") is the lower boundary of the 95 confidence interval for positive log fold changes, and 0 minus the upper boundary of the 95 better than ordering by p-value or by log fold change. See discussion in the package vignette.