Skip to content

Instantly share code, notes, and snippets.

@jdingel

jdingel/FJS_DEX.do

Created May 23, 2020 03:21
Show Gist options
  • Save jdingel/cbca66b5d04ce8fb4fcc7246cf0f4559 to your computer and use it in GitHub Desktop.
Save jdingel/cbca66b5d04ce8fb4fcc7246cf0f4559 to your computer and use it in GitHub Desktop.
Download ZIP
Use DEX data to produce FJS Figure 1
Raw
FJS_DEX.do
import delimited using "https://raw.githubusercontent.com/COVIDExposureIndices/COVIDExposureIndices/master/dex_data/state_dex.csv", clear
keep state date dex num_devices dex_a num_devices_a
gen date_stata = date(date,"YMD")
format %td date_stata
gen baseline_date = mod(date_stata - 21934,7)
bys state baseline_date (date_stata): egen baseline_dex_a = total(dex_a * inrange(date_stata,21974,21982))
by state baseline_date (date_stata): egen baseline_denom = total(inrange(date_stata,21974,21982))
gen dex_a_r = dex_a / (baseline_dex_a / baseline_denom)
collapse (min) dex_a_r_min = dex_a_r (max) dex_a_r_max = dex_a_r (p10) dex_a_r_10 = dex_a_r (p25) dex_a_r_25 = dex_a_r (p50) dex_a_r_50 = dex_a_r (p75) dex_a_r_75 = dex_a_r (p90) dex_a_r_90 = dex_a_r (p5) dex_a_r_05 = dex_a_r (p95) dex_a_r_95 = dex_a_r, by(date_stata)
gen byte ins = (inrange(date_stata, 21934+30,22039-21)==1)
twoway (line dex_a_r_min date_stata, lcol(gray)) (line dex_a_r_10 date_stata, lcol(midblue)) (line dex_a_r_50 date_stata, lcol(blue)) (line dex_a_r_90 date_stata, lcol(midblue)) (line dex_a_r_max date_stata, lcol(gray)) ///
, graphregion(color(white))
twoway (line dex_a_r_min date_stata if ins==1, lcol(gray)) (line dex_a_r_10 date_stata if ins==1, lcol(midblue)) (line dex_a_r_50 date_stata if ins==1, lcol(blue)) (line dex_a_r_90 date_stata if ins==1, lcol(midblue)) (line dex_a_r_max date_stata if ins==1, lcol(gray)) ///
, graphregion(color(white))
graph export "~/Desktop/DEX_FJS.png", replace
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment