Statistical estimation of parameters in a disease transmission model: analysis of a Cryptosporidium outbreak

Population dynamic models, commonly used tools in the study of epidemics and other complex population processes, are implicit non-linear mathematical equations. Inference based on such models can be difficult due to the problems associated with high dimensional parameters that may be non-identified and complex likelihood functions that are difficult to maximize. To address a problem of non-identifiability due to collinearity of parameter estimates in a mathematical model of the 1993 Milwaukee Cryptosporidium parvum outbreak, we examined the utility of a constrained profile likelihood approach. This method was used to study two parameters of interest from the mathematical model: (i). the rate of secondary transmission; (ii). the proportional increase in primary transmission due to water treatment failure. The estimated values of these parameters were shown to depend strongly on poorly understood aspects of Cryptosporidium epidemiology such as asymptomatic proportion and the population immune status. Our analysis demonstrated that the combination of a disease transmission model and a constrained profile likelihood procedure provides an effective approach for inference and estimation of important parameters regulating infectious disease outbreaks.

[1]  Estimation of the contact rate in a seasonal SEIR model: application to chickenpox incidence in France. , 2000, Statistics in medicine.

[2]  Adrian E. Raftery,et al.  Inference from a Deterministic Population Dynamics Model for Bowhead Whales , 1995 .

[3]  M. V. Bhat,et al.  An Efficient Clustering Algorithm , 1976, IEEE Transactions on Systems, Man, and Cybernetics.

[4]  W. M. Mac Kenzie,et al.  Cryptosporidium parvum-specific antibody responses among children residing in Milwaukee during the 1993 waterborne outbreak. , 2001, The Journal of infectious diseases.

[5]  C. Rogier,et al.  Bayesian analysis of an epidemiologic model of Plasmodium falciparum malaria infection in Ndiop, Senegal. , 2000, American journal of epidemiology.

[6]  J. P. Davis,et al.  A massive outbreak in Milwaukee of cryptosporidium infection transmitted through the public water supply. , 1994, The New England journal of medicine.

[7]  Elena N. Naumova,et al.  Did milwaukee experience waterborne cryptosporidiosis before the large documented outbreak in 1993 , 1998 .

[8]  R C Spear,et al.  Dynamic model comparing the bionomics of two isolated Culex tarsalis (Diptera: Culicidae) populations: sensitivity analysis. , 1995, Journal of medical entomology.

[9]  H. Dupont,et al.  Cryptosporidium parvum: intensity of infection and oocyst excretion patterns in healthy volunteers. , 1996, The Journal of infectious diseases.

[10]  D E Peterson,et al.  Massive outbreak of waterborne cryptosporidium infection in Milwaukee, Wisconsin: recurrence of illness and risk of secondary transmission. , 1995, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[11]  R. Spear Eutrophication in peel inlet—II. Identification of critical uncertainties via generalized sensitivity analysis , 1980 .

[12]  R C Spear,et al.  An Analysis of the Milwaukee Cryptosporidiosis Outbreak Based on a Dynamic Model of the Infection Process , 1998, Epidemiology.

[13]  Thomas J. Manetsch,et al.  Time-Varying Distributed Delays and Their Use in Aggregative Models of Large Systems , 1976, IEEE Transactions on Systems, Man, and Cybernetics.

[14]  H. Dupont,et al.  Infectivity of Cryptosporidium parvum in healthy adults with pre-existing anti-C. parvum serum immunoglobulin G. , 1999, The American journal of tropical medicine and hygiene.