# NOT RUN {
params <- list("p_male"= 0.51,
"birth_rate" = stepfun(c(15,40),c(0,0.05,0)),
"death_rate" = gompertz(0.008,0.02))
death_event <- mk_event_individual(type = "death",
intensity_code = "result = death_rate(age(I,t));")
birth_event <- mk_event_individual(type = 'birth',
intensity_code = "if (I.male) result = 0;
else result=birth_rate(age(I,t));",
kernel_code = "newI.male = CUnif(0, 1) < p_male;")
model <- mk_model(characteristics = get_characteristics(population_df),
events = list(death_event,birth_event),
parameters = params)
summary(model)
# }
# NOT RUN {
# }
Run the code above in your browser using DataCamp Workspace