Read district heating prices

.buildlibrary

@@ -1,4 +1,4 @@
-ValidationKey: '585930'
+ValidationKey: '586110'
 AutocreateReadme: yes
 AcceptedWarnings:
 - 'Warning: package ''.*'' was built under R version'

CITATION.cff

@@ -3,7 +3,7 @@ message: If you use this software, please cite it using the metadata from this f
 type: software
 title: 'mredgebuildings: Prepare data to be used by the EDGE-Buildings model'
 version: 0.3.0
-date-released: '2023-06-23'
+date-released: '2023-06-29'
 abstract: Prepare data to be used by the EDGE-Buildings model.
 authors:
 - family-names: Hasse

DESCRIPTION

@@ -2,7 +2,7 @@ Type: Package
 Package: mredgebuildings
 Title: Prepare data to be used by the EDGE-Buildings model
 Version: 0.3.0
-Date: 2023-06-23
+Date: 2023-06-29
 Authors@R: c(
     person("Robin", "Hasse", , "[email protected]", role = c("aut", "cre"),
            comment = c(ORCID = "0000-0003-1818-3186")),

NAMESPACE

@@ -15,6 +15,7 @@ export(convertCensusHub)
 export(convertDaioglou)
 export(convertEEAfloorspace)
 export(convertEUBuildingsDB)
+export(convertEnergiforsk2016)
 export(convertEurObservER)
 export(convertEuropeanCommissionRenovation)
 export(convertEurostatBuildings)
@@ -29,6 +30,7 @@ export(readDaioglou)
 export(readEEAfloorspace)
 export(readEHI)
 export(readEUBuildingsDB)
+export(readEnergiforsk2016)
 export(readEurObservER)
 export(readEuropeanCommissionRenovation)
 export(readEurostatBuildings)
@@ -113,6 +115,7 @@ importFrom(stats,predict)
 importFrom(stats,pweibull)
 importFrom(tidyr,complete)
 importFrom(tidyr,gather)
+importFrom(tidyr,pivot_longer)
 importFrom(tidyr,replace_na)
 importFrom(tidyr,separate)
 importFrom(tidyr,spread)

R/calcBuildingStock.R

@@ -726,7 +726,8 @@ calcBuildingStock <- function(subtype = c("residential", "commercial")) {
         group_by(across(all_of(c("region", "period")))) %>%
         mutate(value = proportions(.data[["value"]]))
       stock <- stock %>%
-        left_join(heatingShares, by = c("region", "period")) %>%
+        left_join(heatingShares, by = c("region", "period"),
+                  relationship = "many-to-many") %>%
         mutate(value = .data[["value.x"]] * .data[["value.y"]]) %>%
         select(-"value.x", -"value.y") %>%
         rename(heating = "carrier")

R/convertEnergiforsk2016.R

@@ -0,0 +1,26 @@
+#' Convert European district heating prices
+#'
+#' We should consider prper extrapolation
+#'
+#' @author Robin Hasse
+#'
+#'@param x MAgPIE object with original district heating prices
+#'
+#' @importFrom madrat toolCountry2isocode toolCountryFill
+#' @importFrom magclass getItems<- getItems collapseDim
+#' @export
+
+convertEnergiforsk2016 <- function(x) {
+
+  data <- x
+
+  # rename regions: ISO2 -> ISO3
+  getItems(data, 1) <- toolCountry2isocode(getItems(data, 1))
+
+  # fill missing countries with NA
+  data <- data %>%
+    collapseDim() %>%
+    toolCountryFill(verbosity = 2)
+
+  return(data)
+}

R/readEHI.R

@@ -1,14 +1,5 @@
 #' Read EHI Heating Market Report
 #'
-#' PDF market reports on Renewables in Europe
-#'
-#' @source https://www.eurobserv-er.org/
-#' @note I download the report and snip a png of the table. Then I extract the
-#' data with [this online tool](https://extracttable.com/) and save it as a
-#' tabular csv in the source folder. The Tool is not fully reliable. This
-#' function throws an error if values don't sum up. You then have to correct the
-#' csv file manually.
-#'
 #' @param subtype <report series>.<variable>
 #' @returns magpie object
 #'

R/readEnergiforsk2016.R

@@ -0,0 +1,21 @@
+#' Read European district heating prices
+#'
+#' @author Robin Hasse
+#'
+#' @importFrom utils read.csv
+#' @importFrom tidyr pivot_longer
+#' @importFrom dplyr select matches mutate rename
+#' @importFrom quitte as.quitte
+#' @importFrom magclass as.magpie
+#' @export
+
+readEnergiforsk2016 <- function() {
+  file <- "european-district-heating-price-series-energiforskrapport-2016-316.csv"
+  read.csv(file, skip = 1) %>%
+    select(-matches("All\\.countries")) %>%
+    pivot_longer(-"Year", names_to = "region", values_to = "value") %>%
+    mutate(region = gsub("\\.", " ", .data[["region"]])) %>%
+    rename(period = "Year") %>%
+    as.quitte() %>%
+    as.magpie()
+}