Abstract
The generation of large metabolomic data sets has created a high demand for software that can fit statistical models to one-metabolite-at-a-time on hundreds of metabolites. We provide the %polynova_2way macro in SAS to identify metabolites differentially expressed in study designs with a two-way factorial treatment and hierarchical design structure. For each metabolite, the macro calculates the least squares means using a linear mixed model with fixed and random effects, runs a 2-way ANOVA, corrects the P-values for the number of metabolites using the false discovery rate or Bonferroni procedure, and calculate the P-value for the least squares mean differences for each metabolite. Finally, the %polynova_2way macro outputs a table in excel format that combines all the results to facilitate the identification of significant metabolites for each factor. The macro code is freely available in the Supporting Information.