R codes and their outputs used for: Evolution of realised niche breadth diversity driven by community dynamics

Why many herbivorous insects are host-plant specialists, with non-negligible exceptions, is a conundrum of evolutionary biology, especially because the host plants are not necessarily optimal larval diets. Here, I present a novel model of host-plant preference evolution of two insect species. Because habitat preference evolution is contingent upon demographic dynamics, I integrate the evolutionary framework with the modern coexistence theory. The results show that the two insect species can evolve into a habitat specialist and generalist when they experience both negative and positive frequency-dependent community dynamics. This happens because the joint action of positive and negative frequency dependence creates multiple (up to nine) eco-evolutionary equilibria. Furthermore, initial condition dependence due to positive frequency dependence allows specialization to poor habitats. Thus, evolved habitat preferences do not necessarily correlate with the performances. The model provides ex..., Individual based (i.e., stochastic) model of eco-evolutionary dynamics. The simulations were run on a 32-cored Win machine. Memory capacity of 32 GB will be sufficient., , # R code used for: Evolution of realised niche breadth diversity driven by community dynamics

The R simulation codes, output files, and R visualisation codes underlying the paper.

## Description of the data and file structure

R code (.R) and output files (.txt and R objects, which should be both loaded as data.frame).

See the 00_Read.Me.txt file for more details.

## Sharing/Access information

NA

## Code/Software

Simularions were run using R software version 4.3.1 for Windows.