r-ggpairs的自定义群均值函数

Question

根据GGally R包中ggpairs((函数的文档，可以指定自定义函数作为"lower"/"upper"参数的输入。对于连续离散变量组合，我想简单地显示分类变量类别中连续变量的平均值(最好使用点，而不是条(，如果可能的话，使用颜色美学通过另一个分类变量进一步分层。

我在以下线程中发现了一些信息：

https://github.com/ggobi/ggally/issues/218

然而，我对ggpairs(和ggplot2(的了解太肤浅，无法从线程中的模板中生成这样的自定义函数(此外，变量名"Species"似乎是硬编码到模板中的，如果可能的话，我宁愿在自定义函数中没有任何硬编码信息(。

如果有人能帮我提供解决方案的模板或草图，我将不胜感激，例如使用以下示例(其中"custom_function"需要用上述函数替换(：

dat <- reshape::tips
pm <- ggpairs(dat,
mapping = aes(color = sex, alpha = 0.3),
columns = c("total_bill", "smoker", "time", "tip"),
showStrips = T,
lower = list(combo = custom_function))
print(pm)

Answer 1

根据@aosmith的评论，我制作了一个自定义函数，它似乎对我的目的足够好，到目前为止还没有进行广泛的测试，但也许它无论如何都有帮助：

library(GGally)
library(ggplot2)
library(ggstance)
gmean_point <- function(data, mapping, ...) {
x <- eval(mapping$x, data)
y <- eval(mapping$y, data)
if(is.numeric(y)) {
p <- ggplot(data) +
geom_blank(mapping) +
stat_summary(mapping,
geom = 'point', fun.y = mean,
position = position_dodge(width = 0.2))
} else {
p <- ggplot(data) +
geom_blank(mapping) +
stat_summaryh(mapping,
geom = 'point', fun.x = mean,
position = position_dodgev(height = 0.2))
}
p
}
pm <- ggpairs(reshape::tips,
mapping = aes(color = sex, alpha = 0.3),
columns = c("total_bill", "smoker", "time", "tip"),
showStrips = T,
lower = list(combo = gmean_point),
upper = list(combo = 'box'))
print(pm)

上方代码生成的图