Как использовать пакет gganimate построить динамическую карту в R? - PullRequest
Новая
       14

Как использовать пакет gganimate построить динамическую карту в R?

2 голосов
/ 07 мая 2019

Как я знаю, gganimate был в версии 1.0.3, мы можем использовать функцию transition_* для построения динамического графика. Но когда я запускаю следующий код, возникает ошибка: 

Error in `$<-.data.frame`(`*tmp*`, "group", value = "") : 
  replacement has 1 row, data has 0

код: 

library(ggmap)
library(gganimate)
world <- map_data("world")
world <- world[world$region!="Antarctica",]
data <- data.frame(state = c("Alabama","Alaska","Alberta","Alberta","Arizona"),
                   lon = c(-86.55,-149.52,-114.05,-113.25,-112.05),
                   lat = c(33.30,61.13,51.05,53.34,33.30)
                   )
ggplot()+
  geom_map(data = world,
           map = world,
           aes(long,lat,map_id = region),
           color = '#333300',
           fill = '#663300') +
  geom_point(data = data,
             aes(x = lon, y = lat),
             size = 2.5) +
  geom_jitter(width = 0.1) +
  transition_states(states = state)

1 Ответ

1 голос
/ 08 мая 2019

У вас нет данных, определенных в строке ggplot() верхнего уровня, поэтому state в transition_* происходит из ниоткуда.

Мне также непонятно, почему у вас есть уровень geom_jitter вваш код.Как и transition_*, он не имеет наследуемых данных / эстетических отображений верхнего уровня, поэтому он также выдал бы ошибку, если бы transition_* сначала не вызвал ошибку.Кроме того, даже если мы добавим сопоставления, учитывая диапазон координат широты / долготы в ваших данных, дрожание на 0,1 вряд ли будет иметь значение визуально.

Вы можете попробовать следующее: 

# put data in top level ggplot()
ggplot(data,
       aes(x = lon, y = lat))+
  geom_map(data = world,
           map = world,
           aes(long,lat,map_id = region),
           color = '#333300', fill = '#663300',
           # lighter background for better visibility
           alpha = 0.5) + 
  geom_point(size = 2.5) +
  # limit coordinates to relevant range
  coord_quickmap(x = c(-180, -50), y = c(25, 85)) +
  transition_states(states = state)

plot

...