old version

analysis/GREENGridModel/GREENGridModel_master_v1.0.Rmd

@@ -46,23 +46,23 @@ spatialec::loadLibraries(myPackages) # this will try to install packages it can'
 projRoot <- here::here()
 
 # --- Project Settings ----
-source(paste0(projRoot, "/_scripts/spatialecSettings.R"))   # loads knitr options & assigns file paths and loads system info
+source(paste0(projRoot, "/packageSettings.R"))   # loads knitr options & assigns file paths and loads system info
 
 # ---- Local Settings ----
-sParams$censusDataSource <- "http://archive.stats.govt.nz/Census/2013-census/data-tables/meshblock-dataset.aspx"
+p_Params$censusDataSource <- "http://archive.stats.govt.nz/Census/2013-census/data-tables/meshblock-dataset.aspx"
 
-sParams$ggPath <- "~/Data/NZ_GREENGrid" # over-rides default for speed
+p_Params$ggPath <- "~/Data/NZ_GREENGrid" # over-rides default for speed
 
 # Print system information
-outputMessage("Running on", sParams$sysName, 
-              "from projLoc =", projRoot)   # sysName is set in spatialecSettings.R
+message("Running on ", p_Params$sysName, 
+              " from projLoc =", projRoot)   # sysName is set in spatialecSettings.R
 ```
 
 # Report Purpose
 
 To use:
 
-* NZ Census 2013 data (from `r sParams$censusDataSource`) and
+* NZ Census 2013 data (from `r p_Params$censusDataSource`) and
 * NZ GREENGrid household power demand data (from [@anderson_new_2018])
 
 To develop initial local area estimates of temporal (half-hourly) power demand in
@@ -78,10 +78,25 @@ This paper will present preliminary results from the development of a neighbourh
 ```{r loadNonPowerData, child="_loadNonPowerData.Rmd"}
 ```
 
-# Create Synthetic census
+# Load IPF weights
 
-```{r includeCensusIPF, child="_createSyntheticCensus.Rmd"}
+## Census 2013
 
+These are created using the GREEN Grid sample, Census 2013 and IPF
+
+```{r loadIPF2013}
+p_Params$ipf2013 <- paste0(p_Params$ggPath, "/safe/ipf/nonZeroWeightsAu2013.csv")
+
+
+file.info(p_Params$ipf2013)
+
+ipfWeights2013DT <- data.table::fread(p_Params$ipf2013)
+
+st <- summary(ipfWeights2013DT)
+
+kableExtra::kable(st, digits = 2,
+                  caption = "Summary of IPF 2013 weights") %>%
+  kable_styling()
 ```