Baqaee, Farhi, Mina, Stock BPEA (2020)
Replication files

Stata code:
NPI_setup_v5.do -- Reads in data and does preliminary transformations, first step in _control and _estn programs
NLSEIRD_PID_v12.do -- SIR model, both for simulation and estimation, called by _control and _estn programs
SEIRD_control_v9.do -- Specifies counterfactuals and does simulations, given estimated model parameters.
SEIRD_estn_v8.do -- Estimation of model parameters
SEIRD_SEs_v3.do -- Standard errors for R-effective and nonparametric smoother estimate of R-effective (produces Fig. 4 in paper)
SEIRD_filecopy.do -- Copies output files into BPEA file folder for figures
SEIRD_generate_figs.do -- Creates figures in BPEA paper

Folder structure
- .do files should be in main folder
- all other folders are subfolders off the main folder for reading in data or writing output files or figures

Output file naming convention:
Example:
fast-SD3_100_s10_Q5_7_W_unrate_fan_SE    NOTE: this is the base second wave case in Figure 1L

Abbreviations:
fast = fast or slow governor
SD3 = scenario. SD3 = baseline fast reopening. Rest of scenario abbreviations are in SEIRD_control_v9.do
_100 = value of coefficient k_theta in paper, times 100 = coefficient on GDP-to-Risk index (times 100), see eq (19) in paper
_s10 = fraction of kids who return to school in the fall * 10 so s10 means 100% of kids return to school, s0 means none return
_Q5 = quarantine rate in percent (Q5 means 5% quarantine)
_7 = 0.7% population IFR
_W_unrate_fan = fan (shaded) chart with weekly deaths W and monthly unemployment rate
_SE = type of figure file:
	_SE = shaded areas are determined by standard errors computed using delta method along the time path. 
	_IF = shaded fans come from varying the IFR from 0.4 to 1.0 instead of from estimation uncertainty.
	_age = age profile plots (BPEA paper Fig. 5)
