Learn R Programming

phenModel (version 1.0)

casement9596: 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")

Arguments

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 <U+00C9>ireann).

Examples

Run this code
# NOT RUN {
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)
# }

Run the code above in your browser using DataLab