# Load data of viral concentration in wastewater
data("ww.data")
# Estimate Rt based on wastewater data
x = estimate_R_ww(
ww.conc = ww.data,
dist.fec = ern::def_dist(
dist = "gamma",
mean = 12.9,
mean_sd = 1.13,
shape = 1.75,
shape_sd = 0.26,
max = 33
),
dist.gi = ern::def_dist(
dist = "gamma",
mean = 6.84,
mean_sd = 0.74,
shape = 2.39,
shape_sd = 0.35,
max = 15
),
silent = TRUE
)
# Diagnostic plot
g = plot_diagnostic_ww(x)
plot(g)
g2 = plot_diagnostic_ww(x, wrap.plots = FALSE, caption = "This is your caption")
plot(g2$wastewater_data)
plot(g2$inferred_incidence)
plot(g2$Rt)
Run the code above in your browser using DataLab