Combine Rt estimates from EpiEstim across multiple runs.
combine_rt_epiestim(rt, samples, q = NULL, rank = TRUE)A list of Rt calculations from
lancelot_rt_trajectories_epiestim() (though any Rt calculation that
confirms to this will work)
A list of samples from mcstate::pmcmc
A vector of quantiles to return values for
A boolean deciding whether to rank trajectories by increasing incidence or not before combining Rt estimates
A list of Rt output in the same structure as the first
element of rt. Rt estimates will be aggregated across
regions (or whatever else you are aggregating on) based on the
parameters in samples.