Abstract
We study infectious disease outbreaks that have evolved in isolation without the influence of one another. If stochastic effects are identified within each outbreak, it is necessary to model the dynamics with stochastic epidemic models. However, the accuracy of the estimated model parameters depends on several factors including the statistical inference methodologies that are used. One approach to making inferences from multiple outbreak data is the use of a Bayesian hierarchical model. This statistical framework allows simultaneous inference for multiple outbreaks and the estimation of model parameters at a group level. A hierarchical model will generally provide improved estimates; however, we show that this is not always true when the variability among model parameter values of the outbreaks is high. We further show that subsets of outbreaks with similar parameter values can be identified prior to a hierarchical analysis using common clustering algorithms such as k-means. When hierarchical analyses are carried out for these pre-identified subsets of outbreaks, parameter estimates are improved compared to those estimated under a hierarchical analysis for the complete set of outbreaks. We have applied our estimation framework within a simulation-based experiment using synthetic data generated from stochastic SIRS models. The framework is generalizable to other biological data.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
Punya Alahakoon is supported by a Melbourne Research Scholarship from the University of Melbourne. P.G. Taylor would like to acknowledge the support of the Australian Research Council via the Centre of Excellence for Mathematical and Statistical Frontiers (ACEMS).
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
All data produced are available online at https://github.com/PunyaAlahakoon/Use_of_Clusetring_in_hierarchical_models.git
https://github.com/PunyaAlahakoon/Use_of_Clusetring_in_hierarchical_models.git