对R中的两组使用散点图使用不同的颜色/形状

Question

我目前正在使用splom函数在R中使用点阵软件包来创建散点图矩阵。 我的数据集中有两个组，分别标记为两个不同的列：

PC1 PC2 PC3 Group1 Group2
1   2   3   A      X
1   2   3   B      X
1   2   3   C      X
1   2   3   D      X
1   2   3   A      Y
1   2   3   B      Y
1   2   3   C      Y
1   2   3   D      Y
1   2   3   A      Z
1   2   3   B      Z
1   2   3   C      Z
1   2   3   D      Z

使用以下代码，我可以使splom函数为组之一使用不同的颜色和形状，但不能同时为两个组使用：

splom(~pcVT[,1:3], data = pcVT, xlab = NULL, groups = Group1, pch = c(1,2,3),
col = super.sym$col[1:3], panel = panel.superpose, 
key = list(points = list(pch = c(1,2,3),col = super.sym$col[1:3]),text = list(mylabels)))

如何获得将这两个组用于着色和形状的方法-即，我希望Group1根据颜色绘制图形，而Group2根据形状绘制图形。 或者，如果splom无法执行此操作，是否有使用gpplot2的不错方法？

非常感谢！

Answer 1

您的示例不可复制，数据集也不太可用。 我又做了一个

这是带有基本图的解决方案：

d <- as.data.frame(princomp(iris[,1:4])$scores)
d$Group1 <- iris$Species
d$Group2 <- factor(sample(c("A","B","C"), 150, replace = TRUE))

mycols <- c("forestgreen", "gold", "dodgerblue")

x11(width = 16/2.54, height = 12/2.54)
pairs(d[,1:4], oma=c(3,3,6,3), 
      col = mycols[as.numeric(d$Group1)], pch = c(1:3)[as.numeric(d$Group2)], gap = 0)
legend("top", col = mycols, legend = levels(d$Group1), pch = 20, 
       xpd = NA, ncol = 3, bty = "n", inset = 0.01, pt.cex = 1.5)
legend("top", pch = 1:3, legend = levels(d$Group2), col = "black",
       xpd = NA, ncol = 3, bty = "n", inset = -0.03)

如果需要ggplot解决方案，请探索GGally:ggpairs可能性（计算速度要慢得多）：

library(GGally)
ggpairs(data=d, mapping = aes(color = Group1, shape = Group2), 
        columns = 1:4, legend = c(2,1))