-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathPandemic.R
144 lines (124 loc) · 3.83 KB
/
Pandemic.R
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
# PDU Projekt nr.3
# 30 maj 2022
# Pandemic Data
library(pander)
library(data.table)
library(tidyverse)
library(ggmap)
library(ggplot2)
april20 <- fread("data/NC/202004-citibike-tripdata.csv")
may20 <- fread("data/NC/202005-citibike-tripdata.csv")
june20 <- fread("data/NC/202006-citibike-tripdata.csv")
april19 <- fread("data/NC/201904-citibike-tripdata.csv")
may19 <- fread("data/NC/201905-citibike-tripdata.csv")
summary(june20)
# Data from 2019
april19_data <- april19 %>%
group_by(`start station id`) %>%
summarise(
start_name = `start station name`[1],
lat = as.numeric(`start station latitude`[1]),
long = as.numeric(`start station longitude`[1]),
n_trips = n() / 31
) %>%
rename(start_id = `start station id`)
may19_data <- may19 %>%
group_by(`start station id`) %>%
summarise(
start_name = `start station name`[1],
lat = as.numeric(`start station latitude`[1]),
long = as.numeric(`start station longitude`[1]),
n_trips = n() / 31
) %>%
rename(start_id = `start station id`)
# Data from 2020
april20_data <- april20 %>%
group_by(`start station id`) %>%
summarise(
start_name = `start station name`[1],
lat = as.numeric(`start station latitude`[1]),
long = as.numeric(`start station longitude`[1]),
n_trips = n() / 31
) %>%
rename(start_id = `start station id`)
may20_data <- may20 %>%
group_by(`start station id`) %>%
summarise(
start_name = `start station name`[1],
lat = as.numeric(`start station latitude`[1]),
long = as.numeric(`start station longitude`[1]),
n_trips = n() / 31
) %>%
rename(start_id = `start station id`)
june20_data <- june20 %>%
group_by(`start station id`) %>%
summarise(
start_name = `start station name`[1],
lat = as.numeric(`start station latitude`[1]),
long = as.numeric(`start station longitude`[1]),
n_trips = n() / 31
) %>%
rename(start_id = `start station id`)
# New York map
nyc_map <- get_map(
location = "Noho, New York",
color = "color",
maptype = "terrain",
source = "google",
zoom = 12
)
# Map for April 2019
ggmap(nyc_map) +
stat_summary_2d(
data = april19_data,
aes(x = long, y = lat, z = n_trips),
fun = mean,
alpha = 0.6,
bins = 30) +
scale_fill_gradient(name = "Number of trips \n per day", low = "green", high = "red", limits = c(0, 300)) +
ggtitle("Average number of trip per day in April 2019")
ggsave("maps/April2019.png")
# Map for May 2019
ggmap(nyc_map) +
stat_summary_2d(
data = may19_data,
aes(x = long, y = lat, z = n_trips),
fun = mean,
alpha = 0.6,
bins = 30) +
scale_fill_gradient(name = "Number of trips \n per day", low = "green", high = "red", limits = c(0, 300)) +
ggtitle("Average number of trip per day in May 2019")
ggsave("maps/Map2019.png")
# Map for April 2020
ggmap(nyc_map) +
stat_summary_2d(
data = april20_data,
aes(x = long, y = lat, z = n_trips),
fun = mean,
alpha = 0.6,
bins = 30) +
scale_fill_gradient(name = "Number of trips \n per day", low = "green", high = "red", limits = c(0, 300)) +
ggtitle("Average number of trip per day in April 2020")
ggsave("maps/April2020.png")
# Map for May 2020
ggmap(nyc_map) +
stat_summary_2d(
data = may20_data,
aes(x = long, y = lat, z = n_trips),
fun = mean,
alpha = 0.6,
bins = 30) +
scale_fill_gradient(name = "Number of trips \n per day", low = "green", high = "red", limits = c(0, 300)) +
ggtitle("Average number of trip per day in May 2020")
ggsave("maps/May2020.png")
# Map for June 2020
ggmap(nyc_map) +
stat_summary_2d(
data = june20_data,
aes(x = long, y = lat, z = n_trips),
fun = mean,
alpha = 0.6,
bins = 30) +
scale_fill_gradient(name = "Number of trips \n per day", low = "green", high = "red", limits = c(0, 300)) +
ggtitle("Average number of trip per day in June 2020")
ggsave("maps/June2020.png")