Data from: Higher-order species interactions cause time-dependent niche and fitness differences: experimental evidence in plant-feeding arthropods

trajectories.csv (the raw data)

id: replicate identifiervariant: co-existence status ("competition" or "monoculture")day: day of experimentspecies: mite species ("CRM" or "WCM")n: population density

model.R

The R script with the GAMM fitted to the trajectory data (the GAMM model is saved as model.RData); also produces simulations from this model (saved as sim.csv).

model.RData

The GAMM for growth rates.

sim.csv (simulations from the GAMM)

day: day of experimentspec_var: combination of co-existence status ("competition" or "monoculture") and species ("CRM" or "WCM")X1:X1000: population densities simulated from the fitted GAMM (on the log scale)

NFD_over_time_monte_carlo_gam.py

Python script to compute niche and fitness differences. Takes sim.csv (densities over time for different instantiations) and model.RData (stores the GAMM from R for the growth rates) as input and generates the file Data_NFD_monte_carlo_multi_c.csv which stores the niche and fitness differences computed for these communities.

figures.R

The R script that produces Figures 2-4.

plot_biotic_model.py

Python code to generate the figures S3 and S4 showing the simulations of a biotic resource competition model. 

plot_abiotic_model.py

Python code to generate the figures S1 and S2 showing simulations of an abiotic resource competition model.