Difference between revisions of "ERFs of environmental pollutants"

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_en5827", "initiate") # Get the latest ovariables from code calculations on page Op_en5918 (ERF of several environmental pollutions).
ERF

ERF 		<- EvalOutput(ERF, N = 3)
ERF@output$ERFResult <- as.numeric(as.character(ERF@output$ERFResult))
ERF@output <- ERF@output[colnames(ERF@output) != "Obs"]

print(xtable(ERF@output, digits = 4), type = 'html')

ggplot(ERF@output[ERF@output$ERF.Parameter == "RR" , ], # Show a bar chart
	aes(x = Pollutant, y = ERFResult, fill = Response.metric)) +
	geom_bar(position="dodge") + 
	theme_grey(base_size = 24) +
	theme(axis.text.x = element_text(angle = 90, hjust = 1))

ggplot(ERF@output[ERF@output$ERF.Parameter == "UR" , ], # Show a bar chart
	aes(x = Health.impact, y = ERFResult, fill = Response.metric)) +
	geom_bar(position="dodge") + 
	theme_grey(base_size = 24) +
	theme(axis.text.x = element_text(angle = 90, hjust = 1))

library(OpasnetUtils)

d <- opbase.data("Op_en5827")
d$Obs <- NULL
colnames(d) <- gsub(" ", "_", colnames(d))
d$Result <- ifelse(d$Result == "", "0", as.character(d$Result))

ERF <- Ovariable("ERF", data = d[d$Observation == "ERF", colnames(d) != "Observation"])

threshold <- Ovariable("Threshold", data = d[d$Observation == "Threshold", colnames(d) != "Observation"])

objects.store(ERF, threshold)
cat("Ovariables ERF, threshold stored. Page: Op_en5827, code: initiate.\n")