Блог

Это следует из моего последнего вопроса. Я потратил около часа, пытаясь понять, как передать переменную, которую я использую для фильтрации моего фрейма данных, в заголовок сгенерированного графика.

Следуя моим предыдущим вопросам.

 library (tidyverse)
library (epitools)


# here's my made up data

DISEASE = c("Marco Polio","Marco Polio","Marco Polio","Marco Polio","Marco Polio",
            "Mumps","Mumps","Mumps","Mumps","Mumps",
            "Chicky Pox","Chicky Pox","Chicky Pox","Chicky Pox","Chicky Pox")
YEAR = c(2011, 2012, 2013, 2014, 2015,
         2011, 2012, 2013, 2014, 2015,
         2011, 2012, 2013, 2014, 2015)
VALUE = c(82,89,79,51,51,
          79,91,69,89,78,
          71,69,95,61,87)
AREA =c("A", "B","C")

DATA = data.frame(DISEASE, YEAR, VALUE,AREA)

DATA<-
DATA %>%
  mutate(POPN = case_when(
    AREA == "A" ~ 2.5,
    AREA == "B" ~ 3,
    AREA == "C" ~ 7,
    TRUE ~ 0)) %>%
  group_by(DISEASE,AREA,POPN) %>%
  count(AREA) %>%
  mutate(res = list(pois.byar(n, POPN))) %>%
  unnest()

DATA%>%filter(DISEASE== "Marco Polio")%>%
  ggplot(aes(x=AREA, y=rate))  geom_point()  
  geom_hline(aes(yintercept=rate[AREA == "A"]), 
             linetype="dashed", color = "red")
  

Я думал, что это

     x_label = "Area!!!"
    y_label = "Rate!!!"
    DATA%>%filter(DISEASE== "Marco Polio")%>%
          ggplot(aes(x=AREA, y=rate))  geom_point()  
          geom_hline(aes(yintercept=rate[AREA == "A"]), 
                     linetype="dashed", color = "red") 
labs(x = x_label,y = y_label) 
ggtitle(DATA$DISEASE)
  

Почему этого не происходит?
Он генерирует диаграмму для Marco Polio, но использует Chicky Pox в качестве заголовка.

Чего я хочу, так это (ложный код) ggtitle == фильтр (болезнь)

Потому что то, что я собираюсь сделать после этого, — это ходить и мурлыкать, чтобы получить каждый график для каждой инфекции, и я хотел бы, чтобы заголовок был автоматически.

Ta.

РЕДАКТИРОВАТЬ: Я попробовал предложение ниже, и оно не совсем работает.

Я пробовал это

 DATA%>%filter(DISEASE== "Mumps")%>%
  ggplot(aes(x=AREA, y=rate))  geom_point()  
  geom_hline(aes(yintercept=rate[AREA == "A"]), 
             linetype="dashed", color = "red") 
  ggtitle(paste(DISEASE))


DATA%>%filter(DISEASE== "Mumps")%>%
  ggplot(aes(x=AREA, y=rate))  geom_point()  
  geom_hline(aes(yintercept=rate[AREA == "A"]), 
             linetype="dashed", color = "red") 
  ggtitle(as.character(DISEASE))
  

и безуспешно.

Имеет ли это какое-то отношение к тому, что БОЛЕЗНЬ становится ФАКТОРОМ при ее группировании?

Похоже, вам нужна функция, в которой вы можете ввести заболевание и создать график.

 disease_plot <- function(disease_of_interest) {
  DATA %>%
    filter(DISEASE == disease_of_interest) %>%
    ggplot(aes(x = AREA, y = rate))   
        geom_point()  
        geom_hline(aes(yintercept = rate[AREA == "A"]),
          linetype = "dashed", color = "red")  
        # labs(x = x_label, y = y_label)  
        ggtitle(disease_of_interest)
}

disease_plot("Marco Polio")
disease_plot("Chicky Pox")
disease_plot("Mumps")
  

Или создать их все сразу…

 map(unique(DATA$DISEASE), disease_plot)
  

В конце концов, я воспользовался советом и помощью Стивена и Марка и совместил это со своим первоначальным планом пройти через это, мурлыкая под нос.

Вот оно:

     walk(unique(DATA$DISEASE), function(disease_of_interest) {
      p <- DATA%>%filter(DISEASE== !!disease_of_interest)%>%
      ggplot(aes(x=AREA, y=rate,y=rate, 
      ymin = rate-lower, ymax = rate upper))  
      geom_point()  
      geom_hline(aes(yintercept=rate[AREA == "A"]), 
                   linetype="dashed", color = "red") 
      labs(x = x_label,y = y_label), 
      ggtitle(paste0("Number of ",disease_of_interest,
      " in 2018")) 
      geom_errorbar(aes(ymin=lower, ymax=upper), width=.1)
      print(p)
      ggsave(paste("drive path",disease_of_interest, "plot.png")) 
      scale_x_discrete(limits=c("C","A","B"))
    })
  

Это потому, что вы используете весь текст DATA$DISEASE в качестве заголовка, и кажется, что он просто захватывает последнее значение в нем column . Гораздо проще dataframe сначала создать отфильтрованный, а затем ввести его в plot , я думаю.

 df <- DATA%>%filter(DISEASE== "Marco Polio")

  ggplot(data = df, aes(x=AREA, y=rate))  geom_point()  
  geom_hline(aes(yintercept=rate[AREA == "A"]), 
             linetype="dashed", color = "red") 
  labs(x = x_label,y = y_label) 
  ggtitle(df$DISEASE)
  

В конце концов, я подумал, что решением было бы вызывать не DATA $, а просто DISEASE
Однако, похоже, это работает не так, как ожидалось, при фильтрации по другому заболеванию. Я думаю, вам пришлось бы также включить DISEASE в ggtitle, или лучше использовать первую функцию, или другой ответ, опубликованный другим пользователем.

РАБОТАЕТ НЕ ТАК, КАК ОЖИДАЛОСЬ:

 DATA%>%filter(DISEASE== "Marco Polio")%>%
  ggplot(aes(x=AREA, y=rate))  geom_point()  
  geom_hline(aes(yintercept=rate[AREA == "A"]), 
             linetype="dashed", color = "red") 
  labs(x = x_label,y = y_label) 
  ggtitle(DISEASE)