Libraries

Read CSV

Filter PA VSOs (NTEECC W30)

Examine Centre County (FIPS 42027)

Filter by VSO Type

Visualize Breakdown of PA VSO Types

Summarize

Map

Libraries

Load Census API Key (Hidden)

ACS Variables

PA Data

Map

library(dplyr)

Summarize by County FIPS

# Merge ACS Vet Estimate with VSO sums by county FIPS
comb_map_data <- merge(pa_acs_B21001_002_map_data, fips_vso_sum, by.x = "GEOID", by.y = "county_fips")

comb_map_data <-
comb_map_data %>%
  mutate(vets_per_vso = estimate / vsos)

comb_map_data

Summarize PA Legions by County FIPS

fips_legion_sum <- 
  pa_vso_legion %>%
    group_by(FIPS) %>%
    tally(sort = TRUE, name = "vsos") %>%
    rename(county_fips = FIPS)

#fips_legion_sum <- drop_na(fips_legion_sum)
fips_legion_sum$county_fips <- gsub("\\..*","",fips_legion_sum$county_fips)
fips_legion_sum$county_fips <- sprintf("%05s", fips_legion_sum$county_fips)

#fips_legion_sum

# Merge ACS Vet Estimate with VSO sums by county FIPS
comb_legion_map_data <- merge(pa_acs_B21001_002_map_data, fips_legion_sum, by.x = "GEOID", by.y = "county_fips")

comb_legion_map_data <-
  comb_legion_map_data %>%
    mutate(vets_per_vso = estimate / vsos)

comb_legion_map_data

Map

mapview(comb_map_data, zcol = "vets_per_vso", layer.name = "All PA WC30 BM VSOs - Total Vets Per VSO")
mapview(comb_legion_map_data, zcol = "vets_per_vso", layer.name = "PA Legions - Total Estimated Vets Per Legion")