@PausalZ@fediscience.org Profile picture
Sep 24, 2020 15 tweets 4 min read Read on X
Herd immunity is a far squishier concept then many seem to be describing in their "shielding" or "stratified herd immunity" plans. Here is the formula for herd immunity threshold for a SIR model Image
where \beta is the effective contact rate, N is the number of individuals, and r is the inverse of the duration

The threshold says if are above that level the disease will disappear / we expect no outbreaks of disease. However, that threshold is neither sufficient nor necessary
To show this, let's talk about a perfect vaccine. If you get this vaccine you are perfectly protected from the infection and thus cannot transmit it (everything also applies to imperfect vaccines but it's messier)

Blue circles are vaccinated individuals and red are unvaccinated
In all the following networks, the elements of R_0 are held constant. The network structure will change; but the \beta look the same despite that

Our R_0 will be 3, meaning the herd immunity threshold is 0.67 (but everything applies to other R_0 > 1)
Network 1: Random Mixing

When someone talks about a threshold for herd immunity, this is the underlying network of what they are generally talking about (setting aside WAIFW for the moment). The threshold calculation applies normally Image
Network 2: Mesh Network

Well it looks like a big ole net. We can guarantee no transmission with only 50% vaccinated. Less than the supposed threshold Image
"Paul, all you have shown is that you can effectively go below the threshold. That doesn't negate the whole concept. We might instead say we need *at least* the threshold"

Bad news my imaginary interlocutor
Network 3: Clustered Network

In a clustered network, we can go above the herd immunity threshold but still have outbreaks. The following network has 75% vaccinated but an outbreak would occur Image
This network is closer to how human interactions are actually structured. This is also closer to WAIFW setups (but WAIFW won't capture the power-law for degree within clusters that the above network has)
Network 4: Household Network

This last network is variation on the previous clustered network. It instead places people into households. Depending on the vaccination strategy used, we can go above or below the 'threshold' Image
If your strategy is to vaccinate the between-household contacts, you can be lower than the threshold Image
If your strategy is instead to vaccinate the persons with household-only contacts then you will have to go above the threshold Image
Ultimately, the (possibly) unobserved network matters. If your model only assumes some random mixing parameter for broad groups of people, your model is probably way too simple to say anything meaningful about what policies we should actually consider
Because I did a bad job of defining WAIFW: it is the Who Acquires Infection From Whom matrix. Rather than random edges between persons, WAIFW generates edges random at \beta_w for connections in the group and \beta_o for outside
Below is a 3 group WAIFW with different Pr of connections in each group (by color) and between groups. Essentially WAIFW is a stratified random mixing model Image

• • •

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

Keep Current with @PausalZ@fediscience.org

@PausalZ@fediscience.org 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 @PausalZ

Nov 17, 2021
a 🧵 on M-Estimation and why I think its a valuable tool that epidemiologist should be using more often
M-Estimation is a general approach of defining an estimator as the solution to estimating equations like the following. Importantly, obs are independent and \psi is a known function that doesn't depend on i or n
I think its a great tool for two reasons: (1) the ability to stack estimating equations together, and (2) the sandwich variance
Read 17 tweets
Jul 7, 2021
Big fan of the "I forced a bot to [...] over 1000" memes. But most of those posts are fake (i.e. human-generated). That's why I decided to make a real one

So I forced a bot to read over 1000 PubMed abstracts in order to generate new abstracts ImageImage
Basically, I pulled a random sample of 5000 abstracts from PubMed using the search terms: (causal inference) AND English[Language]

A random sample of the returned abstracts was used to train a recurrent neural network (RNN)
Basically, a sequence of 40 characters is used to predict the next character. This process can then be repeated with the new character to generate a whole new sentence

So you give the machine a starting point, set a 'creativity dial', and let it go
Read 12 tweets
Sep 20, 2020
8: WHEN CAN I IGNORE THE METHODOLOGISTS
Section 8 discusses when standard analytic approaches are fine (aka time-varying confounding isn't as issue for us). Keeping with the occupation theme, it is presented in the context of when employment history can be ignored Image
First we go through the simpler case of point-exposures (ie only treatment assignment at baseline matters). Note that while we get something similar to the modern definition, I don't think the differentiation from colliders is quite there yet (in the language) ImageImage
Generalization of the point-exposure definition of confounding to time-varying exposures isn't direct Image
Read 13 tweets
Sep 19, 2020
7: MORE ASSUMPTIONS
Section 7 adds some additional a priori assumptions that can allow us to estimate in the context where we don't have all necessary confounders.
We have the beautifully named: A-complete Stage 0 PL-sufficient reduced graph of R CISTG A Image
We start with some rules for reducing graph G_A to a counterpart G_B. Honestly the language in this section isn't clear to me despite reading it several times... ImageImage
I do think the graphs help a bit though. To me it seems we are narrowing the space of the problem. We are going from multiple divisions at t_1 and t_2 to only considering the divisions at t_2 for a single branch. The reduced STG is a single branch ImageImage
Read 6 tweets
Sep 15, 2020
6: NONPARAM TESTS
Section 6 goes through the sharp null hypothesis (that no effect of exposure on any individual). Note that this is weaker than the null of no _average_ effect in the population Image
Another way of thinking about this is if there is no individual causal effect (ICE) then there must be no average causal effect (ACE). The reverse (no ACE then no ICE) is not guaranteed
Robins provides us with the G-null hypothesis as a means of assessing the sharp null (the g-null is that call causal parameters are 0) ImageImageImage
Read 9 tweets
Sep 13, 2020
5: ESTIMATION
After a little hiatus, back to discussing Robins 1986 (with a new keyboard)! Robins starts by reminding us (me) that we are assuming the super-population model for inference Image
If we had a infinite n in our study, we could use NPMLE. However, time-varying exposures have a particular large number of possible intervention plans. We probably don't have anywhere near enough obs to consider all the possible plans Image
Instead we use a parametric projection of the time-varying variables. We hope that the parametric projection is sufficiently flexible to approx the true density function (it is why it is best to include as many splines and interaction terms as feasible)
Read 8 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

Don't want to be a Premium member but still want to support us?

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

Donate via Paypal

Or Donate anonymously using crypto!

Ethereum

0xfe58350B80634f60Fa6Dc149a72b4DFbc17D341E copy

Bitcoin

3ATGMxNzCUFzxpMCHL5sWSt4DVtS8UqXpi copy

Thank you for your support!

Follow Us!

:(