Difference between revisions of "Building model"

### This code is Op_en6289/buildingstest on page [[Building model]].

library(OpasnetUtils)

########## Calculate all building events (constructions, demolitions, renovations)

buildings <- Ovariable("buildings", 
	dependencies = data.frame(Name = c(
		"stockBuildings",
		"changeBuildings",
		"heatingShares", # This can be indexed by building year (Time) or Observation year (Obsyear)
		"efficiencyShares", # Same here
		"renovationShares", 
		"renovationRate", # Fraction of buildings renovated between timepoints
		"obstime" # In the new version obstime is ovariable with column Obsyear.
	)),
	formula = function(...) {

		stock <- stockBuildings * obstime * heatingShares * efficiencyShares

		stock2 <- changeBuildings * obstime # Must contain data about heatings and efficiencies
		stock2@output$Time <- as.factor(stock2@output$Time) # Also stock@output$Time is factor
		buil <- combine(stock, stock2) # Buildings in the whole timeline
		buil <- unkeep(buil, sources = TRUE)
		buil@output$Age <- as.numeric(as.character(buil@output$Obsyear)) - as.numeric(as.character(buil@output$Time))
		buil@marginal <- c(buil@marginal, TRUE)
		buil@output <- buil@output[buil@output$Age >= 0 , ]

		# Note that if stockBuildings and changeBuildings have different marginals (typically policies),
		# there will be NA in these indices. This will be corrected below.
		for(i in colnames(buil@output)[buil@marginal]) {
			if(any(is.na(buil@output[[i]]))) {
				buil@output <- fillna(buil@output, i)
				print(paste("Column", i, "treated with fillna."))
			}
		}
		renovate <- (1 - renovationRate) * (1 - renovationShares) # Assumes that data has row for Renovation: None = 0
		renovate@output <- renovate@output[renovate@output$Renovation == "None" , ]
		renovate <- combine(renovate, renovationRate * renovationShares)
		
		out <- data.frame()
		prevreno <- unkeep(buil * renovate * 0, sources = TRUE) # Previously renovated buildings for the first time point
		prevreno@output <- prevreno@output[prevreno@output$Renovation != "None" , ][1,] # Just take one renovated example row
		marginals <- colnames(prevreno@output)[prevreno@marginal]
		
		for(i in obstime@output$Obsyear) { # Accumulate the building stock
			# Take the building stock
			temp <- buil
			temp@output <- temp@output[temp@output$Obsyear == i , ] # Take buildings of year i
			prevreno@output$Obsyear <- i # update the observation year of the previous renovation
			# Subtrack from the building stock buildings that are already renovated.
			temp <- combine(temp, -1 * oapply(prevreno, cols = "Renovation", FUN = sum)) # Subtract previously renovated
			reno <- unkeep(temp * renovate, sources = TRUE) # Renovate the current non-renovated stock
			# If renovate does not match with this timepoint, mark all buildings here as unrenovated.
			# Is there new renovation in this time point?
			if(sum(result(reno)) == 0) {
				reno <- temp * Ovariable(
					output = data.frame(Renovation = "None", Result = 1), 
					marginal = c(TRUE, FALSE)
				)
			}
			newreno <- combine(prevreno, reno) # Renovated buildings in this time point.
			# Take the previous timepoints, and previously renovated and now possibly renovated together.
			out <- rbind(newreno@output, out) # Add previously and now renovated together
			prevreno <- unkeep(newreno, sources = TRUE)
			prevreno@output <- prevreno@output[prevreno@output$Renovation != "None" , ]
		}
		out <- Ovariable(output = out, marginal = colnames(out) %in% marginals)
		out <- oapply(unkeep(out, sources = TRUE), cols = c("Age", "Time"), FUN = sum)
		colnames(out@output)[colnames(out@output) == "Obsyear"] <- "Time" # From now on, Time means the time of observation
		out@output$Time <- as.numeric(as.character(out@output$Time))
		
		return(out)
	}
)

objects.store(buildings)

cat("Saved ovariable buildings\n")