9. Ro formulation using the Next Generation Matrix (Diekmann et al. 2010)