Minnie Bell - Binge Drinking Choropleth Maps
This analysis will be taking a look at binge drinking in the United States compared to Google Searches for ‘frat parties’ and ‘college alcohol.’
library(tidyverse)
library(broom)
library(plotly)
library(tidycensus) # gets census data that we can use to create maps
library(sf) # helper package for mapping
library(leaflet) # interactive mapping package
library(trendyy)
library(usdata) # this package has a conversion utility for state abbreviations to full names
v19 <- load_variables(2019, "acs5", cache = TRUE)
states <- get_acs(geography = "state", # gets state by state data
variables = "B01003_001", # this is state population
geometry = TRUE, # gets geometry (the maps)
shift_geo = T) # shifts Hawaii and AlaskaGetting data from the 2015-2019 5-year ACS
Using feature geometry obtained from the albersusa package
Please note: Alaska and Hawaii are being shifted and are not to scale.statesFirst off, the words ‘binge drinking’ make me think of the student population in Montana. The maps below show the population of undergrad and graduates students in Montana counties. I retrieved the codes from the most recent ACS survey. The undergrad code represents those enrolled in undergraduate courses and includes up to four races. The grad code was also taken from the same ACS survey, and from the same list.
undergrad_pop_MT <- get_acs(geography = "county",
state = "MT",
variables = "B14007G_017",
geometry = TRUE) Getting data from the 2015-2019 5-year ACS
Downloading feature geometry from the Census website. To cache shapefiles for use in future sessions, set `options(tigris_use_cache = TRUE)`.
Using FIPS code '30' for state 'MT'
|
| | 0%
|
|= | 1%
|
|== | 1%
|
|== | 2%
|
|=== | 2%
|
|=== | 3%
|
|==== | 3%
|
|==== | 4%
|
|===== | 4%
|
|===== | 5%
|
|====== | 5%
|
|======= | 6%
|
|======== | 6%
|
|======== | 7%
|
|========= | 7%
|
|========= | 8%
|
|========== | 8%
|
|========== | 9%
|
|=========== | 9%
|
|=========== | 10%
|
|============ | 10%
|
|============= | 11%
|
|============== | 11%
|
|============== | 12%
|
|=============== | 12%
|
|=============== | 13%
|
|================ | 13%
|
|================ | 14%
|
|================= | 14%
|
|================= | 15%
|
|================== | 15%
|
|=================== | 16%
|
|==================== | 16%
|
|==================== | 17%
|
|===================== | 17%
|
|===================== | 18%
|
|====================== | 18%
|
|====================== | 19%
|
|======================= | 19%
|
|======================= | 20%
|
|======================== | 20%
|
|========================= | 21%
|
|========================== | 22%
|
|=========================== | 22%
|
|=========================== | 23%
|
|============================ | 23%
|
|============================ | 24%
|
|============================= | 24%
|
|============================= | 25%
|
|============================== | 25%
|
|============================== | 26%
|
|=============================== | 26%
|
|================================ | 27%
|
|================================= | 27%
|
|================================= | 28%
|
|================================== | 28%
|
|================================== | 29%
|
|=================================== | 29%
|
|=================================== | 30%
|
|==================================== | 30%
|
|==================================== | 31%
|
|===================================== | 31%
|
|====================================== | 32%
|
|======================================= | 32%
|
|======================================= | 33%
|
|======================================== | 33%
|
|======================================== | 34%
|
|========================================= | 34%
|
|========================================= | 35%
|
|========================================== | 35%
|
|========================================== | 36%
|
|=========================================== | 36%
|
|=========================================== | 37%
|
|============================================ | 37%
|
|============================================= | 37%
|
|============================================= | 38%
|
|============================================== | 38%
|
|============================================== | 39%
|
|=============================================== | 39%
|
|=============================================== | 40%
|
|================================================ | 40%
|
|================================================ | 41%
|
|================================================= | 41%
|
|================================================= | 42%
|
|================================================== | 42%
|
|=================================================== | 43%
|
|==================================================== | 43%
|
|==================================================== | 44%
|
|===================================================== | 44%
|
|===================================================== | 45%
|
|====================================================== | 45%
|
|====================================================== | 46%
|
|======================================================= | 46%
|
|======================================================= | 47%
|
|======================================================== | 47%
|
|========================================================= | 48%
|
|========================================================== | 48%
|
|========================================================== | 49%
|
|=========================================================== | 49%
|
|=========================================================== | 50%
|
|============================================================ | 50%
|
|============================================================ | 51%
|
|============================================================= | 51%
|
|============================================================= | 52%
|
|============================================================== | 52%
|
|=============================================================== | 53%
|
|================================================================ | 53%
|
|================================================================ | 54%
|
|================================================================= | 54%
|
|================================================================= | 55%
|
|================================================================== | 55%
|
|================================================================== | 56%
|
|=================================================================== | 56%
|
|=================================================================== | 57%
|
|==================================================================== | 57%
|
|===================================================================== | 58%
|
|====================================================================== | 58%
|
|====================================================================== | 59%
|
|======================================================================= | 59%
|
|======================================================================= | 60%
|
|======================================================================== | 60%
|
|======================================================================== | 61%
|
|========================================================================= | 61%
|
|========================================================================= | 62%
|
|========================================================================== | 62%
|
|========================================================================== | 63%
|
|=========================================================================== | 63%
|
|============================================================================ | 63%
|
|============================================================================ | 64%
|
|============================================================================= | 64%
|
|============================================================================= | 65%
|
|============================================================================== | 65%
|
|============================================================================== | 66%
|
|=============================================================================== | 66%
|
|=============================================================================== | 67%
|
|================================================================================ | 67%
|
|================================================================================ | 68%
|
|================================================================================= | 68%
|
|================================================================================== | 69%
|
|=================================================================================== | 69%
|
|=================================================================================== | 70%
|
|==================================================================================== | 70%
|
|==================================================================================== | 71%
|
|===================================================================================== | 71%
|
|===================================================================================== | 72%
|
|====================================================================================== | 72%
|
|====================================================================================== | 73%
|
|======================================================================================= | 73%
|
|======================================================================================== | 74%
|
|========================================================================================= | 74%
|
|========================================================================================= | 75%
|
|========================================================================================== | 75%
|
|========================================================================================== | 76%
|
|=========================================================================================== | 76%
|
|=========================================================================================== | 77%
|
|============================================================================================ | 77%
|
|============================================================================================ | 78%
|
|============================================================================================= | 78%
|
|============================================================================================== | 79%
|
|=============================================================================================== | 80%
|
|================================================================================================ | 80%
|
|================================================================================================ | 81%
|
|================================================================================================= | 81%
|
|================================================================================================= | 82%
|
|================================================================================================== | 82%
|
|================================================================================================== | 83%
|
|=================================================================================================== | 83%
|
|=================================================================================================== | 84%
|
|==================================================================================================== | 84%
|
|===================================================================================================== | 85%
|
|====================================================================================================== | 85%
|
|====================================================================================================== | 86%
|
|======================================================================================================= | 86%
|
|======================================================================================================= | 87%
|
|======================================================================================================== | 87%
|
|======================================================================================================== | 88%
|
|========================================================================================================= | 88%
|
|========================================================================================================= | 89%
|
|========================================================================================================== | 89%
|
|=========================================================================================================== | 90%
|
|============================================================================================================ | 90%
|
|============================================================================================================ | 91%
|
|============================================================================================================= | 91%
|
|============================================================================================================= | 92%
|
|============================================================================================================== | 92%
|
|============================================================================================================== | 93%
|
|=============================================================================================================== | 93%
|
|=============================================================================================================== | 94%
|
|================================================================================================================ | 94%
|
|================================================================================================================= | 95%
|
|================================================================================================================== | 95%
|
|================================================================================================================== | 96%
|
|=================================================================================================================== | 96%
|
|=================================================================================================================== | 97%
|
|==================================================================================================================== | 97%
|
|==================================================================================================================== | 98%
|
|===================================================================================================================== | 98%
|
|===================================================================================================================== | 99%
|
|====================================================================================================================== | 99%
|
|====================================================================================================================== | 100%
|
|=======================================================================================================================| 100%undergrad_pop_MT %>%
ggplot() +
geom_sf(aes(fill = estimate), color = NA) +
coord_sf(datum = NA) +
theme_minimal() +
scale_fill_viridis_c() +
labs(title = "UnderGrad Student Population in Montana by County")
grad_pop_MT <- get_acs(geography = "county",
state = "MT",
variables = "B14007G_018",
geometry = TRUE) Getting data from the 2015-2019 5-year ACS
Downloading feature geometry from the Census website. To cache shapefiles for use in future sessions, set `options(tigris_use_cache = TRUE)`.
Using FIPS code '30' for state 'MT'grad_pop_MT %>%
ggplot() +
geom_sf(aes(fill = estimate), color = NA) +
coord_sf(datum = NA) +
theme_minimal() +
scale_fill_viridis_c() +
labs(title = "Grad Student Population in Montana by County")
Both of the maps above indicate that Missoula County has the most grads and undergrads. The University of Montana is located in Missoula which is in Missoula County.
Next, I searched for words related to binge drinking. ‘Alcohol’ came to mind first however, it did not seem specific enough. ‘College binge drinking’ did not have enough hits, which was surprising, and ‘frat parties’ did not have many more than that, which was even more surprising! College alcohol came back with quite a few more hits than the others.
alcohol <- trendy("alcohol",
geo = "US",
from = "2019-01-01", to = "2020-01-01")
alcohol_states <- alcohol %>%
get_interest_region()
alcohol_statesNAcollegebinging <- trendy("college binge drinking",
geo = "US",
from = "2019-01-01", to = "2020-01-01")
collegebinging_states <- collegebinging %>%
get_interest_region()
collegebinging_statesNAfrats <- trendy("frat parties",
geo = "US",
from = "2019-01-01", to = "2020-01-01")
frats_states <- frats %>%
get_interest_region()
frats_statesNAcollege_alcohol <- trendy("college alcohol",
geo = "US",
from = "2019-01-01", to = "2020-01-01")
college_alcohol_states <- college_alcohol %>%
get_interest_region()
college_alcohol_statesNAstates %>%
rename(location = NAME) %>%
inner_join(college_alcohol_states)Joining, by = "location"Simple feature collection with 51 features and 9 fields
Geometry type: MULTIPOLYGON
Dimension: XY
Bounding box: xmin: -2100000 ymin: -2500000 xmax: 2516374 ymax: 732103.3
CRS: +proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs
First 10 features:
GEOID location variable estimate moe hits keyword geo gprop geometry
1 04 Arizona B01003_001 7050299 NA 34 college alcohol US web MULTIPOLYGON (((-1111066 -8...
2 05 Arkansas B01003_001 2999370 NA 24 college alcohol US web MULTIPOLYGON (((557903.1 -1...
3 06 California B01003_001 39283497 NA 49 college alcohol US web MULTIPOLYGON (((-1853480 -9...
4 08 Colorado B01003_001 5610349 NA 23 college alcohol US web MULTIPOLYGON (((-613452.9 -...
5 09 Connecticut B01003_001 3575074 NA 100 college alcohol US web MULTIPOLYGON (((2226838 519...
6 11 District of Columbia B01003_001 692683 NA NA college alcohol US web MULTIPOLYGON (((1960720 -41...
7 13 Georgia B01003_001 10403847 NA 24 college alcohol US web MULTIPOLYGON (((1379893 -98...
8 17 Illinois B01003_001 12770631 NA 23 college alcohol US web MULTIPOLYGON (((868942.5 -2...
9 18 Indiana B01003_001 6665703 NA 32 college alcohol US web MULTIPOLYGON (((1279733 -39...
10 22 Louisiana B01003_001 4664362 NA 25 college alcohol US web MULTIPOLYGON (((1080885 -16...undergrad_leaflet <- get_acs(geography = "state", # gets state by state data
variables = "B14007G_017", # this is state undergrads
geometry = TRUE) # gets geometry (the maps)Getting data from the 2015-2019 5-year ACS
Downloading feature geometry from the Census website. To cache shapefiles for use in future sessions, set `options(tigris_use_cache = TRUE)`.
|
| | 0%
|
|= | 1%
|
|== | 1%
|
|== | 2%
|
|=== | 2%
|
|=== | 3%
|
|==== | 3%
|
|===== | 4%
|
|====== | 5%
|
|======= | 6%
|
|======== | 6%
|
|======== | 7%
|
|========= | 8%
|
|========== | 8%
|
|=========== | 9%
|
|=========== | 10%
|
|============ | 10%
|
|============= | 11%
|
|============== | 12%
|
|=============== | 12%
|
|=============== | 13%
|
|================ | 13%
|
|================ | 14%
|
|================= | 14%
|
|================== | 15%
|
|================== | 16%
|
|=================== | 16%
|
|==================== | 17%
|
|===================== | 18%
|
|====================== | 18%
|
|====================== | 19%
|
|======================= | 19%
|
|======================= | 20%
|
|======================== | 20%
|
|========================= | 21%
|
|========================== | 21%
|
|========================== | 22%
|
|=========================== | 23%
|
|============================ | 24%
|
|============================== | 25%
|
|=============================== | 26%
|
|================================ | 27%
|
|================================= | 28%
|
|================================== | 29%
|
|=================================== | 29%
|
|=================================== | 30%
|
|==================================== | 30%
|
|==================================== | 31%
|
|===================================== | 31%
|
|====================================== | 32%
|
|======================================= | 33%
|
|======================================== | 33%
|
|======================================== | 34%
|
|========================================= | 34%
|
|========================================= | 35%
|
|========================================== | 35%
|
|========================================== | 36%
|
|=========================================== | 36%
|
|=========================================== | 37%
|
|============================================ | 37%
|
|============================================= | 38%
|
|============================================== | 39%
|
|=============================================== | 39%
|
|=============================================== | 40%
|
|================================================ | 40%
|
|================================================ | 41%
|
|================================================= | 41%
|
|================================================== | 42%
|
|=================================================== | 42%
|
|=================================================== | 43%
|
|==================================================== | 43%
|
|==================================================== | 44%
|
|===================================================== | 44%
|
|===================================================== | 45%
|
|====================================================== | 45%
|
|======================================================= | 46%
|
|======================================================== | 47%
|
|========================================================= | 48%
|
|========================================================== | 48%
|
|========================================================== | 49%
|
|=========================================================== | 49%
|
|=========================================================== | 50%
|
|============================================================ | 50%
|
|============================================================ | 51%
|
|============================================================= | 51%
|
|============================================================== | 52%
|
|=============================================================== | 53%
|
|================================================================ | 53%
|
|================================================================ | 54%
|
|================================================================= | 54%
|
|================================================================= | 55%
|
|================================================================== | 55%
|
|================================================================== | 56%
|
|=================================================================== | 56%
|
|=================================================================== | 57%
|
|==================================================================== | 57%
|
|===================================================================== | 58%
|
|====================================================================== | 59%
|
|======================================================================= | 59%
|
|======================================================================= | 60%
|
|======================================================================== | 60%
|
|======================================================================== | 61%
|
|========================================================================= | 61%
|
|========================================================================== | 62%
|
|========================================================================== | 63%
|
|=========================================================================== | 63%
|
|============================================================================ | 64%
|
|============================================================================= | 65%
|
|============================================================================== | 65%
|
|============================================================================== | 66%
|
|=============================================================================== | 66%
|
|================================================================================ | 67%
|
|================================================================================= | 68%
|
|================================================================================== | 69%
|
|=================================================================================== | 70%
|
|==================================================================================== | 70%
|
|==================================================================================== | 71%
|
|===================================================================================== | 71%
|
|===================================================================================== | 72%
|
|====================================================================================== | 72%
|
|======================================================================================= | 73%
|
|======================================================================================== | 74%
|
|========================================================================================= | 74%
|
|========================================================================================= | 75%
|
|========================================================================================== | 76%
|
|=========================================================================================== | 76%
|
|=========================================================================================== | 77%
|
|============================================================================================ | 77%
|
|============================================================================================ | 78%
|
|============================================================================================= | 78%
|
|============================================================================================== | 79%
|
|=============================================================================================== | 80%
|
|================================================================================================ | 80%
|
|================================================================================================ | 81%
|
|================================================================================================= | 82%
|
|================================================================================================== | 82%
|
|=================================================================================================== | 83%
|
|=================================================================================================== | 84%
|
|==================================================================================================== | 84%
|
|===================================================================================================== | 85%
|
|====================================================================================================== | 86%
|
|======================================================================================================= | 86%
|
|======================================================================================================= | 87%
|
|======================================================================================================== | 87%
|
|========================================================================================================= | 88%
|
|========================================================================================================== | 89%
|
|=========================================================================================================== | 90%
|
|============================================================================================================ | 91%
|
|============================================================================================================== | 92%
|
|=============================================================================================================== | 93%
|
|================================================================================================================ | 94%
|
|================================================================================================================= | 95%
|
|================================================================================================================== | 96%
|
|=================================================================================================================== | 96%
|
|=================================================================================================================== | 97%
|
|==================================================================================================================== | 97%
|
|==================================================================================================================== | 98%
|
|===================================================================================================================== | 98%
|
|===================================================================================================================== | 99%
|
|====================================================================================================================== | 99%
|
|=======================================================================================================================| 100% # shift_geo = T # shifts Hawaii and Alaska
grad_leaflet <- get_acs(geography = "state", # gets state by state data
variables = "B14007G_018", # this is state grads
geometry = TRUE) # gets geometry (the maps)Getting data from the 2015-2019 5-year ACS
Downloading feature geometry from the Census website. To cache shapefiles for use in future sessions, set `options(tigris_use_cache = TRUE)`. # shift_geo = T # shifts Hawaii and Alaska
Binge_drinking_percent <- read_csv("Binge_drinking_percent.csv")
-- Column specification ---------------------------------------------------------------------------------------------------------
cols(
State = col_character(),
Percentage = col_double()
)glimpse(Binge_drinking_percent)Rows: 51
Columns: 2
$ State <chr> "AL", "AK", "AZ", "AR", "CA", "CO", "CT", "DE", "DC", "FL", "GA", "HI", "ID", "IL", "IN", "IA", "KS", "KY", ~
$ Percentage <dbl> 12.2, 20.0, 15.0, 15.2, 16.7, 18.1, 18.3, 15.8, 24.4, 17.2, 15.8, 19.8, 14.8, 20.8, 16.6, 21.3, 16.5, 16.1, ~glimpse(college_alcohol_states)Rows: 51
Columns: 5
$ location <chr> "Connecticut", "Pennsylvania", "Massachusetts", "Alabama", "Maryland", "Minnesota", "Oklahoma", "California", ~
$ hits <int> 100, 65, 58, 57, 55, 53, 50, 49, 49, 48, 47, 46, 45, 43, 43, 41, 40, 39, 38, 37, 36, 36, 35, 34, 34, 34, 33, 3~
$ keyword <chr> "college alcohol", "college alcohol", "college alcohol", "college alcohol", "college alcohol", "college alcoho~
$ geo <chr> "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "U~
$ gprop <chr> "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web"~Now a look at how these two search terms relate to the binge drinking percentage of the states. The following map shows the binge drinking percentage by state.
college_colors <- colorNumeric(palette = "viridis", domain = college_alcohol_data$Percentage)
undergrad_leaflet %>%
rename(location = NAME) %>%
inner_join(college_alcohol_data) %>%
leaflet() %>%
addTiles() %>%
addPolygons(weight = 1,
fillColor = ~college_colors(Percentage),
label = ~paste0(location, ", Binge Drinking by Percentage by State = ", Percentage),
highlight = highlightOptions(weight = 2)) %>%
setView(-95, 40, zoom = 4) %>%
addLegend(pal = college_colors, values = ~Percentage)Joining, by = "location"
sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
Need '+proj=longlat +datum=WGS84'In the next map, we will see Google Search hits for college alcohol by state.
college_colors <- colorNumeric(palette = "viridis", domain = college_alcohol_data$hits)
undergrad_leaflet %>%
rename(location = NAME) %>%
inner_join(college_alcohol_data) %>%
leaflet() %>%
addTiles() %>%
addPolygons(weight = 1,
fillColor = ~college_colors(hits),
label = ~paste0(location, ", College Alcohol Hits from Google Search = ", hits),
highlight = highlightOptions(weight = 2)) %>%
setView(-95, 40, zoom = 4) %>%
addLegend(pal = college_colors, values = ~hits)Joining, by = "location"
sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
Need '+proj=longlat +datum=WGS84'Next we will see if there is a correlation between college alcohol and the binge drinking percentage.
college_alcohol_data <- college_alcohol_states %>%
mutate(State = state2abbr(location)) %>%
inner_join(Binge_drinking_percent)
college_alcohol_states %>%
mutate(State = state2abbr(location)) %>%
inner_join(Binge_drinking_percent)college_model <- lm(Percentage ~ hits, data = college_alcohol_data)
glance(college_model)tidy(college_model)sqrt(.1457)[1] 0.3817067The following line graph shows a good correlation between college alcohol search hits and the binge drinking percentage by state, as the p-value also suggests.
college_alcohol_data %>%
drop_na() %>%
plot_ly(x = ~hits,
y = ~Percentage,
hoverinfo = "text",
text = ~paste("State: ", location, "<br>", "'College Alcohol' search rate: ", hits, "<br>", "Percentage: ", Percentage)) %>%
add_markers(showlegend = F) %>%
add_lines(y = ~fitted(college_model)) %>%
layout(title = "Relationship between google searches for 'College Alcohol' and Binge Drinking Rates, by State",
xaxis = list(title = "Google search volume for 'College Alcohol'"),
yaxis = list(title = "State Binge Drinking Rate, per capita"))NA
NANow we will look at Google search hits for frat parties as compared to the binge drinking percentage to see if there is a good correlation there as well.
states %>%
rename(location = NAME) %>%
inner_join(frats_states)Joining, by = "location"Simple feature collection with 51 features and 9 fields
Geometry type: MULTIPOLYGON
Dimension: XY
Bounding box: xmin: -2100000 ymin: -2500000 xmax: 2516374 ymax: 732103.3
CRS: +proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs
First 10 features:
GEOID location variable estimate moe hits keyword geo gprop geometry
1 04 Arizona B01003_001 7050299 NA NA frat parties US web MULTIPOLYGON (((-1111066 -8...
2 05 Arkansas B01003_001 2999370 NA NA frat parties US web MULTIPOLYGON (((557903.1 -1...
3 06 California B01003_001 39283497 NA 47 frat parties US web MULTIPOLYGON (((-1853480 -9...
4 08 Colorado B01003_001 5610349 NA NA frat parties US web MULTIPOLYGON (((-613452.9 -...
5 09 Connecticut B01003_001 3575074 NA NA frat parties US web MULTIPOLYGON (((2226838 519...
6 11 District of Columbia B01003_001 692683 NA NA frat parties US web MULTIPOLYGON (((1960720 -41...
7 13 Georgia B01003_001 10403847 NA NA frat parties US web MULTIPOLYGON (((1379893 -98...
8 17 Illinois B01003_001 12770631 NA 63 frat parties US web MULTIPOLYGON (((868942.5 -2...
9 18 Indiana B01003_001 6665703 NA NA frat parties US web MULTIPOLYGON (((1279733 -39...
10 22 Louisiana B01003_001 4664362 NA NA frat parties US web MULTIPOLYGON (((1080885 -16...glimpse(Binge_drinking_percent)Rows: 51
Columns: 2
$ State <chr> "AL", "AK", "AZ", "AR", "CA", "CO", "CT", "DE", "DC", "FL", "GA", "HI", "ID", "IL", "IN", "IA", "KS", "KY", "LA", "ME~
$ Percentage <dbl> 12.2, 20.0, 15.0, 15.2, 16.7, 18.1, 18.3, 15.8, 24.4, 17.2, 15.8, 19.8, 14.8, 20.8, 16.6, 21.3, 16.5, 16.1, 18.0, 20.~glimpse(frats_states)Rows: 51
Columns: 5
$ location <chr> "New Jersey", "New York", "Pennsylvania", "Virginia", "Illinois", "North Carolina", "Massachusetts", "Florida", "Califo~
$ hits <int> 100, 82, 70, 69, 63, 62, 55, 51, 47, 7, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,~
$ keyword <chr> "frat parties", "frat parties", "frat parties", "frat parties", "frat parties", "frat parties", "frat parties", "frat p~
$ geo <chr> "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US", "US",~
$ gprop <chr> "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", "web", ~frats_data <- frats_states %>%
mutate(State = state2abbr(location)) %>%
inner_join(Binge_drinking_percent)Joining, by = "State"
frats_colors <- colorNumeric(palette = "viridis", domain = frats_data$hits)
undergrad_leaflet %>%
rename(location = NAME) %>%
inner_join(frats_data) %>%
leaflet() %>%
addTiles() %>%
addPolygons(weight = 1,
fillColor = ~frats_colors(hits),
label = ~paste0(location, ", Frat Parties Hits from Google Search = ", hits),
highlight = highlightOptions(weight = 2)) %>%
setView(-95, 40, zoom = 4) %>%
addLegend(pal = frats_colors, values = ~hits)Joining, by = "location"
sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
Need '+proj=longlat +datum=WGS84'The above map shows that there were not many hits for ‘frat parties’ on Google Search.
frats_model <- lm(Percentage ~ hits, data = frats_data)
glance(frats_model)tidy(frats_model)sqrt(.0406)[1] 0.2014944frats_data %>%
drop_na() %>%
plot_ly(x = ~hits,
y = ~Percentage,
hoverinfo = "text",
text = ~paste("State: ", location, "<br>", "'Frat Parties' search rate: ", hits, "<br>", "Percentage: ", Percentage)) %>%
add_markers(showlegend = F) %>%
add_lines(y = ~fitted(frats_model)) %>%
layout(title = "Relationship between google searches for 'Frat Parties' and Binge Drinking Rates, by State",
xaxis = list(title = "Google search volume for 'Frat Parties'"),
yaxis = list(title = "State Binge Drinking Rate, per capita"))NA
NAHowever, the p-values for ‘frat parties’ show a statistically significant relationship to binge drinking percentage by state as well, as shown in the line graph above.
LS0tDQp0aXRsZTogIk1pbm5pZSBCZWxsIC0gQmluZ2UgRHJpbmtpbmcgQ2hvcm9wbGV0aCBNYXBzIg0Kb3V0cHV0OiBodG1sX25vdGVib29rDQotLS0NCg0KVGhpcyBhbmFseXNpcyB3aWxsIGJlIHRha2luZyBhIGxvb2sgYXQgYmluZ2UgZHJpbmtpbmcgaW4gdGhlIFVuaXRlZCBTdGF0ZXMgY29tcGFyZWQgdG8gR29vZ2xlIFNlYXJjaGVzIGZvciAnZnJhdCBwYXJ0aWVzJyBhbmQgJ2NvbGxlZ2UgYWxjb2hvbC4nDQoNCmBgYHtyfQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpsaWJyYXJ5KGJyb29tKQ0KbGlicmFyeShwbG90bHkpDQpsaWJyYXJ5KHRpZHljZW5zdXMpICAgICAgIyBnZXRzIGNlbnN1cyBkYXRhIHRoYXQgd2UgY2FuIHVzZSB0byBjcmVhdGUgbWFwcw0KbGlicmFyeShzZikgICAgICAgICAgICAgICMgaGVscGVyIHBhY2thZ2UgZm9yIG1hcHBpbmcNCmxpYnJhcnkobGVhZmxldCkgICAgICAgICAjIGludGVyYWN0aXZlIG1hcHBpbmcgcGFja2FnZQ0KbGlicmFyeSh0cmVuZHl5KQ0KbGlicmFyeSh1c2RhdGEpICAgICAgICAgICMgdGhpcyBwYWNrYWdlIGhhcyBhIGNvbnZlcnNpb24gdXRpbGl0eSBmb3Igc3RhdGUgYWJicmV2aWF0aW9ucyB0byBmdWxsIG5hbWVzDQoNCmBgYA0KDQoNCmBgYHtyfQ0KdjE5IDwtIGxvYWRfdmFyaWFibGVzKDIwMTksICJhY3M1IiwgY2FjaGUgPSBUUlVFKQ0KYGBgDQoNCg0KDQpgYGB7cn0NCg0Kc3RhdGVzIDwtIGdldF9hY3MoZ2VvZ3JhcGh5ID0gInN0YXRlIiwgICAgICAgICAgICAgICAjIGdldHMgc3RhdGUgYnkgc3RhdGUgZGF0YQ0KICAgICAgICAgICAgICAgICAgdmFyaWFibGVzID0gIkIwMTAwM18wMDEiLCAgICAgICAgICAjIHRoaXMgaXMgc3RhdGUgcG9wdWxhdGlvbg0KICAgICAgICAgICAgICAgICAgZ2VvbWV0cnkgPSBUUlVFLCAgICAgICAgICAgICAgICAgICAjIGdldHMgZ2VvbWV0cnkgKHRoZSBtYXBzKQ0KICAgICAgICAgICAgICAgICAgc2hpZnRfZ2VvID0gVCkgICAgICAgICAgICAgICAgICAgICAjIHNoaWZ0cyBIYXdhaWkgYW5kIEFsYXNrYQ0KDQpgYGANCg0KYGBge3J9DQpzdGF0ZXMNCmBgYA0KDQoNCkZpcnN0IG9mZiwgdGhlIHdvcmRzICdiaW5nZSBkcmlua2luZycgbWFrZSBtZSB0aGluayBvZiB0aGUgc3R1ZGVudCBwb3B1bGF0aW9uIGluIE1vbnRhbmEuIFRoZSBtYXBzIGJlbG93IHNob3cgdGhlIHBvcHVsYXRpb24gb2YgdW5kZXJncmFkIGFuZCBncmFkdWF0ZXMgc3R1ZGVudHMgaW4gTW9udGFuYSBjb3VudGllcy4gSSByZXRyaWV2ZWQgdGhlIGNvZGVzIGZyb20gdGhlIG1vc3QgcmVjZW50IEFDUyBzdXJ2ZXkuIFRoZSB1bmRlcmdyYWQgY29kZSByZXByZXNlbnRzIHRob3NlIGVucm9sbGVkIGluIHVuZGVyZ3JhZHVhdGUgY291cnNlcyBhbmQgaW5jbHVkZXMgdXAgdG8gZm91ciByYWNlcy4gVGhlIGdyYWQgY29kZSB3YXMgYWxzbyB0YWtlbiBmcm9tIHRoZSBzYW1lIEFDUyBzdXJ2ZXksIGFuZCBmcm9tIHRoZSBzYW1lIGxpc3QuDQoNCg0KYGBge3J9DQoNCnVuZGVyZ3JhZF9wb3BfTVQgPC0gZ2V0X2FjcyhnZW9ncmFwaHkgPSAiY291bnR5IiwgDQogICAgICAgICAgICAgICAgICAgICAgICAgc3RhdGUgPSAiTVQiLCANCiAgICAgICAgICAgICAgICAgICAgICAgICB2YXJpYWJsZXMgPSAiQjE0MDA3R18wMTciLCANCiAgICAgICAgICAgICAgICAgICAgICAgICBnZW9tZXRyeSA9IFRSVUUpIA0KDQp1bmRlcmdyYWRfcG9wX01UICU+JSANCiAgZ2dwbG90KCkgKyANCiAgZ2VvbV9zZihhZXMoZmlsbCA9IGVzdGltYXRlKSwgY29sb3IgPSBOQSkgKyANCiAgY29vcmRfc2YoZGF0dW0gPSBOQSkgKyANCiAgdGhlbWVfbWluaW1hbCgpICsgDQogIHNjYWxlX2ZpbGxfdmlyaWRpc19jKCkgKw0KICBsYWJzKHRpdGxlID0gIlVuZGVyR3JhZCBTdHVkZW50IFBvcHVsYXRpb24gaW4gTW9udGFuYSBieSBDb3VudHkiKQ0KDQpgYGANCg0KDQoNCmBgYHtyfQ0KDQpncmFkX3BvcF9NVCA8LSBnZXRfYWNzKGdlb2dyYXBoeSA9ICJjb3VudHkiLCANCiAgICAgICAgICAgICAgICAgICAgICAgICBzdGF0ZSA9ICJNVCIsIA0KICAgICAgICAgICAgICAgICAgICAgICAgIHZhcmlhYmxlcyA9ICJCMTQwMDdHXzAxOCIsIA0KICAgICAgICAgICAgICAgICAgICAgICAgIGdlb21ldHJ5ID0gVFJVRSkgDQoNCmdyYWRfcG9wX01UICU+JSANCiAgZ2dwbG90KCkgKyANCiAgZ2VvbV9zZihhZXMoZmlsbCA9IGVzdGltYXRlKSwgY29sb3IgPSBOQSkgKyANCiAgY29vcmRfc2YoZGF0dW0gPSBOQSkgKyANCiAgdGhlbWVfbWluaW1hbCgpICsgDQogIHNjYWxlX2ZpbGxfdmlyaWRpc19jKCkgKw0KICBsYWJzKHRpdGxlID0gIkdyYWQgU3R1ZGVudCBQb3B1bGF0aW9uIGluIE1vbnRhbmEgYnkgQ291bnR5IikNCg0KYGBgDQoNCg0KQm90aCBvZiB0aGUgbWFwcyBhYm92ZSBpbmRpY2F0ZSB0aGF0IE1pc3NvdWxhIENvdW50eSBoYXMgdGhlIG1vc3QgZ3JhZHMgYW5kIHVuZGVyZ3JhZHMuIFRoZSBVbml2ZXJzaXR5IG9mIE1vbnRhbmEgaXMgbG9jYXRlZCBpbiBNaXNzb3VsYSB3aGljaCBpcyBpbiBNaXNzb3VsYSBDb3VudHkuDQoNCk5leHQsIEkgc2VhcmNoZWQgZm9yIHdvcmRzIHJlbGF0ZWQgdG8gYmluZ2UgZHJpbmtpbmcuICdBbGNvaG9sJyBjYW1lIHRvIG1pbmQgZmlyc3QgaG93ZXZlciwgaXQgZGlkIG5vdCBzZWVtIHNwZWNpZmljIGVub3VnaC4gJ0NvbGxlZ2UgYmluZ2UgZHJpbmtpbmcnIGRpZCBub3QgaGF2ZSBlbm91Z2ggaGl0cywgd2hpY2ggd2FzIHN1cnByaXNpbmcsIGFuZCAnZnJhdCBwYXJ0aWVzJyBkaWQgbm90IGhhdmUgbWFueSBtb3JlIHRoYW4gdGhhdCwgd2hpY2ggd2FzIGV2ZW4gbW9yZSBzdXJwcmlzaW5nISBDb2xsZWdlIGFsY29ob2wgY2FtZSBiYWNrIHdpdGggcXVpdGUgYSBmZXcgbW9yZSBoaXRzIHRoYW4gdGhlIG90aGVycy4NCg0KYGBge3J9DQphbGNvaG9sIDwtIHRyZW5keSgiYWxjb2hvbCIsIA0KICAgICAgICAgICAgICAgICAgIGdlbyA9ICJVUyIsIA0KICAgICAgICAgICAgICAgICAgIGZyb20gPSAiMjAxOS0wMS0wMSIsIHRvID0gIjIwMjAtMDEtMDEiKQ0KDQoNCmFsY29ob2xfc3RhdGVzIDwtIGFsY29ob2wgJT4lDQogIGdldF9pbnRlcmVzdF9yZWdpb24oKQ0KDQphbGNvaG9sX3N0YXRlcw0KDQpgYGANCg0KYGBge3J9DQpjb2xsZWdlYmluZ2luZyA8LSB0cmVuZHkoImNvbGxlZ2UgYmluZ2UgZHJpbmtpbmciLCANCiAgICAgICAgICAgICAgICAgICBnZW8gPSAiVVMiLCANCiAgICAgICAgICAgICAgICAgICBmcm9tID0gIjIwMTktMDEtMDEiLCB0byA9ICIyMDIwLTAxLTAxIikNCg0KDQpjb2xsZWdlYmluZ2luZ19zdGF0ZXMgPC0gY29sbGVnZWJpbmdpbmcgJT4lDQogIGdldF9pbnRlcmVzdF9yZWdpb24oKQ0KDQpjb2xsZWdlYmluZ2luZ19zdGF0ZXMNCg0KYGBgDQoNCg0KDQoNCmBgYHtyfQ0KZnJhdHMgPC0gdHJlbmR5KCJmcmF0IHBhcnRpZXMiLCANCiAgICAgICAgICAgICAgICAgICBnZW8gPSAiVVMiLCANCiAgICAgICAgICAgICAgICAgICBmcm9tID0gIjIwMTktMDEtMDEiLCB0byA9ICIyMDIwLTAxLTAxIikNCg0KDQpmcmF0c19zdGF0ZXMgPC0gZnJhdHMgJT4lDQogIGdldF9pbnRlcmVzdF9yZWdpb24oKQ0KDQpmcmF0c19zdGF0ZXMNCg0KYGBgDQoNCmBgYHtyfQ0KY29sbGVnZV9hbGNvaG9sIDwtIHRyZW5keSgiY29sbGVnZSBhbGNvaG9sIiwgDQogICAgICAgICAgICAgICAgICAgZ2VvID0gIlVTIiwgDQogICAgICAgICAgICAgICAgICAgZnJvbSA9ICIyMDE5LTAxLTAxIiwgdG8gPSAiMjAyMC0wMS0wMSIpDQoNCg0KY29sbGVnZV9hbGNvaG9sX3N0YXRlcyA8LSBjb2xsZWdlX2FsY29ob2wgJT4lDQogIGdldF9pbnRlcmVzdF9yZWdpb24oKQ0KDQpjb2xsZWdlX2FsY29ob2xfc3RhdGVzDQoNCmBgYA0KDQoNCmBgYHtyfQ0Kc3RhdGVzICU+JSANCiAgcmVuYW1lKGxvY2F0aW9uID0gTkFNRSkgJT4lIA0KICBpbm5lcl9qb2luKGNvbGxlZ2VfYWxjb2hvbF9zdGF0ZXMpDQpgYGANCg0KDQoNCmBgYHtyfQ0KdW5kZXJncmFkX2xlYWZsZXQgPC0gZ2V0X2FjcyhnZW9ncmFwaHkgPSAic3RhdGUiLCAgICAgICAjIGdldHMgc3RhdGUgYnkgc3RhdGUgZGF0YQ0KICAgICAgICAgICAgICAgICAgdmFyaWFibGVzID0gIkIxNDAwN0dfMDE3IiwgICAgICAgICAgIyB0aGlzIGlzIHN0YXRlIHVuZGVyZ3JhZHMNCiAgICAgICAgICAgICAgICAgIGdlb21ldHJ5ID0gVFJVRSkgICAgICAgICAgICAgICAgICAgIyBnZXRzIGdlb21ldHJ5ICh0aGUgbWFwcykNCiAgICAgICAgICAgICAgICAgICMgc2hpZnRfZ2VvID0gVCAgICAgICAgICAgICAgICAgICAgIyBzaGlmdHMgSGF3YWlpIGFuZCBBbGFza2ENCg0KYGBgDQoNCg0KDQpgYGB7cn0NCmdyYWRfbGVhZmxldCA8LSBnZXRfYWNzKGdlb2dyYXBoeSA9ICJzdGF0ZSIsICAgICAgICMgZ2V0cyBzdGF0ZSBieSBzdGF0ZSBkYXRhDQogICAgICAgICAgICAgICAgICB2YXJpYWJsZXMgPSAiQjE0MDA3R18wMTgiLCAgICAgICAgICAjIHRoaXMgaXMgc3RhdGUgZ3JhZHMNCiAgICAgICAgICAgICAgICAgIGdlb21ldHJ5ID0gVFJVRSkgICAgICAgICAgICAgICAgICAgIyBnZXRzIGdlb21ldHJ5ICh0aGUgbWFwcykNCiAgICAgICAgICAgICAgICAgICMgc2hpZnRfZ2VvID0gVCAgICAgICAgICAgICAgICAgICAgIyBzaGlmdHMgSGF3YWlpIGFuZCBBbGFza2ENCg0KYGBgDQoNCg0KDQpgYGB7cn0NCkJpbmdlX2RyaW5raW5nX3BlcmNlbnQgPC0gcmVhZF9jc3YoIkJpbmdlX2RyaW5raW5nX3BlcmNlbnQuY3N2IikNCg0KYGBgDQoNCmBgYHtyfQ0KZ2xpbXBzZShCaW5nZV9kcmlua2luZ19wZXJjZW50KQ0KZ2xpbXBzZShjb2xsZWdlX2FsY29ob2xfc3RhdGVzKQ0KYGBgDQoNCk5vdyBhIGxvb2sgYXQgaG93IHRoZXNlIHR3byBzZWFyY2ggdGVybXMgcmVsYXRlIHRvIHRoZSBiaW5nZSBkcmlua2luZyBwZXJjZW50YWdlIG9mIHRoZSBzdGF0ZXMuIFRoZSBmb2xsb3dpbmcgbWFwIHNob3dzIHRoZSBiaW5nZSBkcmlua2luZyBwZXJjZW50YWdlIGJ5IHN0YXRlLg0KDQpgYGB7cn0NCg0KY29sbGVnZV9jb2xvcnMgPC0gY29sb3JOdW1lcmljKHBhbGV0dGUgPSAidmlyaWRpcyIsIGRvbWFpbiA9IGNvbGxlZ2VfYWxjb2hvbF9kYXRhJFBlcmNlbnRhZ2UpDQoNCnVuZGVyZ3JhZF9sZWFmbGV0ICU+JSANCiAgcmVuYW1lKGxvY2F0aW9uID0gTkFNRSkgJT4lIA0KICBpbm5lcl9qb2luKGNvbGxlZ2VfYWxjb2hvbF9kYXRhKSAlPiUgDQogIGxlYWZsZXQoKSAlPiUgDQogIGFkZFRpbGVzKCkgJT4lDQogIGFkZFBvbHlnb25zKHdlaWdodCA9IDEsDQogICAgICAgICAgICAgIGZpbGxDb2xvciA9IH5jb2xsZWdlX2NvbG9ycyhQZXJjZW50YWdlKSwNCiAgICAgICAgICAgICAgbGFiZWwgPSB+cGFzdGUwKGxvY2F0aW9uLCAiLCBCaW5nZSBEcmlua2luZyBieSBQZXJjZW50YWdlIGJ5IFN0YXRlID0gIiwgUGVyY2VudGFnZSksDQogICAgICAgICAgICAgIGhpZ2hsaWdodCA9IGhpZ2hsaWdodE9wdGlvbnMod2VpZ2h0ID0gMikpICU+JSANCiAgc2V0VmlldygtOTUsIDQwLCB6b29tID0gNCkgJT4lIA0KICBhZGRMZWdlbmQocGFsID0gY29sbGVnZV9jb2xvcnMsIHZhbHVlcyA9IH5QZXJjZW50YWdlKQ0KYGBgDQoNCkluIHRoZSBuZXh0IG1hcCwgd2Ugd2lsbCBzZWUgR29vZ2xlIFNlYXJjaCBoaXRzIGZvciBjb2xsZWdlIGFsY29ob2wgYnkgc3RhdGUuDQoNCg0KYGBge3J9DQoNCmNvbGxlZ2VfY29sb3JzIDwtIGNvbG9yTnVtZXJpYyhwYWxldHRlID0gInZpcmlkaXMiLCBkb21haW4gPSBjb2xsZWdlX2FsY29ob2xfZGF0YSRoaXRzKQ0KDQp1bmRlcmdyYWRfbGVhZmxldCAlPiUgDQogIHJlbmFtZShsb2NhdGlvbiA9IE5BTUUpICU+JSANCiAgaW5uZXJfam9pbihjb2xsZWdlX2FsY29ob2xfZGF0YSkgJT4lIA0KICBsZWFmbGV0KCkgJT4lIA0KICBhZGRUaWxlcygpICU+JQ0KICBhZGRQb2x5Z29ucyh3ZWlnaHQgPSAxLA0KICAgICAgICAgICAgICBmaWxsQ29sb3IgPSB+Y29sbGVnZV9jb2xvcnMoaGl0cyksDQogICAgICAgICAgICAgIGxhYmVsID0gfnBhc3RlMChsb2NhdGlvbiwgIiwgQ29sbGVnZSBBbGNvaG9sIEhpdHMgZnJvbSBHb29nbGUgU2VhcmNoID0gIiwgaGl0cyksDQogICAgICAgICAgICAgIGhpZ2hsaWdodCA9IGhpZ2hsaWdodE9wdGlvbnMod2VpZ2h0ID0gMikpICU+JSANCiAgc2V0VmlldygtOTUsIDQwLCB6b29tID0gNCkgJT4lIA0KICBhZGRMZWdlbmQocGFsID0gY29sbGVnZV9jb2xvcnMsIHZhbHVlcyA9IH5oaXRzKQ0KYGBgDQoNCk5leHQgd2Ugd2lsbCBzZWUgaWYgdGhlcmUgaXMgYSBjb3JyZWxhdGlvbiBiZXR3ZWVuIGNvbGxlZ2UgYWxjb2hvbCBhbmQgdGhlIGJpbmdlIGRyaW5raW5nIHBlcmNlbnRhZ2UuDQoNCmBgYHtyfQ0KY29sbGVnZV9hbGNvaG9sX2RhdGEgPC0gY29sbGVnZV9hbGNvaG9sX3N0YXRlcyAlPiUgDQogIG11dGF0ZShTdGF0ZSA9IHN0YXRlMmFiYnIobG9jYXRpb24pKSAlPiUgDQogIGlubmVyX2pvaW4oQmluZ2VfZHJpbmtpbmdfcGVyY2VudCkNCiANCg0KYGBgDQoNCmBgYHtyfQ0KY29sbGVnZV9hbGNvaG9sX3N0YXRlcyAlPiUgDQogIG11dGF0ZShTdGF0ZSA9IHN0YXRlMmFiYnIobG9jYXRpb24pKSAlPiUgDQogIGlubmVyX2pvaW4oQmluZ2VfZHJpbmtpbmdfcGVyY2VudCkNCmBgYA0KDQoNCmBgYHtyfQ0KY29sbGVnZV9tb2RlbCA8LSBsbShQZXJjZW50YWdlIH4gaGl0cywgZGF0YSA9IGNvbGxlZ2VfYWxjb2hvbF9kYXRhKQ0KDQpnbGFuY2UoY29sbGVnZV9tb2RlbCkNCnRpZHkoY29sbGVnZV9tb2RlbCkNCmBgYA0KDQoNCmBgYHtyfQ0Kc3FydCguMTQ1NykNCg0KYGBgDQpUaGUgZm9sbG93aW5nIGxpbmUgZ3JhcGggc2hvd3MgYSBnb29kIGNvcnJlbGF0aW9uIGJldHdlZW4gY29sbGVnZSBhbGNvaG9sIHNlYXJjaCBoaXRzIGFuZCB0aGUgYmluZ2UgZHJpbmtpbmcgcGVyY2VudGFnZSBieSBzdGF0ZSwgYXMgdGhlIHAtdmFsdWUgYWxzbyBzdWdnZXN0cy4NCg0KYGBge3J9DQpjb2xsZWdlX2FsY29ob2xfZGF0YSAlPiUgDQogIGRyb3BfbmEoKSAlPiUgDQogIHBsb3RfbHkoeCA9IH5oaXRzLCANCiAgICAgICAgICB5ID0gflBlcmNlbnRhZ2UsDQogICAgICAgICAgaG92ZXJpbmZvID0gInRleHQiLCANCiAgICAgICAgICB0ZXh0ID0gfnBhc3RlKCJTdGF0ZTogIiwgbG9jYXRpb24sICI8YnI+IiwgIidDb2xsZWdlIEFsY29ob2wnIHNlYXJjaCByYXRlOiAiLCBoaXRzLCAiPGJyPiIsICJQZXJjZW50YWdlOiAiLCBQZXJjZW50YWdlKSkgJT4lIA0KICBhZGRfbWFya2VycyhzaG93bGVnZW5kID0gRikgJT4lIA0KICBhZGRfbGluZXMoeSA9IH5maXR0ZWQoY29sbGVnZV9tb2RlbCkpICU+JSANCiAgbGF5b3V0KHRpdGxlID0gIlJlbGF0aW9uc2hpcCBiZXR3ZWVuIGdvb2dsZSBzZWFyY2hlcyBmb3IgJ0NvbGxlZ2UgQWxjb2hvbCcgYW5kIEJpbmdlIERyaW5raW5nIFJhdGVzLCBieSBTdGF0ZSIsDQogICB4YXhpcyA9IGxpc3QodGl0bGUgPSAiR29vZ2xlIHNlYXJjaCB2b2x1bWUgZm9yICdDb2xsZWdlIEFsY29ob2wnIiksDQogICB5YXhpcyA9IGxpc3QodGl0bGUgPSAiU3RhdGUgQmluZ2UgRHJpbmtpbmcgUmF0ZSwgcGVyIGNhcGl0YSIpKQ0KDQoNCmBgYA0KDQoNCk5vdyB3ZSB3aWxsIGxvb2sgYXQgR29vZ2xlIHNlYXJjaCBoaXRzIGZvciBmcmF0IHBhcnRpZXMgYXMgY29tcGFyZWQgdG8gdGhlIGJpbmdlIGRyaW5raW5nIHBlcmNlbnRhZ2UgdG8gc2VlIGlmIHRoZXJlIGlzIGEgZ29vZCBjb3JyZWxhdGlvbiB0aGVyZSBhcyB3ZWxsLg0KDQoNCmBgYHtyfQ0Kc3RhdGVzICU+JSANCiAgcmVuYW1lKGxvY2F0aW9uID0gTkFNRSkgJT4lIA0KICBpbm5lcl9qb2luKGZyYXRzX3N0YXRlcykNCmBgYA0KDQpgYGB7cn0NCmdsaW1wc2UoQmluZ2VfZHJpbmtpbmdfcGVyY2VudCkNCmdsaW1wc2UoZnJhdHNfc3RhdGVzKQ0KYGBgDQoNCmBgYHtyfQ0KZnJhdHNfZGF0YSA8LSBmcmF0c19zdGF0ZXMgJT4lIA0KICBtdXRhdGUoU3RhdGUgPSBzdGF0ZTJhYmJyKGxvY2F0aW9uKSkgJT4lIA0KICBpbm5lcl9qb2luKEJpbmdlX2RyaW5raW5nX3BlcmNlbnQpDQogDQoNCmBgYA0KDQoNCmBgYHtyfQ0KDQpmcmF0c19jb2xvcnMgPC0gY29sb3JOdW1lcmljKHBhbGV0dGUgPSAidmlyaWRpcyIsIGRvbWFpbiA9IGZyYXRzX2RhdGEkaGl0cykNCg0KdW5kZXJncmFkX2xlYWZsZXQgJT4lIA0KICByZW5hbWUobG9jYXRpb24gPSBOQU1FKSAlPiUgDQogIGlubmVyX2pvaW4oZnJhdHNfZGF0YSkgJT4lIA0KICBsZWFmbGV0KCkgJT4lIA0KICBhZGRUaWxlcygpICU+JQ0KICBhZGRQb2x5Z29ucyh3ZWlnaHQgPSAxLA0KICAgICAgICAgICAgICBmaWxsQ29sb3IgPSB+ZnJhdHNfY29sb3JzKGhpdHMpLA0KICAgICAgICAgICAgICBsYWJlbCA9IH5wYXN0ZTAobG9jYXRpb24sICIsIEZyYXQgUGFydGllcyBIaXRzIGZyb20gR29vZ2xlIFNlYXJjaCA9ICIsIGhpdHMpLA0KICAgICAgICAgICAgICBoaWdobGlnaHQgPSBoaWdobGlnaHRPcHRpb25zKHdlaWdodCA9IDIpKSAlPiUgDQogIHNldFZpZXcoLTk1LCA0MCwgem9vbSA9IDQpICU+JSANCiAgYWRkTGVnZW5kKHBhbCA9IGZyYXRzX2NvbG9ycywgdmFsdWVzID0gfmhpdHMpDQpgYGANCg0KVGhlIGFib3ZlIG1hcCBzaG93cyB0aGF0IHRoZXJlIHdlcmUgbm90IG1hbnkgaGl0cyBmb3IgJ2ZyYXQgcGFydGllcycgb24gR29vZ2xlIFNlYXJjaC4NCg0KYGBge3J9DQpmcmF0c19tb2RlbCA8LSBsbShQZXJjZW50YWdlIH4gaGl0cywgZGF0YSA9IGZyYXRzX2RhdGEpDQoNCmdsYW5jZShmcmF0c19tb2RlbCkNCnRpZHkoZnJhdHNfbW9kZWwpDQpgYGANCg0KYGBge3J9DQpzcXJ0KC4wNDA2KQ0KDQpgYGANCg0KDQoNCmBgYHtyfQ0KZnJhdHNfZGF0YSAlPiUgDQogIGRyb3BfbmEoKSAlPiUgDQogIHBsb3RfbHkoeCA9IH5oaXRzLCANCiAgICAgICAgICB5ID0gflBlcmNlbnRhZ2UsDQogICAgICAgICAgaG92ZXJpbmZvID0gInRleHQiLCANCiAgICAgICAgICB0ZXh0ID0gfnBhc3RlKCJTdGF0ZTogIiwgbG9jYXRpb24sICI8YnI+IiwgIidGcmF0IFBhcnRpZXMnIHNlYXJjaCByYXRlOiAiLCBoaXRzLCAiPGJyPiIsICJQZXJjZW50YWdlOiAiLCBQZXJjZW50YWdlKSkgJT4lIA0KICBhZGRfbWFya2VycyhzaG93bGVnZW5kID0gRikgJT4lIA0KICBhZGRfbGluZXMoeSA9IH5maXR0ZWQoZnJhdHNfbW9kZWwpKSAlPiUgDQogIGxheW91dCh0aXRsZSA9ICJSZWxhdGlvbnNoaXAgYmV0d2VlbiBnb29nbGUgc2VhcmNoZXMgZm9yICdGcmF0IFBhcnRpZXMnIGFuZCBCaW5nZSBEcmlua2luZyBSYXRlcywgYnkgU3RhdGUiLA0KICAgeGF4aXMgPSBsaXN0KHRpdGxlID0gIkdvb2dsZSBzZWFyY2ggdm9sdW1lIGZvciAnRnJhdCBQYXJ0aWVzJyIpLA0KICAgeWF4aXMgPSBsaXN0KHRpdGxlID0gIlN0YXRlIEJpbmdlIERyaW5raW5nIFJhdGUsIHBlciBjYXBpdGEiKSkNCg0KDQpgYGANCg0KDQpIb3dldmVyLCB0aGUgcC12YWx1ZXMgZm9yICdmcmF0IHBhcnRpZXMnIHNob3cgYSBzdGF0aXN0aWNhbGx5IHNpZ25pZmljYW50IHJlbGF0aW9uc2hpcCB0byBiaW5nZSBkcmlua2luZyBwZXJjZW50YWdlIGJ5IHN0YXRlIGFzIHdlbGwsIGFzIHNob3duIGluIHRoZSBsaW5lIGdyYXBoIGFib3ZlLg0KDQoNCg0K