-
Notifications
You must be signed in to change notification settings - Fork 0
/
2021_05_18_tidy_tuesday.Rmd
157 lines (124 loc) · 3.69 KB
/
2021_05_18_tidy_tuesday.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
---
title: "TidyTemplate"
date: 2021-05-20
output: html_output
---
# TidyTuesday
Join the R4DS Online Learning Community in the weekly #TidyTuesday event!
Every week we post a raw dataset, a chart or article related to that dataset, and ask you to explore the data.
While the dataset will be “tamed”, it will not always be tidy! As such you might need to apply various R for Data Science techniques to wrangle the data into a true tidy format.
The goal of TidyTuesday is to apply your R skills, get feedback, explore other’s work, and connect with the greater #RStats community!
As such we encourage everyone of all skills to participate!
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
library(tidyverse)
library(tidytuesdayR)
library(scales)
library(countrycode)
library(maps)
library(ggthemes)
```
# Load the weekly Data
Dowload the weekly data and make available in the `tt` object.
```{r Load}
tt <- tt_load("2021-05-04")
tt
```
tidying up
```{r date}
library(lubridate)
water <- tt$water %>%
mutate(report_date = mdy(report_date)) %>%
rename(lat = lat_deg,
lon = lon_deg,
country = country_name) %>%
separate(water_tech, c ("water_tech", "brand"), sep = " - ",
fill = "right") %>%
mutate(install_year = ifelse(install_year > 2021, NA_real_, install_year)) %>%
filter(!country %in% c("Peru", "Dominican Republic", "Timor-Leste"),
!is.na(country)) %>%
filter(between(lat, -30, 20),
between(lon, -40, 60))
water
```
taking a look at the data
```{r}
water %>%
count(status_id)
water %>%
count(water_source, sort = TRUE)
water %>%
count(water_tech, sort = TRUE)
water %>%
count(water_source, water_tech, sort = TRUE)
water %>%
filter(install_year > 2021)
water %>%
filter(install_year > 1990) %>%
count(install_year, sort = TRUE) %>%
ggplot(aes(install_year, n)) +
geom_col()
water %>%
count(installer, sort = TRUE) %>%
view()
water %>%
count(pay, sort = TRUE) %>%
view()
```
"pay" is not clean a free text column
```{r}
water %>%
count(status_id, status, sort = TRUE) %>%
view()
```
### Making a map
```{r maps}
countries <- unique(water$country)
africa_map_data <- map_data("world") %>%
as_tibble() %>%
mutate(continent = countrycode(region, "country.name", "continent")) %>%
filter(continent == "Africa")
water %>%
#sample_n(200000) %>%
filter(country == "Nigeria",
between(lat, -5, 15),
between(lon, 0, 15)) %>%
ggplot(aes(lon, lat)) +
geom_polygon(aes(long, lat, group = group),
color = "dark gray",
fill = NA,
data = africa_map_data) +
geom_point(size = .3, alpha = .01) +
theme(legend.position = "none") +
labs(x = NULL,
y = NULL,
title = "Water Sources in Central Africa",
subtitle = "2021-05-04",
caption = "From the Water Point Data Exchange (WPDx-Basic)")
```
WHere are water points in Nigeria
```{r Nigeria Map}
water_nigeria <- water %>%
filter(country == "Nigeria",
between(lat, -15, 15),
between(lon, -10, 15)) %>%
ggplot() +
geom_point(aes(lat, lon, color = status_id),
size = .1, alpha = .05) +
borders("Nigeria") +
coord_fixed(ratio = 1) +
theme_map() +
scale_color_discrete(guide = guide_legend(override.aes = list(size = 2, alpha = 1))) +
labs(x = NULL,
y = NULL,
title = "Water Sources in Nigeria",
subtitle = "2021-05-04",
caption = "From the Water Point Data Exchange (WPDx-Basic)")
water_nigeria
?borders
?get_map
```
```{r}
??get_map
get_map("Nigeria", source = "osm")
```