phen_model_gen2: Phenology Model Evaluation for Further Generations
In phenModel: Insect Phenology Model Evaluation Based on Daily Temperatures

Description Usage Arguments Details Value Author(s) See Also Examples

Evaluates the phenology model for specified population quantiles, based on daily temperatures for further generations.

phen_model_gen2(Tmean, day_length,
                bud_pars, weib_pars, lactin_pars,
                sex_mat_day, pop_quantiles = 0.5,
                max_day_length = 14.92, data)

`Tmean`	name of the column in the data representing the mean temperature
`day_length`	name of the column in the data representing the day length
`bud_pars`	vector of parameters for the budburst model
`weib_pars`	matrix of parameters for the Weibull model for each developmental stage
`lactin_pars`	matrix of parameters for the lactin model for each developmental stage
`sex_mat_day`	the day individuals from the first generation reached the sexual maturation stage
`pop_quantiles`	vector of quantiles of the population to be evaluated by the phenology model
`max_day_length`	critical day length threshold
`data`	a `data.frame`

See phen_model

A data.frame containing model evaluations for each life cycle (columns) and population quantile (rows).

Rafael de Andrade Moral (rafael.deandrademoral@mu.ie) and Rowan Fealy

plot_phen_model

data(casement9596)

casement1995 <- data.frame(casement9596, 
                           with(casement9596, 
                                get_total_info(Tmax, Tmin, Tbase = 7, Lat, day)))

data(pvulg_budburst_parameters)
data(pvulg_weibull_parameters)
data(pvulg_lactin_parameters)

phen95 <- phen_model(Tmean = "Tmean",
                     thermal_units = "thermal_units",
                     chill_days = "chill_days",
                     day_length = "day_length",
                     bud_pars = pvulg_budburst_parameters,
                     weib_pars = pvulg_weibull_parameters,
                     lactin_pars = pvulg_lactin_parameters,
                     max_day_length = 14.92,
                     data = casement1995)
           
phen95_gen2 <- phen_model_gen2(Tmean = "Tmean",
                               day_length = "day_length",
                               bud_pars = pvulg_budburst_parameters,
                               weib_pars = pvulg_weibull_parameters,
                               lactin_pars = pvulg_lactin_parameters,
                               sex_mat_day = phen95["0.05","sexual_maturation_day"],
                               pop_quantiles = .05,
                               data = casement1995)