Capture MGCV error

[njtierney]

library(conmat)
library(readr)
library(dplyr)

oz_pop <- abs_pop_age_lga_2020 %>%
group_by(age_group) %>%
summarise(
population = sum(population)
) %>%
mutate(
lower.age.limit = parse_number(as.character(age_group))
) %>%
mutate(
country = "Australia"
)
oz_pop

fit_setting_contacts(
contact_data_list = get_polymod_setting_data(),
population = oz_pop
)
#> Error in if (sum(uconv.ind) == 0) {: missing value where TRUE/FALSE needed

Thanks to @goldingn for working out when this error typically occurs:

Because the interpolated population function was predicting 0 population for some ages, which then led to a non-finite offset

The goal here it to either try and capture that error exactly, and give a more informative error, replacing it with suggestions.

Alternatively, identify this ahead of time somehow, before model fitting, and describe how to fix that issue.