R-language code to model the current and future geographic distribution (under different climate change scenarios) of the functional diversity of European bats. Species distribution models are performed using Bayesian Additive Regression Trees with optimized hyperparameters.