# Run multiple simulations
multi_result <- ssir(N = 10000, T = 300, alpha = 0.3, inf_period = 4,
inf_init = 32, report = 0.02, lag = 7, rep = 5)
# Plot the results
plots <- plot(multi_result)
# Display the first simulation's plots
plots[[1]]
Run the code above in your browser using DataLab