我有一组经度/纬度值,其中第三个变量可以为-1或1。(示例数据:https://pastebin.com/WPq90ka1)
我正在使用stat_summary_hex将这些数据映射到十六进制垃圾箱中:
single <- read.csv(file="mysamplefile.csv",head=TRUE)
pdxMap <- get_stamenmap(bbox = c(left = -122.844435, bottom = 45.42, right =
-122.43, top = 45.665316), zoom = 11, maptype = "toner")
ggmap(pdxMap) +
coord_equal() + theme(aspect.ratio = 1) +
stat_summary_hex(data=single, aes(x=LONGITUDE,y=LATITUDE,z = Value),
fun = function(z) mean(z),
alpha = .9, bins=50) +
scale_fill_gradientn(colors=c('#4393c3', '#ffffff', '#d6604d'))
但是,我宁愿不计算均值(z),而是计算值为1的百分比,即:(值之和== 1)/(所有abs(值)的总和)
分母不是问题,但我不知道如何管理分子。任何指导将不胜感激!