Upload files to ''

This commit is contained in:
Paweł Lewicki 2023-05-30 22:08:48 +02:00
parent 9d98bb75f1
commit 0d28556d38
1 changed files with 89 additions and 0 deletions

89
projekt_2_housing_new.R Normal file
View File

@ -0,0 +1,89 @@
library(tidyverse)
library(eurostat)
# install.packages("ggthemes")
library(leaflet)
library(sf)
library(scales)
library(cowplot)
library(ggthemes)
library(RColorBrewer)
# Load results from project_1
map_df = read.csv(".//data//compound_interest_housing.csv")
map_df
# Load map/polygon data from eurostat
SHP_0 <- get_eurostat_geospatial(resolution = 10,
nuts_level = 0,
year = 2016)
SHP_0 %>%
ggplot() +
geom_sf()
EU28 <- eu_countries %>%
select(geo = code, name)
SHP_28 <- SHP_0 %>%
select(geo = NUTS_ID, geometry) %>%
inner_join(EU28, by = "geo") %>%
arrange(geo) %>%
st_as_sf()
SHP_28 %>%
ggplot() +
geom_sf() +
scale_x_continuous(limits = c(-10, 35)) +
scale_y_continuous(limits = c(35, 65))
# Join datasets
mapdata_new <- left_join(SHP_28, map_df, by="geo", relationship = "many-to-many")
# Delete Greece (null values)
mapdata_new <- mapdata_new[-9,]
mapdata_new$wiki <- paste0( "https://en.wikipedia.org/wiki/", mapdata_new$name )
# Create a continuous palette function
qpal <- colorQuantile("YlOrBr", mapdata_new$substr_house_prices_wo_hicp, n = 5)
qpal_colors <- unique(qpal(sort(mapdata_new$substr_house_prices_wo_hicp))) # hex codes
qpal_labs <- quantile(mapdata_new$substr_house_prices_wo_hicp, seq(0, 1, .2)) # depends on n from pal
qpal_labs <- paste(lag(qpal_labs), qpal_labs, sep = " - ")[-1] # first lag is NA
popup_content <- paste0( "<div style='text-align: center'>"
, "_______________________________<br>"
, "Country: ", mapdata_new$name, "<br><br>"
, "</div>"
, "House prices index: ", mapdata_new$substr_house_prices_wo_hicp, "<br>"
, "Compound interest equivalent: ", mapdata_new$compound_interest, "%<br>"
, "<div style='text-align: center'>"
, "<a href='", mapdata_new$wiki, "' target='_blank'>", "<br>"
, "Click here to view wiki</a>", "<br>"
, "</div>" )
map <- leaflet() %>%
addProviderTiles(providers$CartoDB.Positron) %>%
addPolygons(data=mapdata_new,
fillOpacity = 0.6, # Przezroczystość
stroke = TRUE, # Borders visible
color = "grey", # Border color
weight = 1,
fillColor = ~qpal(mapdata_new$substr_house_prices_wo_hicp),
popup = popup_content,
popupOptions = popupOptions(maxWidth ="100%", closeOnClick = TRUE),
) %>%
setView( lat = 51, lng = 14, zoom = 4) %>%
addLegend("bottomright", colors = qpal_colors,
title = "<span style='white-space: pre-line;'> Real house prices \n index (2022) </span>",
labels = qpal_labs,
opacity = 1)
map