School of Mathematics and Statistics - Research Publications

Permanent URI for this collection

Search Results

Now showing 1 - 10 of 12
  • Item
    No Preview Available
    A multi-criteria framework for disease surveillance site selection: case study for Plasmodium knowlesi malaria in Indonesia
    Harrison, LE ; Flegg, JA ; Tobin, R ; Lubis, IND ; Noviyanti, R ; Grigg, MJ ; Shearer, FM ; Price, DJ (ROYAL SOC, 2024-01-10)
    Disease surveillance aims to collect data at different times or locations, to assist public health authorities to respond appropriately. Surveillance of the simian malaria parasite, Plasmodium knowlesi, is sparse in some endemic areas and the spatial extent of transmission is uncertain. Zoonotic transmission of Plasmodium knowlesi has been demonstrated throughout Southeast Asia and represents a major hurdle to regional malaria elimination efforts. Given an arbitrary spatial prediction of relative disease risk, we develop a flexible framework for surveillance site selection, drawing on principles from multi-criteria decision-making. To demonstrate the utility of our framework, we apply it to the case study of Plasmodium knowlesi malaria surveillance site selection in western Indonesia. We demonstrate how statistical predictions of relative disease risk can be quantitatively incorporated into public health decision-making, with specific application to active human surveillance of zoonotic malaria. This approach can be used in other contexts to extend the utility of modelling outputs.
  • Item
    Thumbnail Image
    Updating estimates of Plasmodium knowlesi malaria risk in response to changing land use patterns across Southeast Asia
    Tobin, RJ ; Harrison, LE ; Tully, MK ; Lubis, IND ; Noviyanti, R ; Anstey, NM ; Rajahram, GS ; Grigg, MJ ; Flegg, JA ; Price, DJ ; Shearer, FM ; Reithinger, R (Public Library of Science (PLoS), 2024-01)
    BACKGROUND: Plasmodium knowlesi is a zoonotic parasite that causes malaria in humans. The pathogen has a natural host reservoir in certain macaque species and is transmitted to humans via mosquitoes of the Anopheles Leucosphyrus Group. The risk of human P. knowlesi infection varies across Southeast Asia and is dependent upon environmental factors. Understanding this geographic variation in risk is important both for enabling appropriate diagnosis and treatment of the disease and for improving the planning and evaluation of malaria elimination. However, the data available on P. knowlesi occurrence are biased towards regions with greater surveillance and sampling effort. Predicting the spatial variation in risk of P. knowlesi malaria requires methods that can both incorporate environmental risk factors and account for spatial bias in detection. METHODS & RESULTS: We extend and apply an environmental niche modelling framework as implemented by a previous mapping study of P. knowlesi transmission risk which included data up to 2015. We reviewed the literature from October 2015 through to March 2020 and identified 264 new records of P. knowlesi, with a total of 524 occurrences included in the current study following consolidation with the 2015 study. The modelling framework used in the 2015 study was extended, with changes including the addition of new covariates to capture the effect of deforestation and urbanisation on P. knowlesi transmission. DISCUSSION: Our map of P. knowlesi relative transmission suitability estimates that the risk posed by the pathogen is highest in Malaysia and Indonesia, with localised areas of high risk also predicted in the Greater Mekong Subregion, The Philippines and Northeast India. These results highlight areas of priority for P. knowlesi surveillance and prospective sampling to address the challenge the disease poses to malaria elimination planning.
  • Item
    No Preview Available
    COVID-19 vaccine coverage targets to inform reopening plans in a low incidence setting
    Conway, E ; Walker, CR ; Baker, C ; Lydeamore, MJ ; Ryan, GE ; Campbell, T ; Miller, JC ; Rebuli, N ; Yeung, M ; Kabashima, G ; Geard, N ; Wood, J ; McCaw, JM ; McVernon, J ; Golding, N ; Price, DJ ; Shearer, FM (ROYAL SOC, 2023-08-30)
    Since the emergence of SARS-CoV-2 in 2019 through to mid-2021, much of the Australian population lived in a COVID-19-free environment. This followed the broadly successful implementation of a strong suppression strategy, including international border closures. With the availability of COVID-19 vaccines in early 2021, the national government sought to transition from a state of minimal incidence and strong suppression activities to one of high vaccine coverage and reduced restrictions but with still-manageable transmission. This transition is articulated in the national 're-opening' plan released in July 2021. Here, we report on the dynamic modelling study that directly informed policies within the national re-opening plan including the identification of priority age groups for vaccination, target vaccine coverage thresholds and the anticipated requirements for continued public health measures-assuming circulation of the Delta SARS-CoV-2 variant. Our findings demonstrated that adult vaccine coverage needed to be at least 60% to minimize public health and clinical impacts following the establishment of community transmission. They also supported the need for continued application of test-trace-isolate-quarantine and social measures during the vaccine roll-out phase and beyond.
  • Item
    No Preview Available
    Forecasting COVID-19 activity in Australia to support pandemic response: May to October 2020
    Moss, R ; Price, DJ ; Golding, N ; Dawson, P ; McVernon, J ; Hyndman, RJ ; Shearer, FM ; McCaw, JM (NATURE PORTFOLIO, 2023-05-30)
    As of January 2021, Australia had effectively controlled local transmission of COVID-19 despite a steady influx of imported cases and several local, but contained, outbreaks in 2020. Throughout 2020, state and territory public health responses were informed by weekly situational reports that included an ensemble forecast of daily COVID-19 cases for each jurisdiction. We present here an analysis of one forecasting model included in this ensemble across the variety of scenarios experienced by each jurisdiction from May to October 2020. We examine how successfully the forecasts characterised future case incidence, subject to variations in data timeliness and completeness, showcase how we adapted these forecasts to support decisions of public health priority in rapidly-evolving situations, evaluate the impact of key model features on forecast skill, and demonstrate how to assess forecast skill in real-time before the ground truth is known. Conditioning the model on the most recent, but incomplete, data improved the forecast skill, emphasising the importance of developing strong quantitative models of surveillance system characteristics, such as ascertainment delay distributions. Forecast skill was highest when there were at least 10 reported cases per day, the circumstances in which authorities were most in need of forecasts to aid in planning and response.
  • Item
    No Preview Available
    A model for malaria treatment evaluation in the presence of multiple species
    Walker, CR ; Hickson, RI ; Chang, E ; Ngor, P ; Sovannaroth, S ; Simpson, JA ; Price, DJ ; McCaw, JM ; Price, RN ; Flegg, JA ; Devine, A (ELSEVIER, 2023-09)
    Plasmodium falciparum and P. vivax are the two most common causes of malaria. While the majority of deaths and severe morbidity are due to P. falciparum, P. vivax poses a greater challenge to eliminating malaria outside of Africa due to its ability to form latent liver stage parasites (hypnozoites), which can cause relapsing episodes within an individual patient. In areas where P. falciparum and P. vivax are co-endemic, individuals can carry parasites of both species simultaneously. These mixed infections complicate dynamics in several ways: treatment of mixed infections will simultaneously affect both species, P. falciparum can mask the detection of P. vivax, and it has been hypothesised that clearing P. falciparum may trigger a relapse of dormant P. vivax. When mixed infections are treated for only blood-stage parasites, patients are at risk of relapse infections due to P. vivax hypnozoites. We present a stochastic mathematical model that captures interactions between P. falciparum and P. vivax, and incorporates both standard schizonticidal treatment (which targets blood-stage parasites) and radical cure treatment (which additionally targets liver-stage parasites). We apply this model via a hypothetical simulation study to assess the implications of different treatment coverages of radical cure for mixed and P. vivax infections and a "unified radical cure" treatment strategy where P. falciparum, P. vivax, and mixed infections all receive radical cure after screening glucose-6-phosphate dehydrogenase (G6PD) normal. In addition, we investigated the impact of mass drug administration (MDA) of blood-stage treatment. We find that a unified radical cure strategy leads to a substantially lower incidence of malaria cases and deaths overall. MDA with schizonticidal treatment was found to decrease P. falciparum with little effect on P. vivax. We perform a univariate sensitivity analysis to highlight important model parameters.
  • Item
    No Preview Available
    Optimal Interruption of P. vivax Malaria Transmission Using Mass Drug Administration
    Anwar, MN ; Hickson, RI ; Mehra, S ; Price, DJ ; McCaw, JM ; Flegg, MB ; Flegg, JA (SPRINGER, 2023-06)
    Plasmodium vivax is the most geographically widespread malaria-causing parasite resulting in significant associated global morbidity and mortality. One of the factors driving this widespread phenomenon is the ability of the parasites to remain dormant in the liver. Known as 'hypnozoites', they reside in the liver following an initial exposure, before activating later to cause further infections, referred to as 'relapses'. As around 79-96% of infections are attributed to relapses from activating hypnozoites, we expect it will be highly impactful to apply treatment to target the hypnozoite reservoir (i.e. the collection of dormant parasites) to eliminate P. vivax. Treatment with radical cure, for example tafenoquine or primaquine, to target the hypnozoite reservoir is a potential tool to control and/or eliminate P. vivax. We have developed a deterministic multiscale mathematical model as a system of integro-differential equations that captures the complex dynamics of P. vivax hypnozoites and the effect of hypnozoite relapse on disease transmission. Here, we use our multiscale model to study the anticipated effect of radical cure treatment administered via a mass drug administration (MDA) program. We implement multiple rounds of MDA with a fixed interval between rounds, starting from different steady-state disease prevalences. We then construct an optimisation model with three different objective functions motivated on a public health basis to obtain the optimal MDA interval. We also incorporate mosquito seasonality in our model to study its effect on the optimal treatment regime. We find that the effect of MDA interventions is temporary and depends on the pre-intervention disease prevalence (and choice of model parameters) as well as the number of MDA rounds under consideration. The optimal interval between MDA rounds also depends on the objective (combinations of expected intervention outcomes). We find radical cure alone may not be enough to lead to P. vivax elimination under our mathematical model (and choice of model parameters) since the prevalence of infection eventually returns to pre-MDA levels.
  • Item
    Thumbnail Image
    Real-time analysis of hospital length of stay in a mixed SARS-CoV-2 Omicron and Delta epidemic in New South Wales, Australia
    Tobin, RJJ ; Wood, JGG ; Jayasundara, D ; Sara, G ; Walker, CRR ; Martin, GEE ; McCaw, JMM ; Shearer, FMM ; Price, DJJ (BMC, 2023-01-17)
    BACKGROUND: The distribution of the duration that clinical cases of COVID-19 occupy hospital beds (the 'length of stay') is a key factor in determining how incident caseloads translate into health system burden. Robust estimation of length of stay in real-time requires the use of survival methods that can account for right-censoring induced by yet unobserved events in patient progression (e.g. discharge, death). In this study, we estimate in real-time the length of stay distributions of hospitalised COVID-19 cases in New South Wales, Australia, comparing estimates between a period where Delta was the dominant variant and a subsequent period where Omicron was dominant. METHODS: Using data on the hospital stays of 19,574 individuals who tested positive to COVID-19 prior to admission, we performed a competing-risk survival analysis of COVID-19 clinical progression. RESULTS: During the mixed Omicron-Delta epidemic, we found that the mean length of stay for individuals who were discharged directly from ward without an ICU stay was, for age groups 0-39, 40-69 and 70 +, respectively, 2.16 (95% CI: 2.12-2.21), 3.93 (95% CI: 3.78-4.07) and 7.61 days (95% CI: 7.31-8.01), compared to 3.60 (95% CI: 3.48-3.81), 5.78 (95% CI: 5.59-5.99) and 12.31 days (95% CI: 11.75-12.95) across the preceding Delta epidemic (1 July 2021-15 December 2021). We also considered data on the stays of individuals within the Hunter New England Local Health District, where it was reported that Omicron was the only circulating variant, and found mean ward-to-discharge length of stays of 2.05 (95% CI: 1.80-2.30), 2.92 (95% CI: 2.50-3.67) and 6.02 days (95% CI: 4.91-7.01) for the same age groups. CONCLUSIONS: Hospital length of stay was substantially reduced across all clinical pathways during a mixed Omicron-Delta epidemic compared to a prior Delta epidemic, contributing to a lessened health system burden despite a greatly increased infection burden. Our results demonstrate the utility of survival analysis in producing real-time estimates of hospital length of stay for assisting in situational assessment and planning of the COVID-19 response.
  • Item
    Thumbnail Image
    A modelling approach to estimate the transmissibility of SARS-CoV-2 during periods of high, low, and zero case incidence
    Golding, N ; Price, DJ ; Ryan, G ; McVernon, J ; McCaw, JM ; Shearer, FM (eLIFE SCIENCES PUBL LTD, 2023-01-20)
    Against a backdrop of widespread global transmission, a number of countries have successfully brought large outbreaks of COVID-19 under control and maintained near-elimination status. A key element of epidemic response is the tracking of disease transmissibility in near real-time. During major outbreaks, the effective reproduction number can be estimated from a time-series of case, hospitalisation or death counts. In low or zero incidence settings, knowing the potential for the virus to spread is a response priority. Absence of case data means that this potential cannot be estimated directly. We present a semi-mechanistic modelling framework that draws on time-series of both behavioural data and case data (when disease activity is present) to estimate the transmissibility of SARS-CoV-2 from periods of high to low - or zero - case incidence, with a coherent transition in interpretation across the changing epidemiological situations. Of note, during periods of epidemic activity, our analysis recovers the effective reproduction number, while during periods of low - or zero - case incidence, it provides an estimate of transmission risk. This enables tracking and planning of progress towards the control of large outbreaks, maintenance of virus suppression, and monitoring the risk posed by re-introduction of the virus. We demonstrate the value of our methods by reporting on their use throughout 2020 in Australia, where they have become a central component of the national COVID-19 response.
  • Item
    Thumbnail Image
    COVID-19 in low-tolerance border quarantine systems: Impact of the Delta variant of SARS-CoV-2
    Zachreson, C ; Shearer, FM ; Price, DJ ; Lydeamore, MJ ; McVernon, J ; McCaw, J ; Geard, N (AMER ASSOC ADVANCEMENT SCIENCE, 2022-04)
    In controlling transmission of coronavirus disease 2019 (COVID-19), the effectiveness of border quarantine strategies is a key concern for jurisdictions in which the local prevalence of disease and immunity is low. In settings like this such as China, Australia, and New Zealand, rare outbreak events can lead to escalating epidemics and trigger the imposition of large-scale lockdown policies. Here, we develop and apply an individual-based model of COVID-19 to simulate case importation from managed quarantine under various vaccination scenarios. We then use the output of the individual-based model as input to a branching process model to assess community transmission risk. For parameters corresponding to the Delta variant, our results demonstrate that vaccination effectively counteracts the pathogen's increased infectiousness. To prevent outbreaks, heightened vaccination in border quarantine systems must be combined with mass vaccination. The ultimate success of these programs will depend sensitively on the efficacy of vaccines against viral transmission.
  • Item
    No Preview Available
    Development of an influenza pandemic decision support tool linking situational analytics to national response policy.
    Shearer, FM ; Moss, R ; Price, DJ ; Zarebski, AE ; Ballard, PG ; McVernon, J ; Ross, JV ; McCaw, JM (Elsevier, 2021-06-19)
    National influenza pandemic plans have evolved substantially over recent decades, as has the scientific research that underpins the advice contained within them. While the knowledge generated by many research activities has been directly incorporated into the current generation of pandemic plans, scientists and policymakers are yet to capitalise fully on the potential for near real-time analytics to formally contribute to epidemic decision-making. Theoretical studies demonstrate that it is now possible to make robust estimates of pandemic impact in the earliest stages of a pandemic using first few hundred household cohort (FFX) studies and algorithms designed specifically for analysing FFX data. Pandemic plans already recognise the importance of both situational awareness i.e., knowing pandemic impact and its key drivers, and the need for pandemic special studies and related analytic methods for estimating these drivers. An important next step is considering how information from these situational assessment activities can be integrated into the decision-making processes articulated in pandemic planning documents. Here we introduce a decision support tool that directly uses outputs from FFX algorithms to present recommendations on response options, including a quantification of uncertainty, to decision makers. We illustrate this approach using response information from within the Australian influenza pandemic plan.