dem.sub <- dem[dem[, "wbcode2"] <= 100, ]
dem.sub.panel <- PanelData(dem.sub, "wbcode2", "year", "dem", "y")
# create subset of data for simplicity
PM.results <- PanelMatch(panel.data = dem.sub.panel, lag = 4,
refinement.method = "ps.match",
match.missing = TRUE,
covs.formula = ~ tradewb,
size.match = 5, qoi = "att",
lead = 0:4,
forbid.treatment.reversal = FALSE)
PE.results <- PanelEstimate(sets = PM.results,
panel.data = dem.sub.panel,
se.method = "unconditional")
summary(PE.results)
summary(PE.results, confidence.level = .9)
Run the code above in your browser using DataLab