Difference between revisions of "Energy balance in Stuttgart"

library(OpasnetUtils)
library(ggplot2)

N <- 10

objects.latest("Op_en5141", code_name = "initiate")

balance <- Ovariable("balance", ddata = "Op_en5815.equations")

balance@data$Equation <- as.character(levels(balance@data$Equation)[balance@data$Equation]) # Change factor into character.
balance@data$Policy[balance@data$Policy == ""] <- NA # Prepare indices for fillna. This should be part of fillna.
balance@data$CHPcapacity[balance@data$CHPcapacity == ""] <- NA
balance@data <- fillna(balance@data, marginals = 1:2) # Fill empty slots in indices.

nonlinearity <- Ovariable("nonlinearity", ddata = "Op_en5815", subset = "Nonlinearity parameters")

energy.balance <- EvalOutput(energy.balance) 

energy.balance@marginal[colnames(energy.balance@output) == "energy.balanceVar"] <- TRUE # This has to be done manually because CheckMarginals does not notice this as a marginal.
energy.balance@marginal[colnames(energy.balance@output) %in% nonlinearity@output$critIndex] <- FALSE # Nonlinear indices are demarginalised because only one of the two equations apply.

oprint(summary(energy.balance))

ggplot(energy.balance@output, aes_string(x = "energy.balanceVar", y = "energy.balanceResult", fill = "Policy")) + 
	geom_boxplot() +
	theme_grey(base_size = 24) + 
	theme(axis.text.x = element_text(angle = 90, hjust = 1))

# objects.store(energy.balance)