Matlab and R files to create the plots in Ekrem et al. 2025 "Temporal niche differentiation often leads to priority effects rather than coexistence: lessons from a marine midge."

Matlab R2024a

R version 2024.04.2+764)

Packages in R:

R.matlab version 3.6.2

ggtern version 3.3.0

grid version 4.0.2

‘RunModFig….m’: generates data (‘DataFig…_variables.mat’) for figures plotted in R (‘ProduceFigure….R’).

The model uses the following files in Matlab:

‘TrackDynamics.m’ tracks the changes in frequencies of 3 strains over one or multiple seasons consisting of 44 weeks. If the input parameter seasonality_option is zero, then the code runs one season and it always starts with the same moon phase (used to create all the figures in the main text). If seasonality_option=n, where n is some integer larger than 0, then the code runs n consecutive seasons and every season starts with a random moon phase.  The three different pure strains that in the text is referred to as E_M1, E_M2 and E_F, are called S1,S2 and A (or As), respectively, or FM, NM and SL, respectively. The three other genotypes are hybrids of the three pure strains.

‘matingprob.m’ calculates mating success for females during emergence events, using the function in file ‘beach.m’.

‘computetriangle.m’ stores the correct data for producing the ternary plots.