Add regionally-summed variables to community_climate_vulnerability

data-raw/get_community_climate_vulnerability.R

@@ -9,12 +9,14 @@ raw.dir <- here::here("data-raw")
 
 # Define input files
 comm_clim_vuln_xlsx <- "NE_CCCVI_00_22.xlsx"
+comm_clim_subregion_xlsx <- "SOE SubRegions_CCCRI_Clean.xlsx"
 
 get_community_climate_vulnerability <- function(save_clean = F){
   # Create data frame to define State-EPU
   df<- data.frame(State = c("ME", "MA", "RI", "CT", "DE", "NH", "NY", "NJ", "MD", "VA", "NC"),
                   EPU = c("NE", "NE", "NE", "NE", "MAB", "NE","MAB","MAB","MAB","MAB","MAB"))
-
+
+  # Transform vulnerability input file into SOE format
   community_climate_vulnerability <- readxl::read_excel(file.path(raw.dir,comm_clim_vuln_xlsx),
                                                         sheet = "NE_CCCVI_Final00_22") |>
     dplyr::mutate(Units = NA_character_) |>
@@ -41,8 +43,63 @@ get_community_climate_vulnerability <- function(save_clean = F){
     left_join(df, by = "State") |>
     tidyr::unite("PortState", c(Port, State), sep = ", ") |>
     tidyr::unite("Var", c(PortState,Var), sep = "-") |>
-    dplyr::select(Time, Var, Value, EPU, Units)  
-
+    dplyr::select(Time, Var, Value, EPU, Units)
+
+  # Transform NE subregion input file into SOE format
+  community_climate_newengland <- readxl::read_excel(file.path(raw.dir,comm_clim_subregion_xlsx),
+                                                        sheet = "NewEngland SOE_CCCRI_Clean") |>
+    dplyr::select("YEAR","TempSum","STCKSum","OASum","SENSSum","VULNSum","TempSumLB","STCKSumLB","OASumLB","SENSSumLB","VULNSumLB") |>
+    tidyr::pivot_longer(cols = c("TempSum",
+                                 "STCKSum",
+                                 "OASum",
+                                 "SENSSum",
+                                 "VULNSum",
+                                 "TempSumLB",
+                                 "STCKSumLB",
+                                 "OASumLB",
+                                 "SENSSumLB",
+                                 "VULNSumLB"),
+                        names_to = "Var", values_to = "Value") |>
+    dplyr::mutate(EPU = "NE",
+                  Var2 = "Regional",
+                  Units = "value") |>
+    dplyr::mutate(Units = case_when(stringr::str_detect(Var, "LB") ~ "pounds",
+                                    stringr::str_detect(Var, "LB", negate = T) ~ "value")) |>
+    dplyr::rename("Time" = "YEAR") |>
+    tidyr::unite("Var", c(Var2,Var), sep = "-") |>
+    dplyr::select(Time, Var, Value, EPU, Units)
+
+  # Transform MAB subregion input file into SOE format
+  community_climate_midatlantic <- readxl::read_excel(file.path(raw.dir,comm_clim_subregion_xlsx),
+                                                     sheet = "MidAtlantic SOE_CCCRI_Clean") |>
+    dplyr::select("YEAR","TempSum","STCKSum","OASum","SENSSum","VULNSum","TempSumLB","STCKSumLB","OASumLB","SENSSumLB","VULNSumLB") |>
+    tidyr::pivot_longer(cols = c("TempSum",
+                                 "STCKSum",
+                                 "OASum",
+                                 "SENSSum",
+                                 "VULNSum",
+                                 "TempSumLB",
+                                 "STCKSumLB",
+                                 "OASumLB",
+                                 "SENSSumLB",
+                                 "VULNSumLB"),
+                        names_to = "Var", values_to = "Value") |>
+    dplyr::mutate(EPU = "MAB",
+                  Var2 = "Regional",
+                  Units = "value") |>
+    dplyr::mutate(Units = case_when(stringr::str_detect(Var, "LB") ~ "pounds",
+                                    stringr::str_detect(Var, "LB", negate = T) ~ "value")) |>
+    dplyr::rename("Time" = "YEAR") |>
+    tidyr::unite("Var", c(Var2,Var), sep = "-") |>
+    dplyr::select(Time, Var, Value, EPU, Units)
+
+  community_climate_vulnerability <- rbind(community_climate_vulnerability,
+                                           community_climate_newengland,
+                                           community_climate_midatlantic)
+
+
+
+
     if (save_clean){
       usethis::use_data(community_climate_vulnerability, overwrite = T)
     } else {