Data analysis
We applied separate linear mixed-effects models to each group of
communities defined by their functional similarities and for each type
of climatic event (moderate and extreme dry and wet events). Resilience
and resistance were treated as response variables and were modeled as a
function of a single fixed effect of biodiversity descriptors (taxonomic
richness or diversity, functional redundancy or functional response
diversity), and the season of the climatic event as a random effect
(summer, spring, autumn, winter). The fixed effects were standardized to
allow comparison of the effects among biodiversity descriptors. We used
the packages lme4 (Bates et al. 2015), lmerTest (Kuznetsovaet al. 2017) and MuMIn (BartoĊ 2020) to run all models.
For data treatment and manipulation, we used tidyverse package (Wickhamet al. 2019). The figures were created using ggplot2, from
tidyverse package. All data treatment and analyses were performed in R
(R Core Team 2019).