r - 更改刻面环绕中刻面的位置，按年份排序

Question

我想更改我的 ggplot Boxplot 的顺序和实验室。这是我到目前为止所拥有的：

ggplot(total_nr_obj_week_corr.3, aes(x=distance, y=n)) +
  stat_boxplot(geom = "errorbar")+
  geom_boxplot()+
  facet_wrap(~year_month, strip.position = "bottom",scales = "free", nrow=2)+
  scale_y_continuous(expand = c(0,0),limits = c(0, 70), breaks = seq(0, 80, by = 10)) +
  scale_x_discrete(limits = c("0", "5","10", "15","20"), breaks = c("0", "5","10", "15","20"))+
  theme(strip.placement="outside", 
        strip.background = element_rect(color="black",
fill="white",size = rel(1)),
        axis.title.x = element_text(margin = margin(t = 10, r = 0 , b = 0, l = 0)),
        axis.title.y = element_text(margin = margin(t = 0, r = 10 , b = 0, l = 0)))  

这是我的阴谋

我想改变不同方面的位置。第一排 2017 05-08 第二排 2018 05-09

我在 facet_wrap() 中尝试了 as.table = F，但仍然没有得到正确的顺序。

其次，我想更改标签。在 x 轴下方，我想要月份（05-08/09），在大标签上的每一行上方，每行带有年份（2017/2018）。

我现在用谷歌搜索了几个小时，但找不到解决方案。我也试过fact_grid，但不是我想要的结果

谢谢！！

标签： rggplot2boxplotfacetfacet-wrap

Answer 1

您需要更改 year_month 的级别，好的，首先我们创建一个函数来生成您的绘图，请注意我包含了一个 facet_order，我使用它重新调整您的 year_month：

func=function(DA,facet_order){

ggplot(data=DA, aes(x=distance, y=n)) +
  stat_boxplot(geom = "errorbar")+
  geom_boxplot()+
  facet_wrap(~factor(year_month,levels=facet_order), strip.position = "bottom",scales = "free", nrow=2)+
  scale_y_continuous(expand = c(0,0),limits = c(0, 70), breaks = seq(0, 80, by = 10)) +
  scale_x_discrete(limits = c("0", "5","10", "15","20"), breaks = c("0", "5","10", "15","20"))+
  theme(strip.placement="outside", 
        strip.background = element_rect(color="black",
fill="white",size = rel(1)),
        axis.title.x = element_text(margin = margin(t = 10, r = 0 , b = 0, l = 0)),
        axis.title.y = element_text(margin = margin(t = 0, r = 10 , b = 0, l = 0))) 
}

我模拟了一些看起来像你的数据的东西：

library(ggplot2)

DATES = c("2017 05","2017 06","2017 07","2017 08",
 "2018 05","2018 06","2018 07","2018 08","2018 09")

total_nr_obj_week_corr.3 = lapply(DATES,function(i){
  data.frame(
  distance=factor(rep(seq(0,20,by=5),each=10)),
  n = rnbinom(50,mu=20,siz=1),
  year_month = i
)
})

total_nr_obj_week_corr.3 = do.call(rbind,total_nr_obj_week_corr.3)

现在，如果我们使用原始关卡，我们会得到与您的情节类似的东西：

lvl=levels(total_nr_obj_week_corr.3$year_month)
func(total_nr_obj_week_corr.3,lvl)+ggtitle("original levels")

一种“转置”关卡的快速方法：

NCOL=5
NROW=2
LAY_FACET = matrix(NA,ncol=NCOL,nrow=NROW)
LAY_FACET[1:length(lvl)]=lvl
new_lvl = na.omit(c(t(LAY_FACET)))

func(total_nr_obj_week_corr.3,new_lvl)+ggtitle("new levels")