`add_trace( y = colQuantiles(extract$symptoms, probs = 0.975), fill = "tonexty", fillcolor = "rgba(0,0,0.5,0.3)", showlegend = TRUE, name = "CI 5%-95%") %>%` in symptom_report.R. the 2.5% - 97.5% quantiles are extracted, but are labelled 5% - 95%