Here’s a great example of how sometimes modelling produces counter-intuitive results: if (under certain conditions) you randomly infect 100,000 people with covid in the next few weeks, you’d save 250 lives. What??? Yes, that’s right, more infections means fewer deaths. 1/
I should be clear before we go any further that I’m not proposing this as a practical suggestion, for reasons that will become clear. It’s more a mathematical curiosity, but it helps to expose a couple of interesting dynamics, so I think it’s worth sharing with you. 2/
So, to start with let’s set up my model with a moderately-sized exit wave, as predicted by the Warwick or Imperial models (NB this is not my central case – I am currently predicting no exit wave, or a very small one – but it’s not impossible). 3/
Now, let’s introduce a “spike” in infections in early May. I do this by pushing R up above its trend for one week, and then pushing it down below trend on week 2 – in order to neutralise the impact of the spike on future cases – and then back to trend from week 3 onwards. 4/
Let’s look at the infection curve in both the original scenario – in blue – and our adjusted scenario (with spike) in orange. It’s a bit hard to see the difference, but perhaps you can just make out the spike in early May, and then a reduction in the height of the exit wave. 5/
It’s a lot clearer if we just look at the difference between the two infection curves: now you can see the spike very clearly, but also a gradual reduction in the size of the exit wave, and then a slight increase at the end (as the wave is pushed slightly later in time). 6/
As it turns out, the number of extra infections in the spike is around 107,000. But the number of infections in the exit wave is reduced by 195,000. So the net effect of adding the spike has been to *reduce* the total number of infections by about 88,000. 7/
But why is it doing that? Well, simply put the exit wave happens because we’re not yet at the herd immunity threshold when we unlock. And when we do unlock, infections will rise until we reach the HIT, when they peak (with R=1) and then start to decline. 8/
If we unlock at low prevalence, the exit wave will be roughly symmetric around the point where we do hit the HIT. So if we add one infection before we unlock, we’re reducing the gap to HIT by one infection. So that’s one less infection in the left-hand side of the exit wave 9/
…and (with a roughly symmetric wave) also one less in the right-hand side. So we’re exchanging one infection prior to unlocking for (roughly) two after unlocking. That’s a good trade – I can now see how adding a spike might reduce total infections, and hence deaths. 10/
But that’s not the whole story. If we look at the difference in deaths for our two scenarios – original and adjusted – we can see again the spike (with 155 additional deaths) and the reduced exit wave (with 405 fewer deaths). Hence a reduction of 250 deaths overall. 11/
But that’s a bit more than we’d expect. It looks like the ratio of deaths to infections (IFR) is a bit lower in the spike than in the exit wave. So when we trade 1 infection in the spike for ~2 in the exit wave, that’s not the only benefit: the infections we are saving 12/
…later on are more deadly than the one we are adding in the spike. So while we’re saving 1.8 infections for every one we add (1.8 = 195 / 107), we’re actually saving 2.6 deaths for every death we cause in the spike (2.6 = 405 / 155). OK, that’s even better news. But why? 13/
Well, if we plot the IFR over time we can see what’s happening. IFR in my model starts out around 0.8%, and falls as we vaccinate the vulnerable population to a low of about 0.15% in late April. But then it starts to drift up over June & July, reaching 0.24% in September. 14/
OK so mathematically I can see what’s happening, but why is the IFR going back up? To answer that we need to remember that the IFR is not an average over the whole population, but over the people who are infected. And to understand that who gets infected changes over time. 15/
The people who get infected are, by definition, those who are susceptible to the disease. So let’s look at that susceptible group. The number of them is reducing over time, but what we’re interested in here is the composition, so we’ll use a 100% stacked area chart: 16/
We can see that as we vaccinate more, the %ge of susceptibles who have some partial immunity (and are therefore protected against severe disease and death) is growing – see the grey area. By definition the IFR in this group is zero, so as it grows, the overall IFR reduces. 17/
The second effect has to do with the mix of the other categories. Broadly speaking, the higher % that the vulnerable groups (orange & yellow) form of the total, the higher IFR is. And we can see that proportion going down fast in Jan/Feb/Mar as we vaccinate those groups. 18/
But then it reaches a low point in Apr/May, and starts to climb back up again in Jun/Jul as we vaccinate the younger groups. This doesn’t mean those vaccinations are a bad idea –it’s still a big positive overall– but it does have the side-effect of meaning IFR goes up a bit. 19/
So this explains why, if we increase infections in April/May, it has a lower effect on deaths than if we did the same in June/July. And hence why “trading” higher infections in April/May for a lower infection rate in June/July is a good deal. 20/
One way of looking at this is we’ve temporarily created a “Great Barrington”-like situation where the vulnerable are mostly protected, but the young are still susceptible. So if you know you need more infections to get over the HIT, now’s the time to take them. 21/
But now I need to remind you why I’m not actually promoting this as a plan. There are two big reasons:

1) This whole analysis depends on there being an exit wave – without it, there would be no reduction in infections or deaths, and you’d just be adding extra infections 22/
…and deaths in the spike. And I’m not expecting there to be an exit wave when we unlock in June.

2) In reality we can’t control R in the way that I have done in the model, to ensure that the ‘spike’ is neutralised and we don’t set off a chain of growing infections. 23/
Without the neutralisation, we give the exit wave more momentum, and we end up killing 2000 extra people, rather than saving 250. So the real-world risks of accidentally making things a lot worse are very high. So it’s really not something I would recommend trying. 24/
But I do think it’s interesting that if you were in govt, convinced by SAGE’s briefing papers that an exit wave was a nailed-on certainty, and someone gave you the choice of pushing infections up for 1 week, or down for 1 week, the correct answer is: make infections go up. /end

• • •

Missing some Tweet in this thread? You can try to force a refresh
 

Keep Current with James Ward

James Ward Profile picture

Stay in touch and get notified when new unrolls are available from this author!

Read all threads

This Thread may be Removed Anytime!

PDF

Twitter may remove this content at anytime! Save it as PDF for later use!

Try unrolling a thread yourself!

how to unroll video
  1. Follow @ThreadReaderApp to mention us!

  2. From a Twitter thread mention us with a keyword "unroll"
@threadreaderapp unroll

Practice here first or read more on our help page!

More from @JamesWard73

13 Apr
So today’s big covid Twitter debate has been the relative contribution of lockdown (or, more broadly, controls / NPIs) and vaccinations to bringing down case rates. It occurs to me that we already have some numbers – or at least assumptions – that we can look at: 1/4 Image
This chart shows the assumptions in my model for the relative reduction in R from controls, and from immunity. As you’d expect, in January it’s nearly all controls – and by June it’s mostly immunity. The cross-over point is in mid-April – in fact, round about now. 2/4
But of course immunity is a composite of vaccination-acquired immunity, and infection-acquired immunity. So let’s split those out. We can see that infections grow slowly, and have already been overtaken by vaccines. But vaccines don’t overtake controls until early June. 3/4 Image
Read 5 tweets
9 Apr
I’ve been getting increasingly confident in recent weeks that we’re on the right path to defeating B.1.1.7 in the UK, and stand a good chance of ‘unlocking’ in May/June with a no exit wave, or a very small one. But there’s one thing still bothering me: new variants. 1/
Up ‘til now, I’ve not been sure how worried to be about these. Many people (including some whose opinions and expertise I trust) tell me not to worry, there won’t be much ‘immunity escape’, and in any event T-cells will ensure that we remain protected vs. severe disease. 2/
But the questions still bothering me are: 1) how much ‘immunity escape’ would be needed to cause a problem? 2) and how likely is that to occur? I can’t answer the second question – that’s one for the virologists and immunologists – but I can have a go at the first. 3/
Read 37 tweets
7 Apr
Having spent a good chunk of yesterday going through the recent SPI-M papers including the latest models from Warwick, Imperial and LSHTM, and then come back up for air, what have I learned? Here’s my summary of ten important messages: 1/14
1. While the headlines have once again been full of doom-laden reports of a huge ‘exit wave’ when we unlock in the summer, the reality is that both Warwick and Imperial models are now predicting a much smaller wave with only 15-20k deaths in their central scenarios. 2/14
2. Looking at the Warwick model in more detail, the change from previous iterations is mostly due to improved (more optimistic – or in my view, realistic) assumptions for vaccine rollout, vaccine take-up, vaccine efficacy vs. severe disease, and efficacy vs. infection. 3/14
Read 15 tweets
2 Apr
After yesterday’s model update () a few people have asked: so if there’s no exit wave, or just a small one, does that mean we could open up earlier? And the answer is: maybe, but it’s complicated. If you’re interested in the details, read on… 1/n
All of yesterday’s modelling assumed that we follow the government’s roadmap to 21st June, when we remove all remaining restrictions and controls, and behaviour returns to normal. (note that point: all controls removed on 21st June, it will become important later on). 2/n
But what happens if we vary this date? I’ve adjusted this in the model, bringing it back 1 week at a time up to 5th April, and also extending it the other way to 26th July, for completeness. As you can see (from blue line below), opening on 21st June looks OK, and you could 3/n
Read 15 tweets
1 Apr
I’ve updated my model for recent news, and after a series of assumption changes that mostly net out, it still predicts a very small exit wave with ~4k additional deaths (reduced from 9k previously) – and with a small added dose of seasonality, no exit wave at all. 1/n
The key changes in the model are shown in the waterfall chart below, and explored in more detail in the thread that follows this summary. 2/n
This conclusion feels more solid than before, as no single-factor sensitivity (within plausible ranges of uncertainty) takes the exit-wave deaths over 10k- see 1st chart. To get significantly higher deaths I have to move multiple assumptions at the same time– see 2nd chart. 3/n
Read 17 tweets
24 Mar
I’ve spent a lot of time over the last 2 days going through the most recent paper thelancet.com/action/showPdf… from the Warwick modelling group that feeds into SAGE, and trying to work out why it predicts a large exit wave, when my model (mostly) doesn’t. My conclusions are: 1/12
1. Their model in fact doesn’t guarantee an exit wave: with fast vaccine rollout (4m per week) and high (85%) transmission blocking, and either high (85-95%) vaccine uptake or high (94%) protection vs. severe disease, there is no material exit wave - see the orange line 2/12
2. Their choice of central assumptions is (imo) consistently on the pessimistic side. In fact, my central assumptions frequently correspond closely to their “optimistic” upside scenario. To avoid a very long thread I’ve built a table that compares their assumptions to mine 3/12
Read 13 tweets

Did Thread Reader help you today?

Support us! We are indie developers!


This site is made by just two indie developers on a laptop doing marketing, support and development! Read more about the story.

Become a Premium Member ($3/month or $30/year) and get exclusive features!

Become Premium

Too expensive? Make a small donation by buying us coffee ($5) or help with server cost ($10)

Donate via Paypal Become our Patreon

Thank you for your support!

Follow Us on Twitter!