/home/wpurdemo/resistancesuffolk/wp-content/mu-plugins
Estimating COVID-19’s Rt in Real-Time
Kevin Systrom – April 17
In any epidemic, Rt is the measure known as the effective reproduction number. It’s the number of people who become infected per infectious person at time t. The most well-known version of this number is the basic reproduction number: Rt when t=0. However, R0 is a single measure that does not adapt with changes in behavior and restrictions.
As a pandemic evolves, increasing restrictions (or potential releasing of restrictions) changes Rt. Knowing the current Rt is essential. When Rt is >1, the pandemic will spread through a large part of the population. If Rt is <1 the pandemic will slow quickly before it has a chance to infect many people. The lower the Rt: the more manageable the situation. In general, any Rt <1 means things are under control.
The value of Rt helps us in two ways. (1) It helps us understand how effective our measures have been controlling an outbreak and (2) it gives us vital information about whether we should increase or reduce restrictions based on our competing goals of economic prosperity and human safety. Well-respected epidemiologists argue that tracking Rt is the only way to manage through this crisis.
Yet, today, we don’t yet use Rt in this way. In fact, the only real-time measure I’ve seen has been for Hong Kong. More importantly, it is not useful to understand Rt at a national level. Instead, to manage this crisis effectively, we need a local (state, county and/or city) granularity of Rt.
What follows is a solution to this problem at the US State level. It’s a modified version of a solution created by Bettencourt & Ribeiro 2008 to estimate real-time Rt using a Bayesian approach. While this paper estimates a static R value, here we introduce a process model with Gaussian noise to estimate a time-varying Rt.
If you have questions, comments, or improvements feel free to get in touch: hello@systrom.com. And if it’s not entirely clear, I’m not an epidemiologist. At the same time, data is data, and statistics are statistics and this is based on work by well-known epidemiologists so you can calibrate your beliefs as you wish. In the meantime, I hope you can learn something new as I did by reading through this example. Feel free to take this work and apply it elsewhere – internationally or to counties in the United States.
Additionally, a huge thanks to Frank Dellaert who suggested the addition of the process and to Adam Lerer who implemented the changes. Not only did I learn something new, it made the model much more responsive.
If you want more math, check this out: https://github.com/k-sys/covid-19/blob/master/Realtime%20R0.ipynb
]]>