Data and code from: Adaptive harvesting of two trophic levels stabilises predator-prey dynamics. Simulations with Eurasian lynx and European roe deer

This dataset originates from the article "Adaptive harvesting of two trophic levels stabilises predator-prey dynamics. Simulations with Eurasian lynx and European roe deer," which investigates the population dynamics of roe deer (Capreolus capreolus) and lynx (Lynx lynx) under different harvest management scenarios. This dataset contains R scripts with functions used to implement a predator-prey matrix population model, parameter files used in the model to set vital rates and their stochasticity, model output containing simulated population trajectories of roe deer and lynx, and the R script used to produce the figures. The model explores the impact of multiple harvest management scenarios on the predator-prey dynamics over a period of 25 years. These files allow full reproduction of the model results and figures presented in the article.