Developing next generation matrices from existing data
13th May 2020
The estimates of transmission rates in each sector (i.e., home, school, work and other locations) and the overall reproduction number are derived from baseline estimates of the daily, age-specific contact rates between individuals of different age groups. These contact rates are provided by the analysis in Prem et al.(1) where data from population-based contact diaries in eight European countries were projected to generate contact intensities for 144 other countries using Bayesian modelling techniques. The inferred values ð‘ð‘Žð‘Ž′ give the number of (pre-COVID-19) typical daily contacts an individual of age ð‘Ž′ makes with an individual of age ð‘Ž. In the dataset, age bands are separated into 5 year age groups and contacts are further divided into four locations: work, home, school and other.
To estimate the transmission capacity associated with these contacts we convert the contact intensity matrices to next-generation matrices, ð¾, whose elements, ð‘˜ð‘Žð‘Ž′, give the number of new infections of age ð‘Ž generated by individuals of age ð‘Ž′. As a first step, we compute an unscaled next-generation matrix ð¾Ì… by weighting the elements of the contact matrix ð‘ð‘Žð‘Ž′ by the age-dependent relative susceptibility (ðœŽð‘Ž) and infectivity (ð›½ð‘Ž) of individuals in the population and the distribution of susceptible (ð‘ ð‘Ž) and total (ð‘›ð‘Ž) individuals in each age group. In particular, the elements, ð‘˜Ì…
ð‘Žð‘Ž′, of the unscaled next-generation matrix (NGM), ð¾Ì…, are given by
ð‘˜Ì…ð‘Žð‘Ž′=ðœŽð‘Žð‘ ð‘Žð‘ð‘Žð‘Ž′ð›½ð‘Ž′ð‘›ð‘Ž′.
Here ðœŽð‘Ž is the relative susceptibility to infection for an individual in age group ð‘Ž and ð›½ð‘Ž is their corresponding transmissibility once infected. Since the population is entirely susceptible upon first introduction of the infection such that ð‘ ð‘Ž=ð‘›ð‘Ž.
For symmetry, we assume that the age-dependent susceptibility and transmissibility profiles are equal equivalent, i.e., ðœŽð‘Ž=ð›½ð‘Ž, and are given by the following parametric equation:
ðœŽð‘Ž=1−ðœŽrel2tanh(ð‘(ð‘Ž−ð‘))+1+ðœŽrel2
where ðœŽrel is approximately equal to the relative susceptibility between individuals in the youngest (<5) and those in the oldest (>80) age groups. In the following analysis we assume baseline values of ðœŽmin = 0.1, ð‘ = 0.3 and ð‘ = 27.
We choose values to match the proportion of each age group infected in China (the country used to calibrate the model) and then applied the calibrated values to Australian mixing matrices.
Back to top