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

上一篇: Custom group mean function for ggpairs

下一篇: Coloring points based on variable with R ggpairs