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

View source: R/phen_model_forecast.R

A forecast based on the phenology model for multiple years, based on daily temperatures.

1 2 3 | ```
phen_model_forecast(Tmean, thermal_units, chill_days, day_length,
bud_pars, weib_pars, lactin_pars,
max_day_length = 14.92, pop_quantile, year, data)
``` |

`Tmean` |
name of the column in the data representing the mean temperature |

`thermal_units` |
name of the column in the data representing the thermal units |

`chill_days` |
name of the column in the data representing the number of chill days |

`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 |

`max_day_length` |
critical day length threshold |

`pop_quantile` |
quantile of the population to be evaluated by the phenology model |

`year` |
a vector of indices representing each yearly cycle (usually the output from |

`data` |
a |

See `phen_model`

A `data.frame`

with the number of rows equal to the total number of yearly cycles, and each column representing the developmental stage of the insect.

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

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 | ```
data(valentia)
valentia$cycle <- get_cycle(data = valentia, day = "day", month = "month",
day.end = 31, month.end = 10)
valentia_split <- split(valentia, valentia$cycle)
valentia_split <- lapply(valentia_split,
function(x) data.frame(x,
with(x, get_total_info(Tmax = x$max_temp,
Tmin = x$min_temp,
Tbase = 7,
lat = x$Lat,
x$day_of_year))))
valentia <- do.call("rbind", valentia_split)
data(pvulg_budburst_parameters)
data(pvulg_weibull_parameters)
data(pvulg_lactin_parameters)
phen_forecast <- phen_model_forecast(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,
pop_quantile = .5,
year = "cycle",
data = valentia)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.