寻找 plot 成对散点图矩阵的方法，其中变量分为两组

Question

寻找一种方法来 plot 成对散点图矩阵，其中变量位于两组中，即图需要限制为组之间的对，不包括组内对。 例如，碱基pairs()产生这个

    pairs(iris[1:4])

例如，我希望仅将萼片变量与花瓣变量相关联，并且对萼片（或花瓣）长度与宽度的相关性不感兴趣。 所以，我需要一个 2x2 矩阵，而不是 4x4 矩阵。

解决方案看起来像这样

pairs()似乎支持公式参数，但是pairs(Sepal.Length + Sepal.Width ~ Petal.Length + Petal.Width, data=iris)看起来不正确。 它给出了一个 3x3 矩阵的正确方向，但由于某种原因，这两个 sepal 变量被分组。

最好在ggplot2或GGally::ggpairs()中寻找解决方案。

Answer 1

这是一个带有 tidyr 预处理的香草 ggplot 方法。

首先，我们将重塑数据，使宽度和长度变量具有相同的列。

library(ggplot2)

df <- tidyr::pivot_longer(iris, c(Petal.Length, Petal.Width), 
                          names_to = "petal_metric", values_to = "petal_value")
df <- tidyr::pivot_longer(df, c(Sepal.Length, Sepal.Width),
                          names_to = "sepal_metric", values_to = "sepal_value")

接下来，我们可以简单地分析花瓣/萼片指标。

p <- ggplot(df, aes(sepal_value, petal_value, colour = Species)) +
  geom_point()
  
p + facet_grid(petal_metric ~ sepal_metric)

我们可以更好地装饰 plot 以将条带用作轴标题。

p + facet_grid(petal_metric ~ sepal_metric, switch = "both") +
  theme(strip.placement = "outside",
        strip.background = element_blank(),
        axis.title = element_blank())

^{由代表 package (v2.0.1) 于 2022 年 1 月 11 日创建}