I would like to save the results of `did_imputation` to be used later when generating coefficient plots (by running the command `event_plot`). To be more specific, I would like to run the coefficient plots once I decide what scale I would like to fix the Y-axis at. Because I will only decide this after I have estimated results for various outcomes, I would like to run `event_plot` later on instead of right after each one of the estimations.

I am not sure how I can store the output of `did_imputation` to then call `event_plot`. I tried using `estimates save` but it does not seem to work with `did_imputation` (nothing is returned when calling the stored results with `estimates use`).

Code:
. did_imputation Y i t Ei, allhorizons minn(15) pretrend(5)

                                                         Number of obs = 3,000
------------------------------------------------------------------------------
           Y | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
        tau0 |    4.87038   .1956373    24.89   0.000     4.486938    5.253822
        tau1 |   3.827304   .2102286    18.21   0.000     3.415263    4.239344
        tau2 |    3.26428   .2490712    13.11   0.000      2.77611    3.752451
        tau3 |   2.074517   .3117521     6.65   0.000     1.463494     2.68554
        tau4 |   .7844959   .3371969     2.33   0.020      .123602     1.44539
        tau5 |   1.030274   .5526556     1.86   0.062    -.0529109    2.113459
        pre1 |   -.131531   .3122927    -0.42   0.674    -.7436135    .4805514
        pre2 |  -.2804478        .27    -1.04   0.299     -.809638    .2487424
        pre3 |  -.1544596   .2480803    -0.62   0.534    -.6406881     .331769
        pre4 |  -.4179183   .2060002    -2.03   0.042    -.8216713   -.0141654
        pre5 |   .2538401   .2042816     1.24   0.214    -.1465444    .6542246
------------------------------------------------------------------------------

. return list

matrices:
              r(table) :  9 x 11

. event_plot

. estimates save "saving_output", replace
file saving_output.ster saved

. clear

. estimates use "saving_output"

. return list

.