在本文中，我将创建纽约市机动车碰撞的市镇级热图。
数据来自纽约市数据。特别是，我将从镇级碰撞到街道级碰撞。
可下载资源
我使用下面的三个函数来处理我的数据。boro（）函数子集用于与指定自治市中的街道名称。accident_freq（）函数计算每条街道的碰撞频率，然后将这些数字合并。Assign_col（）函数获取特定市镇的碰撞级别数据集（使用accident_freq（）函数创建），并为每条街道分配从白色到指定颜色（例如绿色，红色等）的颜色。碰撞更多的街道将更暗。
man_col=assign_col(man_freq,'dodgerblue')
bronx_col=assign_col(bronx_freq,'darkred')
brook_col=assign_col(brook_freq,'violet')
si_col=assign_col(si_freq,'darkgreen')
q_col=assign_col(q_freq,'darkgoldenrod4')最后，我使用ggmap的get_map（）函数获取NYC的样式图并添加geom_path图层。每个市只有一个geom_path（）层。Geom_path（）使用直线或“路径”连接同一条街道上的所有经度和纬度点。组中的所有坐标都已连接。然后使用col =参数为每行赋予一个由assign_col（）确定的颜色。
ny_plot+
geom_path(data=man,size=1,aes(x=man$long, y=man$lat,group=man$ON.STREET.NAME),col=man_col[man_freq$freqPerc])+
geom_path(data=bronx,size=1,aes(x=bronx$long, y=bronx$lat,group=bronx$ON.STREET.NAME),col=bronx_col[bronx_freq$freqPerc])+
geom_path(data=brook,size=1,aes(x=brook$long, y=brook$lat,group=brook$ON.STREET.NAME),col=brook_col[brook_freq$freqPerc])+
geom_path(data=si,size=1,aes(x=si$long, y=si$lat,group=si$ON.STREET.NAME),col=si_col[si_freq$freqPerc])+
geom_path(data=q,size=1,aes(x=q$long, y=q$lat,group=q$ON.STREET.NAME),col=q_col[q_freq$freqPerc])+
ggtitle('Street-Level NYC Vehicle Accidents by Borough')+
xlab(" ")+ylab(" ")
在下面，我加载ggmap包和数据。
library(ggmap)
comm=regexpr(',',d_clean$LOCATION)
# create year variable
d_clean$year=substr(d_clean$DATE,7,10)
【视频】因子分析简介及R语言应用实例：对地区经济研究分析重庆市经济指标
R语言宏观经济学：IS-LM曲线可视化货币市场均衡
R语言代做编程辅导和解答GLM Coursework
R语言分析上海空气质量数据：kmean聚类、层次聚类、时间序列分析：arima模型、指数平滑法
