We integrated this assumption into our estimations, and fitted to the entire dataset model M DT a , with heterogeneous time scales for the parameters (cf. model Table 1 and Methods), where a was kept constant for the whole epidemic history, while the rates b t and e t changed with frequency DT.