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).