Estimating and Simulating a SIRD Model of COVID-19 for Many - - PowerPoint PPT Presentation

estimating and simulating a sird model of covid 19 for
SMART_READER_LITE
LIVE PREVIEW

Estimating and Simulating a SIRD Model of COVID-19 for Many - - PowerPoint PPT Presentation

Estimating and Simulating a SIRD Model of COVID-19 for Many Countries, States, and Cities Jes us Fern andez-Villaverde and Chad Jones August 28, 2020 0 xkcd: Everyones an Epidemiologist Help macroeconomists with data and model 1


slide-1
SLIDE 1

Estimating and Simulating a SIRD Model of COVID-19 for Many Countries, States, and Cities

Jes´ us Fern´ andez-Villaverde and Chad Jones

August 28, 2020

slide-2
SLIDE 2

xkcd: Everyone’s an Epidemiologist Help macroeconomists with data and model

1

slide-3
SLIDE 3

Outline

  • Setup
  • SIR model with a time-varying R0t
  • Recover R0t as the “Solow residual” of SIR to fit deaths
  • Estimation and simulation
  • Different countries, U.S. states, and global cities
  • “Forecasts” from each of the last 7 days
  • Re-opening and herd immunity
  • How much can we relax social distancing?

Our dashboard contains 30+ pages of results for each of 100 cities, states, and countries

2

slide-4
SLIDE 4

Basic Model

3

slide-5
SLIDE 5

Notation

  • Number of people who are (stocks):

St = Susceptible It = Infectious Rt = Resolving Dt = Dead Ct = ReCovered

  • Constant population size is N

St + It + Rt + Dt + Ct = N

4

slide-6
SLIDE 6

SIRD Model: Overview

  • Susceptible get infected at rate βtIt/N

New infections = βtIt/N · St

  • Fraction γ of Infectious resolve each day, so the average number
  • f days that a person is infectious is 1/γ so γ = .2 ⇒ 5 days
  • Fraction θ of Post-infectious cases resolve each day. E.g.

θ = .1 ⇒ 10 days

  • Resolution happens in one of two ways:
  • Death: fraction δ
  • Recovery: fraction 1 − δ

5

slide-7
SLIDE 7

SIRD Model: Laws of Motion ∆St+1 = −βtStIt/N

  • new infections

∆It+1 = βtStIt/N

  • new infections

− γIt

  • resolving infectious

∆Rt+1 = γIt

  • resolving infectious

− θRt

  • cases that resolve

∆Dt+1 = δθRt

  • die

∆Ct+1 = (1 − δ)θRt

  • reCovered

6

slide-8
SLIDE 8

Recyled notation R0: Initial infection rate

  • Initial reproduction number R0 ≡ β/γ

R0 = β × 1/γ

# of infections from one sick person # of “adequate” contacts per day # of days contacts are infectious

  • R0 = expected number of infections via the first sick person
  • R0 > 1 ⇒ disease initially grows
  • R0 < 1 ⇒ disease dies out: infectious generate less than 1

new infection

  • If 1/γ = 5, then easy to have R0 >> 1

7

slide-9
SLIDE 9

Basic Properties of Differential System (Hethcote 2000)

  • Continuous time, constant β
  • Initial exponential growth rate of infections is

β − γ = γ (R0 − 1)

  • Let st ≡ St/N = fraction susceptible
  • Infectious grow at β − γ = γ (R0st − 1)
  • If R0st > 1, the virus spreads, otherwise declines
  • As t → ∞, the total fraction of people ever infected, e∗, solves

(assuming s0 ≈ 1) e∗ = − 1 R0 log(1 − e∗) Long run is pinned down by R0 (and death rate), γ and θ affect timing

8

slide-10
SLIDE 10

Social Distancing

  • What about a time-varying infection rate βt?
  • Disease characteristics – fixed, homogeneous
  • Regional factors (NYC vs Montana) – fixed, heterogeneous
  • Social distancing – varies over time and space
  • Reasons why βt may change over time
  • Policy changes on social distancing
  • Individuals voluntarily change behavior to protect

themselves and others

  • Masks, superspreading events

9

slide-11
SLIDE 11

Recovering βt and R0t

  • Recover βt, a latent variable, from the data:
  • Like the Solow Residual of the SIRD model!
  • Notation
  • Dt+1: stock of deaths as of the end of date t + 1
  • ∆Dt+1 ≡ dt+1: number of people who die on date t + 1
  • With algebra, “invert” the SIRD model to obtain:

βt = N St

  • γ +

1 θ∆∆dt+3 + ∆dt+2 1 θ∆dt+2 + dt+1

  • St+1 = St
  • 1 − βt

1 δγN 1 θ ∆dt+2 + dt+1

  • 10
slide-12
SLIDE 12

Recovering βt and R0t (continued) βt = N St

  • γ +

1 θ∆∆dt+3 + ∆dt+2 1 θ∆dt+2 + dt+1

  • St+1 = St
  • 1 − βt

1 δγN 1 θ ∆dt+2 + dt+1

  • Use data on dt, and initial condition S0/N ≈ 1,
  • Iterate forward in time and recover βt and St+1
  • Uses future deaths over the next 3 days to tell us about βt today
  • More general point about SIRD models
  • State-space representation that we can exploit
  • Richer structure possible (heterogeneity, general functions)

11

slide-13
SLIDE 13

An endogenous R0t when simulating future outcomes

  • Individuals react endogenously to risk
  • Much of the reaction is not even government-mandated
  • Could solve a complex dynamic programming problem
  • Instead, Cochrane (2020) suggests:

R0t = Constant · e−αdt where dt is daily deaths per million people.

  • We estimate α from our data on R0t

12

slide-14
SLIDE 14

Estimates and Simulations

13

slide-15
SLIDE 15

Parameters assumed fixed and homogeneous

  • γ = 0.2: average duration is 5 days (or γ = 0.1)
  • θ = 0.1: average duration post-infectious is 10 days.
  • average case takes 10+5 = 15 days to resolve.
  • Long tail for exponential distribution
  • α = 0.05: estimate αi for each location i.
  • Tremendous heterogeneity across locations
  • R0t falls by 5 percent with each daily death
  • We report results with α = 0 and α = .05.

14

slide-16
SLIDE 16

Mortality rate (IFR): δ = 1.0%

  • Evidence from seroepidemiological national survey in Spain:
  • Stratified random sample of 61,000 people
  • δ in Spain is between 1% and 1.1%.
  • Correction by demographics to other countries
  • Most countries cluster around 1%.
  • U.S.: 0.76% without correcting for life expectancy and 1.05%

correcting by it.

15

slide-17
SLIDE 17

Heterogeneity in Mortality Rates by Age

  • Mortality rates vary substantially by age
  • IFR for ages 65-69 in Spain = 1%
  • Gompertz Law: mortality rate grows exponentially with age
  • COVID-19: doubles for every 5-year age group (?)
  • 70 vs 20 year olds: 50 years = 10 doublings ⇒ 1000-fold
  • 2 in 100 versus 2 in 100,000
  • Our estimation does not feature this heterogeneity — lack of data

May underestimate herd immunity if many young people are increasingly infected

16

slide-18
SLIDE 18

Estimation based entirely on death data

  • Johns Hopkins University CSSE data
  • Excess death issue
  • Currently no correction, just using the JHU/CSSE data
  • (previously adjusted upward by 33%)
  • We use 7-day moving averages (centered)
  • Otherwise, very serious “weekend effects” in which deaths

are underreported

  • Even zero sometimes, followed by a large spike
  • Further smoothing: HP-filter with smoothing parameter 800

after taking moving average

17

slide-19
SLIDE 19

New York City: Estimates of R0t = βt/γ

Mar Apr May Jun Jul Aug Sep 2020 0.5 1 1.5 2 2.5 3 3.5

R0(t) New York City (only) = 0.010 =0.10 =0.20

18

slide-20
SLIDE 20

New York City: Daily Deaths and HP Filter

Feb Mar Apr May Jun Jul Aug 2020 10 20 30 40 50 60 70 80 90 100

New York City (only): Daily deaths, d = 0.010 =0.10 =0.20

19

slide-21
SLIDE 21

New York City: Change in Smoothed Daily Deaths

Feb Mar Apr May Jun Jul Aug 2020

  • 4
  • 3
  • 2
  • 1

1 2 3 4 5

New York City (only): Delta d = 0.010 =0.10 =0.20

20

slide-22
SLIDE 22

New York City: Change in (Change in Smoothed Daily Deaths)

Feb Mar Apr May Jun Jul Aug 2020

  • 0.6
  • 0.4
  • 0.2

0.2 0.4 0.6 0.8 1

New York City (only): Delta (Delta d) = 0.010 =0.10 =0.20

21

slide-23
SLIDE 23

New York City: Percent of the Population Currently Infectious

Mar Apr May Jun Jul Aug Sep 2020 1 2 3 4 5 6

Percent currently infectious, I/N (percent) New York City (only) Peak I/N = 5.11% Final I/N = 0.02% = 0.010 =0.10 =0.20

22

slide-24
SLIDE 24

Estimates of R0t = βt/γ

Mar Apr May Jun Jul Aug Sep 2020 0.5 1 1.5 2 2.5 3

R0(t) = 0.010 =0.10 =0.20

New York City (only) Paris, Ile-de-France Miami Brazil SF Bay Area

23

slide-25
SLIDE 25

Percent of the Population Currently Infectious

Feb Mar Apr May Jun Jul Aug Sep 2020 1 2 3 4 5 6

Percent Infectious, I/N

New York City (only) Lombardy, Italy Boston+Middlesex District of Columbia Brazil Arizona Florida

24

slide-26
SLIDE 26

Daily Growth Rate of Daily Deaths, Past Week

Mar Apr May Jun Jul Aug Sep 2020

  • 10
  • 5

5 10 15 20 25 30 35 40

Growth Rate of Daily Deaths, Past Week (percent) = 0.010 =0.10 =0.20

New York City (only) Paris, Ile-de-France Miami Kentucky Arizona

25

slide-27
SLIDE 27

Dashboard Table (link)

Total (pm) Growth — R0 — % Ever % Infectious Deaths, t rate initial today infected peak today New York City (only) 2836

  • 2.36

0.94 28.6% 5.11% 0.02% Lombardy, Italy 1675

  • 2.20

0.20 16.8% 2.36% 0.01% Stockholm, Sweden 1465

  • 2.20

0.20 14.7% 1.68% 0.03% Madrid, Spain 1289

  • 2.22

0.20 12.9% 2.37% 0.04% Boston+Middlesex 1292 6.5% 2.08 1.90 13.1% 1.75% 0.12% District of Columbia 853 2.7% 1.82 1.21 8.6% 0.88% 0.08% Paris, France 837

  • 2.13

0.72 8.4% 1.33% 0.02% Miami 811

  • 1.71

1.04 8.9% 0.69% 0.51% London, U.K. 652 2.3% 2.11 1.10 6.6% 1.21% 0.00% Arizona 644

  • 1.24

0.82 6.8% 0.60% 0.22% United States 532

  • 1.80

1.02 5.5% 0.39% 0.15% Brazil 532

  • 1.33

1.06 5.6% 0.25% 0.23% Texas 492

  • 1.28

0.99 5.5% 0.52% 0.41% Mexico 461

  • 1.23

0.91 4.9% 0.28% 0.20% California 304

  • 1.35

1.01 3.2% 0.18% 0.16% Kentucky 251 4.3% 1.51 1.27 2.6% 0.15% 0.14% SF Bay Area 143

  • 1.16

0.96 1.5% 0.08% 0.05% Germany 111

  • 1.50

1.04 1.1% 0.15% 0.00% Israel 93 1.1% 1.17 1.06 1.0% 0.08% 0.08% Norway 49

  • 1.33

0.20 0.4% 0.08% 0.02%

26

slide-28
SLIDE 28

Dashboard Table (link)

Total (pm) Growth — R0 — % Ever % Infectious Deaths, t rate initial today infected peak today New York City (only) 2836

  • 2.36

0.94 28.6% 5.11% 0.02% Lombardy, Italy 1675

  • 2.20

0.20 16.8% 2.36% 0.01% Stockholm, Sweden 1465

  • 2.20

0.20 14.7% 1.68% 0.03% Madrid, Spain 1289

  • 2.22

0.20 12.9% 2.37% 0.04% Boston+Middlesex 1292 6.5% 2.08 1.90 13.1% 1.75% 0.12% District of Columbia 853 2.7% 1.82 1.21 8.6% 0.88% 0.08% Paris, France 837

  • 2.13

0.72 8.4% 1.33% 0.02% Miami 811

  • 1.71

1.04 8.9% 0.69% 0.51% London, U.K. 652 2.3% 2.11 1.10 6.6% 1.21% 0.00% Arizona 644

  • 1.24

0.82 6.8% 0.60% 0.22% United States 532

  • 1.80

1.02 5.5% 0.39% 0.15% Brazil 532

  • 1.33

1.06 5.6% 0.25% 0.23% Texas 492

  • 1.28

0.99 5.5% 0.52% 0.41% Mexico 461

  • 1.23

0.91 4.9% 0.28% 0.20% California 304

  • 1.35

1.01 3.2% 0.18% 0.16% Kentucky 251 4.3% 1.51 1.27 2.6% 0.15% 0.14% SF Bay Area 143

  • 1.16

0.96 1.5% 0.08% 0.05% Germany 111

  • 1.50

1.04 1.1% 0.15% 0.00% Israel 93 1.1% 1.17 1.06 1.0% 0.08% 0.08% Norway 49

  • 1.33

0.20 0.4% 0.08% 0.02%

27

slide-29
SLIDE 29

Dashboard Table (link)

Total (pm) Growth — R0 — % Ever % Infectious Deaths, t rate initial today infected peak today New York City (only) 2836

  • 2.36

0.94 28.6% 5.11% 0.02% Lombardy, Italy 1675

  • 2.20

0.20 16.8% 2.36% 0.01% Stockholm, Sweden 1465

  • 2.20

0.20 14.7% 1.68% 0.03% Madrid, Spain 1289

  • 2.22

0.20 12.9% 2.37% 0.04% Boston+Middlesex 1292 6.5% 2.08 1.90 13.1% 1.75% 0.12% District of Columbia 853 2.7% 1.82 1.21 8.6% 0.88% 0.08% Paris, France 837

  • 2.13

0.72 8.4% 1.33% 0.02% Miami 811

  • 1.71

1.04 8.9% 0.69% 0.51% London, U.K. 652 2.3% 2.11 1.10 6.6% 1.21% 0.00% Arizona 644

  • 1.24

0.82 6.8% 0.60% 0.22% United States 532

  • 1.80

1.02 5.5% 0.39% 0.15% Brazil 532

  • 1.33

1.06 5.6% 0.25% 0.23% Texas 492

  • 1.28

0.99 5.5% 0.52% 0.41% Mexico 461

  • 1.23

0.91 4.9% 0.28% 0.20% California 304

  • 1.35

1.01 3.2% 0.18% 0.16% Kentucky 251 4.3% 1.51 1.27 2.6% 0.15% 0.14% SF Bay Area 143

  • 1.16

0.96 1.5% 0.08% 0.05% Germany 111

  • 1.50

1.04 1.1% 0.15% 0.00% Israel 93 1.1% 1.17 1.06 1.0% 0.08% 0.08% Norway 49

  • 1.33

0.20 0.4% 0.08% 0.02%

28

slide-30
SLIDE 30

Dashboard Table (link)

Total (pm) Growth — R0 — % Ever % Infectious Deaths, t rate initial today infected peak today New York City (only) 2836

  • 2.36

0.94 28.6% 5.11% 0.02% Lombardy, Italy 1675

  • 2.20

0.20 16.8% 2.36% 0.01% Stockholm, Sweden 1465

  • 2.20

0.20 14.7% 1.68% 0.03% Madrid, Spain 1289

  • 2.22

0.20 12.9% 2.37% 0.04% Boston+Middlesex 1292 6.5% 2.08 1.90 13.1% 1.75% 0.12% District of Columbia 853 2.7% 1.82 1.21 8.6% 0.88% 0.08% Paris, France 837

  • 2.13

0.72 8.4% 1.33% 0.02% Miami 811

  • 1.71

1.04 8.9% 0.69% 0.51% London, U.K. 652 2.3% 2.11 1.10 6.6% 1.21% 0.00% Arizona 644

  • 1.24

0.82 6.8% 0.60% 0.22% United States 532

  • 1.80

1.02 5.5% 0.39% 0.15% Brazil 532

  • 1.33

1.06 5.6% 0.25% 0.23% Texas 492

  • 1.28

0.99 5.5% 0.52% 0.41% Mexico 461

  • 1.23

0.91 4.9% 0.28% 0.20% California 304

  • 1.35

1.01 3.2% 0.18% 0.16% Kentucky 251 4.3% 1.51 1.27 2.6% 0.15% 0.14% SF Bay Area 143

  • 1.16

0.96 1.5% 0.08% 0.05% Germany 111

  • 1.50

1.04 1.1% 0.15% 0.00% Israel 93 1.1% 1.17 1.06 1.0% 0.08% 0.08% Norway 49

  • 1.33

0.20 0.4% 0.08% 0.02%

29

slide-31
SLIDE 31

Dashboard Table (link)

Total (pm) Growth — R0 — % Ever % Infectious Deaths, t rate initial today infected peak today New York City (only) 2836

  • 2.36

0.94 28.6% 5.11% 0.02% Lombardy, Italy 1675

  • 2.20

0.20 16.8% 2.36% 0.01% Stockholm, Sweden 1465

  • 2.20

0.20 14.7% 1.68% 0.03% Madrid, Spain 1289

  • 2.22

0.20 12.9% 2.37% 0.04% Boston+Middlesex 1292 6.5% 2.08 1.90 13.1% 1.75% 0.12% District of Columbia 853 2.7% 1.82 1.21 8.6% 0.88% 0.08% Paris, France 837

  • 2.13

0.72 8.4% 1.33% 0.02% Miami 811

  • 1.71

1.04 8.9% 0.69% 0.51% London, U.K. 652 2.3% 2.11 1.10 6.6% 1.21% 0.00% Arizona 644

  • 1.24

0.82 6.8% 0.60% 0.22% United States 532

  • 1.80

1.02 5.5% 0.39% 0.15% Brazil 532

  • 1.33

1.06 5.6% 0.25% 0.23% Texas 492

  • 1.28

0.99 5.5% 0.52% 0.41% Mexico 461

  • 1.23

0.91 4.9% 0.28% 0.20% California 304

  • 1.35

1.01 3.2% 0.18% 0.16% Kentucky 251 4.3% 1.51 1.27 2.6% 0.15% 0.14% SF Bay Area 143

  • 1.16

0.96 1.5% 0.08% 0.05% Germany 111

  • 1.50

1.04 1.1% 0.15% 0.00% Israel 93 1.1% 1.17 1.06 1.0% 0.08% 0.08% Norway 49

  • 1.33

0.20 0.4% 0.08% 0.02%

30

slide-32
SLIDE 32

Repeated “Forecasts” from the past 7 days of data

– After peak, forecasts settle down. – Before that, very noisy!

31

slide-33
SLIDE 33

Guide to Graphs

  • 7 days of forecasts: Rainbow color order!

ROY-G-BIV (old to new, low to high)

  • Black=current
  • Red = oldest, Orange = second oldest, Yellow =third oldest...
  • Violet (purple) = one day earlier
  • For robustness graphs, same idea
  • Black = baseline (e.g. δ = 1.0%)
  • Red = lowest parameter value (e.g. δ = 0.8%)
  • Green = highest parameter value (e.g. δ = 1.2%)

32

slide-34
SLIDE 34

Guide to Graphs (continued)

  • R0 in subtitle:
  • Initial / Today / Final
  • “%Infect”
  • Today / t+30 / Final
  • This is the percent ever infected
  • (so fraction δ will eventually be deaths)

33

slide-35
SLIDE 35

New York City (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec 2020 10 20 30 40 50 60 70 80

Daily deaths per million people New York City (plus) R0=2.3/1.1/1.1 = 0.010 =0.05 =0.1 %Infect=24/24/24

DATA THROUGH 24-AUG-2020

34

slide-36
SLIDE 36

New York City (7 days): Cumulative Deaths per Million (Future)

Mar Apr May Jun Jul Aug Sep Oct Nov Dec 2020 500 1000 1500 2000 2500 3000

Cumulative deaths per million people New York City (plus) R0=2.3/1.1/1.1 = 0.010 =0.05 =0.1 %Infect=24/24/24

DATA THROUGH 24-AUG-2020

35

slide-37
SLIDE 37

New York excl NYC (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 1 2 3 4 5 6 7 8

Daily deaths per million people New York excluding NYC R0=1.8/0.3/0.3 = 0.010 =0.05 =0.1 %Infect= 4/ 4/ 4

DATA THROUGH 24-AUG-2020

36

slide-38
SLIDE 38

New York excl NYC (7 days): Cumulative Deaths per Million (Future)

Mar Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 1 2 4 8 16 32 64 128 256 512 1024 2048 4096

Cumulative deaths per million people New York excluding NYC R0=1.8/0.3/0.3 = 0.010 =0.05 =0.1 %Infect= 4/ 4/ 4

New York City Italy 37

slide-39
SLIDE 39

Stylized Experiences to Keep in Mind Total (pm) Peak Daily Deaths, t Deaths (pm) New York City (only) 2800 80 Paris = London = Washington DC 800 20 NY excl NYC = Atlanta 400 5

38

slide-40
SLIDE 40

Lombardy (7 days): Daily Deaths per Million People

Mar Apr May Jun Jul Aug Sep Oct Nov Dec 2020 5 10 15 20 25 30 35 40 45

Daily deaths per million people Lombardy, Italy R0=2.2/0.2/0.2 = 0.010 =0.05 =0.1 %Infect=17/17/17

DATA THROUGH 24-AUG-2020

39

slide-41
SLIDE 41

Italy (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec 2020 2 4 6 8 10 12 14

Daily deaths per million people Italy R0=2.0/1.0/1.0 = 0.010 =0.05 =0.1 %Infect= 6/ 6/ 6

DATA THROUGH 24-AUG-2020

40

slide-42
SLIDE 42

Spain (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec 2020 2 4 6 8 10 12 14 16 18 20

Daily deaths per million people Spain R0=2.1/0.2/0.2 = 0.010 =0.05 =0.1 %Infect= 6/ 6/ 6

DATA THROUGH 24-AUG-2020

41

slide-43
SLIDE 43

Paris (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec 2020 2 4 6 8 10 12 14 16 18

Daily deaths per million people Paris, Ile-de-France R0=2.1/1.3/1.1 = 0.010 =0.05 =0.1 %Infect= 6/ 7/10

DATA THROUGH 24-AUG-2020

42

slide-44
SLIDE 44

London (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec 2020 5 10 15 20 25

Daily deaths per million people London, U.K. R0=2.1/1.1/1.1 = 0.010 =0.05 =0.1 %Infect= 7/ 7/ 7

DATA THROUGH 24-AUG-2020

43

slide-45
SLIDE 45

U.K. (7 days): Daily Deaths per Million

Apr May Jun Jul Aug Sep Oct Nov Dec 2020 2 4 6 8 10 12 14

Daily deaths per million people United Kingdom R0=2.0/0.8/0.8 = 0.010 =0.05 =0.1 %Infect= 6/ 6/ 6

DATA THROUGH 24-AUG-2020

44

slide-46
SLIDE 46

Stockholm (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec 2020 5 10 15 20 25 30 35 40

Daily deaths per million people Stockholm, Sweden R0=2.2/0.2/0.2 = 0.010 =0.05 =0.1 %Infect=15/15/15

DATA THROUGH 24-AUG-2020

45

slide-47
SLIDE 47

Boston (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 5 10 15 20 25 30 35 40 45

Daily deaths per million people Boston+Middlesex R0=2.1/1.8/1.3 = 0.010 =0.05 =0.1 %Infect=14/18/25

DATA THROUGH 24-AUG-2020

46

slide-48
SLIDE 48

District of Columbia (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 2 4 6 8 10 12 14 16 18

Daily deaths per million people District of Columbia R0=1.8/1.2/1.1 = 0.010 =0.05 =0.1 %Infect= 9/10/12

DATA THROUGH 24-AUG-2020

47

slide-49
SLIDE 49

Miami (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 2 4 6 8 10 12 14 16 18

Daily deaths per million people Miami R0=1.7/1.1/1.2 = 0.010 =0.05 =0.1 %Infect=10/12/19

DATA THROUGH 24-AUG-2020

48

slide-50
SLIDE 50

Miami (7 days): Cumulative Deaths per Million (Future)

Mar Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 1 2 4 8 16 32 64 128 256 512 1024 2048 4096

Cumulative deaths per million people Miami R0=1.7/1.1/1.2 = 0.010 =0.05 =0.1 %Infect=10/12/19

New York City Italy 49

slide-51
SLIDE 51

Arizona (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 2 4 6 8 10 12

Daily deaths per million people Arizona R0=1.2/0.9/1.1 = 0.010 =0.05 =0.1 %Infect= 7/ 8/ 9

DATA THROUGH 24-AUG-2020

50

slide-52
SLIDE 52

Arizona (7 days): Daily Deaths with α = 0

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 2 4 6 8 10 12

Daily deaths per million people Arizona R0=1.2/0.8/0.8 = 0.010 =0.00 =0.1 %Infect= 7/ 7/ 8

DATA THROUGH 24-AUG-2020

51

slide-53
SLIDE 53

Mexico (7 days): Daily Deaths per Million People

May 2020 Jun 2020 Jul 2020 Aug 2020 Sep 2020 Oct 2020 Nov 2020 Dec 2020 Jan 2021 Feb 2021 1 2 3 4 5 6 7

Daily deaths per million people Mexico R0=1.2/0.9/1.1 = 0.010 =0.05 =0.1 %Infect= 5/ 6/ 8

DATA THROUGH 24-AUG-2020

52

slide-54
SLIDE 54

Mexico (7 days): Cumulative Deaths per Million (Future)

Mar 2020 Apr 2020 May 2020 Jun 2020 Jul 2020 Aug 2020 Sep 2020 Oct 2020 Nov 2020 Dec 2020 Jan 2021 Feb 2021 1 2 4 8 16 32 64 128 256 512 1024 2048 4096

Cumulative deaths per million people Mexico R0=1.2/0.9/1.1 = 0.010 =0.05 =0.1 %Infect= 5/ 6/ 8

New York City Italy 53

slide-55
SLIDE 55

Peru (7 days): Daily Deaths per Million People

May Jun Jul Aug Sep Oct Nov Dec Jan 2020 2 4 6 8 10 12

Daily deaths per million people Peru R0=1.5/1.1/1.2 = 0.010 =0.05 =0.1 %Infect= 9/10/14

DATA THROUGH 24-AUG-2020

54

slide-56
SLIDE 56

Brazil (7 days): Daily Deaths per Million People

May Jun Jul Aug Sep Oct Nov Dec Jan 2020 1 2 3 4 5 6

Daily deaths per million people Brazil R0=1.3/1.1/1.1 = 0.010 =0.05 =0.1 %Infect= 6/ 7/11

DATA THROUGH 24-AUG-2020

55

slide-57
SLIDE 57

Texas (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 2 4 6 8 10 12

Daily deaths per million people Texas R0=1.3/1.0/1.1 = 0.010 =0.05 =0.1 %Infect= 6/ 8/13

DATA THROUGH 24-AUG-2020

56

slide-58
SLIDE 58

SF Bay Area (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6

Daily deaths per million people SF Bay Area R0=1.2/1.0/1.0 = 0.010 =0.05 =0.1 %Infect= 2/ 2/ 2

DATA THROUGH 24-AUG-2020

57

slide-59
SLIDE 59

Los Angeles (7 days): Daily Deaths per Million People

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

Daily deaths per million people Los Angeles R0=1.5/1.1/1.1 = 0.010 =0.05 =0.1 %Infect= 5/ 7/10

DATA THROUGH 24-AUG-2020

58

slide-60
SLIDE 60

Reopening and Herd Immunity

59

slide-61
SLIDE 61

Percent Ever Infected would be very informative

— Percent Ever Infected (today) — δ = 0.5% δ = 1.0% δ = 1.2% New York City (only) 57 29 24 Lombardy, Italy 34 17 14 Stockholm, Sweden 29 15 12 Madrid, Spain 26 13 11 Boston+Middlesex 27 14 11 Philadelphia 23 11 9 Belgium 18 9 7 District of Columbia 18 9 7 Paris, France 17 8 7 Miami 19 10 8 London, U.K. 13 7 5 Spain 12 6 5 Arizona 14 7 6 Italy 12 6 5 United States 11 6 5 Texas 12 6 5 Mexico 10 5 4 New York excluding NYC 8 4 3 Houston (Harris Co.) 10 5 4 Kentucky 6 3 2 SF Bay Area 3 2 1 Israel 2 1 1 Norway 1

60

slide-62
SLIDE 62

Herd Immunity

  • How far can we relax social distancing?
  • Let s(t) = S(t)/N = the fraction still susceptible
  • The disease will die out as long as

R0(t)s(t) < 1

  • That is, if the “new” R0 is smaller than 1/s(t)
  • Today’s infected people infect fewer than 1 person on

average

  • We can relax social distancing to raise R0(t) to 1/s(t)

61

slide-63
SLIDE 63

Herd Immunity and Opening the Economy? δ = 1.0%

Percent R0(t+30) Percent Susceptible with no way back R0 R0(t) t+30

  • utbreak

to normal New York City (only) 2.4 1.0 71.4 1.4 31.8 Lombardy, Italy 2.2 0.2 83.2 1.2 49.9 Stockholm, Sweden 2.2 0.2 85.3 1.2 48.2 Madrid, Spain 2.2 0.2 87.1 1.1 46.0 Chicago 1.9 1.3 88.9 1.1

  • 40.0

Belgium 2.1 1.1 91.0 1.1 2.2 District of Columbia 1.8 1.2 90.5 1.1

  • 13.3

Paris, France 2.1 0.7 91.5 1.1 26.0 Miami 1.7 1.1 87.7 1.1 11.4 London, U.K. 2.1 1.1 93.4 1.1

  • 2.7

United Kingdom 2.0 0.8 93.9 1.1 19.2 Italy 2.0 1.0 94.1 1.1 1.8 Sweden 1.8 1.3 93.8 1.1

  • 61.5

United States 1.8 1.0 93.5 1.1 5.6 Brazil 1.3 1.1 92.7 1.1 6.8 France 2.0 1.2 95.2 1.1

  • 18.4

Mexico 1.2 0.9 94.1 1.1 42.2 California 1.4 1.0 95.6 1.0 8.0 Kentucky 1.5 1.2 95.3 1.0

  • 57.9

SF Bay Area 1.2 1.0 98.2 1.0 25.5 Germany 1.5 1.0 98.9 1.0

  • 5.8

Israel 1.2 1.0 98.3 1.0

  • 26.0

Norway 1.3 0.2 99.6 1.0 71.1

62

slide-64
SLIDE 64

Simulations of Re-Opening

  • Begin with the basic estimates shown already
  • Different policies are then adopted starting around July 20
  • Black: assumes R0(today) remains in place forever
  • Red: assumes R0(suppress)= 1/s(today)
  • Green: we move 25% of the way from R0(today) back to

initial R0 = “normal”

  • Purple: we move 50% of the way from R0(today) back to

initial R0 = “normal”

  • We assume these R0 values adjust to daily deaths via α
  • Each daily death reduces R0 by 5 percent

63

slide-65
SLIDE 65

New York City: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people New York City (plus) R0(t)=1.1, R0(suppress)=1.3, R0(25/50)=1.4/1.7, = 0.010, =0.05

64

slide-66
SLIDE 66

New York excluding NYC: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people New York excluding NYC R0(t)=0.3, R0(suppress)=1.0, R0(25/50)=0.7/1.1, = 0.010, =0.05

65

slide-67
SLIDE 67

California: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people California R0(t)=1.0, R0(suppress)=1.0, R0(25/50)=1.3/1.5, = 0.010, =0.05

66

slide-68
SLIDE 68

Miami: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people Miami R0(t)=1.1, R0(suppress)=1.1, R0(25/50)=1.3/1.5, = 0.010, =0.05

67

slide-69
SLIDE 69

Washington, DC: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people District of Columbia R0(t)=1.2, R0(suppress)=1.1, R0(25/50)=1.4/1.6, = 0.010, =0.05

68

slide-70
SLIDE 70

Mexico: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people Mexico R0(t)=0.9, R0(suppress)=1.1, R0(25/50)=1.2/1.5, = 0.010, =0.05

69

slide-71
SLIDE 71

Stockholm: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people Stockholm, Sweden R0(t)=0.2, R0(suppress)=1.2, R0(25/50)=0.7/1.2, = 0.010, =0.05

70

slide-72
SLIDE 72

Brazil: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people Brazil R0(t)=1.1, R0(suppress)=1.1, R0(25/50)=1.3/1.5, = 0.010, =0.05

71

slide-73
SLIDE 73

Houston: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people Houston (Harris Co.) R0(t)=0.9, R0(suppress)=1.1, R0(25/50)=1.2/1.5, = 0.010, =0.05

72

slide-74
SLIDE 74

Arizona: Re-Opening

Apr May Jun Jul Aug Sep Oct Nov Dec Jan 2020 10 20 30 40 50 60 70 80 90 100

Daily deaths per million people Arizona R0(t)=0.9, R0(suppress)=1.1, R0(25/50)=1.2/1.4, = 0.010, =0.05

73

slide-75
SLIDE 75

Macroeconomic Outcomes

74

slide-76
SLIDE 76

Economic Policy Trade Off, Holding Health Policy and Luck Constant

COVID DEATHS PER MILLION PEOPLE GDP LOSS (PERCENT)

Shut down economy Keep economy open

75

slide-77
SLIDE 77

Health Policy Decisions and Luck Can Shift the Tradeoff

COVID DEATHS PER MILLION PEOPLE GDP LOSS (PERCENT)

Good health policy

  • r good luck

Bad health policy

  • r bad luck

76

slide-78
SLIDE 78

Putting together...

COVID DEATHS PER MILLION PEOPLE GDP LOSS (PERCENT)

Shut down economy Keep economy open Good health policy

  • r good luck

Bad health policy

  • r bad luck

77

slide-79
SLIDE 79

International Covid Deaths and Lost GDP, 2020

100 200 300 400 500 600 700 800 900

  • 5

5 10 15 20 25 30

United States Austria Belgium Chile Czechia Estonia France Germany Indonesia Italy Japan Korea, South Mexico Netherlands Philippines Portugal Singapore Slovakia Spain Sweden Taiwan United Kingdom

COVID DEATHS PER MILLION PEOPLE GDP LOSS (PERCENT)

78

slide-80
SLIDE 80

Google Activity Tracker: International Evidence

Feb Mar Apr May Jun Jul Aug Sep 2020

  • 100
  • 80
  • 60
  • 40
  • 20

20

Percent change relative to baseline

Italy U.S. Spain U.K. Germany

79

slide-81
SLIDE 81

Google Activity Tracker for Key Global Cities

Feb Mar Apr May Jun Jul Aug Sep 2020

  • 100
  • 80
  • 60
  • 40
  • 20

20

Percent change relative to baseline

Lombardy NYC Madrid London Stockholm Seoul Tokyo Paris

80

slide-82
SLIDE 82

Cumulative Google Activity and Lost GDP

50 100 150 200 250 300

  • 5

5 10 15 20 25 30

United States Austria Belgium Chile Czechia Estonia France Germany Indonesia Italy Japan Korea, South Mexico Netherlands Portugal Singapore Slovakia Spain Sweden Taiwan United Kingdom

OLS Slope = 0.079

  • Std. Err. = 0.017

R2 = 0.52

GOOGLE CUMULATIVE REDUCED ACTIVITY GDP LOSS (PERCENT)

81

slide-83
SLIDE 83

Covid Deaths and Cumulative Google Activity

100 200 300 400 500 600 700 800 900 50 100 150 200 250 300

United States Italy Germany United Kingdom Sweden Norway Japan Korea, South Indonesia Taiwan Mexico Portugal Austria Belgium Spain France Denmark Switzerland

COVID DEATHS PER MILLION PEOPLE CUMULATIVE REDUCED ACTIVITY (GOOGLE)

82

slide-84
SLIDE 84

Global Cities: Covid Deaths and Cumulative Google Activity

500 1000 1500 2000 2500 3000 50 100 150 200 250 300 350 400

New York City Lombardy, Italy London Paris Madrid, Spain Stockholm Tokyo, Japan Seoul, Korea Boston Los Angeles SF Bay Area Miami Chicago Houston

COVID DEATHS PER MILLION PEOPLE CUMULATIVE REDUCED ACTIVITY (GOOGLE) 83

slide-85
SLIDE 85

U.S. States: Covid Deaths and Cumulative Google Activity

200 400 600 800 1000 1200 1400 1600 1800 130 140 150 160 170 180 190 200 210 220 230

CA TX FL NY PA IL OH GA NC NJ MI VA WA AZ MA

COVID DEATHS PER MILLION PEOPLE CUMULATIVE REDUCED ACTIVITY (GOOGLE)

84

slide-86
SLIDE 86

Summary

COVID DEATHS GDP LOSS

California

[lucky? too tight?]

New York City Lombardy United Kingdom Madrid

[unlucky? bad policy?]

Sweden

[unlucky? too loose?]

Germany Norway Japan, S. Korea Kentucky

[lucky? good policy?]

85

slide-87
SLIDE 87

Global Cities: Monthly Evolution from March to August

1 2 4 8 16 32 64 128 256 512 1024 2048 10 20 30 40 50 60 70 80 90

New York City Lombardy, Italy London Paris Madrid, Spain Stockholm Tokyo, Japan Seoul, Korea

COVID DEATHS PER MILLION PEOPLE REDUCED ACTIVITY (GOOGLE) (PERCENT) 86

slide-88
SLIDE 88

Global Cities: Monthly Evolution from March to August

1 2 4 8 16 32 64 128 256 512 1024 2048 10 20 30 40 50 60 70 80

Boston Los Angeles SF Bay Area Miami Chicago Houston

COVID DEATHS PER MILLION PEOPLE REDUCED ACTIVITY (GOOGLE) (PERCENT)

87

slide-89
SLIDE 89

Conclusions and Questions Speculations based on model, we are not epidemiologists

  • Things I wonder about
  • What is the distribution of outcomes for young people?

(hospitalizations? long-term effects?)

  • What if every at risk person had an N95 mask?
  • What did Tokyo and Seoul do that we should learn from?
  • Next wave?
  • Paper/saliva tests

Our dashboard contains 30+ pages of results for each of 100 cities, states, and countries

88