As you can imagine, I’ve been trying for the last few days to work out what’s going on in the case data, and why it is that our models haven’t predicted the take-off in cases over the last week. I think I’m starting to get my head around it, but it’s a work-in-progress. 1/n
In essence I think we’re still thinking too much in terms of national averages, and not quite recognising that this is a wave composed of a number of local/regional outbreaks, which are concentrated in specific age groups. This has a number of implications for modelling. 2/n
In particular I want to focus on the age-group dynamic, because I think this has been under-recognised, and has led to a couple of specific problems with the models. Let’s take a quick look at the growth rate by age group in the most recent case data: 3/n
This chart shows a sort of pseudo-R for each 5-year age group. We can see the very fast growth in the 20-30s, with a doubling time of 5-7 days (R getting close to 2). Looking at the raw data for the 20-24s makes the recent growth appear even more dramatic: 4/n
Meanwhile the growth in school-age children has taken a short break over half term (due in part to reduced testing, and also hopefully to reduced infections), but looks ready to increase again, with a spike of new cases in the last few days: 5/n
This got me thinking. And to explain my thinking I’d like to borrow @BristOliver’s analogy of the cities of Novaxia and Bigpharmia, and the hats. If you’re not familiar with this, it’s essential reading (and well worth a refresher if you are): 6/n
Now let’s imagine that the population of Bigpharmia is composed of two tribes: the Overs and the Unders. The Overs are at higher risk from the disease, and have managed to ensure that they have first claim on the magic hats – in fact most of them have got two hats already. 7/n
The Unders, by contrast, mostly haven’t been offered a hat at all yet. Thankfully, they’re also at less risk from the disease, but have put up with all sorts of restrictions on their daily lives over the last year, mostly for the benefit of protecting the Overs from dying. 8/n
Now the King of Bigpharmia has announced that, with disease rates relatively low, and with the (vulnerable) Overs nearly all having at least one hat, it’s time to get back to normal life. And in particular he has re-opened the pubs, taverns and ale-houses for indoor service. 9/n
Before doing this, the wizards and sages of BigPharmia advised the King that this would probably be OK, because on average there’s a lot of hats around that will stop the disease from spreading, and not that many people go to ale-houses each day so the impact will be small. 10/n
Unfortunately, this advice failed to take account of some important differences between Overs and Unders 1) Nearly all the hats are in the possession of Overs; very few Unders have a hat 2) Unders are a more sociable tribe than Overs – they tend to meet up with more people 11/n
...and specifically they like meeting up *with other Unders*. So most of the people they meet don’t have a hat either. 3) Unders also *love* going to ale-houses. So while the average impact across the whole population of opening up the ale-houses might be small, 12/n
if we look at the impact just on Unders, it’s much larger. 4) Unders are a generally responsible crowd, and they’ve put up with a lot over the last year, mostly for the benefit of Overs. But some of the more boisterous Unders may have taken the opening of the ale-houses 13/n
...as a bit of a signal that life should return to normal in other ways. And, being at low risk from the disease themselves, they may not always be strictly following all the other rules that the King has still in force to control the spread of the disease. 14/n
All of that might have been fine were it not for the fact that, in the very same week that the ale-houses were opened, a new form of the disease arrived from Novaxia, which spread faster than the old disease and made people more likely to need treatment by the leech doctor. 15/n
I think you know what happened next. The new disease quickly spread amongst the Unders, who were:
-Mostly hatless
-Highly sociable with other hatless Unders
-Making good use of the newly-reopened ale-houses
-Not always staying two sword-lengths away from each other. 16/n
The wizards and sages were slightly perplexed by this, because they hadn’t quite anticipated the way in which the Unders would mix with other Unders, and the disproportionate impact that opening the ale-houses would have on the ability of that group to socialise. 17/n
But they were smart wizards and quickly worked out what they needed to do. They knew that if the disease continued to spread like wildfire in the Unders, then eventually some Overs would catch it too (those without a hat, or with just one hat). 18/n
And then there was a good chance of the doctor becoming overwhelmed, and probably running out of leeches. So they advised the King to issue hats immediately to as many Unders as possible. (as well as to those Overs who didn’t yet have a second hat). 19/n
Did it work? Well, I’m honestly not sure yet. The sociable nature of Unders with each other works in our favour, because the King can quickly give hats to quite a large proportion of the Unders, and in theory that should have a big impact on the spread of the disease. 20/n
But on the other hand, the new form of the disease makes the hats work less well – and in particular, means that you really need 2 hats, not one, if you want to stop the disease in its tracks. So I’m not sure yet whether giving 1 hat to each Under will be enough. 21/n
Dropping out of the analogy, I feel I need to build a model just for the under-30s (who are, obviously, the equivalent of the Unders tribe in our real-life version of Bigpharmia). This needs to consider age-group-specific contact rates (with each other and other groups), 22/n
…the impact of the Steps specifically on their lives, age-specific vaccine take-up and the changing levels of immunity in that group, as well as the impact of Delta more generally. If we can bring the epidemic in that group under control, then it feels like we’ve won. 23/n
If on the other hand it continues to grow, then it will inevitably leak into other groups, and we’ll be worrying about stocks of leeches (ICU beds?) again. I’m not sure my existing model, or my intuition, is much help in deciding which of those two futures is more likely. 24/n
As I said at the start, this is a work-in-progress, and I’d welcome additions, challenges and corrections to this line of argument. Big thanks to @vaneylenv@pmholling@rufusSG@doug_no1 and others who have, knowingly or unknowingly, contributed parts of the jigsaw. 25/n
All errors, bad jokes and stretched analogies remain mine, of course. (except where stolen from @BristOliver). And to be clear: I’m not blaming anyone for this situation, just observing where the data is leading us, and trying to ensure we fully understand the dynamics. /end
PS the implied criticism of wizards above is mostly intended as self-critique – I feel we should have been able to spot this in advance, and I didn’t. To be fair to the real sages on SPI-M, their models do contain a good dose of age-stratification, and it’s possible this dynamic
…is a partial explanation for why those models were more pessimistic about the potential for a rise in cases at Step 3, compared to national-average models such as mine. But there isn’t much evidence in the published papers to date of them using the age-stratified outputs
...of the models to inform policy-making, nor afaik were they modelling the impact of the Steps with different factors for each age group – which may turn out with hindsight to be a missed opportunity. /endPS
@kallmemeg has kindly nominated me to undertake some vaccine efficacy (VE) estimates from the data in Table 6 of the latest PHE Tech Briefing. I know the number of deaths in the “double vaxxed” column has been causing concern, but I think it’s OK. 1/n
My overall conclusions are:
1)The vaccine’s efficacy vs. disease is hard to deduce from this table, as it’s sensitive to the exposure risk assumption, but a rough estimate is consistent with the PHE’s figures in Table 18 (i.e. 33% after 1 dose and 80% after 2 doses) 2/n
2)The VE vs. hospitalisation and death appears robust to different assumptions of exposure risk, and implies that the vaccines are maintaining good protection (~80%) after 1 dose, and very strong protection {>95%) after 2 doses. 3/n
There’s something a bit odd in that CO-CIN data I tweeted on Monday, showing shorter stays in hospital more recently. The problem is, if I use the full dataset to predict the number of inpatients, I get more than there actually were in wave2. 1/n
Let’s try that to start with, using a model scenario (I’ll bring in real data in a bit). I’m taking the worst scenario from my recent model update, which has very similar peak hospitalisations to the actual peak in Jan 2021. 2/n
I can run that admissions curve through a simple model which combines it with the CO-CIN data for length of stay, and predicts what the inpatients curve would have been. In fact I get two curves: one using the “recent” data (yellow), and one for the whole dataset (blue): 3/n
thanks to the person who pointed me towards this report: assets.publishing.service.gov.uk/government/upl… (you know who you are!) - lots of interesting data. In particular it does support the anecdotal evidence that average length of stay in hospital is shorter in recent admissions: see this chart 1/3
...and compare to the older data in the chart below. So after 10 days, about 75% have been discharged in the recent data, compared to 40% in the older data. (but note deaths are lower in the recent data, which slightly offsets the benefits for hospital bed occupancy). 2/3
what's not entirely clear to me is whether this is just the result of a lower average age of patients being admitted - see below- (and would reverse if more older cases arrived), or whether it says something more fundamental about the progression of covid disease with Delta. 3/3
Model update klaxon!! I’ve spent a while this weekend refitting my model to current data and trends, and as a result I’ve changed my mind on at least 3 things – and also got closer to a policy recommendation for Step 4 re-opening scheduled for 21st June. Thread follows: 1/n
First to note a few assumptions:
-Baseline controls and cautious behaviour (e.g. continued WFH) are assumed to continue after 21st June, until the end of 2021, and to reduce R by ~25%
-I’ve adjusted my vaccination schedules to more precisely match current supply estimates 2/n
– with thanks as always to @PaulMainwood for his insights on this.
-In particular, now that we have MHRA approval for 12-15 year olds to receive Pfizer, I’ve assumed that 12-17s get their first doses in September (& second doses in November) 3/n
I think the age dynamics of the current case growth are interesting, particularly as they have shifted in the last week or two. Growth is now very clearly dominated by the 20-30s, while growth in school-age children has tailed off: 1/6
If we look more closely e.g. at 5-9 year olds, we can see that tailing off of growth. I’m not 100% sure why this has happened, it seems slightly early for a half term effect, so maybe local PH efforts and bubble closures are having an impact? &/or parents being vaxxed?? 2/6
Meanwhile the growth in 20-24s is pretty dramatic, doubling in just under a week, and now back to levels not seen since February: 3/6
A few thoughts on yesterday’s PHE data: overall, the report is a bit of a mixed bag. Of the three significant bits of news, one is positive, one negative, and the other neutral. But you might not get that impression from the reaction here on Twitter or in the media. 1/n
Let’s take the negative first: we now have data suggesting that the new variant (Delta) is more likely to lead to people being admitted to hospital – perhaps 2.5 times more likely than when infected with the old Alpha variant. 2/n
That’s clearly bad news, and will have an impact on model projections for Step 4. But it’s not necessarily a complete disaster: if we can control the spread of the virus, then it doesn’t matter what the hospitalisation ratio is, because very few people will be catching it. 3/n