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.stock', 'energy.per.area', 'heating.fraction'),
		Ident = c('Op_en5932/initiate')
	),
	formula = function(...){

		building.stock@output <- building.stock@output[building.stock@output$Observation == "Area" , ]
		colnames(building.stock@output)[colnames(building.stock@output) == "Observation"] <- "building.stockObservation"

		energy.per.area@output <- energy.per.area@output[energy.per.area@output$Observation == observation , ]
		colnames(energy.per.area@output)[colnames(energy.per.area@output) == paste("energy.per.area", observation, sep = "")] <- "energy.per.areaResult"

		out <- building.stock * heating.fraction / 100 * energy.per.area

		return(out)
	} 
)

test <- EvalOutput(test)

oprint(summary(test))

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

ograph(test, x = "Heating", fill = "Type")

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

library(OpasnetUtils)

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

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

heating.fraction <- Ovariable(
	name = 'heating.fraction',
	ddata = 'Op_en5932',
	subset = 'Fractions of detached houses according heating type',
	getddata = FALSE
)

energy.per.volume <- Ovariable(
	name = 'energy.per.volume',
	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.stock, 
	energy.per.area,
	energy.per.volume,
	heating.fraction,
	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))