Блог

Я использую ggmap для построения тепловой карты продолжительности поездок на такси в Нью-Йорке. Я хотел бы использовать продолжительность поездки в качестве своей цветовой шкалы. У меня есть большое количество наблюдений с широтой / длиной захвата, которые я могу легко нанести на карту Нью-Йорка. Продолжительность отключения указана в секундах и доступна для каждого наблюдения.

Тем не менее, я пытаюсь найти наилучший способ использовать что-то вроде stat_density2d построения контуров вокруг определенных областей, которые имеют более высокую продолжительность поездки (возможно, можно использовать среднее значение в определенной области?). Я также группирую свои графики по месяцам.

Как я могу это сделать? Мой приведенный ниже код отображает это, но я считаю, что мой масштаб выражен в терминах # observations, а не числовой переменной trip_duration .

 library(ggmap)
nyc_map <- get_map(location = c(lon = median(nyc_train2$pickup_longitude), lat = median(nyc_train2$pickup_latitude)), maptype = "terrain", zoom = 13)

pickup_trip_duration_plot <- ggmap(nyc_map)   
                  geom_density2d(data = nyc_train2, aes(x = pickup_longitude, y = pickup_latitude, fill = trip_duration), size = 0.3)   
                  stat_density2d(data = nyc_train2, aes(x = pickup_longitude, y = pickup_latitude, fill = ..level.., alpha = ..level..), size = 0.01, bins = 16, geom = "polygon")   
                  facet_wrap(~ pickup_month_name, ncol = 3)  
                  scale_fill_gradient(low = "yellow", high = "red")   
                  ggtitle("Trip Duration Heatmap by Month for Pickup Locations")  
                  scale_alpha(range = c(0, 0.3), guide = FALSE)
pickup_trip_duration_plot

# I believe this may have something to do with why I am having plotting issues
# with a numeric value?
Warning message:
Ignoring unknown aesthetics: fill