The metabolic data were analyzed by one-way analysis of variance (ANOVA) in SPSS 21.0 software (SPSS Inc., Chicago, IL). When sampling time was included, data were analyzed using a linear mixed model with treatment and treatment by sampling time interaction as fixed effects, and sampling time as a repeated measures variable. The microbial community composition data was analyzed using Wilcoxon rank-sum test in the JMP Pro software (JMP Pro version 13.2.1, SAS Institute Inc. SAS Institute, Cary, NC, USA). The genomic ranks of attributes were evaluated by Correlation, ReliefF, Symmetrical Uncert, and multi-cluster feature selection (MCFS) methods in software of Waikato Environment for Knowledge Analysis (WEKA) (version 3.8.4, Hamilton, New Zealand) [32 ], and further comprehensively analyzed by RobustRankAggreg (RRA) R package [33 (link)]. All p values were adjusted for False Discovery Rate (FDR) using the Benjamini-Hochberg method. Statistical significance was declared at p ≤ 0.05 and tendencies at 0.05 < p ≤ 0.10.
Free full text: Click here