library(tidyverse)
library(broom)
library(plotly)
library(tidycensus) # gets census data that we can use to create maps
library(sf) # helper package for mappingLinking to GEOS 3.9.1, GDAL 3.4.0, PROJ 8.1.1; sf_use_s2() is TRUElibrary(leaflet) # interactive mapping package
library(trendyy)
library(usdata) # this package has a conversion utility for state abbreviations to full namesAssignment
census_api_key("686d80e4f7ee63c3023fcefac14d068f56cc4d19", install = TRUE)Error: A CENSUS_API_KEY already exists. You can overwrite it with the argument overwrite=TRUEI first entered my census api key.
smokers <- read_csv("Cig_smoking_percent.csv")Rows: 52 Columns: 5
── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): State
dbl (3): Cig_percent, Low_Confidence_Limit, High_Confidence_Limit
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.I then read in the excel data for Cigarette Smokers
states_leaflet <- get_acs(geography = "state",
variables = "cig_percent",
geometry = TRUE) Getting data from the 2016-2020 5-year ACS
Downloading feature geometry from the Census website. To cache shapefiles for use in future sessions, set `options(tigris_use_cache = TRUE)`.
Error: Your API call has errors. The API message returned is error: error: unknown variable 'cig_percentE'.I used the code above to obtain state population and map data from the census bureau.
nicorette_colors <- colorNumeric(palette = "viridis", domain = nicorette_data$Cig_percent)
states_leaflet %>%
rename(location = NAME) %>%
inner_join(nicorette_data) %>%
leaflet() %>%
addTiles() %>%
addPolygons(weight = 1,
fillColor = ~nicorette_colors(Cig_percent),
label = ~paste0(location, ", Cig Percent = ", Cig_percent),
highlight = highlightOptions(weight = 2)) %>%
setView(-95, 40, zoom = 4) %>%
addLegend(pal = nicorette_colors, values = ~Cig_percent)Joining, by = "location"Warning: sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
Need '+proj=longlat +datum=WGS84'Lastly, I created a chloropleth map using leaflet to show the percent of cigarette smokers among US states that you can hover over to display the information.
nicorette <- trendy("nicorette",
geo = "US",
from = "2020-01-01", to = "2021-01-01")
nicorette_states <- nicorette %>%
get_interest_region()
nicorette_statesI chose the Google search term nicorette because I felt it would be closely related to populations of smokers as most would probably like to quit the habit at some point or another.
nicorette_colors <- colorNumeric(palette = "viridis", domain = nicorette_states$hits)
states_leaflet %>%
rename(location = NAME) %>%
inner_join(nicorette_states) %>%
leaflet() %>%
addTiles() %>%
addPolygons(weight = 1,
fillColor = ~nicorette_colors(hits),
label = ~paste0(location, ", Search volume for nicorette = ", hits),
highlight = highlightOptions(weight = 2)) %>%
setView(-95, 40, zoom = 4) %>%
addLegend(pal = nicorette_colors, values = ~hits)Joining, by = "location"Warning: sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
Need '+proj=longlat +datum=WGS84'