ggpairs的自定义组平均函数
根据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)
基于@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)
由上面的代码产生的剧情
链接地址: http://www.djcxy.com/p/30921.html