Significance of group classification, batch, enrichment, and of the covariable age at entry was determined using the following linear mixed effects model implemented in R studio (2015) statistical software:
where Υ is the trait (metabolite concentration); Group is the group classification of pigs (resilient, middle, susceptible and dead); Batch is the fixed batch effect (N = 1 to 15); NurEnrich is the fixed effect of environmental enrichment (barren or enriched); Age is the covariate of age when the pig entered the quarantine nursery; Pen is the random effect of pen by batch in the quarantine nursery; Litter is the random effect of the litter environmental effect, and e is the error. The residuals of the model were plotted and visually inspected for the presence of outliers. To test differences between group classification and enrichment conditions, the Least Square Means (LSMs) for each pair-wise comparison were estimated and their P values were adjusted using false discovery ratio correction (FDR). Differences were considered significant if P ≤ 0.05 and as a tendency towards significance if P > 0.05 and ≤ 0.10.