Skip to content

Commit

Permalink
fix weightings on figure S1
Browse files Browse the repository at this point in the history
  • Loading branch information
Rob committed Apr 4, 2020
1 parent 46cf42c commit e0e6517
Show file tree
Hide file tree
Showing 3 changed files with 21 additions and 0 deletions.
Binary file modified FigureS1.pdf
View file
Binary file not shown.
Binary file modified FigureS1.png
View file
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
21 changes: 21 additions & 0 deletions R/nyt_county.R
View file
Original file line number Diff line number Diff line change
Expand Up @@ -408,6 +408,7 @@ ua1 = ggplot(cv_ex %>% filter(date == last_day) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_tavg_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1) +
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_tavg_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.25)) +
xlab('Average Temperature (C)') +
ylab('Density') +
ggtitle(last_day)
Expand All @@ -418,6 +419,7 @@ ua2 = ggplot(cv_ex %>% filter(date == last_day - 14) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_tavg_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_tavg_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.25)) +
xlab('Average Temperature (C)') +
ylab('Density') +
ggtitle(last_day - 14)
Expand All @@ -427,6 +429,7 @@ ua3 = ggplot(cv_ex %>% filter(date == last_day - 28) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_tavg_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_tavg_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.25)) +
xlab('Average Temperature (C)') +
ylab('Density') +
ggtitle(last_day - 28)
Expand All @@ -435,27 +438,31 @@ ub1 = ggplot(cv_ex %>% filter(date == last_day) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_tmin_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_tmin_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.25)) +
xlab('Min Temperature (C)') +
ylab('Density')

ub2 = ggplot(cv_ex %>% filter(date == last_day - 14) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_tmin_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_tmin_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.25)) +
xlab('Min Temperature (C)') +
ylab('Density')

ub3 = ggplot(cv_ex %>% filter(date == last_day - 28) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_tmin_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_tmin_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.25)) +
xlab('Min Temperature (C)') +
ylab('Density')

ubb1 = ggplot(cv_ex %>% filter(date == last_day) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_tmax_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_tmax_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.25)) +
xlab('Max Temparature (C)') +
ylab('Density')

Expand All @@ -465,13 +472,15 @@ ubb2 = ggplot(cv_ex %>% filter(date == last_day - 14) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_tmax_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_tmax_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.25)) +
xlab('Max Temperature (C)') +
ylab('Density')

ubb3 = ggplot(cv_ex %>% filter(date == last_day - 28) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_tmax_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_tmax_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.25)) +
xlab('Max Temperature (C)') +
ylab('Density')

Expand All @@ -480,6 +489,7 @@ uc1 = ggplot(cv_ex %>% filter(date == last_day) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_srad_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_srad_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.0008)) +
xlab('Solar Radiation (kJ m^-2 / day)') +
ylab('Density')

Expand All @@ -489,20 +499,23 @@ uc2 = ggplot(cv_ex %>% filter(date == last_day - 14) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_srad_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_srad_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.0008)) +
xlab('Solar Radiation (kJ m^-2 / day)') +
ylab('Density')

uc3 = ggplot(cv_ex %>% filter(date == last_day - 28) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_srad_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_srad_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.0008)) +
xlab('Solar Radiation (kJ m^-2 / day)') +
ylab('Density')

ucc1 = ggplot(cv_ex %>% filter(date == last_day) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_prec_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_prec_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.03)) +
xlab('Precipitation (mm)') +
ylab('Density')

Expand All @@ -512,13 +525,15 @@ ucc2 = ggplot(cv_ex %>% filter(date == last_day - 14) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_prec_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_prec_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.03)) +
xlab('Precipitation (mm)') +
ylab('Density')

ucc3 = ggplot(cv_ex %>% filter(date == last_day - 28) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_prec_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_prec_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 0.03)) +
xlab('Precipitation (mm)') +
ylab('Density')

Expand All @@ -527,6 +542,7 @@ ucd1 = ggplot(cv_ex %>% filter(date == last_day) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_wind_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_wind_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 1)) +
xlab('Wind Speed (m/s)') +
ylab('Density')

Expand All @@ -536,13 +552,15 @@ ucd2 = ggplot(cv_ex %>% filter(date == last_day - 14) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_wind_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_wind_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 1)) +
xlab('Wind Speed (m/s)') +
ylab('Density')

ucd3 = ggplot(cv_ex %>% filter(date == last_day - 28) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_wind_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_wind_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 1)) +
xlab('Wind Speed (m/s)') +
ylab('Density')

Expand All @@ -551,6 +569,7 @@ ud1 = ggplot(cv_ex %>% filter(date == last_day) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_vapr_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_vapr_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 5.5)) +
xlab('Water Vapor Pressure (kPa)') +
ylab('Density')

Expand All @@ -560,13 +579,15 @@ ud2 = ggplot(cv_ex %>% filter(date == last_day - 14) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_vapr_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_vapr_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 5.5)) +
xlab('Water Vapor Pressure (kPa)') +
ylab('Density')

ud3 = ggplot(cv_ex %>% filter(date == last_day - 28) %>% filter(!is.na(pop))) +
geom_density(aes(x=wc2.1_2.5m_vapr_03, weight=cases/sum(cases, na.rm=T)), colour='darkred', fill='darkred', alpha=0.1)+
geom_density(data=all_ex2, aes(x=wc2.1_2.5m_vapr_03, weight=pop/sum(pop, na.rm=T)), colour='darkblue', fill='darkblue', alpha=0.1) +
theme_minimal() +
ylim(c(0, 5.5)) +
xlab('Water Vapor Pressure (kPa)') +
ylab('Density')

Expand Down

0 comments on commit e0e6517

Please sign in to comment.