Dear all,
I'm trying to plot a density on a worldmap using ggplot2. I have three
variables: Freq, Latitude and Longitude. However most of my points are around
the same coordinates, which is why my density is displayed red there. Is there
a way that the displayed density dependts on the variable Freq instead of the
amount of points in the area?
My current code:
rel = as.data.frame(matrix(0, nrow = 6, ncol = 3))
rel[,1] = c(1,0,0,0.73,0.94,1)
rel[,2] = c(34.55040,34.16807,36.02210,36.09917,36.35868,32.12694)
rel[,3] = c(40.88923, 45.10331, 43.40581, 43.32750, 43.15194, 45.23083)
colnames(rel) = c("Freq", "Latitude", "Longitude")
map = get_map(location = c(lon=43, lat=36), zoom = 6)
ggmap(map) +
geom_density2d(aes(x=Longitude, y=Latitude),
data = rel ,size = 0.3) +
stat_density2d(data = rel, aes(x=Longitude, y=Latitude,
fill = ..level..,
alpha = ..level..), size =0.01, bins = 16,
geom = "polygon") +
scale_fill_gradient(low = "green", high = "red") +
scale_alpha(range = c(0, 0.3), guide = FALSE) +
geom_point(aes(x=Longitude, y=Latitude, colour = Freq),alpha = 0.75, size=5,
data = rel)+
scale_colour_gradient(low = "#3afff4", high = "#fc0000")
I would like the density to be displayed in red everywhere where red points are
present.
Thanks,
Nico
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
