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", "calculations") # Get the latest ovariables from code calculations on page Op_en5917.

objects.get("3i7fJFPcIoH6CKiN")

diseaseRisk <- EvalOutput(diseaseRisk)

print(xtable(diseaseRisk@output), type = 'html') # Show a result table

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 = TRUE, save = FALSE)
colnames(diseaseRisk@data)
diseaseRisk@data <- diseaseRisk@data[ , colnames(diseaseRisk@data) != "Obs"]

objects.put(diseaseRisk)

cat("Ovariable diseaseRisk saved. Save the key for surther use.\n") # page_indent = Op_en5917, code_name = calculations\n")

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