Tạo ma trận biểu đồ phân tán (cặp () tương đương) trong ggplot2

Có thể vẽ một ma trận gồm các biểu đồ phân tán ggplot2, bằng cách sử dụng ggplotcác tính năng tốt đẹp như ánh xạ các yếu tố bổ sung với màu sắc, hình dạng, v.v. và thêm mượt mà?

Tôi đang nghĩ về một cái gì đó tương tự như basechức năng pairs.

Bạn có thể muốn thử plotmatrix:

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

đối với tôi mpg (cột đầu tiên trong mtcars) không nên là một yếu tố. Tôi đã không kiểm tra nó, nhưng không có lý do tại sao nó phải là một. Tuy nhiên tôi nhận được một âm mưu phân tán :)

Lưu ý: Để tham khảo trong tương lai, plotmatrix()hàm đã được thay thế bằng ggpairs()hàm từ GGallygói như @ naught101 gợi ý trong một câu trả lời khác bên dưới cho câu hỏi này.

Tôi tiếp tục muốn làm điều này, nhưng plotmatrix là tào lao. Hadley khuyên bạn nên sử dụng gói GGally thay thế. Nó có một hàm, ggpairs là một biểu đồ cặp được cải tiến rất nhiều (cho phép bạn sử dụng các biến không liên tục trong khung dữ liệu của mình). Nó vẽ các ô khác nhau trong mỗi ô vuông, tùy thuộc vào các loại biến:

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

Nếu một người muốn có được một ggplotđối tượng (không phải ggmatrixnhư trong trường hợp của ggpairs()), giải pháp là làm tan chảy dữ liệu hai lần, sau đó ggplotvới khía cạnh. facet_wrapsẽ tốt hơn là facet_gridgiới hạn diện tích được vẽ, với scales = 'free'tham số được cung cấp.

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')
}