data-raw/test.R

In common2016/didWD: DID with staggred treatment

rm(list = ls())
load('/Users/yangnay/elements/RawE/17_HuaDong/research/AI/RCodeEctrics/data-raw/regdata.rdata') # From TFP.R
library(pacman)
p_load(didWD,lmtest,plm, reshape2, devtools, ggplot2, tidyverse, matlab, did)
load_all()

regdata$id <- paste(regdata$indus, regdata$iso2c, sep = '_')
regdata$w <- 0
regdata$w[regdata$operationalstock > 100] <- 1

cn <- regdata[regdata$iso2c %in% c('CN') & regdata$year <= 2014,]

cn$wcont <- log(cn$operationalstock+1)
fit <- didWD(cn, 'id', 'year', y = 'TFP', w = 'w', wcontinuous = 'wcont', pretime = T)
fit <- didWD(cn, 'id', 'year', y = 'TFP', w = 'w', pretime = T)
summary(fit$fit)
ans <- dyneff(fit, se = 'arellano', type = 'sss')

xname <- coef(fit$fit) %>% names()
aggeff(fit$fit, xname)

aggbeta <- xname[str_detect(xname,'d2008')]
aggeff(fit$fit, aggbeta)

ans <- fit$data