Crear una matriz de diagramas de dispersión (pares() equivalentes) en ggplot2

Question

¿Es posible trazar una matriz de diagramas de dispersión con ggplot2, usando las características agradables de ggplot como mapear factores adicionales al color, forma, etc. y agregando suavizado?

Estoy pensando en algo similar a la función basepairs.

Answer 1

Es posible que desee probar plotmatrix:

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

me mpg (primera columna de mtcars) no debe ser un factor. No lo he comprobado, pero no hay ninguna razón por la que debería ser uno. Sin embargo consigo un gráfico de dispersión :)

---

Nota: Para referencia futura, la función plotmatrix() ha sido sustituida por la función ggpairs() del paquete GGally como @ naught101 sugiere in another response below a esta pregunta.

Answer 2

Sigo queriendo hacer esto, pero plotmatrix es una porquería. Hadley recommends utilizando el GGally package en su lugar. Tiene una función, ggpairs que es una parcela de pares muy mejorada (le permite usar variables no continuas en sus marcos de datos). Traza diferentes parcelas en cada cuadrado, dependiendo de los tipos de variables:

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

Answer 3

Si se desea obtener un objeto ggplot (no ggmatrix como en el caso de ggpairs()), la solución consiste en fundir el datos dos veces, luego ggplot con facetado. facet_wrap sería mejor que facet_grid al limitar el área trazada, dado que se proporciona el parámetro 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') 
}