library("unitdid")
library(unitdid)
mdl_base <- base_heterocp |>
unitdid(yname = "y",
iname = "id",
tname = "year",
ename = "cyear",
bname = "byear")
summary(mdl_base, agg = "event_age")# A tibble: 159 × 4
rel_time mean n event_age
<int> <dbl> <dbl> <int>
1 0 -0.0138 16 15
2 0 -0.0937 52 16
3 0 -0.155 82 17
4 0 -0.132 103 18
5 0 -0.0927 127 19
6 0 -0.0674 157 20
7 0 -0.104 189 21
8 0 -0.0709 222 22
9 0 -0.0832 233 23
10 0 -0.0557 279 24
# ℹ 149 more rows