-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy path.Rhistory
512 lines (512 loc) · 29.8 KB
/
.Rhistory
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
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
#pm_ct <- pm_ct + theme(legend.position="top", legend.key = element_blank())
pm_ct <- pm_ct + coord_fixed()
pm_ct <- pm_ct + theme_nothing(legend=TRUE)
pm_ct <- pm_ct + labs(x=NULL, y=NULL, title="Where tow yards target")
pm_ct
pm_ct <- ggmap(hart_map)
pm_ct <- pm_ct + stat_density2d(data = tows_no_na, show.legend=F, aes(x=tow_lon, y=tow_lat, fill=..level.., alpha=..level..), geom="polygon",size=.5,bins=10)
pm_ct <- pm_ct + geom_polygon(data = joined_tracts2, aes(x=long, y=lat, group=group), fill=NA, color = "black", size=0.2)
pm_ct <- pm_ct + geom_polygon(data = town_borders, aes(x=long, y=lat, group=group), fill=NA, color = "black", size=0.4)
#pm_ct <- pm_ct + geom_polygon(data = ct_only, aes(x=long, y=lat, group=group), fill="seagreen2", color = "gray93", size=0.2)
#pm_ct <- pm_ct + gg_circle(r=9, xc=-73, yc=42, color="white", fill=NA, alpha=0.2, size=40)
pm_ct <- pm_ct + scale_fill_gradient(low="deepskyblue2", high="firebrick1", name="Distribution")
#pm_ct <- pm_ct + scale_fill_discrete()
#extra_lat <- c(46.358685, 35.872715)
#extra_lon <- c(-64.209938, -79.735653)
#pm_ct <- pm_ct + theme(legend.position="top", legend.key = element_blank())
pm_ct <- pm_ct + coord_fixed()
pm_ct <- pm_ct + theme_nothing(legend=TRUE)
pm_ct <- pm_ct + labs(x=NULL, y=NULL, title="Where tow yards target")
pm_ct <- pm_ct + facet_wrap(~Tow_Firm)
pm_ct
library(dplyr)
library(lubridate)
library(leaflet)
library(ggmap)
library(knitr)
library(stringr)
library(geosphere)
library(ggplot2)
hartbox <- make_bbox(lon = tows_no_na$tow_lon, lat =tows_no_na$tow_lat, f = .1)
hart_map <- get_map(location = hartbox, maptype = "roadmap", source = "google")
pm_ct <- ggmap(hart_map)
pm_ct <- pm_ct + stat_density2d(data = tows_no_na, show.legend=F, aes(x=tow_lon, y=tow_lat, fill=..level.., alpha=..level..), geom="polygon",size=.5,bins=10)
#pm_ct <- pm_ct + geom_polygon(data = joined_tracts2, aes(x=long, y=lat, group=group), fill=NA, color = "black", size=0.2)
#pm_ct <- pm_ct + geom_polygon(data = town_borders, aes(x=long, y=lat, group=group), fill=NA, color = "black", size=0.4)
#pm_ct <- pm_ct + geom_polygon(data = ct_only, aes(x=long, y=lat, group=group), fill="seagreen2", color = "gray93", size=0.2)
#pm_ct <- pm_ct + gg_circle(r=9, xc=-73, yc=42, color="white", fill=NA, alpha=0.2, size=40)
pm_ct <- pm_ct + scale_fill_gradient(low="purple", high="firebrick1", name="Distribution")
#pm_ct <- pm_ct + scale_fill_discrete()
#extra_lat <- c(46.358685, 35.872715)
#extra_lon <- c(-64.209938, -79.735653)
#pm_ct <- pm_ct + theme(legend.position="top", legend.key = element_blank())
pm_ct <- pm_ct + coord_fixed()
pm_ct <- pm_ct + theme_nothing(legend=TRUE)
pm_ct <- pm_ct + labs(x=NULL, y=NULL, title="Where tow yards target")
#pm_ct <- pm_ct + facet_wrap(~Tow_Firm)
#pm_ct <- pm_ct + theme(text = element_text(size=15), panel.background = element_rect(fill = 'gray93', color=NA))
pm_ct <- pm_ct + theme(plot.title=element_text(face="bold", hjust=.4))
pm_ct <- pm_ct + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
pm_ct <- pm_ct + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
#pm_ct <- pm_ct + theme(legend.key.size = unit(1, "cm"))
pm_ct
tow_hist <- ggplot(tow_hour2, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=paste("Tows by time"),
caption="SOURCE: City of Hartford\nAndrew Ba Tran/TrendCT.org")
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
tow_hist
address_list <- unique(tow_hour2$Tow_From_Address)
for (i in 1:count(address_list)) {
print(i)
}
for (i in 1:length(address_list)) {
print(i)
}
i <- 1
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(tow_hour2, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=paste("Tows by time")
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 5, height = 4, type = "cairo-png")
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(tow_hour2, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 5, height = 4, type = "cairo-png")
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(tow_hour2, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 4, height = 2, type = "cairo-png")
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(histo, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 4, height = 2, type = "cairo-png")
for (i in 1:length(address_list)) {
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(histo, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 5, height = 2.5, type = "cairo-png")
}
tow_hour <- tows %>%
group_by(Tow_From_Address, hour) %>%
summarise(count=n()) %>%
filter(count > 4)
arrange(-count)
tow_hour2 <- spread(tow_hour, hour, count)
tow_hour2 <- gather(tow_hour2, hour, count, 2:25)
# Heatmap by hour
# Histograms of address by hour
address_list <- unique(tow_hour2$Tow_From_Address)
library(dplyr)
library(tidyr)
tow_hour <- tows %>%
group_by(Tow_From_Address, hour) %>%
summarise(count=n()) %>%
filter(count > 4)
arrange(-count)
tow_hour2 <- spread(tow_hour, hour, count)
tow_hour2 <- gather(tow_hour2, hour, count, 2:25)
# Heatmap by hour
# Histograms of address by hour
address_list <- unique(tow_hour2$Tow_From_Address)
for (i in 1:length(address_list)) {
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(histo, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 5, height = 2.5, type = "cairo-png")
}
i <- 1
library(dplyr)
library(lubridate)
library(leaflet)
library(ggmap)
library(knitr)
library(stringr)
library(geosphere)
library(ggplot2)
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(histo, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 5, height = 2.5, type = "cairo-png")
tow_hour <- tows %>%
group_by(Tow_From_Address, hour) %>%
summarise(count=n()) %>%
filter(count > 4)
arrange(-count)
tow_hour2 <- spread(tow_hour, hour, count)
tow_hour2 <- gather(tow_hour2, hour, count, 2:25)
# Heatmap by hour
tows$hour <- hour(tows$Time)
tow_hour <- tows %>%
group_by(Tow_From_Address, hour) %>%
summarise(count=n()) %>%
filter(count > 4) %>%
arrange(-count)
tow_hour2 <- spread(tow_hour, hour, count)
tow_hour2 <- gather(tow_hour2, hour, count, 2:25)
# Heatmap by hour
# Histograms of address by hour
address_list <- unique(tow_hour2$Tow_From_Address)
i <- 1
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(histo, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 5, height = 2.5, type = "cairo-png")
str(tow_hour2)
tow_hour2$hour <- as.numeric(tow_hour2$hour)
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(histo, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 5, height = 2.5, type = "cairo-png")
for (i in 1:length(address_list)) {
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(histo, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 5, height = 2.5, type = "cairo-png")
}
address2$png <- gsub(" ", "", address2$Tow_From_Address)
address2$png <- paste0("<div style='width:200px'><img src='http://projects.ctmirror.org/content/trend/2016/08/towed/hours/", address2$png, ".png</img></div>")
address2$pop <- ifelse(count > 4, paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count), paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count, address2$png))
leaflet(address2) %>% addTiles('http://a.tiles.mapbox.com/v3/borzechowski.gcj2gonc/{z}/{x}/{y}.png', attribution='<a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a>') %>%
setView(-72.690940, 41.751426, zoom = 13) %>%
addCircles(~tow_lon, ~tow_lat, popup=address2$pop, weight = 3, radius=address2$count*1.5,
color="#ffa500", stroke = TRUE, fillOpacity = 0.2) %>%
addLegend("bottomright", colors= "#ffa500", labels="Towed'", title="In Hartford")
address2$png <- gsub(" ", "", address2$Tow_From_Address)
address2$png <- paste0("<div style='width:200px'><img src='http://projects.ctmirror.org/content/trend/2016/08/towed/hours/", address2$png, ".png</img></div>")
address2$pop <- ifelse(address2$count > 4, paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count), paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count, address2$png))
leaflet(address2) %>% addTiles('http://a.tiles.mapbox.com/v3/borzechowski.gcj2gonc/{z}/{x}/{y}.png', attribution='<a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a>') %>%
setView(-72.690940, 41.751426, zoom = 13) %>%
addCircles(~tow_lon, ~tow_lat, popup=address2$pop, weight = 3, radius=address2$count*1.5,
color="#ffa500", stroke = TRUE, fillOpacity = 0.2) %>%
addLegend("bottomright", colors= "#ffa500", labels="Towed'", title="In Hartford")
View(address2)
address2$png <- gsub(" ", "", address2$Tow_From_Address)
address2$png <- paste0("<div style='width:200px'><img src='http://projects.ctmirror.org/content/trend/2016/08/towed/hours/", address2$png, ".png'></img></div>")
address2$pop <- ifelse(address2$count > 4, paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count), paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count, address2$png))
leaflet(address2) %>% addTiles('http://a.tiles.mapbox.com/v3/borzechowski.gcj2gonc/{z}/{x}/{y}.png', attribution='<a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a>') %>%
setView(-72.690940, 41.751426, zoom = 13) %>%
addCircles(~tow_lon, ~tow_lat, popup=address2$pop, weight = 3, radius=address2$count*1.5,
color="#ffa500", stroke = TRUE, fillOpacity = 0.2) %>%
addLegend("bottomright", colors= "#ffa500", labels="Towed'", title="In Hartford")
View(address2)
address2$png <- gsub(" ", "", address2$Tow_From_Address)
address2$png <- paste0("<img src='http://projects.ctmirror.org/content/trend/2016/08/towed/hours/", address2$png, ".png' width='250px'></img>")
address2$pop <- ifelse(address2$count > 4, paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count), paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count, "<br />", address2$png))
leaflet(address2) %>% addTiles('http://a.tiles.mapbox.com/v3/borzechowski.gcj2gonc/{z}/{x}/{y}.png', attribution='<a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a>') %>%
setView(-72.690940, 41.751426, zoom = 13) %>%
addCircles(~tow_lon, ~tow_lat, popup=address2$pop, weight = 3, radius=address2$count*1.5,
color="#ffa500", stroke = TRUE, fillOpacity = 0.2) %>%
addLegend("bottomright", colors= "#ffa500", labels="Towed'", title="In Hartford")
View(address2)
address2$png <- gsub(" ", "", address2$Tow_From_Address)
address2$png <- paste0("<img src='http://projects.ctmirror.org/content/trend/2016/08/towed/hours/", address2$png, ".png' width='250px'></img>")
address2$pop <- ifelse(address2$count < 4, paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count), paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count, "<br />", address2$png))
leaflet(address2) %>% addTiles('http://a.tiles.mapbox.com/v3/borzechowski.gcj2gonc/{z}/{x}/{y}.png', attribution='<a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a>') %>%
setView(-72.690940, 41.751426, zoom = 13) %>%
addCircles(~tow_lon, ~tow_lat, popup=address2$pop, weight = 3, radius=address2$count*1.5,
color="#ffa500", stroke = TRUE, fillOpacity = 0.2) %>%
addLegend("bottomright", colors= "#ffa500", labels="Towed'", title="In Hartford")
address2$png <- gsub(" ", "", address2$Tow_From_Address)
address2$png <- paste0("<img src='http://projects.ctmirror.org/content/trend/2016/08/towed/hours/", address2$png, ".png' width='250px'></img>")
address2$pop <- ifelse(address2$count < 4, paste0("<strong>Tows: </strong>", address2$count), paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count, "<br />", address2$png))
leaflet(address2) %>% addTiles('http://a.tiles.mapbox.com/v3/borzechowski.gcj2gonc/{z}/{x}/{y}.png', attribution='<a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a>') %>%
setView(-72.690940, 41.751426, zoom = 13) %>%
addCircles(~tow_lon, ~tow_lat, popup=address2$pop, weight = 3, radius=address2$count*1.5,
color="#ffa500", stroke = TRUE, fillOpacity = 0.2) %>%
addLegend("bottomright", colors= "#ffa500", labels="Towed'", title="In Hartford")
address2$png <- gsub(" ", "", address2$Tow_From_Address)
address2$png <- paste0("<img src='http://projects.ctmirror.org/content/trend/2016/08/towed/hours/", address2$png, ".png' width='250px'></img>")
address2$pop <- ifelse(address2$count < 4, paste0("<strong>Tows: </strong>", address2$count), paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count, "<br />", address2$png))
leaflet(address2) %>% addTiles('http://a.tiles.mapbox.com/v3/borzechowski.gcj2gonc/{z}/{x}/{y}.png', attribution='<a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a>') %>%
setView(-72.690940, 41.751426, zoom = 13) %>%
addCircles(~tow_lon, ~tow_lat, popup=address2$pop, weight = 3, radius=address2$count*1.5,
color="#ffa500", stroke = TRUE, fillOpacity = 0.2) %>%
addLegend("bottomright", colors= "#ffa500", labels="Towed'", title="In Hartford")
address2$png <- gsub(" ", "", address2$Tow_From_Address)
address2$png <- paste0("<img src='http://projects.ctmirror.org/content/trend/2016/08/towed/hours/", address2$png, ".png' width='250px'></img>")
address2$pop <- ifelse(address2$count < 4, paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count), paste0("<strong>Tows: </strong>", address2$count, "<br />", address2$png))
leaflet(address2) %>% addTiles('http://a.tiles.mapbox.com/v3/borzechowski.gcj2gonc/{z}/{x}/{y}.png', attribution='<a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a>') %>%
setView(-72.690940, 41.751426, zoom = 13) %>%
addCircles(~tow_lon, ~tow_lat, popup=address2$pop, weight = 3, radius=address2$count*1.5,
color="#ffa500", stroke = TRUE, fillOpacity = 0.2) %>%
addLegend("bottomright", colors= "#ffa500", labels="Towed'", title="In Hartford")
View(address2)
tows$hour <- hour(tows$Time)
tow_hour <- tows %>%
group_by(Tow_From_Address, hour) %>%
summarise(count=n()) %>%
filter(count > 4) %>%
arrange(-count)
tow_hour2 <- spread(tow_hour, hour, count)
tow_hour2 <- gather(tow_hour2, hour, count, 2:25)
tow_hour2$hour <- as.numeric(tow_hour2$hour)
# Heatmap by hour
# Histograms of address by hour
address_list <- unique(tow_hour2$Tow_From_Address)
for (i in 1:length(address_list)) {
addy <- address_list[i]
addy_short <- gsub(" ", "", addy)
histo <- subset(tow_hour2, Tow_From_Address==addy)
tow_hist <- ggplot(histo, aes(x=hour, y=count))
tow_hist <- tow_hist + geom_bar(stat="identity")
#tow_hist <- tow_hist + facet_wrap(~Tow_From_Address)
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", hjust=.4))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(l=20)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=12), color="#7a7d7e", hjust=0))
tow_hist <- tow_hist + labs(x="Hour", y="Tows", title=addy)
tow_hist <- tow_hist + scale_x_continuous(limits=c(0,24),
breaks=0:12*2,
labels=c("12 am", paste(1:5*2,"am"),
"12 pm",
paste(7:11*2-12,"pm"),
"12 am"))
tow_hist <- tow_hist + theme_bw(base_family="Calibri")
tow_hist <- tow_hist + theme(panel.border=element_blank())
tow_hist <- tow_hist + theme(text = element_text(size=10))
tow_hist <- tow_hist + theme(plot.title=element_text(face="bold", family="Lato Black", size=22))
tow_hist <- tow_hist + theme(plot.subtitle=element_text(face="italic", size=9, margin=margin(b=12)))
tow_hist <- tow_hist + theme(plot.caption=element_text(size=12, margin=margin(t=10, r=80), color="#7a7d7e"))
tow_hist <- tow_hist + theme(legend.position="none")
file_name <- paste0("hours/",addy_short, ".png")
ggsave(tow_hist, file = file_name, width = 5, height = 2.5, type = "cairo-png")
}
address2$png <- gsub(" ", "", address2$Tow_From_Address)
address2$png <- paste0("<img src='http://projects.ctmirror.org/content/trend/2016/08/towed/hours/", address2$png, ".png' width='250px'></img>")
address2$pop <- ifelse(address2$count < 5, paste0(address2$Tow_From_Address, "<br /><strong>Tows: </strong>", address2$count), paste0("<strong>Tows: </strong>", address2$count, "<br />", address2$png))
leaflet(address2) %>% addTiles('http://a.tiles.mapbox.com/v3/borzechowski.gcj2gonc/{z}/{x}/{y}.png', attribution='<a href="http://www.openstreetmap.org/copyright">OpenStreetMap</a>') %>%
setView(-72.690940, 41.751426, zoom = 13) %>%
addCircles(~tow_lon, ~tow_lat, popup=address2$pop, weight = 3, radius=address2$count*1.5,
color="#ffa500", stroke = TRUE, fillOpacity = 0.2) %>%
addLegend("bottomright", colors= "#ffa500", labels="Towed'", title="In Hartford")
main <- subset(tows, Tow_From_Address=="1000 MAIN ST")
View(main)
main$month <- month(main$Date)
table(main$month)
View(address2)
write.csv(address2, "address2.csv")