Stochastic multi-type SIR epidemics among a population partitioned into households

We consider a stochastic model for the spread of an SIR (susceptible --> infective --> removed) epidemic among a closed, finite population that contains several types of individuals and is partitioned into households. The infection rate between two individuals depends on the types of the transmitting and receiving individuals and also on whether the infection is local (i.e., within a household) or global (i.e., between households). The exact distribution of the final outcome of the epidemic is outlined. A branching process approximation for the early stages of the epidemic is described and made fully rigorous, by considering a sequence of epidemics in which the number of households tends to infinity and using a coupling argument. This leads to a threshold theorem for the epidemic model. A central limit theorem for the final outcome of epidemics which take off is derived, by exploiting an embedding representation.