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)
library(xtable)

diseaseRisk <- Ovariable("diseaseRisk", ddata = "Op_en5917", getddata = FALSE, save = TRUE)

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

print(xtable(EvalOutput(diseaseRisk, N = 1)@output), type = 'html')