Table of TukeyHSD corrected p-values given in form of a heatmap.

ANOVA was calculated for differences between groups for each of the SCFAs, pH and cultured bacterial counts (columns of the heatmap). ANOVA was followed by TukeyHSD honesty test to calculate p-values (shown in each cell of the heatmap) for each to each group-to-group comparison (rows). The heatmap was coloured by p-values according to the key below the graph in order to emphasise significant differences (in dark red). For example the lower right corner cell shows that groups E+FM+Cp and CTRL group (labelled E+FM+Cp – CTRL on the graph) are showing slight differences (TukeyHSD corrected P-value = 0.054184) in concentration of propionic acid. The clustering of heatmap rows and columns was done to identify similarities in SCFA changes between conditions and was based on Euclidean distance and complete agglomeration method in R package ggplot. This identified a hot-spot in the lower right corner of the heatmap of sample comparisons showing significant differences in isobutyric, isovaleric and propionic acid.