  BEF studies aim to understand how ecosystems respond to a gradient of species
diversity. Diversity-Interactions (DI) models are suitable for analysing the
BEF relationship. These models relate an ecosystem function response of a
community to the identity of the species in the community, their evenness
(proportions) and interactions. The number of species in the community
(richness) is also implicitly modelled through this approach. It is common in
BEF studies to model an ecosystem function as a function of richness; while
this can uncover trends in the BEF relationship, by definition, species
diversity is much broader than richness alone, and important patterns in the
BEF relationship may remain hidden. In this paper, we introduce the DImodels R
package for implementing DI models. We also compare DI models to traditional
modelling approaches to highlight the advantages of using a multi-dimensional
definition of species diversity. We show that using DI models can lead to
considerably improved model fit over other methods; it does this by
incorporating variation due to the multiple facets of species diversity.
Predicting from a DI model is not limited to the study design points, the model
can extrapolate to predict for any species composition and proportions
(assuming there is sufficient coverage of this space in the study design).
Expressing the BEF relationship as a function of richness alone can be useful
to capture overall trends. However, collapsing the multiple dimensions of
species diversity to a single dimension (such as richness) can result in
valuable ecological information being lost. DI modelling provides a framework
to test the multiple components of species diversity in the BEF relationship.
It facilitates uncovering a deeper ecological understanding of the BEF
relationship and can lead to enhanced inference.
