Estimating the within-household infection rate in emerging SIR epidemics among a community of households

This paper is concerned with estimation of the within-household infection rate λL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _L$$\end{document} for a susceptible →\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rightarrow $$\end{document} infective →\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rightarrow $$\end{document} recovered epidemic among a population of households, from observation of the early, exponentially growing phase of an epidemic. Specifically, it is assumed that an estimate of the exponential growth rate is available from general data on an emerging epidemic and more-detailed, household-level data are available in a sample of households. Estimates of λL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _L$$\end{document} obtained using the final size distribution of single-household epidemics are usually biased owing to the emerging nature of the epidemic. A new method, which accounts correctly for the emerging nature of the epidemic, is developed by exploiting the asymptotic theory of supercritical branching processes and proved to yield a strongly consistent estimator of λL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _L$$\end{document} as the population and sampled households both tend to infinity in an appropriate fashion. The theory is illustrated by simulations which demonstrate that the new method is feasible for finite populations and numerical studies are used to explore how changes to the parameters governing the spread of an epidemic affect the bias of estimates based on single-household final size distributions.