Studying the recovery procedure for the time-dependent transmission rate(s) in epidemic models

Determining the time-dependent transmission function that exactly reproduces disease incidence data can yield useful information about disease outbreaks, including a range potential values for the recovery rate of the disease and could offer a method to test the "school year" hypothesis (seasonality) for disease transmission. Recently two procedures have been developed to recover the time-dependent transmission function, β(t), for classical disease models given the disease incidence data. We first review the β(t) recovery procedures and give the resulting formulas, using both methods, for the susceptible-infected-recovered (SIR) and susceptible-exposed-infected-recovered (SEIR) models. We present a modification of one procedure, which is then shown to be identical to the other. Second, we explore several technical issues that appear when implementing the procedure for the SIR model; these are important when generating the time-dependent transmission function for real-world disease data. Third, we extend the recovery method to heterogeneous populations modeled with a certain SIR-type model with multiple time-dependent transmission functions. Finally, we apply the β(t) recovery procedure to data from the 2002-2003 influenza season and for the six seasons from 2002-2003 through 2007-2008, for both one population class and for two age classes. We discuss the consequences of the technical conditions of the procedure applied to the influenza data. We show that the method is robust in the heterogeneous cases, producing comparable results under two different hypotheses. We perform a frequency analysis, which shows a dominant 1-year period for the multi-year influenza transmission function(s).