phenModel (version 1.0)

valentia: Valentia Data

Description

Daily minimum and maximum temperature data obtained at the Valentia synoptic weather station for the period of 1st November 1939 to 31st October 2000.

Usage

data("valentia")

Format

A data frame with 22281 observations on the following 7 variables.

year

a numeric vector

month

a numeric vector

day

a numeric vector

Lat

a numeric vector

max_temp

a numeric vector

min_temp

a numeric vector

day_of_year

a numeric vector

Details

The data was obtained from the Irish National Meteorological Agency (Met <U+00C9>ireann).

Examples

# NOT RUN {
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)

plot_phen_model(obj = phen_forecast[,1:7], obj_gen2 = phen_forecast[,8:12],
                forecast = TRUE)

plot_phen_forecast(obj = phen_forecast, year_from = 1939)
# }