if (FALSE) {
# Fit fect model with bootstrap
fit <- fect(Y ~ D + X, data = panel_data, keep.sims = TRUE)
# Calculate cumulative effects for all treated units
results <- effect(fit)
# Calculate period-specific effects for specific units
results_specific <- effect(fit,
cumu = FALSE,
id = c("unit1", "unit2"),
period = c(1, 4))
# View results
print(results$est.catt)
}
Run the code above in your browser using DataLab