Updated code to run pfs, compute quantiles, and create plots for SIR model