Create a summary

[RodrigoZepeda]

For the rstan object

[RodrigoZepeda]

So the summary function would be something like this:

library(diseasenowcasting)

data(denguedat)

# Change to more iterations
predictions <- nowcast(denguedat, "onset_week", "report_week",
                       now = as.Date("1990-10-01"),
                       dist = "NegativeBinomial", cores = 4, iter = 1000)

summary <- function(predictions, quantiles = c(0.025, 0.975)){

  #1) Summarize the predictions (rstan object)
  data_predict <- function(predictions, quantiles)

  #2) Join the predicted to the real data
  data_predict |> left_join(real_data)

}

[RodrigoZepeda]

We are first doing #19 before this one here

[MatteoPS]

done with summary_nowcast() #45