Matt Harris 9/8/2021

library(tmap) # mapping, install if you don't have it

This assignment is for you to complete a short version of the lab notes, but you have to complete a number of the steps yourself. You will then knit this to a markdown (not an HTML) and push it to your GitHub repo. Unlike HTML, the RMarkdown knit to github_document can be viewed directly on GitHub. You will them email your lab instructor with a link to your repo.

Steps in this assignment:

  1. Make sure you have successfully read, run, and learned from the MUSA_508_Lab2_sf.Rmd Rmarkdown

  2. Find two new variables from the 2019 ACS data to load. Use vars <- load_variables(2019, "acs5") and View(vars) to see all of the variable from that ACS. Note that you should not pick something really obscure like count_38yo_cabinetmakers because you will get lots of NAs.

  3. Pick a neighborhood of the City to map. You will need to do some googling to figure this out. Use the PHL Track Explorer to get the GEOID10 number from each parcel and add them to the myTracts object below. This is just like what was done in the exercise, but with a different neighborhood of your choice. Remember that all GEOIDs need to be 10-characters long.

  4. In the first code chunk you will do that above and then edit the call-outs in the dplyr pipe sequence to rename and mutate your data.

  5. You will transform the data to WGS84 by adding the correct EPSG code. This is discussed heavily in the exercise.

  6. You will produce a map of one of the variables you picked and highlight the neighborhood you picked. There are call-out within the ggplot code for you to edit.

  7. You can run the code chunks and lines of code as you edit to make sure everything works.

  8. Once you are done, hit the knit button at the top of the script window (little blue knitting ball) and you will see the output. Once it is what you want…

  9. Use the Git tab on the bottom left of right (depending on hour your Rstudio is laid out) and click the check box to stage all of your changes, write a commit note, hit the commit button, and then the Push button to push it to Github.

  10. Check your Github repo to see you work in the cloud.

  11. Email your lab instructor with a link!

  12. Congrats! You made a map in code!

Load data from {tidycensus}

census_api_key("791448772c9a051612b70516247f56b54176cfbf", overwrite = TRUE)

vars <- load_variables(2019, "acs5")


#B08006_048 Estimate!!Total:!!Female:!!Bicycle SEX OF WORKERS BY MEANS OF TRANSPORTATION TO WORK
#B08006_009 Estimate!!Total:!!Public transportation (excluding taxicab):!!Bus SEX OF WORKERS BY MEANS OF TRANSPORTATION TO WORK
#B08006_043 Estimate!!Total:!!Female:!!Public transportation (excluding taxicab):!!Bus SEX OF WORKERS BY MEANS OF TRANSPORTATION TO Work

acs_vars_hw <- c( "B08006_048", "B08006_043", "B08006_009", "B08006_014") 

#Washington Square West
myTracts_hw <- c("42101000600", 

acsTractsPHL.2019.sf <- get_acs(geography = "tract",
                             year = 2019,
                             variables = acs_vars_hw,
                             geometry = TRUE,
                             state  = "PA",
                             county = "Philadelphia",
                             output = "wide") %>%
  dplyr::select (GEOID, NAME, all_of(paste0(acs_vars_hw,"E"))) %>%
  rename(womenbike = B08006_048E,
         totalbike = B08006_014E,
        womenbus = B08006_043E,
        totalbus = B08006_009E) %>%
  mutate(pct_womenbus = (womenbus/totalbus)*100 )%>%
  mutate(pct_womenbike = (womenbike/totalbike)*100)%>%

  mutate(Neighborhood = ifelse(GEOID %in% myTracts_hw,
                               "REST OF PHILADELPHIA"))%>%


Transform to WGS84 with {sf}

acsTractsPHL.2019.sf <- acsTractsPHL.2019.sf %>% 
  st_transform(crs = "EPSG:26918")

Plot with {ggplot2}

## Loading required package: viridisLite