寻找由ggplot 2.2.0打破的gtable_add_grob代码的解决方法 [英] Seeking workaround for gtable_add_grob code broken by ggplot 2.2.0
问题描述
在具有多个方面变量的图中，ggplot2重复外部变量的facet标签，而不是在内部变量的所有级别上都有单个跨越方面条。我已经使用了一些代码，用于使用 不幸的是，这个代码不再适用于ggplot2 2.2.0，因为facet条的grob结构发生了变化。具体来说，在以前的ggplot2版本中，每一个facet标签都有自己的一组grob。然而，在版本2.2.0中，每个垂直堆栈的小平面标签看起来都是单一的grob。这打破了我的代码，我不知道如何解决它。 下面是一个具体的例子，摘自 现在我们将介绍前两名我们没有重复剥离标签： 这就是ggplot2 2.1.0的情况： 啰嗦的介绍，这里是我的问题：我如何创建与ggplot2版本2.2.0跨越的方面带，看起来像我使用ggplot2版本2.1.0创建的 ggplot 2.1.0 gtable_add_grob 从 gtable 使用单个生成的小平面条覆盖重复的外部小平面标签c $ c>包。
pg = ggplotGrob（p）
＃为物种
pos = c（4,11）
添加spanning strip标签（i in 1：2）{
pg < - gtable_add_grob（pg，
list（rectGrob（gp = gpar（col =grey50，fill =grey90）），
textGrob（unique（densityAGRLKA $ species）[i]，
gp = gpar（cex = 0.8））），t = 3,1 = pos [i]，b = 3，r = pos [i] +7，
name = c（a，b））
}
＃为位置
pos = c（4,7,11,15）
添加生成带标签（i in 1：4）{
pg = gtable_add_grob（pg，
list（rectGrob（gp = gpar（col =grey50，fill =grey 90）），
textGrob（rep（unique（densityAGRLKA $ location），2）[i]，
gp = gpar（cex = 0.8））），t = 4，l = pos [i ]，b = 4，r = pos [i] +3，
name = c（c，d））
}
grid.draw（pg ）
gtable_add_grob ？我希望有一个简单的调整，但如果它需要大手术，那么也没关系。
pg
TableGrob（9 x 19）布局：45 grobs
z单元格名称grob
2 1（3- 3，4- 4）strip-top absoluteGrob [strip.absoluteGrob.147]
3 2（4- 4，4- 4）strip-top absoluteGrob [strip.absoluteGrob.195]
4 3（5- 5，4- 4）strip-top absoluteGrob [strip.absoluteGrob.243]
5 4（ 3-6,6-6）strip-top absoluteGrob [strip.absoluteGrob.201]
7 5（4- 4,6-6）strip-top absoluteGrob [strip.absoluteGrob.201]
7 6（5-6,6-6）strip-top absoluteGrob [strip.absoluteGrob.249]
8 7（3-3,8-8）strip-top absoluteGrob [strip.absoluteGrob.159]
9 8（4- 4,8-8）strip-top absoluteGrob [strip.absoluteGrob.207]
10 9（5- 5,8-8）strip-top absoluteGrob [strip.absoluteGrob.255]
11 10（3- 3,10-10）strip-top absoluteGrob [strip.absoluteGrob.16 5]
12 11（4- 4,10-10）strip-top absoluteGrob [strip.absoluteGrob.213]
13 12（5- 5,10-10）strip-top absoluteGrob [strip。绝对Grob.261]
14 13（3- 3,12-12）strip-top absoluteGrob [strip.absoluteGrob.171]
15 14（4- 4,12-12）strip-top absoluteGrob [ strip.absoluteGrob.219]
16 15（5- 5,12-12）strip-top absoluteGrob [strip.absoluteGrob.267]
17 16（3- 3,14-14）strip-top绝对Grob [strip.absoluteGrob.177]
18 17（4- 4,14-14）strip-top absoluteGrob [strip.absoluteGrob.225]
19 18（5- 5,14-14）strip -top absoluteGrob [strip.absoluteGrob.273]
20 19（3- 3,16-16）strip-top absoluteGrob [strip.absoluteGrob.183]
21 20（4- 4,16-16 ）strip-top absoluteGrob [strip.absoluteGrob.231]
22 21（5- 5,16-16）strip-top absoluteGrob [strip.absoluteGrob.279]
23 22（3- 3,18 -18）strip-top absoluteGrob [strip.absoluteGrob.189]
24 23（4- 4,18-18）strip-top absoluteGrob [strip.absoluteGrob.237]
25 24（5- 5 ，18-1 8）strip-top absoluteGrob [strip.absoluteGrob.285]
ggplot2 2.2.0
pg
TableGrob（11 x 21）layout：42 grobs
z单元格名称grob
28 2（6- 6，4- 4）strip-t-1 gtable [strip]
29 2（6- 6,6-6）strip-t-2 gtable [strip]
30 2（6- 6,8-8）strip-t-3 gtable [strip]
31 2（6- 6,10-10）strip-t-4 gtable [strip]
32 2（6- 6,12-12）strip-t-5 gtable [strip]
33 2（6- 6,14-14）strip-t-6 gtable [strip]
34 2 （6- 6,16-16）strip-t-7 gtable [strip]
35 2（6- 6,18-18）strip-t-8 gtable [strip]
library（ggplot2）
library（gtable）
library（grid）
＃您的数据
df =结构（列表（位置=结构（c（1L，1L，1L，1L，1L，1L，
1L，1L，1L，1L，2L ，2L，2L，2L，2L，2L，2L，2L，2L，2L，1L，1L，
1L，1L，1L，1L，1L，1L，1L，1L，2L，2L，2L，2L ，2L，2L，2L，2L，
2L，2L），.Label = c（SF，SS），class =factor），species = structure（c（1L，$ b $ 1L，1L，1L，1L，1L，1L，1L，1L，1L，1L，1L，1L，1L，1L，1L，1L，1L，1L，2L，2L，2L， ，2L，2L，2L，2L，2L，2L，2L，2L，2L，2L，2L，2L，2L，2L，2L，2L），标签= c（AGR，LKA ），class =factor），
position = structure（c（1L，1L，1L，1L，1L，2L，2L，2L，2L，
2L，1L，1L，1L，1L ，1L，2L，2L，2L，2L，2L，2L，1L，1L，1L，1L，1L，1L，1L，2L，2L，2L，2L，2L， ，2L，2L，
2L），.Label = c（top，bottom），class =factor），density = c（0.41，
0.41,0.43,0.33,0.35 ，0.43,0.34,0.46,0.32,0.32， 0.4，
0.4,0.45,0.34,0.39,0.39,0.31,0.38,0.48,0.3,0.42，
0.34,0.35,0.4,0.38,0.42,0.36,0.34,0.46,0.38,0.36，
0.39,0.38,0.39,0.39,0.39,0.36,0.39,0.51,0.38）），.Names = c（location，
species，position，density）， row.names = c（NA，-40L），class =data.frame）
＃你有三个方面的ggplot
p = ggplot（df，aes（，density ））+
geom_boxplot（width = 0.7，position = position_dodge（0.7））+
theme_bw（）+
facet_grid（。 〜物品+位置+位置）+
主题（panel.spacing =单元（0，lines），
strip.background = element_rect（color =grey30，fill =grey90），
panel.border = element_rect（color =grey90），
axis.ticks.x = element_blank（））+
labs（x =）
＃获取ggplot grob
pg = ggplotGrob（p）
＃获取最左边的条带
index = which（pg $ layout $ name ==strip-t- 1）
strip1 = pg $ grobs [[index]]
＃绘制条形图
grid.newpage（）
grid.draw（strip1）
＃检查其布局
strip1 $ layout
gtable_show_layout（strip1）
获得外带标签'跨越'内标签的一种简单方法是从零开始构建带子：
＃从原始图中获得条形图，列表中的条形图为
strip = list（）
（我在1：8中）{
index = which（pg $ layout $ name == paste0 （strip-t - ，i））
strip [[i]] = pg $ grobs [[index]]
}
＃构建gtable以包含新条带
newStrip = gtable（widths = unit（rep（1，8），null），heights = strip [[1]] $ heights）
## Populate gtable
＃Top row
for（i in 1：2）{
newStrip = gtable_add_grob（newStrip，strip [[4 * i-3]] [1]，
t = 1，l = 4 * i-3，r = 4 * i）
}
＃中间行
（i在1：4）{
newStrip = gtable_add_grob（newStrip，strip [[2 * i-1]] [2]，
t = 2，l = 2 * i-1，r = 2 * i）
}
＃底行
（我在1：8）{
newStrip = gtable_add_grob（newStrip，strip [[i]] [3]，
t = 3，l = i ）
}
＃将条带放入图中
＃（可能会更好地移除原始条带。
＃在这种情况下，用彩色背景，它并不重要）
pgNew = gtable_add_grob（pg，newStrip，t = 6，l = 4，r = 18）
＃绘制绘图
grid.newpage（）
grid.draw（pgNew）
$ b $
pg = ggplotGrob（p ）
＃从原始图中获取条带列表
strip = lapply（grep（strip -t，pg $ layout $ name），function（x）{pg $ grobs [[x]]}）
＃构造gtable以包含新条带
newStrip = gtable（widths = unit（rep（1，8），null），heights = strip [[1]] $ heights）
##填充gtable
＃顶行
cols = seq（1，by = 4，length.out = 2）
newStrip = gtable_add_grob（newStrip，lapply（strip [cols]，`[`，1），t = 1，l = cols，r = cols + 3）
＃中间行
cols = seq（1，by = 2，length.out = 4）
newStrip = gtable_add_grob（newStrip，lapply（strip [cols]，`[`，2），t = 2，l = cols，r = cols + 1）
＃底行
newStrip = gtable_add_grob（newStrip，lapply（strip，```，3），t = 3，l = 1：8）
＃将条带放入图中
pgNew = gtable_add_grob（pg，newStrip，t = 6，l = 4，r = 18）
＃绘制图
grid.newpage（）
grid.draw（pgNew）
In plots with multiple facet variables, ggplot2 repeats the facet label for the "outer" variable, rather than having a single spanning facet strip across all the levels of the "inner" variable. I have some code that I've been using to cover the repeated outer facet labels with a single spanning facet strip using gtable_add_grob from the gtable package.
Unfortunately, this code no longer works with ggplot2 2.2.0 due to changes in the grob structure of the facet strips. Specifically, in previous versions of ggplot2, each row of facet labels got their own set of grobs. However, in version 2.2.0 it looks like each vertical stack of facet labels is a single grob. This breaks my code and I'm not sure how to fix it.
Here's a concrete example, taken from an SO question I answered a few months ago:
# Data
df = structure(list(location = structure(c(1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L), .Label = c("SF", "SS"), class = "factor"), species = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("AGR", "LKA"), class = "factor"),
position = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L), .Label = c("top", "bottom"), class = "factor"), density = c(0.41,
0.41, 0.43, 0.33, 0.35, 0.43, 0.34, 0.46, 0.32, 0.32, 0.4,
0.4, 0.45, 0.34, 0.39, 0.39, 0.31, 0.38, 0.48, 0.3, 0.42,
0.34, 0.35, 0.4, 0.38, 0.42, 0.36, 0.34, 0.46, 0.38, 0.36,
0.39, 0.38, 0.39, 0.39, 0.39, 0.36, 0.39, 0.51, 0.38)), .Names = c("location",
"species", "position", "density"), row.names = c(NA, -40L), class = "data.frame")
# Begin with a regular ggplot with three facet levels
p=ggplot(df, aes("", density)) +
geom_boxplot(width=0.7, position=position_dodge(0.7)) +
theme_bw() +
facet_grid(. ~ species + location + position) +
theme(panel.margin=unit(0,"lines"),
strip.background=element_rect(color="grey30", fill="grey90"),
panel.border=element_rect(color="grey90"),
axis.ticks.x=element_blank()) +
labs(x="")
We start with a plot that has three levels of facets.
Now we'll cover the top two facet strips with spanning strips so that we don't have repeated strip labels:
pg = ggplotGrob(p)
# Add spanning strip labels for species
pos = c(4,11)
for (i in 1:2) {
pg <- gtable_add_grob(pg,
list(rectGrob(gp=gpar(col="grey50", fill="grey90")),
textGrob(unique(densityAGRLKA$species)[i],
gp=gpar(cex=0.8))), t=3,l=pos[i],b=3,r=pos[i]+7,
name=c("a","b"))
}
# Add spanning strip labels for location
pos=c(4,7,11,15)
for (i in 1:4) {
pg = gtable_add_grob(pg,
list(rectGrob(gp = gpar(col="grey50", fill="grey90")),
textGrob(rep(unique(densityAGRLKA$location),2)[i],
gp=gpar(cex=0.8))), t=4,l=pos[i],b=4,r=pos[i]+3,
name = c("c","d"))
}
grid.draw(pg)
This is what this plot looks like with ggplot2 2.1.0:
However, if I try the same code with ggplot2 2.2.0, I get the original plot back, with no changes to the strip labels. A look at the grob structure of the original plot p suggests why this is happening. I've pasted in the grob tables at the bottom of this question. In order to save space, I've included only the rows related to the facet strips.
Looking at the cells column, note that in the 2.1.0 version of the plot the first two numbers in each row are either 3, 4, or 5, indicating the vertical position of the grob relative to the other grobs in the plot. In the code above, the t and l arguments to gtable_add_grob are set to values of 3 or 4 because those are the facet strip rows that I wanted to cover with spanning strips.
Now look at the cells column in the 2.2.0 version of the plot: Note that the first two numbers are always 6. Also note that the facet strips are comprised of only 8 grobs instead of 24 in version 2.1.0. In version 2.2.0, it seems that each stack of three facet labels is now a single grob instead of three separate grobs. So even if I change the t and b arguments in gtable_add_grob to 6, all three facet strips are covered. Here's an example:
pg = ggplotGrob(p)
# Add spanning strip labels for species
pos = c(4,11)
for (i in 1:2) {
pg <- gtable_add_grob(pg,
list(rectGrob(gp=gpar(col="grey50", fill="grey90")),
textGrob(unique(densityAGRLKA$species)[i],
gp=gpar(cex=0.8))), t=6,l=pos[i],b=6,r=pos[i]+7,
name=c("a","b"))
}
So, after that very long-winded introduction, here's my question: How can I create spanning facet strips with ggplot2 version 2.2.0 that look like the ones I created using gtable_add_grob with ggplot2 version 2.1.0? I'm hoping there's a simple tweak, but if it requires major surgery, well, that's okay too.
ggplot 2.1.0
pg
TableGrob (9 x 19) "layout": 45 grobs z cells name grob 2 1 ( 3- 3, 4- 4) strip-top absoluteGrob[strip.absoluteGrob.147] 3 2 ( 4- 4, 4- 4) strip-top absoluteGrob[strip.absoluteGrob.195] 4 3 ( 5- 5, 4- 4) strip-top absoluteGrob[strip.absoluteGrob.243] 5 4 ( 3- 3, 6- 6) strip-top absoluteGrob[strip.absoluteGrob.153] 6 5 ( 4- 4, 6- 6) strip-top absoluteGrob[strip.absoluteGrob.201] 7 6 ( 5- 5, 6- 6) strip-top absoluteGrob[strip.absoluteGrob.249] 8 7 ( 3- 3, 8- 8) strip-top absoluteGrob[strip.absoluteGrob.159] 9 8 ( 4- 4, 8- 8) strip-top absoluteGrob[strip.absoluteGrob.207] 10 9 ( 5- 5, 8- 8) strip-top absoluteGrob[strip.absoluteGrob.255] 11 10 ( 3- 3,10-10) strip-top absoluteGrob[strip.absoluteGrob.165] 12 11 ( 4- 4,10-10) strip-top absoluteGrob[strip.absoluteGrob.213] 13 12 ( 5- 5,10-10) strip-top absoluteGrob[strip.absoluteGrob.261] 14 13 ( 3- 3,12-12) strip-top absoluteGrob[strip.absoluteGrob.171] 15 14 ( 4- 4,12-12) strip-top absoluteGrob[strip.absoluteGrob.219] 16 15 ( 5- 5,12-12) strip-top absoluteGrob[strip.absoluteGrob.267] 17 16 ( 3- 3,14-14) strip-top absoluteGrob[strip.absoluteGrob.177] 18 17 ( 4- 4,14-14) strip-top absoluteGrob[strip.absoluteGrob.225] 19 18 ( 5- 5,14-14) strip-top absoluteGrob[strip.absoluteGrob.273] 20 19 ( 3- 3,16-16) strip-top absoluteGrob[strip.absoluteGrob.183] 21 20 ( 4- 4,16-16) strip-top absoluteGrob[strip.absoluteGrob.231] 22 21 ( 5- 5,16-16) strip-top absoluteGrob[strip.absoluteGrob.279] 23 22 ( 3- 3,18-18) strip-top absoluteGrob[strip.absoluteGrob.189] 24 23 ( 4- 4,18-18) strip-top absoluteGrob[strip.absoluteGrob.237] 25 24 ( 5- 5,18-18) strip-top absoluteGrob[strip.absoluteGrob.285]
ggplot2 2.2.0
pg
TableGrob (11 x 21) "layout": 42 grobs z cells name grob 28 2 ( 6- 6, 4- 4) strip-t-1 gtable[strip] 29 2 ( 6- 6, 6- 6) strip-t-2 gtable[strip] 30 2 ( 6- 6, 8- 8) strip-t-3 gtable[strip] 31 2 ( 6- 6,10-10) strip-t-4 gtable[strip] 32 2 ( 6- 6,12-12) strip-t-5 gtable[strip] 33 2 ( 6- 6,14-14) strip-t-6 gtable[strip] 34 2 ( 6- 6,16-16) strip-t-7 gtable[strip] 35 2 ( 6- 6,18-18) strip-t-8 gtable[strip]
Indeed, ggplot2 v2.2.0 constructs complex strips column by column, with each column a single grob. This can be checked by extracting one strip, then examining its structure. Using your plot:
library(ggplot2)
library(gtable)
library(grid)
# Your data
df = structure(list(location = structure(c(1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L), .Label = c("SF", "SS"), class = "factor"), species = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("AGR", "LKA"), class = "factor"),
position = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L), .Label = c("top", "bottom"), class = "factor"), density = c(0.41,
0.41, 0.43, 0.33, 0.35, 0.43, 0.34, 0.46, 0.32, 0.32, 0.4,
0.4, 0.45, 0.34, 0.39, 0.39, 0.31, 0.38, 0.48, 0.3, 0.42,
0.34, 0.35, 0.4, 0.38, 0.42, 0.36, 0.34, 0.46, 0.38, 0.36,
0.39, 0.38, 0.39, 0.39, 0.39, 0.36, 0.39, 0.51, 0.38)), .Names = c("location",
"species", "position", "density"), row.names = c(NA, -40L), class = "data.frame")
# Your ggplot with three facet levels
p=ggplot(df, aes("", density)) +
geom_boxplot(width=0.7, position=position_dodge(0.7)) +
theme_bw() +
facet_grid(. ~ species + location + position) +
theme(panel.spacing=unit(0,"lines"),
strip.background=element_rect(color="grey30", fill="grey90"),
panel.border=element_rect(color="grey90"),
axis.ticks.x=element_blank()) +
labs(x="")
# Get the ggplot grob
pg = ggplotGrob(p)
# Get the left most strip
index = which(pg$layout$name == "strip-t-1")
strip1 = pg$grobs[[index]]
# Draw the strip
grid.newpage()
grid.draw(strip1)
# Examine its layout
strip1$layout
gtable_show_layout(strip1)
One crude way to get outer strip labels 'spanning' inner labels is to construct the strip from scratch:
# Get the strips, as a list, from the original plot
strip = list()
for(i in 1:8) {
index = which(pg$layout$name == paste0("strip-t-",i))
strip[[i]] = pg$grobs[[index]]
}
# Construct gtable to contain the new strip
newStrip = gtable(widths = unit(rep(1, 8), "null"), heights = strip[[1]]$heights)
## Populate the gtable
# Top row
for(i in 1:2) {
newStrip = gtable_add_grob(newStrip, strip[[4*i-3]][1],
t = 1, l = 4*i-3, r = 4*i)
}
# Middle row
for(i in 1:4){
newStrip = gtable_add_grob(newStrip, strip[[2*i-1]][2],
t = 2, l = 2*i-1, r = 2*i)
}
# Bottom row
for(i in 1:8) {
newStrip = gtable_add_grob(newStrip, strip[[i]][3],
t = 3, l = i)
}
# Put the strip into the plot
# (It could be better to remove the original strip.
# In this case, with a coloured background, it doesn't matter)
pgNew = gtable_add_grob(pg, newStrip, t = 6, l = 4, r = 18)
# Draw the plot
grid.newpage()
grid.draw(pgNew)
OR using vectorised gtable_add_grob (see the comments):
pg = ggplotGrob(p)
# Get a list of strips from the original plot
strip = lapply(grep("strip-t", pg$layout$name), function(x) {pg$grobs[[x]]})
# Construct gtable to contain the new strip
newStrip = gtable(widths = unit(rep(1, 8), "null"), heights = strip[[1]]$heights)
## Populate the gtable
# Top row
cols = seq(1, by = 4, length.out = 2)
newStrip = gtable_add_grob(newStrip, lapply(strip[cols], `[`, 1), t = 1, l = cols, r = cols + 3)
# Middle row
cols = seq(1, by = 2, length.out = 4)
newStrip = gtable_add_grob(newStrip, lapply(strip[cols], `[`, 2), t = 2, l = cols, r = cols + 1)
# Bottom row
newStrip = gtable_add_grob(newStrip, lapply(strip, `[`, 3), t = 3, l = 1:8)
# Put the strip into the plot
pgNew = gtable_add_grob(pg, newStrip, t = 6, l = 4, r = 18)
# Draw the plot
grid.newpage()
grid.draw(pgNew)
这篇关于寻找由ggplot 2.2.0打破的gtable_add_grob代码的解决方法的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！
