nuts_use_case.R

#
rm(list=ls())
library(fadnUtils)
requiredPackages = c('fadnUtils','data.table', 'devtools','jsonlite', 'ggplot2', 'tidyverse')
for(p in requiredPackages){
  if(!require(p,character.only = TRUE)) install.packages(p)
  library(p,character.only = TRUE)
}

CurrentProjectDirectory = "D:/data/fadn/lieferung_20210414/yang/fadn_work_space"

set.data.dir(CurrentProjectDirectory)

# load str data
str.dir <-  "str_dir"
str_data <- readRDS(paste0(get.data.dir(),"/rds/",str.dir,"/fadn.str.all.rds"))

fadnUtils::nuts.heatmap.group(str_data$info, "NUTS1", countries = "DEU", onepage = FALSE)
fadnUtils::nuts.heatmap.group(str_data$info, "NUTS2", countries = "DEU", onepage = FALSE)
fadnUtils::nuts.heatmap.group(str_data$info, "NUTS3")
fadnUtils::nuts.heatmap.group(str_data$info, "REGION", onepage = TRUE)

# todo ----
#  Konvertierung

nuts1.trans



str_data$info %>%  left_join(nuts2.trans, by = c( "NUTS2", "COUNTRY")) %>%
  mutate(regional_surface_nuts2 = if_else(is.na(regional_surface_nuts2),1,regional_surface_nuts2),
         change_nuts2 = if_else(is.na(change_nuts2),"no change", change_nuts2),
         NUTS2_new = if_else(is.na(NUTS2_new), NUTS2, NUTS2_new )) %>%
  mutate_at(vars(UAA,UAAOWNED), list( new = ~. * regional_surface_nuts2) ) %>%
  # mutate_at(vars(UAA,UAAOWNED), list( ~. * regional_surface_nuts2) ) %>%
  select(-c(TF8:REGION,SIZEUR,NUTS1,NUTS3,ALTITUDE,NAT2000:INCOME.TAXES)) %>%
  filter (change_nuts2 != "no change")