[hr]


tab time, gen(time_dummy)

       time |      Freq.     Percent        Cum.
------------+-----------------------------------
       2020 |         46       20.00       20.00
       2021 |         46       20.00       40.00
       2022 |         46       20.00       60.00
       2023 |         46       20.00       80.00
       2024 |         46       20.00      100.00
------------+-----------------------------------
      Total |        230      100.00

.  gen base =(time == 2022)

. foreach var of varlist time_dummy*{
  2. gen treated_`var' = treated*`var'
  3.
.
.
.
.
. }

. reg outcome treated_* if base==1
note: treated_time_dummy1 omitted because of collinearity.
note: treated_time_dummy2 omitted because of collinearity.
note: treated_time_dummy4 omitted because of collinearity.
note: treated_time_dummy5 omitted because of collinearity.

      Source |       SS           df       MS      Number of obs   =        4
> 6
-------------+----------------------------------   F(1, 44)        =     14.2
> 9
       Model |  10398.0728         1  10398.0728   Prob > F        =    0.000
> 5
    Residual |  32019.5735        44  727.717579   R-squared       =    0.245
> 1
-------------+----------------------------------   Adj R-squared   =    0.228
> 0
       Total |  42417.6463        45  942.614362   Root MSE        =    26.97
> 6

-----------------------------------------------------------------------------
> -
     outcome | Coefficient  Std. err.      t    P>|t|     [95% conf. interval
> ]
-------------+---------------------------------------------------------------
> -
treated_ti~1 |          0  (omitted)
treated_ti~2 |          0  (omitted)
treated_ti~3 |   31.14812   8.240184     3.78   0.000     14.54112    47.7551
> 2
treated_ti~4 |          0  (omitted)
treated_ti~5 |          0  (omitted)
       _cons |   6.671198   6.542699     1.02   0.313    -6.514747    19.8571
> 4
-----------------------------------------------------------------------------
> -