Hello from Australia!

  • MUA Department of Econometrics and Business Statistics is the largest general statistics group in Victoria
  • Business Analytics is broad-ranging, including fisheries, data visualisation, extreme weather events, and more theory-based research.

My roadmap

  • Bachelor of Mathematical Sciences at the University of Adelaide (completed 2012)
  • Masters of Philosophy at the University of Adelaide (completed 2015)
  • PhD at the University of Melbourne (completed 2019)
  • Joined Monash University in mid-2019
  • Including honorary position at the SaferCare Victoria, Department of Health & Human Services
  • Recruited into the COVID-19 response in late February/early March
  • Appointed “Modelling & Forecasting Lead” in late March
  • Appointed Manager of Analytics in September
  • Rejoined Monash in 2021, including the Doherty COVID-19 Modelling Consortium
  • Joined EBS in September 2021

Burden of Healthcare Acquired Infections in Australia

Coauthors: Brett Mitchell, Tracey Bucknall, Allen Cheng, Phil Russo, Andrew Stewardson

Healthcare Associated Infections

  • Healthcare associated infections (HAIs) associated with increased morbidity and mortality
  • Five of the most common HAIs are:
    • Clostridiodes difficile – causes severe damage to the colon, can be fatal
    • Bloodstream infection (sepsis) – estimated mortality rate of 15-30%
    • Urinary track infection – low mortality but associated with multi-drug resistance and significantly longer hospital stays
    • Healthcare acquired pneumonia – mortality rate of 40-70%, increasing dramatically with age
    • Surgical site infection – significantly increases length of stay

In Australia, HAIs are not notifiable => We have no robust way to track whether their prevalence is increasing or decreasing

Europe tracks these closely

HAIs are actively monitored across Europe through the ECDC

  • In 2016 (based on 2012 data), 2,609,911 new HAIs are estimated to have occurred.

The data were obtained in a point prevalence survey on an enormous scale

  • 273, 753 patients
  • 1,149 hospitals

Point prevalence survey

A point prevalence survey counts the number of people with a condition on a given day

For the Australian PPS:

  • Adults in 19 large acute care public hospitals were sampled
  • All acute care wards were include, non-acute, paediatric, neonatal ICUs, rehab and emergency departments were excluded.

The hospitals sampled make up approximately 60% of all overnight separations in Australia

Point prevalence survey

  • 2767 patients were sampled between 6 Aug and 29 Nov 2018
  • Median age: 67 (range 18-104)
  • 52.9% male, 46.6% female, 0.5% unknown/other
  • 85.7% patients in major city hospitals

Estimation methodology

Step 1: Hospital prevalence

Hospital prevalence, P, estimated as

P = r \times Beta(n_{\text{obs}}, N-n_{\text{obs}}+1) + (1-r) \times Beta(n_{\text{obs}}+1, N-n_{\text{obs}})

Estimation methodology

Step 2: Estimation of hospital incidence

Hospital incidence, I, calculated as

I = P \frac{\text{LA}}{\text{LOI}}

where

  • P is the hospital prevalence from step 1,
  • \text{LA} is the mean length of stay and
  • \text{LOI} is the length of infection.

\text{LA} is taken from the AIHW 2018 statistics on all public hospitals.

Estimation methodology

Step 2a: Estimation of length of infection

No data on length of infection, only LOI_\text{pps}, the length of stay until the date of survey.

We can calculate

P(LOI_\text{pps} = 1),

the probability a patient is in the first day of their HAI. Then, approximate LOI with

E[LOI] = 1 / P(LOI_\text{pps} = 1).

For small sample sizes, results has shown this is biased, so we use a mixture of this estimator and the empirical mean.

Estimation methodology

Step 3: Estimation of population incidence

Calculate population incidence as

I_{\text{pop}} = I \times N_{\text{discharges}}.

For us, N_{\text{discharges}} = 3,713.513, which is 60% of the total admissions in the year.

This is one of the few pieces of data Australia has that the ECDC doesn’t!

Estimation methodology

Step 4: Stratification by age and sex

Use a multinomial likelihood with a Dirichlet prior, with weights taken from the number of cases in each age/sex category.

A psuedocount is added to each strata (0.001\sum \text{weights}) to ensure likelihood can be calculated with empty strata

This psuedocount almost surely induces bias, and there are better techniques out there.

Estimation methodology

Step 5: Adjustment for life expectancy

Use the McCabe score, which gives the life expectancy according to severity of disease. Patients are categorised as:

  • non-fatal
  • fatal (life expectancy 3 years)
  • rapidly fatal (average life expectancy of 0.5 years)

These scores, combined with disease outcome trees, give DALYs and deaths.

Disease outcome trees

Key results

Number of HAIs
(95% CI)
Deaths
(95% CI)
DALYs
(95% CI)
SSI 44,238
(31,176 - 73,797)
876
(617 - 1,263)
13,197
(9,298 - 19,001)
UTI 42,408
(25,200 - 68,735)
729
(259 - 1,772)
16,087
(5,939 - 37,218)
CDI 5,125
(2,360 - 10,740)
262
(13 - 836)
2,757
(241 - 8,655)
HAP 51,499
(31,343 - 82,877)
1,904
(462 - 4,430)
39,276
(17,608 - 77,915)
BSI 23,979
(15,658 - 36,245)
3,512
(1,874 - 6,075)
46,773
(26,205 - 79,104)
All 170,574
(135,779 - 213,898)
7,583
(4,941 - 11,135)
122,376
(85,136 - 172,784)

That’s 1 in 20 admissions resulting in an avoidable infection!

Key results

Key results

Key results

Novelty

  • First estimate of HAI burden in Australia using (relatively) robust survey data in an established framework

  • Based on first point prevalence survey since 1984

  • There is no routine surveillance of HAIs in Australia

  • Point prevalence surveys remain the only way to understand the burden of these conditions

Summary

  • 498 DALYs per 100,000 is a large amount
    • Motor vehicles: 180 DALYs
    • Infectious diseases: 370 DALYs
    • Respiratory diseases: 1380 DALYs

This work has informed guidance on HAI surveillance in Australia, including new funding schemes to better understand these conditions.

And all this based on just 2767 patients from 19 hospitals…

HAIs are largely preventable.

They represent great opportunity for improvement, and we have a long way to go to prevent them entirely.

Carbapenemase-producing enterobacterales

D. Wu; T. Donker; B. Cooper; M. Easton; N. Geard; C. Gorrie; D. Hennessy; B. Howden; A. Peleg; A. Turner; A. Wilson; A. Stewardson on behalf of the ECHIDNA study group

Carbanememase-producing Enterobacterales

  • Acquired primarily in hospital
  • Resistant to last-line antibiotics
  • Mortality, given clinical symptoms, in the order of 30%
  • Carriage for >= 12 months
  • Testing via rectal swab
  • Can carry with no symptoms
  • Unknown if asymptomatic people can transmit

Incidence in Victoria, Australia

Notification system

Notifications are required to be sent to the Department of Health within 24 hours of a positive result.

However, this is a one-way system. Hospitals can not look up a patient to see if they have tested positive previously at a different health service.

Notification system

Hospital outbreaks

The immediate question:

When an outbreak is detected, what other hospitals could be affected?

Natural to consider patients in a hospital as a network, that has the nodes as the hospitals and edges between them related to the probability of transfer.

Current methods can grossly over-estimate the rate of spread of disease between hospitals

Current state of the field

Most methods are based on line-listed admissions data:

patient_id admission_date separation_date location diagnosis_codes
1 2020-01-01 2020-01-02 H1 A B C
2 2020-01-01 2020-01-02 H2 A B C
1 2020-01-03 2020-01-04 H2 A B C

This is interpreted as one transfer from H1 to H2, with a one-day gap.

Our dataset is also linked to notifiable diseases, genomic information and some ward-level information.

Data summary

The standard network

The standard network on a map

The standard network on a map

But movement changes all the time

Different ways to collapse the data

Some notation:

  • u, v are locations
  • s, t are times, t \geq s
  • Hazard for patients returning home, never readmitted: \zeta_u(t)
  • w_{uv}(s,t) patients discharged from u at s, then readmitted to v at t

Different ways to collapse the data

Different ways to collapse the data

Then, we have the following set of governing equations:

  • Hospital to home: \lambda(u(s) \to \emptyset) = \zeta_u(s)
  • Direct transfers: \lambda(u(s) \to v(t)) = d_{uv}(s, t)
  • Indirect transfers:
    • Hospital to home: \lambda(u(s) \to z_{uv}(t)) = \eta_{uv}(s, t)
    • Home to hospital: \lambda(z_{uv}(t) \to v(t)) = \rho_{uv}(s, t)

Model 1: Naive Static Method

d_{uv}(s,t) = \frac{\sum_{s,t} w_{uv}(s,t)}{T_\Sigma N_u}

Tip

All movements are instanteous, and the rate of movement is the mean over the entire observation period.

Model 2: Improved Static Method

Choose a threshold value, \omega, to separate direct and indirect transfers

d_{uv}(s,t) = d_{uv} = \frac{\sum_{s,t: (t-s) < \omega} w_{uv}(s,t)}{T_\Sigma N_u} \eta_{uv}(s,t) = \eta_{uv} = \frac{\sum_{s,t: (t-s) \geq \omega} w_{uv}(s,t)}{T_\Sigma N_u} \rho_{uv}(s,t) = \rho_{uv} = \left[ \frac{\sum_{s,t: (t-s) \geq \omega} (t-s) w_{uv}(s,t)}{\sum_{s,t: (t-s) \geq \omega} w_{uv}(s,t)} \right]^{-1}

Model 3: Snapshot Model

Choose a snapshot duration \omega. This defines the threshold duration of an indirect transfer.

\begin{aligned} \lambda(u(s) \to {z'}_{uv}(t)) &= \eta_{uv}(s, t)\\ \lambda({z'}_{uv}(t) \to z_{uv}(t)) &= \delta(t \; \mathrm{ mod } \; \omega)\\ \end{aligned}

z' contains the individuals that would enter z within the duration of a given snapshot [t, t+\omega], so that they do not immediately readmit to their next hospital.

Model 4: Temporal Model

Effecitvely, choose a very small snapshot window, \omega. Individuals are not re-admitted to z', instead they are re-admitted to v uniformly between t and t+\omega. This gives,

d_{uv}(s, t) = \frac{\sum_{t: (t-s) < \omega} w_{uv}(s, t)}{\omega N_u} \eta_{uv}(s, t) = \frac{\sum_{t: (t-s) \geq \omega} w_{uv}(s, t)}{\omega N_u} \rho_{uv}(s,t) = \frac{1}{\lceil t \rceil_\omega - t}

Putting them together

Current state-of-the-art: Patients move instantly, no concept of home, using average return rates over the entire data period

Our models: from least to most ‘accurate’

  1. State-of-the-art
  2. State-of-the-art + patients return home
  3. Return rates change over time, patients held at home in “blocks”
  4. Return rates change over time, patients held at home for actual time

What difference does it make?

What difference does it make?

Takeaway

Time is critical when considering an outbreak

  • Still need to consider surveillance algorithms on these weighted, directed, temporal networks
  • Patients going home slows down infection spread
  • Most hospital transfer patterns aren’t stationary
  • Using the wrong collapsing process means you will probably miss infections
  • May declare a hospital “disease-free” when it truly hasn’t arrived yet

Hospital outbreaks

The immediate question:

When an outbreak is detected, what other hospitals could be affected?

Changes over time, but it’s not as fast as you think it is

Other ongoing work

  • Machine learning methods to identify risk factors for hospital infections
    Leong Zhuan Kee, Vis; Ewilly Liew
  • Estimating rates of contact in the absence of local data, conmat
    Nick Tierney; Nick Golding
  • Understanding how non-linear dimension reduction warps your data
    Jayani Lakshika; Di Cook
  • Predicting extreme weather events using short- and long-term climate drivers
    Kate Saunders; Jarryd Chapman
  • Approximating patient movement patterns with piecewise constant networks
    David Wu; Andrew Stewardson
  • Developing non-network based network layout algorithms
    Krisanat Anukarnsakulchularp; Di Cook

…and more.

If your work is broadly related to infectious diseases, or biosciences, I would love to talk more.

Bonus project: CPE Isolation Guidelines

Patients colonised with CPE have historically required isolation on every hospital admission post-diagnosis.

This means:

  • Contact precautions for staff and
  • Single-patient room for patients

independent of colonisation or future testing status.

Changing the guidelines

These rules have put a large strain on single-patient rooms, and has negative patient experiences in hospital.

Recently, this has been changed where individuals can be released from contact precautions if they have not tested positive for CPE in the last 12 months

Testing the impact of this change

We wanted to determine how many extra beds may be available under the new guidelines.

We assume contact precautions start on the ‘onset date’ of CPE for a patient.

  • Under the old guidelines, patients are isolated at all future admissions.
  • Under the new guidelines, after 12 months, patients are no longer isolated.

Result

Does free up beds (roughly 16 in 2019), but within 12 months this gain has been eroded.