covid19hr_reporting_delay_simulation

library(tidyverse)
library(jsonlite)
set.seed(666)
data <-
  "https://www.koronavirus.hr/json/?action=po_osobama" %>%
  fromJSON() %>% 
  as_tibble() %>%
  transmute(report = as.Date(Datum)) %>%
  mutate(
    test_delay = rgamma(n(), 2, 2),
    exam_delay = rgamma(n(), 3, 2),
    onset_delay = rgamma(n(), 3, .4),
    infection_delay = rlnorm(n(), log(6), log(1.5)),
    total_delay = infection_delay + onset_delay + exam_delay + test_delay,
    infection = report - floor(total_delay)
    ) %>%
    select(report, infection, total_delay, everything())