Weight Control SDID

Function Plot

NM <- read_excel("V_NAME_All.xlsx")

with.overlay = function(est, s) { attr(est,'overlay') = s; est }

estimators = function(sdid,sc,s) {
  estimator.list = list(with.overlay(sdid, s), sc, sdid)
  names(estimator.list)=c('SDID+fixed effect', 'synth. control', 'synth. diff-in-diff') #some of this need to remove but i still dont get it
  estimator.list
}

plot.estimators = function(ests, alpha.multiplier) {
  p = synthdid_plot(ests, se.method='none',
                    alpha.multiplier=alpha.multiplier, facet=rep(1,length(ests)),
                    trajectory.linetype = 1, effect.curvature=-.4,
                    trajectory.alpha=.5, effect.alpha=.5, diagram.alpha=1)
  suppressMessages(p + scale_alpha(range=c(0,1), guide='none'))
}

TNF

TNF <- read_excel("TNF_State.xls", sheet = "TNF")

#DATE CLEANS
TNF$DATE<-as.Date(TNF$DATE)
TNF<-arrange(TNF,STATE)
TNF<-filter(TNF,DATE>="2000-01-01")

#VariablTNF Clean
TNF<-select(TNF,STATE,TNF,DATE)

dumy_date=seq.Date(from = as.Date("2000-01-01"),to = as.Date("2022-01-01"),by = "month")
Fix_df<-data_frame(STATE=rep(unique(TNF$STATE),length(dumy_date)))
Fix_df<-arrange(Fix_df,STATE)
Fix_df$DATE<-rep(dumy_date,length(unique(TNF$STATE)))

TNF<-left_join(Fix_df,TNF,by=c("STATE","DATE"))

TNF<-TNF %>% # Remove All NA STATE
  group_by(STATE) %>%  
  filter(DATE<="2008-01-01",!any(is.na(TNF)))

TNF$Treat<-ifelse(between(TNF$DATE,as.Date("2005-08-01"),as.Date("2022-06-01"))
                        &TNF$STATE=="LA",1,0)
TNF<-as.data.frame(TNF)
TNF$STATE<-factor(TNF$STATE)
TNF$L_TNF<-log(TNF$TNF)

# SDID PAN
matx<-panel.matrices(TNF,unit = 1,time = 2,outcome = 5,treatment = 4)

sdid.TNF<-synthdid_estimate(matx$Y, matx$N0, matx$T0)
sc.TNF = sc_estimate(matx$Y, matx$N0, matx$T0)

wg<-summary(sdid.TNF)$controls
wg<-data.frame(State=rownames(wg),Controls=wg[,1])
rownames(wg)<-NULL

wg<-left_join(statepop,wg,by=c("abbr"="State"))
wg

## # A tibble: 51 × 5
##    fips  abbr  full                 pop_2015 Controls
##    <chr> <chr> <chr>                   <dbl>    <dbl>
##  1 01    AL    Alabama               4858979    0.023
##  2 02    AK    Alaska                 738432   NA    
##  3 04    AZ    Arizona               6828065   NA    
##  4 05    AR    Arkansas              2978204    0.035
##  5 06    CA    California           39144818   NA    
##  6 08    CO    Colorado              5456574   NA    
##  7 09    CT    Connecticut           3590886   NA    
##  8 10    DE    Delaware               945934   NA    
##  9 11    DC    District of Columbia   672228   NA    
## 10 12    FL    Florida              20271272   NA    
## # ℹ 41 more rows

p4 = plot.estimators(estimators(sdid = sdid.TNF,sc = sc.TNF,1),   alpha.multiplier=c(1, 1,1))

plot.theme = theme(legend.position="none", legend.direction='horizontal',
                   legend.key=element_blank(), legend.background=element_blank(), 
                   plot.title = element_text(hjust = 0.5,size = 9),axis.text = element_text(size=7),axis.title = element_text(size=7))

#tit=ggtitle(label = NM[NM$VAR=="TNF",]$Variable) #Change this

pl.TNF<-p4 + plot.theme+ #xlab("Year") + 
  ylab("Total\n Nonfarm")+
  scale_x_continuous(labels = c("2000-02","2002-11","2005-08","2008-05"))

pl.TNF

plot_usmap(regions = "state",data = wg,values = "Controls")+
  scale_fill_gradient2(na.value = "white")+ggtitle("Total Nonfarm Control State Map For LA")