# 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