Under circumstances (e.g. low counts, high variability in R) it is advisable not to show the quantiles vs. time of different simulation runs (replicates) but, if still readable, the individual time series. See https://arxiv.org/abs/2007.05035, https://twitter.com/CT_Bergstrom/status/1282583203316817920
(See also further visualization suggestions, e.g. probability to reach a given count vs. time.)
One solution would be to add a parameter in the plot function, e.g. type='quantiles' or type='timeseries', the latter being the default.
Under circumstances (e.g. low counts, high variability in R) it is advisable not to show the quantiles vs. time of different simulation runs (replicates) but, if still readable, the individual time series. See https://arxiv.org/abs/2007.05035, https://twitter.com/CT_Bergstrom/status/1282583203316817920
(See also further visualization suggestions, e.g. probability to reach a given count vs. time.)
One solution would be to add a parameter in the
plotfunction, e.g.type='quantiles'ortype='timeseries', the latter being the default.