R/gibbs_draws.R
compute_si_cutoff.RdAcross a matrix of discretised probability distributions
(see estimate_advantage
this function returns the largest index
(across all columns) such that the
cumulative probability mass before index is
1 - miss_at_most.
compute_si_cutoff(si_distr, miss_at_most = 0.05)a matrix with two columns, each containing the probability mass function for the discrete serial interval for each of the two pathogen/strain/variants, starting with the probability mass function for day 0 in the first row, which should be 0. each column in the matrix should sum to 1
numeric. probability mass in the tail of the SI distribution
integer