使用线和多因素连接ggplot箱图

Question

我正在尝试将ggplot2箱线图与geom_lines连接以考虑多种因素。 到目前为止，我已经能够完成将所有箱线图与线连接的操作，请参见附件图片。 但我希望通过相应的因素来连接唯一的箱形图。

例如，对于我的变量FL，我只想连接这两个箱形图，而不将它们与其余变量连接。 类似地，对于变量RW，将这两个性别盒图连接起来，而没有剩下的。

library("MASS")  
data(crabs)  
melt_crabs <- melt(crabs,id.var=c("sp","sex","index"))   
ggplot(melt_crabs, aes(x = variable, y = value)) +   geom_line(aes(group = index), size = 0.05, alpha = 0.7) +   geom_boxplot(aes(fill = sp), alpha = 0.5) + facet_grid(sex~.)

有谁知道如何实现这一目标？ 我希望我能以最清晰的方式向自己解释。

非常感谢和最良好的祝愿，

Answer 1

我不知道如何将点精确连接到您生成的图中。 但我可以向您展示如何做类似的事情。

困难在于，属于一对箱形图的所有点共享相同的x坐标（即， variable值相同）。 因此，我已将interaction(sex, variable)用作x坐标，以便每个箱线图都有自己的x值。 但是，这意味着成对的箱形图不太明显。 另一方面，连接线沿另一个方向工作。

为了绘制线，将它们按interaction(index, variable)分组，这意味着当数据点共享相同的index和variable值时，数据点已连接。

聊够了，下面是代码：

ggplot(melt_crabs, aes(x = interaction(sex, variable), y = value)) +
    geom_boxplot(aes(fill = sex), alpha = 0.5) +
    geom_line(aes(group = interaction(index, variable)),
              alpha = 0.5, colour = "darkgrey") +
    facet_grid(sp~.) +
    scale_x_discrete(labels = rep(levels(melt_crabs$variable), each = 2))

Answer 2

@Stibu的答案绝对是获得所需输出的最快，最干净的方法。 建立他的答案，因为对于给定的变量，曲线不再相邻，因此您可能希望将每个曲线放在各自的曲线中。 只需对@Stibu的代码中的facet_grid调用进行较小的调整即可轻松实现：

ggplot(melt_crabs, aes(x = interaction(sex, variable), y = value)) +
  geom_boxplot(aes(fill = sex), alpha = 0.5) +
  geom_line(aes(group = interaction(index, variable)),
            alpha = 0.5, colour = "darkgrey") +
  facet_grid(sp~variable,scales="free_x") +
  scale_x_discrete(labels = "")

另外，如果您想将顶部的条带移动到底部，我也做了一些挖掘，这并不困难。 从以下代码改编代码： 刻面时如何在图下方显示条形标签？

我们得到：

p<-ggplot(melt_crabs, aes(x = interaction(sex, variable), y = value)) +
  geom_boxplot(aes(fill = sex), alpha = 0.5) +
  geom_line(aes(group = interaction(index, variable)),
            alpha = 0.5, colour = "darkgrey") +
  facet_grid(sp~variable,scales="free_x") +
  scale_x_discrete(labels = "")

# Convert the plot to a grob
gt <- ggplotGrob(p)

# Get the positions of the panels in the layout: t = top, l = left, ...
panels <-c(subset(gt$layout, name == "panel", select = t:r))

# Get the strip grob & x-axis label grob
stripGrob = gtable_filter(gt, "strip-top")

#Replace x-axis ticks with strip
gt = gtable_add_grob(gt, stripGrob, t = max(panels$b)+1, l = min(panels$l), r = max(panels$r))

# remove the old strip
gt = gt[-(min(panels$t)-1), ]

grid.newpage()
grid.draw(gt)