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

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

ERF <- Ovariable("ERF", ddata = "Op_en5827", getddata = FALSE, save = TRUE)
ERF@data <- ERF@data[colnames(ERF@data) != "Obs"]

cat("Ovariable ERF saved. Page: Op_en5827, code: calculations.\n")

ERF <- EvalOutput(ERF, N = 3)

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