我正在以较大的置信区间绘制带有mean_cl_boot
的一些X值
如何导出每个组中fun.y = mean
和fun.data = mean_cl_boot
值的文本?
我在mean_cl_boot
中有一个值的间隔,我想绘制它们并导出它们。
ggplot(iris, aes(x = Species, y = Petal.Length)) +
geom_jitter(width = 0.5) + stat_summary(fun.y = mean, geom = "point", color = "red") +
stat_summary(fun.data = mean_cl_boot, fun.args=(conf.int=0.9999), geom = "errorbar", width = 0.4)
我必须绘制平均值(fun.y = mean
),并带有:
stat_summary(fun.y=mean, geom="text", aes(label=sprintf("%1.1f", ..y..)),size=3, show.legend=FALSE
但是我不能与mean_cl_boot
相同。
您可以通过stat_summary
访问ggplot_build
的数据。
首先,将ggplot调用存储在对象中:
g <- ggplot(iris, aes(x = Species, y = Petal.Length)) +
geom_jitter(width = 0.5) +
stat_summary(fun.y = mean, geom = "point", color = "red") +
stat_summary(fun.data = mean_cl_boot, fun.args=(conf.int=0.9999), geom = "errorbar", width = 0.4)
然后,加上:
ggplot_build(g)$data[[3]]
您得到用mean_cl_boot
计算的值:
x group y ymin ymax PANEL xmin xmax colour size linetype width alpha 1 1 1 1.462 1.386000 1.543501 1 0.8 1.2 black 0.5 1 0.4 NA 2 2 2 4.260 4.024899 4.462202 1 1.8 2.2 black 0.5 1 0.4 NA 3 3 3 5.552 5.337199 5.798202 1 2.8 3.2 black 0.5 1 0.4 NA
为了正确设置标签,您可以这样做:
# extract the data
mcb <- ggplot_build(g)$data[[3]]
# add the labels to the plot
g + geom_text(data = mcb,
aes(x = group, y = ymin, label = round(ymin,2)),
color = "blue",
vjust = 1)
结果:
但是可能更好的选择是使用ggrepel包:
library(ggrepel)
g + geom_label_repel(data = mcb,
aes(x = group, y = ymin, label = round(ymin,2)),
color = "blue",
nudge_x = 0.2,
nudge_y = -0.2)
其结果: