We used the R (v.3.4.1; R Core Team 2017) package {compositions} (v.1.40‐1; van den Boogaart, Tolosana, & Bren, 2014) for our compositional analysis in acomp geometry (van den Boogaart & Tolosana‐Delgado, 2008). We first examined annual diets, where we reported compositional descriptive statistics and variance for annual diets as the closed geometric mean and variance matrix of our centered log‐ratio (clr) transformed data set. The compositional geometric mean better represents the center of compositional data points than the arithmetic mean, and dispersion of compositional data is summarized using a variance matrix of pairwise log‐ratios (Aitchison, 2003). Conventional univariate measures of dispersion (e.g., SD of the arithmetic mean proportion of protein) are not considered informative for multivariate compositional data. For comparison, however, we report both conventional arithmetic and geometric means. We plotted predicted 2‐sigma and 3‐sigma region ellipsoids around the geometric mean.