Data and R code for: "Nineteenth-century land use shape the current occurrence of some plant species, but weakly affects richness and total composition of Central European grasslands"'

IndVal.all.habitats.csv: the results of the IndVal statistics (De Cáceres et al. 2013) for 1,498 species for the historical land use categories calculated across the entire dataset; IndVal.separate.habitats.csv: the results of the IndVal statistics for 1,498 species for the historical land use categories calculated for each habitat type (dry grasslands, mesic grasslands, wet grasslands) separately; ecological.and.disturbance.values.csv: the original Ellenberg-type and disturbance indicator values, and the varimax-rotated components (‘RC’) used in the analysis (data obtained from Tichý et al. 2023 and Midolo et al. 2023; accessible at the FloraVeg.eu website https://floraveg.eu/download/); R code and data for reproducibility. The R code is for illustration purposes only and is based on a subset of 1,184 mesic grassland vegetation plots located in the Czech Republic and in the study area. This is part of the Czech National Phytosociological Database (Chytrý & Rafajová 2003) and the European Vegetation Archive (Chytrý et al. 2016). The data includes the following: data folder: i. indicator.values.csv: the original indicator values for 831 species; ii. plot.data.csv: data for each of the 1,184 vegetation plots, including their historical land use, plot size, bioclimatic variables (‘bio’; Karger et al. 2017), and soil pH (Hengl et al. 2017); iii. species.matrix.csv: community matrix reporting the relative abundance of species (columns) and plot sites (rows). R scripts for species richness, species composition, and species indicator analyses. R script are also rendered in .html with R Markdown.