ggplot2 geom_pointrange 中具有显着性注释的方面的不同颜色值
Different colored values by facet in ggplot2 geom_pointrange with significance annotation
我改编了我在 Google 组 post 中找到的代码："Annotations to display significant differences" 最初由 "Tom W"
编写
此代码生成以下点范围图：
我希望值（点 + 置信区间带）在方面 A 中为蓝色，在方面 B 中为红色。现在，代码允许我为每个层指定不同的颜色，但不能为每个方面指定不同的颜色。
我发现了相关的堆栈溢出 post，但 none 适用于此数据的结构方式。有人看到解决方案了吗？
提前致谢，
塔拉
library(ggplot2)
library(grid)
meanstable <-
structure(
list(
x_categories = structure(
c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L),
.Label = c("Facet A",
"Facet B"),
class = "factor"),
strata = structure(
c(1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L),
.Label = c("Strata 1", "Strata 2",
"Strata 3", "Strata 4" ),
class = "factor"),
mu = c(0.02, 0.004, 0.07, 0.12,
-0.02, 0.04, 0.04, 0.11),
lo = c(-0.001, -0.01, 0.03, 0.07,
-0.03, 0.02, -0.06, -0.01),
hi = c(0.04, 0.03, 0.16, 0.23,
0.01, 0.07, 0.09, 0.20)),
.Names = c("x_categories", "strata", "mu", "lo", "hi" ),
row.names = c(NA, 8L), class = "data.frame")
segdf <-
structure(
list(
x = c(1, 1, 2, 2, 2, 3, 1, 1, 4),
y = c(0.05, 0.05, 0.05, 0.17, 0.17, 0.17, 0.24, 0.24, 0.24),
xend = c(1, 2, 2, 2, 3, 3, 1, 4, 4),
yend = c(0.045, 0.05, 0.045, 0.165, 0.17, 0.165, 0.235, 0.24, 0.235),
x_categories = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L),
class = "factor",
.Label = "Facet A")),
.Names = c("x", "y", "xend", "yend", "x_categories"),
row.names = c(NA, -9L), class = "data.frame")
segdf2 <-
structure(
list(
x = c(1, 1, 2, 1, 1, 4 ),
y = c(0.08, 0.08, 0.08, 0.21, 0.21, 0.21),
xend = c(1, 2, 2, 1, 4, 4),
yend = c(0.075, 0.08, 0.075, 0.205, 0.21, 0.205),
x_categories = structure(c(1L, 1L, 1L, 1L, 1L, 1L),
class = "factor",
.Label = "Facet B")),
.Names = c("x", "y", "xend", "yend", "x_categories"),
row.names = c(NA, -6), class = "data.frame")
anodf <-
structure(
list(
x = c(1.5, 2.5, 2.5),
y = c(0.055, 0.175, 0.245),
x_categories = structure(c(1L, 1L, 1L),
class = "factor",
.Label = "Facet A")),
.Names = c("x", "y", "x_categories"),
row.names = c(NA, -3L),
class = "data.frame")
anodf2 <-
structure(
list(
x = c(1.5, 2.5),
y = c(0.085, 0.215),
x_categories = structure(c(1L, 1L),
class = "factor",
.Label = "Facet B")),
.Names = c("x", "y", "x_categories"),
row.names = c(NA, -2L),
class = "data.frame")
ggplot(meanstable) +
geom_hline(yintercept=0, linetype="dashed", colour="grey55") +
geom_pointrange(size = 1.2,
aes(x = strata, ymin = lo, ymax = hi, y = mu,
color = strata)) +
facet_wrap(~x_categories, nrow = 1) +
scale_color_manual(values=c("grey55","grey55", "grey55", "grey55")) +
scale_x_discrete("X Axis Label") +
scale_y_continuous("Y Axis Label") +
theme(legend.position = "none",
strip.text.x = element_text(size = rel(1.5)),
axis.title.y = element_text(vjust=1.4, size = rel(1.4)),
axis.title.x = element_text(vjust=-0.2, size = rel(1.4)),
axis.text = element_text(size = rel(1.1)),
plot.margin = unit(c(1, 1, 1, 1), "cm")) +
geom_segment(data = segdf, size = .8,
aes(x=x, y=y, xend=xend, yend=yend, x_categories = x_categories)) +
geom_text(data = anodf, aes(x=x, y=y, x_categories = x_categories),
label=c("*", "**", "*"), size = 8) +
geom_segment(data = segdf2, size = .8,
aes(x=x, y=y, xend=xend, yend=yend, x_categories = x_categories)) +
geom_text(data = anodf2, aes(x=x, y=y, x_categories = x_categories),
label=c("***", "**"), size = 8)
您希望颜色随着小平面而不是层而变化，因此在您想要的元素（点范围）的美学映射中更改它：
geom_pointrange(size = 1.2,
aes(x = strata, ymin = lo, ymax = hi, y = mu,
color = x_categories))
然后更改比例定义中的比例映射：
scale_colour_manual(values = c("Facet A" = "blue", "Facet B" = "red"))
把这些放在一起就可以得到
ggplot(meanstable) +
geom_hline(yintercept=0, linetype="dashed", colour="grey55") +
geom_pointrange(size = 1.2,
aes(x = strata, ymin = lo, ymax = hi, y = mu,
color = x_categories)) +
facet_wrap(~x_categories, nrow = 1) +
scale_colour_manual(values = c("Facet A" = "blue", "Facet B" = "red")) +
scale_x_discrete("X Axis Label") +
scale_y_continuous("Y Axis Label") +
theme(legend.position = "none",
strip.text.x = element_text(size = rel(1.5)),
axis.title.y = element_text(vjust=1.4, size = rel(1.4)),
axis.title.x = element_text(vjust=-0.2, size = rel(1.4)),
axis.text = element_text(size = rel(1.1)),
plot.margin = unit(c(1, 1, 1, 1), "cm")) +
geom_segment(data = segdf, size = .8,
aes(x=x, y=y, xend=xend, yend=yend, x_categories = x_categories)) +
geom_text(data = anodf, aes(x=x, y=y, x_categories = x_categories),
label=c("*", "**", "*"), size = 8) +
geom_segment(data = segdf2, size = .8,
aes(x=x, y=y, xend=xend, yend=yend, x_categories = x_categories)) +
geom_text(data = anodf2, aes(x=x, y=y, x_categories = x_categories),
label=c("***", "**"), size = 8)
我改编了我在 Google 组 post 中找到的代码："Annotations to display significant differences" 最初由 "Tom W"
编写此代码生成以下点范围图：
我希望值（点 + 置信区间带）在方面 A 中为蓝色，在方面 B 中为红色。现在，代码允许我为每个层指定不同的颜色，但不能为每个方面指定不同的颜色。
我发现了相关的堆栈溢出 post，但 none 适用于此数据的结构方式。有人看到解决方案了吗？
提前致谢， 塔拉
library(ggplot2)
library(grid)
meanstable <-
structure(
list(
x_categories = structure(
c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L),
.Label = c("Facet A",
"Facet B"),
class = "factor"),
strata = structure(
c(1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L),
.Label = c("Strata 1", "Strata 2",
"Strata 3", "Strata 4" ),
class = "factor"),
mu = c(0.02, 0.004, 0.07, 0.12,
-0.02, 0.04, 0.04, 0.11),
lo = c(-0.001, -0.01, 0.03, 0.07,
-0.03, 0.02, -0.06, -0.01),
hi = c(0.04, 0.03, 0.16, 0.23,
0.01, 0.07, 0.09, 0.20)),
.Names = c("x_categories", "strata", "mu", "lo", "hi" ),
row.names = c(NA, 8L), class = "data.frame")
segdf <-
structure(
list(
x = c(1, 1, 2, 2, 2, 3, 1, 1, 4),
y = c(0.05, 0.05, 0.05, 0.17, 0.17, 0.17, 0.24, 0.24, 0.24),
xend = c(1, 2, 2, 2, 3, 3, 1, 4, 4),
yend = c(0.045, 0.05, 0.045, 0.165, 0.17, 0.165, 0.235, 0.24, 0.235),
x_categories = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L),
class = "factor",
.Label = "Facet A")),
.Names = c("x", "y", "xend", "yend", "x_categories"),
row.names = c(NA, -9L), class = "data.frame")
segdf2 <-
structure(
list(
x = c(1, 1, 2, 1, 1, 4 ),
y = c(0.08, 0.08, 0.08, 0.21, 0.21, 0.21),
xend = c(1, 2, 2, 1, 4, 4),
yend = c(0.075, 0.08, 0.075, 0.205, 0.21, 0.205),
x_categories = structure(c(1L, 1L, 1L, 1L, 1L, 1L),
class = "factor",
.Label = "Facet B")),
.Names = c("x", "y", "xend", "yend", "x_categories"),
row.names = c(NA, -6), class = "data.frame")
anodf <-
structure(
list(
x = c(1.5, 2.5, 2.5),
y = c(0.055, 0.175, 0.245),
x_categories = structure(c(1L, 1L, 1L),
class = "factor",
.Label = "Facet A")),
.Names = c("x", "y", "x_categories"),
row.names = c(NA, -3L),
class = "data.frame")
anodf2 <-
structure(
list(
x = c(1.5, 2.5),
y = c(0.085, 0.215),
x_categories = structure(c(1L, 1L),
class = "factor",
.Label = "Facet B")),
.Names = c("x", "y", "x_categories"),
row.names = c(NA, -2L),
class = "data.frame")
ggplot(meanstable) +
geom_hline(yintercept=0, linetype="dashed", colour="grey55") +
geom_pointrange(size = 1.2,
aes(x = strata, ymin = lo, ymax = hi, y = mu,
color = strata)) +
facet_wrap(~x_categories, nrow = 1) +
scale_color_manual(values=c("grey55","grey55", "grey55", "grey55")) +
scale_x_discrete("X Axis Label") +
scale_y_continuous("Y Axis Label") +
theme(legend.position = "none",
strip.text.x = element_text(size = rel(1.5)),
axis.title.y = element_text(vjust=1.4, size = rel(1.4)),
axis.title.x = element_text(vjust=-0.2, size = rel(1.4)),
axis.text = element_text(size = rel(1.1)),
plot.margin = unit(c(1, 1, 1, 1), "cm")) +
geom_segment(data = segdf, size = .8,
aes(x=x, y=y, xend=xend, yend=yend, x_categories = x_categories)) +
geom_text(data = anodf, aes(x=x, y=y, x_categories = x_categories),
label=c("*", "**", "*"), size = 8) +
geom_segment(data = segdf2, size = .8,
aes(x=x, y=y, xend=xend, yend=yend, x_categories = x_categories)) +
geom_text(data = anodf2, aes(x=x, y=y, x_categories = x_categories),
label=c("***", "**"), size = 8)
您希望颜色随着小平面而不是层而变化，因此在您想要的元素（点范围）的美学映射中更改它：
geom_pointrange(size = 1.2,
aes(x = strata, ymin = lo, ymax = hi, y = mu,
color = x_categories))
然后更改比例定义中的比例映射：
scale_colour_manual(values = c("Facet A" = "blue", "Facet B" = "red"))
把这些放在一起就可以得到
ggplot(meanstable) +
geom_hline(yintercept=0, linetype="dashed", colour="grey55") +
geom_pointrange(size = 1.2,
aes(x = strata, ymin = lo, ymax = hi, y = mu,
color = x_categories)) +
facet_wrap(~x_categories, nrow = 1) +
scale_colour_manual(values = c("Facet A" = "blue", "Facet B" = "red")) +
scale_x_discrete("X Axis Label") +
scale_y_continuous("Y Axis Label") +
theme(legend.position = "none",
strip.text.x = element_text(size = rel(1.5)),
axis.title.y = element_text(vjust=1.4, size = rel(1.4)),
axis.title.x = element_text(vjust=-0.2, size = rel(1.4)),
axis.text = element_text(size = rel(1.1)),
plot.margin = unit(c(1, 1, 1, 1), "cm")) +
geom_segment(data = segdf, size = .8,
aes(x=x, y=y, xend=xend, yend=yend, x_categories = x_categories)) +
geom_text(data = anodf, aes(x=x, y=y, x_categories = x_categories),
label=c("*", "**", "*"), size = 8) +
geom_segment(data = segdf2, size = .8,
aes(x=x, y=y, xend=xend, yend=yend, x_categories = x_categories)) +
geom_text(data = anodf2, aes(x=x, y=y, x_categories = x_categories),
label=c("***", "**"), size = 8)