R script to simulate the phenology of the pine processionary moth, Thaumetopoea pityocampa, from the egg to the last larval instar (L5).

<p>This folder contains :</p> <p>- ReadMe file with the following explanation</p> <p>- the R script to simulate the phenology of the pine processionary moth (PPM_phenological_model.R),</p> <p>- the flight curve of the pine processionary moth in Orleans (France) in 2019 (Flight_Orleans_2019.csv&nbsp;; col 1 = date, col 2 = number of adult catches, col 3 = cumulated number of catches)&nbsp;</p> <p>- temperature dataset from June 2019 to December 2020 (Temperature.csv; col 1 = mean daily temperature, col 2 = Date, col 3 = Julian day continuously counted from one year to another). These meteorological data were recorded by the agroclimatic station at Orleans (number 45234, north latitude 47.827 &deg;, east longitude 1.909 &deg;), part of the INRAE national agroclimatic network managed by the service unit AgroClim (Avignon, France).</p> <p>The R script was used in R version 3.5.1 (2018-07-02)</p> <p>The R script will call automatically the flight curve and the temperature dataset given that you have placed them in your working directory, it requires the following R packages&nbsp;: lubridate, hydroGOF, ggplot2 (to be installed beforehand). The function coded to simulate the phenology of the pine processionary moth is called &laquo;&nbsp;param&nbsp;&raquo;. Input parameters are&nbsp;:<br> - Rm (Maximal development rate in days^(-1) ),<br> - Tm (Optimum temperature in &deg;C)<br> - To (Spread of curve in &deg;C)<br> for each life stage from the egg to L5.<br> This phenology model accumulates the daily development rate given by the Taylor equation and provides for each life stage (resegg, resL1,&hellip; resL5)&nbsp;: &quot;min First day&quot; (first day with an individual entering at this stage), &quot;max First day&quot; (last day with an individual entering at this stage), &quot;min Last day&quot; (first day with an individual achieving this stage development), &quot;max Last day&quot; (last day with an individual achieving this stage development), &quot;mean duration&quot; (mean duration of this life stage across all individuals). This is a summary of the phenology but &ldquo;tabnb&rdquo; provides the number of individuals present in each stage according to the date.</p> <p>To download the software R, please visit: <a href="https://www.r-project.org/">https://www.r-project.org/</a></p> <p>This work was supported by: the national project called PHENEC (grant from the French Research Agency &ndash; ANR &ndash; number19-CE32-0007; <a href="https://www6.inrae.fr/phenec/Project-description">https://www6.inrae.fr/phenec/Project-description</a>), the SOERE called TEMPO (<a href="https://tempo.pheno.fr/soere-tempo_eng/">https://tempo.pheno.fr/soere-tempo_eng/</a>), and the French region Val de Loire.</p> <p>Associated scientific article:<br> Poitou L, Laparie M, Pincebourde S, Rousselet J, Suppo C, Robinet C (2022) Warming causes atypical phenology in a univoltine moth with differentially sensitive larval stages. Frontiers in Ecology and Evolution, DOI: 10.3389/fevo.2022.825875.</p>