Share

PANDA | Please hear us out, NICD!

0:00
play article
Subscribers can listen to this article

The NICD needs to be more open to different types of modelling in order to fully understand the Covid-19 pandemic in South Africa, writes Pandemic Data and Analytics (PANDA). 


Dr Harry Moultrie of SACMC's article regarding its Epi Model for Covid was well received, suggesting a near-term release of their model's assumptions and code, as we and others have been requesting for some time.

PANDA published its Western Cape and South African models to coincide with the Epi Model update and they continue to perform well, despite having much narrower confidence intervals than the Epi Model.

However, this has resulted in PANDA receiving a slew of requests from doctors, hospital administrators and journalists all over the country who had become frustrated with the NICD's guidance and wanted a steer from us on what they should plan for.

We are a small, not-for-profit business staffed by part-time volunteers, and do not yet have the infrastructure to deal with such requests, as much as we'd like to.

The picture for South Africa as a whole looks like this:

graph panda

As can be seen, the Epi Model has now departed too far from reality to be useful. (See also our website, www.pandata.org.za.)

The NICD have constantly declined to engage with us, but we would like to propose a way out of the morass. First, a few observations: We get the sense from comments on our Twitter feed that SACMC members' objection to our modelling approach has more to do with it being descriptive (based on analysis of international epidemic curves and a study of inter-country mortality variations) as opposed to having a theoretical causal structure, as have the SEIR-type models, on which the Epi Model is based.

Though the "Gompertz" curves that we deploy, as proposed by Nobel laureate Professor Michael Levitt of Stanford University, achieve very good fits to observed experience elsewhere, there is a type of epidemiological model, the so-called "D-model" - introduced to us by Dr Azwinndini Muronga of NMMU - which we have found to fit almost as well as the Gompertz curves.

graph panda
graph panda

The diagram above includes a D-model fit for South Africa, which we have been evaluating since 10 July.

Note its similarity to the Gompertz curves, which form the lower and upper confidence intervals and within which the data for South Africa, Gauteng and the Western Cape have been tracking since the calibration on 10 June.

We note the recent aberrations in death reporting for both the Eastern Cape and Gauteng and consciously describe these even though it must be argued that smoothing would be justifiable.

An advantage of the D-models is its mathematical simplicity. They have closed form solutions and can easily be implemented on a spreadsheet in a few hours - no need for pages of code or numerical solutions to complex differential equations.

A potential sticking point is that D-models assume that recovered individuals are not reinserted into the population. This is equivalent to assuming that the infected population is set up rapidly or that the infection process is not mainly dependent on human-to-human contact. The good fit of the D-model may surprise in the context of Covid-19 given that:

  • The main mode of transmission is thought to be human-to-human;
  • Outbreaks seem to last long enough for some individuals to become symptomatic after others have recovered;
  • It is thought that the incubation or pre-symptomatic period for the disease is quite short.

This calls into question whether the other key assumptions underpinning the D-model are correct. These classically involve a random walk process generating human contact, uniform mortality and binary (as opposed to graduated) susceptibility. If those assumptions are wrong for Covid, this could explain why the SEIR models being so widely applied cannot be made to fit empirical data as closely as Gompertz and D-models.

Nevertheless, we beg the NICD to consider taking up this proposal, so that it can provide useful estimates to our healthcare institutions, who desperately need them.

To do so effectively, it would have to adjust some assumptions and infer ICU and hospital bed demand from death numbers described, for which purpose sufficient data are at hand.

Based on rumour and reverse-engineering, we understand that the Epi Model deploys a susceptibility parameter of approximately 20% and an infection fatality rate of approximately 0.4%. This is similar to the result from the New York State seroprevalence study.

We think the latter is too high in the context of the former, given differences in the age structure, comorbidity prevalence and obesity prevalence in the South African population, and would suggest moderating this to 0.1%.

Our assessment is that SEIR-type models have generated bad fits to the data, no matter what assumptions are used and despite having more parameters.

Note that none of these Gompertz and D-models make any assumptions regarding continuation or suspension of lockdown, because nothing in the international data suggests that lockdown stringency or mobility has any material curve-flattening effect.

This observation comes from detailed statistical analysis of single country curve behaviour following lockdown inception and cessation, and from an inter-country analysis that found no correlation between lockdown stringency and country mortality rates or durations to peak deaths.

We had come to suspect it given the universal fit of the Gompertz models, which would be inconsistent with there being any impact from mitigation measures in general. We have reviewed several papers claiming to find beneficial effects from lockdown, but these have found small effects, which can safely be ignored, and are methodologically dubious in any event, and therefore unlikely to survive replication.

Also note that the model does not make any downward adjustment for "deaths with" Covid. These will not place the same kind of burden on resources as "deaths from" Covid, so it is important to have an estimate for them.

Our observation of the US data indicates that the proportion of "deaths with" may be significant. We have heard many accounts of the same thing happening in South Africa, despite the excellent guidance note of Dr SSS Buthelezi (Director General of Health). This may go some way to explaining the rest of the Epi Model's 10x overestimate for Western Cape ICU demand, which would not be totally removed by the re-parameterisation suggested above.

We remain willing to engage and share our models with the NICD and sincerely hope its position towards us moderates.

We have further ideas around a scaling approach that might be adopted to reduce the provincial models to single-institution forecasts, which hospitals are likely to find useful.

Such an approach would relieve us from an administrative burden we are incapable of meeting and put an end to what has become an unseemly public exchange between PANDA and the NICD.

- written by Pandemic Data and Analytics (PANDA)

References

Amaro, J.E., J´er´emie, D., Orce, J.N., 15 May 2020. Global analysis of the COVID-19 pandemic using simple epidemiological models. https://arxiv.org/pdf/2005.06742.pdf

Muronga, A. 2020. Epidemiological Models for Covid-19. Nelson Mandela University, Faculty of Science.


*Want to respond to the columnist? Send your letter or article to opinions@news24.com with your name, profile picture, contact details and location. We encourage a diversity of voices and views in our readers' submissions and reserve the right not to publish any and all submissions received.

Disclaimer: News24 encourages freedom of speech and the expression of diverse views. The views of columnists published on News24 are therefore their own and do not necessarily represent the views of News24.

We live in a world where facts and fiction get blurred
Who we choose to trust can have a profound impact on our lives. Join thousands of devoted South Africans who look to News24 to bring them news they can trust every day. As we celebrate 25 years, become a News24 subscriber as we strive to keep you informed, inspired and empowered.
Join News24 today
heading
description
username
Show Comments ()
Voting Booth
Do you think corruption-accused National Assembly Speaker Nosiviwe Mapisa-Nqakula will survive a motion of no confidence against her?
Please select an option Oops! Something went wrong, please try again later.
Results
No, her days are numbered
41% - 502 votes
Yes, the ANC caucus will protect her
59% - 717 votes
Vote
Rand - Dollar
18.90
+0.2%
Rand - Pound
23.86
+0.2%
Rand - Euro
20.37
+0.3%
Rand - Aus dollar
12.31
+0.2%
Rand - Yen
0.12
+0.2%
Platinum
908.05
0.0%
Palladium
1,014.94
0.0%
Gold
2,232.75
-0.0%
Silver
24.95
-0.1%
Brent Crude
87.00
+1.8%
Top 40
68,346
0.0%
All Share
74,536
0.0%
Resource 10
57,251
0.0%
Industrial 25
103,936
0.0%
Financial 15
16,502
0.0%
All JSE data delayed by at least 15 minutes Iress logo
Editorial feedback and complaints

Contact the public editor with feedback for our journalists, complaints, queries or suggestions about articles on News24.

LEARN MORE