Difference between revisions of "Building stock in Kuopio"

library(OpasnetUtils)
library(xtable)       # A package with table formats
library(ggplot2)      # A package with fancy graph formats

objects.latest("Op_en5932", "initiate")

test <- Ovariable(
	name = 'test',
	dependencies = data.frame(
		Name = c('building.registry.data', 'baseline.energy.consumption.per.area.unit'),
		Ident = c('Op_en5932/initiate', 'Op_en5932/initiate')
	),
	formula = function(dependencies, ...){
		ComputeDependencies(dependencies, ...)

		colnames(building.registry.data@output)[colnames(building.registry.data@output) == "Observation"] <- "building.registry.dataObservation"

		out <- building.registry.data * baseline.energy.consumption.per.area.unit

		out@output <- out@output[
			!is.na(out@output$Result) & 
			!out@output$building.registry.dataObservation %in% c("Number", "Volume") &
			out@output$Observation == observation , 
		]

		return(out)
	} 
)

test <- EvalOutput(test)

oprint(test@output[test@output$Iter==1, ])

ograph(test, x = "Constructed", fill = "Type", type = geom_bar())

#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))

library(OpasnetUtils)

building.registry.data <- Ovariable(
	name = 'building.registry.data',
	ddata = 'Op_en5932',
	subset = 'Building registry data',
	getddata = FALSE
)

baseline.energy.consumption.per.area.unit <- Ovariable(
	name = 'baseline.energy.consumption.per.area.unit',
	ddata = 'Op_en5932',
	subset = 'Baseline energy consumption per area unit',
	getddata = FALSE
)


baseline.energy.consumption.per.volume.unit <- Ovariable(
	name = 'baseline.energy.consumption.per.area.unit',
	ddata = 'Op_en5932',
	subset = 'Baseline energy consumption per volume unit',
	getddata = FALSE
)


fraction.of.houses.renovated.per.year <- Ovariable(
	name = 'fraction.of.houses.renovated.per.year',
	ddata = 'Op_en5932',
	subset = 'Fraction of houses renovated per year',
	getddata = FALSE
)


floor.area.of.new.houses.and.additional.construction.per.year <- Ovariable(
	name = 'floor.area.of.new.houses.and.additional.construction.per.year',
	ddata = 'Op_en5932',
	subset = 'Floor area of new houses and additional construction per year',
	getddata = FALSE
)


Maximum.allowed.energy.consumption.per.unit <- Ovariable(
	name = 'Maximum.allowed.energy.consumption.per.unit',
	ddata = 'Op_en5932',
	subset = 'Maximum allowed energy consumption per unit (= E-value)',
	getddata = FALSE
)

energy.saving.potential.of.different.renovations <- Ovariable(
	name = 'energy.saving.potential.of.different.renovations',
	ddata = 'Op_en5932',
	subset = 'Energy saving potential of different renovations',
	getddata = FALSE
)

energy.saving.potential.of.different.renovations <- Ovariable(
	name = 'energy.saving.potential.of.different.renovations',
	ddata = 'Op_en5932',
	subset = 'Energy saving potential of different renovations',
	getddata = FALSE
)


energy.classes.of.buildings <- Ovariable(
	name = 'energy.classes.of.buildings',
	ddata = 'Op_en5932',
	subset = 'Energy classes of buildings',
	getddata = FALSE
)
objects.store(
	building.registry.data, 
	baseline.energy.consumption.per.area.unit,
	baseline.energy.consumption.per.volume.unit,
	fraction.of.houses.renovated.per.year,
	floor.area.of.new.houses.and.additional.construction.per.year,
	energy.saving.potential.of.different.renovations,
	energy.classes.of.buildings
)

cat("Objects initiated!\n")

library(OpasnetUtils)
library(ggplot2)

if(secret) {
	objects.get("ozJoC2EaIChFI0uC") 
	building.stock <- objects.decode(etable, key)
} else {
	building.stock <- Ovariable("building.stock", ddata = "Op_en5932", getddata = FALSE)@output
}

colnames(building.stock) <- as.vector(t(building.stock[1, ]))
building.stock <- building.stock[2:nrow(building.stock) , ]

for(i in c(1,3,4,5,7,8)) {
	building.stock[[i]] <- as.numeric(gsub(",", "", building.stock[[i]]))
}

building.stock <- dropall(building.stock)

cat("Mean floor area of buildings in Kuopio by heating type.\n")

oprint(as.data.frame(as.table(tapply(
	building.stock$kokonaisala, 
	building.stock[c("paaasiallinen_lammitystapa")], 
	mean, 
	na.rm = TRUE
))))

ggplot(building.stock, aes(x = paaasiallinen_lammitystapa, weight = kokonaisala, fill = kayttotarkoitus)) + 
	geom_bar() +
	theme_grey(base_size = 24) + 
	theme(axis.text.x = element_text(angle = 90, hjust = 1))