Using Multilevel Models to Estimate Variation in Foraging Returns

Distributions of human foraging success across age have implications for many aspects of human evolution. Estimating the distribution of foraging returns is complicated by (1) the zero-inflated nature of hunting returns, as many if not most trips fail, and (2) the substantial variation among hunters, independent of age. We develop a multilevel mixture analysis of human foraging data to address these difficulties. Using a previously published 20-year record of hunts by 147 individual Aché hunters in eastern Paraguay, we estimate returns-by-age functions for both hunting failures and the size of harvests, while also estimating the heterogeneity among hunters. Consistent with previous analyses, we find that most hunters peak around 40 years of age. We can also show, however, that much more of the variation among Aché hunters arises from heterogeneity in failure rates (zero returns), not harvest sizes. We also introduce a new R package, glmer2stan, to assist in defining and fitting similar multilevel mixture models.