У меня есть набор значений lat / long с третьей переменной, которая может быть либо -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) / (сумма всех абс (значений)))
Знаменатель не проблема, но я не могу понять, как управлять числителем. Любое руководство будет высоко ценится!