Output

This command is implemented as an r-class command, which leaves behind the
following results saved in r():

r(maxinfect) - maximum size of the infected subpopulation throughout
the simulation timeframe. This will be measured in percent if the
option percent was specified, otherwise it is an absolute number
of infected.

r(d_maxinfect) - this is the step of simulation on which the maximum
infected size was reached, and can take values 0, 1, 2, .... If number of steps of simulation per day is 1, then this is also the day-number when the maximum infected size was reached.

r(t_maxinfext) - this is the same as r(d_maxinfect) if no option
day0() was specified and simulation steps per day is 1. But if the
day0() option was specified, then it takes the calendar date
(in Stata's date) format for the date where the maximum infected size
was reached. For example, a value 22025 will correspond to
April 20, 2020. If the number of steps of simulation in more than 1 per day, expect this value to be non-integer.

r(o_maxinfext) - this is the observation number, in which maximum
infected size was reached.

The simulation of the model steps is also saved into a matrix with the name
of the model: r(sir) or r(seir). The columns of this
matrix area the components of the model (time, susceptible, etc), and each
row corresponds to one day of the simulation. The number of rows will be 1
plus the value of the option days(). The time variable in the matrix
is always the day of the model simulation, not the date, even when the
option day0() was specified.

The command also produces output in the data area in memory (which can be
cleared with the option clear if it already contains anything). The
generated variables will have names corresponding to the model components
(for example S for Susceptible, I for Infected, etc).
Model components will be named in upper case if they represent absolute
numbers, and lower case if they represent shares (if option percent
was specified). The time variable t is always lower case.

The model simulation command will build a graph plotting the
trajectories of all the model componets over the simulation time frame.
This graph can be suppressed by specifying an option nograph or it
can be adjusted by specifying common twoway graphing options.

Finally, a report can be created as a file if the pdfreport() is specified.
The report file is always replaced.