casement9596 {phenModel}

Casement 1995-1996 Data

Description

Daily minimum and maximum temperature data obtained at the Casement Aerodrome weather station for the period of 1st November 1995 to 31st October 1996.

Usage

data("casement9596")

Format

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

Year

a numeric vector

Month

a numeric vector

Day

a numeric vector

Lat

a numeric vector

Tmin

a numeric vector

Tmax

a numeric vector

day

a numeric vector

Details

The data was obtained from the Irish National Meteorological Agency (Met Éireann).

Source

https://www.met.ie/climate/available-data/historical-data

Examples

data(casement9596)

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

## evaluating model
quantiles_seq <- seq(.001, .95, .001)

## loading fitted parameters
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,
                     pop_quantiles = quantiles_seq,
                     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)

plot_phen_model(obj = phen95, obj_gen2 = phen95_gen2)