BPEA Appendix Homebase Worker’s Hours

print(paste("Started at", Sys.time()))

## [1] "Started at 2020-07-19 12:29:18"

date_min <- as.Date("2019-10-01")

#-------------------------------------------------------------------------------
# Share of hours from workers in base period

# Multiple Years
fout <- "worker_hours_shours_base_worker"
dt_share <- rbind(readRDS(paste0(dir_exp,fout,".rds"))[,year:=2020],
                  readRDS(paste0(dir_exp,fout,"_2019.rds"))[,year:=2019],
                  readRDS(paste0(dir_exp,fout,"_2018.rds"))[,year:=2018])
dt_share[year==2019, date:=date+as.integer(as.Date("2020-01-19")-as.Date("2019-01-20"))]
dt_share[year==2018, date:=date+as.integer(as.Date("2020-01-19")-as.Date("2018-01-21"))]
ggplot(dt_share[user_base==1 & date>=as.Date("2020-01-19") & date<=max(dt_share[year==2020,date])],aes(x=date,y=s_hours,color=factor(year))) + ggtheme +
  geom_line(size=1.2) +
  scale_color_brewer(palette = "Set1") +
  coord_cartesian(ylim=c(0.5,1)) +
  labs(x="Date",y="Share of Hours from Workers\nin Base Period",color="Year")

# ggsave(paste0(dir_exp, fout,"_181920.png"), scale=1, width = 9, height = 6, dpi = 300, units = "in")
ggsave(paste0(dir_ofig,"PA_worker_hrs_base_worker.png"), scale=1, width = 9, height = 6, dpi = 300, units = "in")
fwrite(dt_share[user_base==1 & date>=as.Date("2020-01-19") & date<=max(dt_share[year==2020,date])],
       paste0(dir_ofigd,"PA_worker_hrs_base_worker.csv"))

fout <- "worker_hours_shours_base_worker"
dt_share <- rbind(readRDS(paste0(dir_exp,fout,".rds"))[,year:="Worker"],readRDS(paste0(dir_exp,fout,"f.rds"))[,year:="Worker-Firm Pair"])
# dt_share[year==2019, date:=date+as.integer(as.Date("2020-01-19")-as.Date("2019-01-20"))]
ggplot(dt_share[user_base==1 & date>=as.Date("2020-01-19")],aes(x=date,y=s_hours,color=factor(year))) + ggtheme +
  geom_line(size=1.2) +
  scale_color_brewer(palette = "Set1") +
  coord_cartesian(ylim=c(0.5,1)) +
  labs(x="Date",y="Share of Hours by Wokers in Base Period",color="Definition")

# ggsave(paste0(dir_exp, fout,"f.png"), scale=1, width = 9, height = 6, dpi = 300, units = "in")

#-------------------------------------------------------------------------------
# Distribution of hours worked for those with positive hours

fout <- "worker_hours_hours_dist_active_worker"
dt_dist <- readRDS(paste0(dir_exp,fout,".rds"))
ggplot(dt_dist[week>=f_wk(date_min)],aes(x=weekd)) + ggtheme +
  geom_line(aes(y=hours,color=factor(stats)), size=1.2) +
  geom_point(aes(y=hours,color=factor(stats),shape=factor(stats)), size=1.8) +
  scale_color_brewer(palette = "Set2") +
  labs(x="Date",y="Hours",color="Percentile",shape="Percentile") +
  theme(legend.position = "bottom")

# ggsave(paste0(dir_exp, fout,".png"), scale=1, width = 9, height = 6, dpi = 300, units = "in")
ggsave(paste0(dir_ofig,"PA_worker_hrs_dist.png"), scale=1, width = 9, height = 6, dpi = 300, units = "in")
fwrite(dt_dist[week>=f_wk(date_min)],paste0(dir_ofigd,"PA_worker_hrs_dist.csv"))