Calculate a contact matrix for age groups based on Polymod contact survey data

Usage

contactMatrixPolymod(agelims, agepops = NULL)

Arguments

agelims: minimum age in years for each age group. The maximum valid age limit is 90, as the socialmixr contact_matrix function supports ages up to 90. Age limits greater than 90 will be replaced with 90 and a warning will be issued.
agepops: population size of each group, defaulting to demography of Polymod survey population. If provided, must match the length of the age groups defined by agelims (after any adjustments for exceeding the 90-year limit).

Value

A symmetric contact matrix with row and column names indicating the age groups.

Details

The socialmixr contact_matrix function supports age limits up to 90. Any age limits above 90 will be adjusted to 90 with a warning, and the corresponding populations will be aggregated into a single "90+" group.

Examples

#Default population distribution uses population data from POLYMOD survey locations:
contactMatrixPolymod(agelims = c(0, 5, 18))
#>          contact.age.group
#> age.group    under5    5to17       18+
#>    under5 2.2370031 1.504903  6.086132
#>    5to17  0.5148409 9.846792  7.272662
#>    18+    0.3754707 1.311486 10.828758
#Specifying the age distribution will lead to an adjusted version:
contactMatrixPolymod(agelims = c(0, 5, 18), agepops = c(500, 1300, 8200))
#>          contact.age.group
#> age.group    under5    5to17       18+
#>    under5 2.2518104 1.347447  6.198479
#>    5to17  0.5182488 8.816532  7.406912
#>    18+    0.3779560 1.174266 11.028652