Как сделать хороший сюжет по этому сценарию

Question

У меня есть этот код:

new_flts.df %>%
  mutate(DELAY_CAT = fct_collapse(DELAY_CAT,
          OnTime = c("On Time", "Short Delay"),
          Delayed = c("Medium Delay", "Long Delay", "Very Long Delay"),
          Canceled =("CNL")
          ))%>%
  filter(DIR=="D") %>%
  group_by(AIRL, ORGDEST, DELAY_CAT)%>%
  tally()%>%mutate(perc = percent(n/sum(n),accuracy = 1))

Он производит следующий вывод (часть вывода):

# A tibble: 260 x 5
# Groups:   AIRL, ORGDEST [94]
    AIRL  ORGDEST DELAY_CAT     n perc 
    <fct> <fct>   <fct>     <int> <chr>
  1 ADR   DUS     OnTime       43 49%  
  2 ADR   DUS     Delayed      31 36%  
  3 ADR   DUS     Canceled     13 15%  
  4 ADR   FRA     OnTime     1673 86%  
  5 ADR   FRA     Delayed     237 12%  
  6 ADR   FRA     Canceled     26 1%   
 24 ART   DUS     OnTime       39 57%  
 25 ART   DUS     Delayed      26 38%  
 26 ART   DUS     Canceled      4 6%   
 27 ART   MLH     OnTime       77 53%  
 28 ART   MLH     Delayed      52 36%  
 29 ART   MLH     Canceled     17 12%  
 30 ART   MUC     OnTime       22 61%  
 31 ART   MUC     Delayed      11 31%  
 32 ART   MUC     Canceled      3 8%   
 44 BUR   DUS     OnTime       53 50%  
 45 BUR   DUS     Delayed      45 43%  
 46 BUR   DUS     Canceled      7 7%   
 47 BUR   HAJ     OnTime       17 45%  
 48 BUR   HAJ     Delayed      20 53%  
 49 BUR   HAJ     Canceled      1 3%   
 50 BUR   MLH     OnTime      158 57%  
 51 BUR   MLH     Delayed     101 37%  
 52 BUR   MLH     Canceled     16 6%   
 53 BUR   OTP     OnTime       30 45%  
 54 BUR   OTP     Delayed      30 45%  
 55 BUR   OTP     Canceled      6 9%   

После построения графика я получил лучшее (это изображение ниже). Но я не доволен этим. Я хочу показать хороший график с процентами производительности для каждой АВИАКОМПАНИИ и НАЗНАЧЕНИЯ. Любая помощь в создании хорошего информационного сюжета!

код участка:

ggplot(data, aes(AIRL, perc))+
  geom_col(aes(fill=DELAY_CAT))+facet_wrap(.~ORGDEST, scales = "free")

Answer 1

Моя интуиция заключается в том, что вам было бы полезно получить всю информацию на одном участке и не использовать facet_wrap в этом случае? Вы можете объединить две переменные, которые вас интересуют, используя paste. Смотрите этот пример.

library(tidyverse)

df <- mtcars

df <- df %>% 
  mutate(CylVs = paste(cyl,vs))

ggplot(df,aes(x = carb,fill = CylVs)) + 
  geom_histogram()

Удачи!