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.
, , # Simulations of predator-prey dynamics with Eurasian lynx and roe deer

This dataset contains R scripts and simulation outputs associated with the article:
**Adaptive harvesting of two trophic levels stabilises predator–prey dynamics: Simulations with Eurasian lynx and European roe deer.**

The study implements two stochastic stage-structured matrix population models integrating European roe deer (*Capreolus capreolus*) and Eurasian lynx (*Lynx lynx*) linked in a coupled predator–prey system. The model evaluates how adaptive harvesting strategies applied to one or both trophic levels influence predator–prey dynamics over a 25-year simulation period.

Simulations were replicated 500 times for each harvest management scenario and combination of initial population densities. The provided .rds files correspond to the output of the simulations used to generate the results and figures presented in the associated publication.

Because demographic stochasticity is included in vital rates, reru..., ,