Удаление одного уровня/группы из Facet_wrap ggplot2 в R

Мой tbl_df:

> str(p2p_dt_SKILL_A)
Classes ‘tbl_dt’, ‘tbl’, ‘data.table’ and 'data.frame': 693 obs. of  35 variables:
 $ Patch          : Factor w/ 7 levels "BVG1","BVG11",..: 1 2 3 4 5 6 7 1 2 3 ...
 $ Skill          : Factor w/ 15 levels "A","BROADBAND",..: 1 1 1 1 1 1 1 1 1 1 ...
 $ Date           : Date, format: "2015-12-04" "2015-12-04" "2015-12-04" "2015-12-04" ...
 $ SOD_FWIH_A     : num  1.09 1.14 1.09 1.1 1.09 1.07 1.09 1.07 1.12 1.07 ...
 $ Prod_MWF       : num  5.06 5.06 5.44 5.34 4.22 4.72 4.89 4.68 4.68 5.22 ...
 $ Prod_MA        : num  4.41 4.41 4.81 4.41 4.68 4.15 4.26 4.44 4.44 4.63 ...
 $ Prod_DL        : num  3.43 3.49 3.45 3.57 3.43 3.47 3.28 3.44 3.42 3.72 ...

создал следующий линейный график, используя ggplot2:

ggplot(p2p_dt_SKILL_A,aes(x=Date,y=Prod_DL)) +
        geom_line(data = p2p_dt_SKILL_A[p2p_dt_SKILL_A$Patch %in% c("BVG1"),],aes(colour="red"),lwd=1.3) +
        geom_smooth() +
        scale_x_date(labels=date_format("%b-%y"),breaks ="1 month")+
        geom_vline(xintercept = as.numeric(p2p_dt_SKILL_A$Date[p2p_dt_SKILL_A$Date=="2015-09-18"]))+
        geom_vline(xintercept = as.numeric(p2p_dt_SKILL_A$Date[p2p_dt_SKILL_A$Date=="2015-10-02"]))+
        geom_vline(xintercept = as.numeric(p2p_dt_SKILL_A$Date[p2p_dt_SKILL_A$Date=="2015-10-23"]))+
        theme(axis.title.y = element_text(size = 15,face="bold",color="red"),
              plot.title = element_text(size = 15,lineheight = .8,face="bold",color="red"),
              axis.title.x = element_blank(),
              legend.position="none")

Теперь я пытаюсь использовать facet_wrap, но без "BVG1" Patch. Но facet_wrap также поддерживает патч BVG1:

ggplot(p2p_dt_SKILL_A,aes(x=Date,y=Prod_DL)) +
        geom_line(data = p2p_dt_SKILL_A[p2p_dt_SKILL_A$Patch != "BVG1",],aes(colour="red"),lwd=1.3) +
        scale_colour_discrete(drop = FALSE)+
        geom_smooth() +
        geom_line(stat = "hline", yintercept = "mean")+
        scale_x_date(labels=date_format("%b-%y"),breaks ="2 month")+
        geom_vline(xintercept = as.numeric(p2p_dt_SKILL_A$Date[p2p_dt_SKILL_A$Date=="2015-09-18"]))+

         geom_vline(xintercept = as.numeric(p2p_dt_SKILL_A$Date[p2p_dt_SKILL_A$Date=="2015-10-02"]))+
         geom_vline(xintercept = as.numeric(p2p_dt_SKILL_A$Date[p2p_dt_SKILL_A$Date=="2015-10-23"]))+
        ylab("DL Prod for All Skills")+
        ggtitle("BVG1 DL Prod for All Skills 2014-2015")+
        theme(axis.title.y = element_text(size = 15,face="bold",color="red"),
              plot.title = element_text(size = 15,lineheight = .8,face="bold",color="red"),
              axis.title.x = element_blank(),
              legend.position="none")+
        facet_wrap(~Patch, ncol = 3)

Как заставить facet_wrap просто печатать от BVG11 до BVG16 и исключить группу BVG1?

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

som_data <- p2p_dt %>%
                        filter(Patch != "BVG1" & Skill %in% c("A"))

> str(som_data)
Classes ‘tbl_dt’, ‘tbl’, ‘data.table’ and 'data.frame': 594 obs. of  39 variables:
 $ Patch          : Factor w/ 6 levels "BVG11","BVG12",..: 1 2 3 4 5 6 1 2 3 4 ...
 $ Skill          : Factor w/ 1 level "A": 1 1 1 1 1 1 1 1 1 1 ...
 $ Date           : Date, format: "2015-12-04" "2015-11-27" "2015-11-20" "2015-11-13" ...
 $ SOD_FWIH_A     : num  1.14 1.12 1.15 1.09 1.13 1.14 1.15 1.16 1.14 1.11 ...
 $ Prod_MWF       : num  5.06 4.68 4.69 4.76 4.63 4.98 4.45 5.12 4.44 4.54 ...
 $ Prod_MA        : num  4.41 4.44 4.14 5.63 4.86 4.97 4.65 4.99 4.52 4.9 ...
 $ Prod_DL        : num  3.49 3.42 3.6 3.29 3.19 3.1 3.18 3.28 3.26 3.21 ...

используя вновь созданный кадр данных выше, он создает следующий график, и это также неверно

ggplot(som_data,aes(x=Date,y=Prod_DL)) +
        geom_line(aes(colour="red"),lwd=1.3) +
        geom_smooth() +
        geom_line(stat = "hline", yintercept = "mean")+
        scale_x_date(labels=date_format("%b-%y"),breaks ="2 month")+
        geom_vline(xintercept = as.numeric(som_data$Date[som_data$Date=="2015-09-18"]))+
        geom_vline(xintercept = as.numeric(som_data$Date[som_data$Date=="2015-10-02"]))+
        geom_vline(xintercept = as.numeric(som_data$Date[som_data$Date=="2015-10-23"]))+
        ylab("DL Prod for All Skills")+
        ggtitle("BVG1 SOMs DL Prod for All Skills 2014-2015")+
        theme(axis.title.y = element_text(size = 15,face="bold",color="red"),
              plot.title = element_text(size = 15,lineheight = .8,face="bold",color="red"),
              axis.title.x = element_blank(),
              legend.position="none")+
        facet_wrap(~Patch)

Любые идеи???

r time-series ggplot2 facet-wrap

Shery 13.12.2015 источник

comment

Возможно, вы можете использовать ?droplevels - talat 13.12.2015

comment

Почему бы не использовать som_data <- p2p_dt_SKILL_A %>% filter(Patch != "BVG1"), а затем, возможно, p2p_dt_SKILL_A <- droplevels(p2p_dt_SKILL_A). - Axeman 13.12.2015

comment

в вашем первом наборе данных Skill имеет 15 уровней, во втором - только 1 уровень - MLavoie 13.12.2015

comment

@docendodiscimus Я сбросил уровни, и у него есть только один навык и оставшиеся патчи (уверяю вас в этом), но все равно не повезло - Shery 13.12.2015

comment

Первый кадр данных также фильтруется на Skill=="A" - Shery 13.12.2015

comment

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

comment

@Axeman обновил ответ с помощью dput - Shery 13.12.2015

Ответы (1)

arrow_upward
4
arrow_downward

Я не вижу, где проблема. Вот минимальный пример с предоставленными вами данными.

library(dplyr)
# I first need to get rid of the data.table
sample_data <- as.data.frame(sample_data)
# I don't want all those other variables you gave me either..
sample_data <- select(sample_data, Date, Prod_DL, Patch)

# Now filter out the Patch that we don't want
sample_data2 <- filter(sample_data, Patch != "BVG1")

# Create a minimal plot
ggplot(sample_data2, aes(x = Date, y = Prod_DL)) +
  geom_line(colour="red", lwd = 1.3) +
  geom_smooth() +
  facet_wrap(~Patch, ncol = 3)

Не то чтобы отличный сюжет (данные ограничены), но всего шесть панелей.

Axeman 13.12.2015

comment

это смешно, но я решил проблему. Я не делал ничего особенного, просто перезапустил команду, и она сработала! но спасибо за ваши усилия в любом случае :) - Shery; 13.12.2015

comment

Пожалуйста, в следующий раз сделайте минимальный воспроизводимый пример, чтобы вы могли найти любые ошибки. - Axeman; 13.12.2015

Удаление одного уровня/группы из Facet_wrap ggplot2 в R

Ответы (1)

Вопросы по теме