Projecting social contact matrices to populations stratified by binary attributes with known homophily

Contact networks are heterogeneous. People with similar characteristics are more likely to interact, a phenomenon called assortative mixing or homophily. Empirical age-stratified social contact matrices have been derived by extensive survey work. We lack however similar empirical studies that provide social contact matrices for a population stratified by attributes beyond age, such as gender, sexual orientation, or ethnicity. Accounting for heterogeneities with respect to these attributes can have a profound effect on model dynamics. Here, we introduce a new method, which uses linear algebra and non-linear optimization, to expand a given contact matrix to populations stratified by binary attributes with a known level of homophily. Using a standard epidemiological model, we highlight the effect homophily can have on model dynamics, and conclude by briefly describing more complicated extensions. The available Python source code enables any modeler to account for the presence of homophily with respect to binary attributes in contact patterns, ultimately yielding more accurate predictive models.