Создайте матрицу диаграмм рассеяния (пар() эквивалентных) в ggplot2

Можно ли построить матрицу графиков рассеяния с ggplot2, используя ggplot приятные функции, такие как отображение дополнительных факторов в цвет, форму и т.д. и добавление более плавного?

Я думаю о чем-то подобном функции base pairs.

Ответ 1

Возможно, вы захотите попробовать plotmatrix:

  library(ggplot2)
  data(mtcars)
  plotmatrix(mtcars[,1:3])

мне mpg (первый столбец в mtcars) не должен быть фактором. Я не проверял его, но нет причин, почему он должен быть одним. Однако я получаю график рассеяния:)

Примечание.. Для справки в будущем функция plotmatrix() была заменена функцией ggpairs() из пакета GGally, поскольку @naught101 предлагает в другом ответе ниже на этот вопрос.

Ответ 2

Я все время хочу этого делать, но сюжетница - это дерьмо. Hadley рекомендует вместо пакет GGally. Он имеет функцию ggpairs, которая представляет собой значительно улучшенный парный сюжет (позволяет использовать в ваших кадрах данные непрерывные переменные). Он отображает разные графики в каждом квадрате, в зависимости от типов переменных:

library(GGally)
ggpairs(iris, aes(colour = Species, alpha = 0.4))

Ответ 3

Если вы хотите получить объект ggplot (не ggmatrix, как в случае ggpairs()), решение состоит в том, чтобы дважды расплавить данные, а затем ggplot с фасеткой. facet_wrap было бы лучше, чем facet_grid в ограничении области графика, если задан параметр scales = 'free'.

require(ggplot2) 
require(dplyr)
require(tidyr)

gatherpairs <- function(data, ..., 
                        xkey = '.xkey', xvalue = '.xvalue',
                        ykey = '.ykey', yvalue = '.yvalue',
                        na.rm = FALSE, convert = FALSE, factor_key = FALSE) {
  vars <- quos(...)
  xkey <- enquo(xkey)
  xvalue <- enquo(xvalue)
  ykey <- enquo(ykey)
  yvalue <- enquo(yvalue)

  data %>% {
    cbind(gather(., key = !!xkey, value = !!xvalue, !!!vars,
                 na.rm = na.rm, convert = convert, factor_key = factor_key),
          select(., !!!vars)) 
  } %>% gather(., key = !!ykey, value = !!yvalue, !!!vars,
               na.rm = na.rm, convert = convert, factor_key = factor_key)
}

iris %>% 
  gatherpairs(Sepal.Length, Sepal.Width, Petal.Length, Petal.Width) %>% {
  ggplot(., aes(x = .xvalue, y = .yvalue, color = Species)) +
      geom_point() + 
      geom_smooth(method = 'lm') +
      facet_wrap(.xkey ~ .ykey, ncol = length(unique(.$.ykey)), scales = 'free', labeller = label_both) +
      scale_color_brewer(type = 'qual')
}