## Not run:
#
# ##= load toxoplasmosis example ============================
# setDALYexample(2)
#
# ##= perform DALY calculation for different scenarios ======
# toxo_00 <- getDALY(aw = FALSE, dr = 0)
# toxo_03 <- getDALY(aw = FALSE, dr = 0.03)
# toxo_13 <- getDALY(aw = TRUE, dr = 0.03)
#
# ##= combine scenarios in 'DALY_list' ======================
# toxo <- DALY_list(toxo_00, toxo_03, toxo_13)
#
# ##= plot YLL and YLL for toxo_00 ==========================
# par(mar = c(4, 4, 1, 1) + .5)
# par(mfrow = c(1, 2))
# scatterplot(toxo_00, plot = "YLL",
# outcomes = FALSE, legend = FALSE)
# scatterplot(toxo_00, plot = "YLD",
# outcomes = FALSE, legend = FALSE)
#
# ##= plot all three scenarios ==============================
# par(mfrow = c(1, 1))
# scatterplot(toxo,
# legend = c("DALY[0,0]",
# "DALY[0,0.03]",
# "DALY[1,0.03]"),
# legend_pos = "topleft",
# log = "xy",
# main = "Scenario analysis")
# ## End(Not run)
Run the code above in your browser using DataLab