Community ecology aims to understand what factors determine the assembly and dynamics of species assemblages at different spatiotemporal scales. To facilitate the integration between conceptual and statistical approaches in community ecology, we propose Hierarchical Modelling of Species Communities (HMSC) as a general, flexible framework for modern analysis of community data. While non‐manipulative data allow for only correlative and not causal inference, this framework facilitates the formulation of data‐driven hypotheses regarding the processes that structure communities. We model environmental filtering by variation and covariation in the responses of individual species to the characteristics of their environment, with potential contingencies on species traits and phylogenetic relationships. We capture biotic assembly rules by species‐to‐species association matrices, which may be estimated at multiple spatial or temporal scales. We operationalise theHMSCframework as a hierarchical Bayesian joint species distribution model, and implement it as R‐ and Matlab‐packages which enable computationally efficient analyses of large data sets. Armed with this tool, community ecologists can make sense of many types of data, including spatially explicit data and time‐series data. We illustrate the use of this framework through a series of diverse ecological examples.

A large array of species distribution model (SDM) approaches has been developed for explaining and predicting the occurrences of individual species or species assemblages. Given the wealth of existing models, it is unclear which models perform best for interpolation or extrapolation of existing data sets, particularly when one is concerned with species assemblages. We compared the predictive performance of 33 variants of 15 widely applied and recently emerged SDMs in the context of multispecies data, including both joint SDMs that model multiple species together, and stacked SDMs that model each species individually combining the predictions afterward. We offer a comprehensive evaluation of these SDM approaches by examining their performance in predicting withheld empirical validation data of different sizes representing five different taxonomic groups, and for prediction tasks related to both interpolation and extrapolation. We measure predictive performance by 12 measures of accuracy, discrimination power, calibration, and precision of predictions, for the biological levels of species occurrence, species richness, and community composition. Our results show large variation among the models in their predictive performance, especially for communities comprising many species that are rare. The results do not reveal any major trade‐offs among measures of model performance; the same models performed generally well in terms of accuracy, discrimination, and calibration, and for the biological levels of individual species, species richness, and community composition. In contrast, the models that gave the most precise predictions were not well calibrated, suggesting that poorly performing models can make overconfident predictions. However, none of the models performed well for all prediction tasks. As a general strategy, we therefore propose that researchers fit a small set of models showing complementary performance, and then apply a cross‐validation procedure involving separate data to establish which of these models performs best for the goal of the study.

Functional diversity (FD) is an important component of biodiversity that quantifies the difference in functional traits between organisms. However, FD studies are often limited by the availability of trait data and FD indices are sensitive to data gaps. The distribution of species abundance and trait data, and its transformation, may further affect the accuracy of indices when data is incomplete. Using an existing approach, we simulated the effects of missing trait data by gradually removing data from a plant, an ant and a bird community dataset (12, 59, and 8 plots containing 62, 297 and 238 species respectively). We ranked plots by FD values calculated from full datasets and then from our increasingly incomplete datasets and compared the ranking between the original and virtually reduced datasets to assess the accuracy of FD indices when used on datasets with increasingly missing data. Finally, we tested the accuracy of FD indices with and without data transformation, and the effect of missing trait data per plot or per the whole pool of species. FD indices became less accurate as the amount of missing data increased, with the loss of accuracy depending on the index. But, where transformation improved the normality of the trait data, FD values from incomplete datasets were more accurate than before transformation. The distribution of data and its transformation are therefore as important as data completeness and can even mitigate the effect of missing data. Since the effect of missing trait values pool-wise or plot-wise depends on the data distribution, the method should be decided case by case. Data distribution and data transformation should be given more careful consideration when designing, analysing and interpreting FD studies, especially where trait data are missing. To this end, we provide the R package “traitor” to facilitate assessments of missing trait data.

Detecting interaction between species is notoriously difficult, and disentangling species associations in host‐related gut communities is especially challenging. Nevertheless, due to contemporary methods, including metabarcoding and 16S sequencing, collecting observational data on community composition has become easier and much more common. We studied the previously collected datasets of intestinal bacterial microbiota and parasite compositions within longitudinally followed mouse lemurs by analysing the potential interactions with diversity metrics and novel joint species distribution modelling. Both methods showed statistical association between certain parasite species and bacterial microbiota composition. Unicellular Eimeria sp. had an effect on diversity of gut microbiota. The cestode Hymenolepis diminuta had negative associations with several bacterial orders, whereas closely related species Hymenolepis nana had positive associations with several bacterial orders. Our results reveal potential interactions between some, but not all, intestinal parasites and gut bacterial microbiota. Host variables contributed over half of the total variation explained with the model, and sex was the most important single host variable; especially with microbiota, there were sex‐related differences in the community composition. This study shows how joint species distribution modelling can incorporate both within‐host dynamics of several taxa and host characteristics to model potential interactions in intestinal community. These results provide new hypothesis for interactions between and among parasites and bacterial microbiota to be tested further with experimental studies.

Our current economic and political structures have an increasingly devastating impact on the Earth’s climate and ecosystems: we are facing a biospheric emergency, with catastrophic consequences for both humans and the natural world on which we depend. Life scientists – including biologists, medical scientists, psychologists and public health experts – have had a crucial role in documenting the impacts of this emergency, but they have failed to drive governments to take action in order to prevent the situation from getting worse. Here we, as members of the movement Scientist Rebellion, call on life scientists to re-embrace advocacy and activism – which were once hallmarks of academia – in order to highlight the urgency and necessity of systemic change across our societies. We particularly emphasise the need for scientists to engage in nonviolent civil resistance, a form of public engagement which has proven to be highly effective in social struggles throughout history.