Moisture content and water activity relations in honey: A Bayesian multilevel meta-analysis

Abstract

Research objectives were to investigate i) linear moisture content (m) - water activity (aw) relations in honeys, ii) quantify variability due to geographical origin, iii) prediction of aw from m using this quantified variability. Bayesian multilevel modeling was used to reach these objectives using thirty literature data sets. The rate of change in aw as function of m (the slope) was identical per origin, whereas actual values of aw at the same m (intercepts) differed. Geographic variation characterized by multilevel modeling (partial pooling) was quantified and compared to single-level modeling with all data pooled (complete pooling) and single-level modeling for each origin (no-pooling). Multilevel modeling results in an overall slope and intercept at the population level but also in individual intercepts per origin, hence variation is characterized at two levels. The obtained multilevel population parameters predicted almost exactly the relation between aw and m of glucose-fructose solutions resembling honey, confirming that aw-m relations in honeys are determined by glucose/fructose but not the actual aw-values themselves. Multilevel modelling, a compromise between over- and underfitting, gives the best prediction of aw of honeys from m including origin variability. The applied procedure is a blueprint to characterize variation in all types of food.