Add seasonal patterns

Adds seasonal forcing to the model using a Fourier function of rainfall.

set_seasonality(params, season_params, floor = 0.001)

Arguments

params

Other malariasimple parameters

season_params

List of Fourier parameters for manual seasonal input

floor

Minimum value of seasonal forcing

Value

Updates the input parameter list to include seasonal parameters

Examples

season_params <- list(g0 = 0.28,
g = c(0.2, -0.07, -0.001),
h = c(0.2, -0.07, -0.1),
theta_c = 1)
params <- get_parameters() |>
          set_seasonality(season_params = season_params, floor = 0.005) |>
          set_equilibrium(init_EIR = 5)