library(ggplot2)
library(ggthemes)
library(socviz)
library(maps)
library(mapproj)
library(viridis)
library(rsconnect)
library(rmarkdown)Lab7
Subways In Arkansas
library(tidyverse)
library(tidycensus)
library(leaflet)
library(stringr)
census_api_key("91d75f2026ef71780c02dbff95f245e8c42f2f15")
WY_pop <-
get_acs(geography = "county",
variables = "B01003_001",
state = "WY",
geometry = TRUE)
|
| | 0%
|
| | 1%
|
|= | 1%
|
|= | 2%
|
|== | 2%
|
|== | 3%
|
|=== | 4%
|
|=== | 5%
|
|==== | 5%
|
|==== | 6%
|
|===== | 7%
|
|===== | 8%
|
|====== | 8%
|
|====== | 9%
|
|======= | 10%
|
|======== | 11%
|
|======== | 12%
|
|========= | 13%
|
|========== | 14%
|
|========== | 15%
|
|=========== | 16%
|
|============ | 17%
|
|============ | 18%
|
|============= | 19%
|
|============== | 19%
|
|============== | 20%
|
|=============== | 21%
|
|=============== | 22%
|
|================ | 22%
|
|================ | 23%
|
|================= | 24%
|
|================= | 25%
|
|================== | 26%
|
|=================== | 27%
|
|=================== | 28%
|
|==================== | 29%
|
|===================== | 29%
|
|===================== | 30%
|
|====================== | 31%
|
|====================== | 32%
|
|======================= | 32%
|
|======================= | 33%
|
|======================== | 34%
|
|======================== | 35%
|
|========================= | 35%
|
|========================= | 36%
|
|========================== | 37%
|
|========================== | 38%
|
|=========================== | 38%
|
|=========================== | 39%
|
|============================ | 40%
|
|============================= | 41%
|
|============================== | 42%
|
|============================== | 43%
|
|=============================== | 44%
|
|=============================== | 45%
|
|================================ | 46%
|
|================================= | 46%
|
|================================= | 47%
|
|================================== | 48%
|
|================================== | 49%
|
|=================================== | 50%
|
|==================================== | 51%
|
|==================================== | 52%
|
|===================================== | 52%
|
|===================================== | 53%
|
|====================================== | 54%
|
|====================================== | 55%
|
|======================================= | 55%
|
|======================================= | 56%
|
|======================================== | 57%
|
|======================================== | 58%
|
|========================================= | 58%
|
|========================================= | 59%
|
|========================================== | 60%
|
|========================================== | 61%
|
|=========================================== | 61%
|
|============================================ | 62%
|
|============================================ | 63%
|
|============================================= | 64%
|
|============================================== | 65%
|
|============================================== | 66%
|
|=============================================== | 67%
|
|=============================================== | 68%
|
|================================================ | 68%
|
|================================================ | 69%
|
|========================================================= | 81%
|
|========================================================== | 83%
|
|=========================================================== | 84%
|
|=========================================================== | 85%
|
|============================================================ | 86%
|
|============================================================= | 87%
|
|============================================================= | 88%
|
|============================================================== | 89%
|
|=============================================================== | 89%
|
|================================================================ | 91%
|
|================================================================ | 92%
|
|================================================================= | 92%
|
|================================================================= | 93%
|
|================================================================== | 94%
|
|================================================================== | 95%
|
|=================================================================== | 95%
|
|=================================================================== | 96%
|
|==================================================================== | 97%
|
|==================================================================== | 98%
|
|===================================================================== | 98%
|
|===================================================================== | 99%
|
|======================================================================| 100%head(as_tibble(WY_pop))# A tibble: 6 × 6
GEOID NAME variable estimate moe geometry
<chr> <chr> <chr> <dbl> <dbl> <MULTIPOLYGON [°]>
1 56045 Weston County, Wyoming B01003_… 6870 NA (((-105.0808 43.96622, -…
2 56029 Park County, Wyoming B01003_… 29878 NA (((-111.0546 45.00095, -…
3 56013 Fremont County, Wyomi… B01003_… 39402 NA (((-110.0532 44.00802, -…
4 56003 Big Horn County, Wyom… B01003_… 11690 NA (((-108.6218 44.99613, -…
5 56007 Carbon County, Wyoming B01003_… 14609 NA (((-107.93 41.49996, -10…
6 56035 Sublette County, Wyom… B01003_… 8801 NA (((-110.6234 42.7792, -1…MapPalette <- colorQuantile(palette = "magma", domain = WY_pop$estimate, n = 20)
library(sf)
WY_pop %>%
st_transform(crs = "+proj=longlat +datum=WGS84") %>%
leaflet(width = "100%", height = 500) %>%
addProviderTiles(provider = "Esri.WorldPhysical") %>%
addPolygons(popup = ~NAME,
stroke = FALSE,
smoothFactor = 0,
fillOpacity = 0.7,
color = ~ MapPalette(estimate)) %>%
addLegend("bottomright",
pal = MapPalette,
values = ~ estimate,
title = "Population Percentiles",
opacity = 1) Northwest Arkansas and Central Arkansas have the biggiest population percentiles. It seems the outer right hand portion of Arkansas also has a high population percentile. As you move down the map, the population percentiles get lower and lower.
subway_locations_in_us <- read.csv("subway_locations_in_us.csv")
subbar <- subway_locations_in_us %>% filter(state=="WY")
subbar %>% leaflet(width = "100%") %>%
addTiles() %>%
setView(-108, 43, zoom = 6) %>%
addMarkers(lat = ~latitude,
lng = ~longitude,
popup = subbar$Riverton)Northwest Arkansas, FortSmith Area, and Central Arkansas have the biggest cluster of subways. The other subways are spread out and mainly on big highways. Overall, a lot of subways are found throughout Arkansas.
MapPalette <- colorQuantile(palette = "viridis", domain = WY_pop$estimate, n = 20)
WY_pop %>%
st_transform(crs = "+proj=longlat +datum=WGS84") %>%
leaflet(width = "100%", height = 500) %>%
addProviderTiles(provider = "Esri.WorldPhysical") %>%
addPolygons(popup = ~NAME,
stroke = FALSE,
smoothFactor = 0,
fillOpacity = 0.7,
color = ~ MapPalette(estimate)) %>%
addLegend("bottomright",
pal = MapPalette,
values = ~ estimate,
title = "Population Percentiles",
opacity = 1) %>%
addCircleMarkers(data = subbar,
lat = subbar$latitude,
lng = subbar$longitude,
popup = subbar$name,
weight = 1,
radius=3,
color = "white",
opacity = 9)