It is now possible to extend the
formulation developed for the riskstructured SIS models to an
agestructured SIR model. This involves three major changes:
1) Given that we are now dealling with an SIR model, (at least) two
compartments need to be modelled per age group. (In fact we
should model S, I and R in each agegroup and then use
the fact that the sum of all the population compartments equals 1 to
eliminate one of the variables. However, by assuming that the size each
of the ageclasses (n_{C} and n_{A}) is at equilbrium,
we can simplify the dynamics and just consider two variable per
agegroup.)
2) Demography needs to be included, this essentially means that
individuals need to be born into the youngest agegroup. In addtion,
natural mortality can be agestructured.
3) Finally, individuals move through the ageclasses as they mature. In
the two ageclass example below individuals leave the childhood class
and move into the adult class, taking their infectious state with them.
In this set of programs we consider the simplest case where just two
agegroups (children and adults) are models. Again we use a
frequentist, such that S_{C}, I_{C} and R_{C}
refer to the proportion of the entire population that are
susceptible, infectious or recovered and also in the childhood
ageclass. Putting these elements together, we
arrive at the following differential equations:
Parameters
β 
is the matrix of transmission
rates and incorporates the encounter rate between susceptible and
infectious individuals together with the probability of transmission. 
γ 
is
called the removal
or recovery rate, though often we are more interested in its reciprocal
(1/γ) which determines the average infectious period. 
l_{C} 
is
the rate at which children mature and move into the adult age class

n_{C} 
is
the
proportion of the population that are in the childhood group. This is
not defined explicitly but calculated as: n_{C} = μ_{A}/(l_{C}+μ_{A})

μ_{C} 
is
the death rate in the childhood group

μ_{A} 
is
the death rate in the adult group

ν 
is
the birth rate into the childhood class. This is not defined explicitly
but calculated as: ν = (l_{C}+μ_{A}) n_{C} 
S_{C}(0) 
is
the initial
proportion of the population that are both susceptible and in the
childhood group. 
I_{C}(0) 
is
the initial
proportion of the population that are both infectious and in the
childhood group.

S_{A}(0) 
is
the initial
proportion of the population that are both susceptible and in the adult
group 
I_{A}(0) 
is
the initial
proportion of the population that are both infectious and in the adult
group. 
All rates are
specified
in years.
Requirements.
All parameters must be positive,
S_{C}(0) + I_{C}(0)≤
n_{C}, and S_{A}(0)
+ I_{A}(0)≤
n_{A }= 1 n_{C}_{
}
Files
C++ Program, Python Program, Fortran Program, Parameters, MATLAB Code.
