Difference between revisions of "Disease risk"

library(OpasnetUtils) # A package with Opasnet functionalities
library(xtable)       # A package with table formats
library(ggplot2)      # A package with fancy graph formats

objects.latest("Op_en5917", "initiate") # Get the latest ovariables from code calculations on page Op_en5917.


diseaseRisk <- EvalOutput(diseaseRisk)

oprint(summary(diseaseRisk))

ggplot(diseaseRisk@output, # Show a bar chart
	aes(x = Disease, y = diseaseRiskResult, fill = Response.metric)) +
	geom_bar(position="dodge") + 
	theme_grey(base_size = 24) +
	theme(axis.text.x = element_text(angle = 90, hjust = 1)) +
	scale_y_log10() # limits = c(1, 1000000))

library(OpasnetUtils)

d <- opbase.data("Op_en5917")
d$Obs <- NULL
colnames(d) <- gsub(" ", "_", colnames(d))
d$Result[d$Unit == "cases /100000py"] <- d$Result[d$Unit == "cases /100000py"] / 100000 # Units to cases/py
levels(d$Unit)[levels(d$Unit) == "cases /100000py"] <- "cases /py"

disincidence <- Ovariable("disincidence", data = d, save = TRUE)

cat("Ovariable disincidence saved. page_indent = Op_en5917, code_name = initiate\n")