Code used in https://arxiv.org/abs/2006.05532
Project description
Masks and COVID-19: a causal framework for imputing value to public-health interventions
Code to reproduce Masks and COVID-19.
This is a refactored version of the original code.
Install
pip install babino2020masks
How to use
Gather data
ny = API(api_settings['NYS'][:2])
df = ny.get_all_data_statewide()
ax = plot_data_and_fit(df, 'Date', 'Odds', None, None, None, figsize=(10, 7))
ax.set_title(f'{df.tail(1).Date[0]:%B %d, %Y}, Positivity Odds:{df.tail(1).Odds[0]:2.3}');
Fit the model
sdf = df.loc[df.Date<='15-05-2020'].copy()
lics = LassoICSelector(sdf['Odds'], 'bic')
lics.fit_best_alpha()
Positivity Odds in NYS
sdf['Fit'], sdf['Odds_l'], sdf['Odds_u'] = lics.odds_hat_l_u()
ax = plot_data_and_fit(sdf, 'Date', 'Odds', 'Fit', 'Odds_l', 'Odds_u', figsize=(10, 7))
instantaneous reproduction number, $R_t$
sdf['R'], sdf['Rl'], sdf['Ru'] = lics.rt()
ax = plot_data_and_fit(sdf, 'Date', None, 'R', 'Rl', 'Ru', figsize=(10, 7), logy=False, palette=[palette[1],palette[1]])
Counterfactual Scenario without Masks
sdf['Cf. Odds'], sdf['cf_odds_l'], sdf['cf_odds_u'] = lics.counterfactual()
ax = plot_data_and_fit(sdf, 'Date', 'Odds', 'Fit', 'Odds_l', 'Odds_u', figsize=(10, 7))
plot_data_and_fit(sdf, 'Date', None, 'Cf. Odds', 'cf_odds_l', 'cf_odds_u', palette=[palette[2],palette[2]], ax=ax);
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
babino2020masks-0.0.5.tar.gz
(13.2 kB
view hashes)
Built Distribution
Close
Hashes for babino2020masks-0.0.5-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 01412386767c9c70a2c0d55aa477f95c2c93baff395ca7d918889fd87f76ba21 |
|
MD5 | 6a50acb368be58c50c2a62690cf0731a |
|
BLAKE2b-256 | ec0c4daf406720e40494260cb7bd746d4a083cf4fc52ba96b9917a3e76c942a4 |