calculate solution to ODEs for each pathway, and extract median, 95\
calc_sol_chain(filename, model, solving_time_max = 72)
\itemfilename.RData file produced by lazymcmc
\itemmodelstring: "full", "endosomal_only", "tmprss2_only"
\itemsolving_time_maxtime to solve until. Default = 72 hpi.
tibble with columns model, t, V, T1... and probs, which can take the values 0.025, 0.5, 0.975 or "max_LL"
calculate solution to ODEs for each pathway, and extract median, 95\