ggplot2：3路互动叠加条形图的分组条形图

Question

我尝试在面板中绘制叠加条形图+ 2路互动，其中包括4个实验的相同图表。但是，我不能躲避酒吧取决于自变量之一。以下是我的数据。

首先我通过下面的代码读取数据。

a<-read.table(file.choose(), header=T, dec=",") 

Exp. \t Gest \t lag \t Sint12 \t Rev12 \t c12 \t t1pi \t t2pi \t t1i \t t2i \t IntWeak \t inc \t Total 
 
1 \t 1 \t 1 \t 15,88 \t 3,28 \t 22,52 \t 11,76 \t 4,08 \t 2,28 \t 16,76 \t 3,24 \t 20,2 \t 100 
 
1 \t 1 \t 3 \t 0,88 \t 1,2 \t 61,36 \t 11,84 \t 8,4 \t 1,84 \t 2,32 \t 0,8 \t 11,36 \t 100 
 
1 \t 1 \t 8 \t 0,24 \t 0,24 \t 65,2 \t 10,24 \t 9,2 \t 1,84 \t 2,4 \t 0,48 \t 10,16 \t 100 
 
1 \t 2 \t 1 \t 14,96 \t 4 \t 25,28 \t 15,12 \t 1,92 \t 0,68 \t 16,8 \t 1,56 \t 19,68 \t 100 
 
1 \t 2 \t 3 \t 1,2 \t 0,72 \t 79,36 \t 8,64 \t 2,88 \t 0,64 \t 0,64 \t 0,64 \t 5,28 \t 100 
 
1 \t 2 \t 8 \t 0,16 \t 0,16 \t 86,72 \t 5,36 \t 3,2 \t 0,08 \t 0,48 \t 0,64 \t 3,2 \t 100 
 
2 \t 1 \t 1 \t 30,6 \t 2,2 \t 24,48 \t 4,56 \t 1,32 \t 0,4 \t 17,8 \t 1 \t 17,64 \t 100 
 
2 \t 1 \t 3 \t 0,96 \t 1,04 \t 87,2 \t 5,04 \t 2,16 \t 0,16 \t 0,4 \t 0,8 \t 2,24 \t 100 
 
2 \t 1 \t 8 \t 0,88 \t 0,24 \t 91,92 \t 3,28 \t 1,52 \t 0 \t 0,32 \t 0,88 \t 0,96 \t 100 
 
2 \t 2 \t 1 \t 20,16 \t 2,32 \t 16,52 \t 14,24 \t 0,72 \t 0,44 \t 15,96 \t 1,76 \t 27,88 \t 100 
 
2 \t 2 \t 3 \t 1,04 \t 0,64 \t 83,84 \t 5,84 \t 2 \t 0,08 \t 0,72 \t 1,12 \t 4,72 \t 100 
 
2 \t 2 \t 8 \t 0,24 \t 0 \t 91,04 \t 4,16 \t 1,52 \t 0,08 \t 0 \t 0,72 \t 2,24 \t 100 
 
3A \t 1 \t 1 \t 35,83 \t 3,92 \t 27,42 \t 2,42 \t 2,08 \t 0,25 \t 7,42 \t 3,63 \t 17,04 \t 100,01 
 
3A \t 1 \t 3 \t 1,58 \t 1 \t 81 \t 4,5 \t 3,33 \t 0,25 \t 0,33 \t 1,08 \t 6,92 \t 99,99 
 
3A \t 1 \t 8 \t 1 \t 0 \t 86,92 \t 3,17 \t 1,75 \t 0,08 \t 0,42 \t 0,33 \t 6,33 \t 100 
 
3A \t 2 \t 1 \t 43,46 \t 2,38 \t 21,29 \t 1,88 \t 1,17 \t 0,17 \t 5,46 \t 4,21 \t 20 \t 100,02 
 
3A \t 2 \t 3 \t 2 \t 0,75 \t 78,67 \t 3,75 \t 3,25 \t 0,17 \t 0,83 \t 0,92 \t 9,67 \t 100,01 
 
3A \t 2 \t 8 \t 1,33 \t 0,33 \t 83,25 \t 3 \t 2,17 \t 0 \t 0,67 \t 0,83 \t 8,42 \t 100 
 
3B \t 1 \t 1 \t 35,5 \t 2,54 \t 29,33 \t 3,04 \t 1,88 \t 0,54 \t 7,46 \t 7,46 \t 12,25 \t 100 
 
3B \t 1 \t 3 \t 1,58 \t 0,67 \t 79,42 \t 4,58 \t 2,83 \t 0,42 \t 0,67 \t 2,75 \t 7,08 \t 100 
 
3B \t 1 \t 8 \t 0,83 \t 0,17 \t 88,83 \t 3,17 \t 2,83 \t 0,08 \t 0,42 \t 0,5 \t 3,17 \t 100 
 
3B \t 2 \t 1 \t 32,33 \t 1,75 \t 17,21 \t 4,5 \t 2,21 \t 0,42 \t 13,21 \t 4,96 \t 23,42 \t 100,01 
 
3B \t 2 \t 3 \t 2,5 \t 0,25 \t 67,58 \t 8,42 \t 4,25 \t 0,5 \t 1 \t 4,58 \t 10,92 \t 100 
 
3B \t 2 \t 8 \t 1 \t 0,08 \t 76,83 \t 6,25 \t 4,5 \t 0,08 \t 0,33 \t 3 \t 7,92 \t 99,99

二我用下面的代码变换它宽以长格式。

b <- reshape(a, 
     varying = c("Sint12", "Rev12", "c12", "t1pi", "t2pi", "t1i", "t2i", "IntWeak", "inc"), 
     v.names = "score", 
     timevar = "variable", 
     times = c("Sint12", "Rev12", "c12", "t1pi", "t2pi", "t1i", "t2i", "IntWeak", "inc"), 
     new.row.names = 1:1000, 
     direction = "long") 

和改造后的数据看起来象下面这样：

Exp. Gest lag Total variable score id 
 
1  1 1 1 100.00 Sint12 15.88 1 
 
2  1 1 3 100.00 Sint12 0.88 2 
 
3  1 1 8 100.00 Sint12 0.24 3 
 
4  1 2 1 100.00 Sint12 14.96 4 
 
5  1 2 3 100.00 Sint12 1.20 5 
 
6  1 2 8 100.00 Sint12 0.16 6 
 
7  2 1 1 100.00 Sint12 30.60 7 
 
8  2 1 3 100.00 Sint12 0.96 8 
 
9  2 1 8 100.00 Sint12 0.88 9 
 
10  2 2 1 100.00 Sint12 20.16 10 
 
11  2 2 3 100.00 Sint12 1.04 11 
 
12  2 2 8 100.00 Sint12 0.24 12 
 
13 3A 1 1 100.01 Sint12 35.83 13 
 
14 3A 1 3 99.99 Sint12 1.58 14 
 
15 3A 1 8 100.00 Sint12 1.00 15 
 
16 3A 2 1 100.02 Sint12 43.46 16 
 
17 3A 2 3 100.01 Sint12 2.00 17 
 
18 3A 2 8 100.00 Sint12 1.33 18 
 
19 3B 1 1 100.00 Sint12 35.50 19 
 
20 3B 1 3 100.00 Sint12 1.58 20 
 
21 3B 1 8 100.00 Sint12 0.83 21 
 
22 3B 2 1 100.01 Sint12 32.33 22 
 
23 3B 2 3 100.00 Sint12 2.50 23 
 
24 3B 2 8 99.99 Sint12 1.00 24

我想要的是什么; 1。 4个绘图（对于每个实验），2.通过Gest和滞后来制作交互图。第3;用变量的颜色填充堆栈。

为了做到这一点，我使用了下面的代码。

ggplot（data = b，aes（x = interaction（Gest，lag），y = score，fill = variable，））+ geom_bar（stat =“identity”）+ facet_wrap（〜Exp。，ncol = 2 ）

Plot

现在，剧情已准备就绪。但是，当我将position = dodge参数传递给geom_bar时;这是行不通的。我想有一个情节，在1.1 & 2.1之间没有差距; 1.3 & 2.3和1.8 & 2.8（X轴标签）。另外，我想指定.1-.3和.8之间的差距。

在此先感谢。

Answer 1

你是什么意思的“它不工作”？当我补充说，我得到了正确闪避酒吧的阴谋。

如果您只是想要修改间距，可以通过向沿x绘制的因子添加附加（空）级别来实现。请注意，在每个附加水平的不同数量的空间：

ggplot(data=b 
     , aes(x= factor(interaction(Gest, lag) 
         , levels = c(1.1,2.1," ",1.3,2.3," ",1.8,2.8)) 
      ,y=score, fill = variable)) + 
    geom_bar(stat="identity", position = "dodge") + 
    facet_wrap(~Exp., ncol=2) + 
    scale_x_discrete(drop = FALSE) 

您可以进一步走几步，如果你使用cowplot包。这里的优点是，如果分别制作每个实验图，则可以在lag变量上进行切面。然后，你可以将它们缝合在一起。在这里，我压制了个人传奇，并在底部添加了一个共同的常见传奇。

sepPlots <- lapply(unique(b$Exp.), function(thisExp){ 
    b %>% 
    filter(Exp. == thisExp) %>% 
    ggplot(aes(x = as.factor(Gest) 
       , y = score 
       , fill = variable)) + 
    geom_bar(stat="identity", position = "dodge") + 
    facet_wrap(~lag, nrow = 1 
       , labeller = label_both 
       , switch = "x") + 
    xlab("Gest") + 
    ggtitle(paste("Experiment:", thisExp)) 

}) 


expPlots <- 
    plot_grid(plotlist = lapply(sepPlots, function(x){x + guides(fill = "none")})) 

plot_grid(expPlots 
      , get_legend(sepPlots[[1]] + theme(legend.direction = "horizontal")) 
      , nrow = 2 
      , rel_heights = c(1, 0.1))