For two of the four special cases (a = '), and for the case considered by Müller et al [7] (a = 1, t inf = ', t lat = 0), Y(p c ) can be calculated numerically in MathematicaH. For all other cases, including the real infections, stochastic simulations were needed.