Measurability of the epidemic reproduction number in data-driven contact networks

Significance The analysis of real epidemiological data has raised issues of the adequacy of the classic homogeneous modeling framework and quantities, such as the basic reproduction number in real-world situations. Based on high-quality sociodemographic data, here we generate a multiplex network describing the contact pattern of the Italian and Dutch populations. By using a microsimulation approach, we show that, for epidemics spreading on realistic contact networks, it is not possible to define a steady exponential growth phase and a basic reproduction number. We show the operational use of the instantaneous reproduction rate as a good descriptor of the transmission dynamics. The basic reproduction number is one of the conceptual cornerstones of mathematical epidemiology. Its classical definition as the number of secondary cases generated by a typical infected individual in a fully susceptible population finds a clear analytical expression in homogeneous and stratified mixing models. Along with the generation time (the interval between primary and secondary cases), the reproduction number allows for the characterization of the dynamics of an epidemic. A clear-cut theoretical picture, however, is hardly found in real data. Here, we infer from highly detailed sociodemographic data two multiplex contact networks representative of a subset of the Italian and Dutch populations. We then simulate an infection transmission process on these networks accounting for the natural history of influenza and calibrated on empirical epidemiological data. We explicitly measure the reproduction number and generation time, recording all individual-level transmission events. We find that the classical concept of the basic reproduction number is untenable in realistic populations, and it does not provide any conceptual understanding of the epidemic evolution. This departure from the classical theoretical picture is not due to behavioral changes and other exogenous epidemiological determinants. Rather, it can be simply explained by the (clustered) contact structure of the population. Finally, we provide evidence that methodologies aimed at estimating the instantaneous reproduction number can operationally be used to characterize the correct epidemic dynamics from incidence data.