Skip to content
Snippets Groups Projects
makeFile.R 2.06 KiB
Newer Older
# loads the data and runs the Rmd render

# Packages ----
library(data.table)
library(here)

# Functions ----
Ben Anderson's avatar
Ben Anderson committed
makeReport <- function(f){
  # default = html
Ben Anderson's avatar
Ben Anderson committed
  rmarkdown::render(input = paste0(here::here("itsTheCatsStupid", f), ".Rmd"),
Ben Anderson's avatar
Ben Anderson committed
                    params = list(title = title,
                                  subtitle = subtitle,
                                  authors = authors),
                    output_file = paste0(here::here("docs/"), f, ".html")
  )
}

dp <- "~/Dropbox/data/"

#> define yaml ----
Ben Anderson's avatar
Ben Anderson committed
rmdFile <- "itsTheCatsStupid" # not the full path
Ben Anderson's avatar
Ben Anderson committed
title = "#backOfaFagPacket: Its the Cats, stupid"
Ben Anderson's avatar
Ben Anderson committed
subtitle = "Does cat ownership correlate with home energy demand?"
authors = "Ben Anderson"

#> load the postcode data here (slow)
Ben Anderson's avatar
Ben Anderson committed

postcodes_elec_dt <- data.table::fread(paste0(dp, "beis/subnationalElec/Postcode_level_all_meters_electricity_2015.csv"))
postcodes_elec_dt[, pcd_sector := data.table::tstrsplit(POSTCODE, " ", keep = c(1))]
pc_sector_elec_dt <- postcodes_elec_dt[, .(nPostcodes = .N, 
                                           total_elec_kWh = sum(`Consumption (kWh)`),
                                           nElecMeters = sum(`Number of meters`)
                                           ), keyby = .(pcd_sector)]
nrow(pc_sector_elec_dt)

postcodes_gas_dt <- data.table::fread(paste0(dp, "beis/subnationalGas/Experimental_Gas_Postcode_Statistics_2015.csv"))
postcodes_gas_dt[, pcd_sector := data.table::tstrsplit(POSTCODE, " ", keep = c(1))]
pc_sector_gas_dt <- postcodes_gas_dt[, .(total_gas_kWh = sum(`Consumption (kWh)`),
                                           nGasMeters = sum(`Number of meters`)), keyby = .(pcd_sector)]
nrow(pc_sector_gas_dt)

setkey(pc_sector_elec_dt, pcd_sector)
setkey(pc_sector_gas_dt, pcd_sector)

pc_sector_energy_dt <- pc_sector_gas_dt[pc_sector_elec_dt]

pc_sector_region_dt <- data.table::fread(here::here("data", "postcode_sectors_dt.csv"))
setkey(pc_sector_region_dt, pcd_sector)

pc_sector_energy_dt <- pc_sector_region_dt[pc_sector_energy_dt]

#> re-run report here ----
Ben Anderson's avatar
Ben Anderson committed
makeReport(rmdFile)