将矩阵与 qq-plots 配对

Question

能够查看您的数据很有帮助。当您有多个变量时，您可以形成散点图矩阵，例如，pairs()。散点图矩阵为您提供一组数据的 2D 边缘投影。

set.seed(8092)
X <- matrix(rnorm(80), ncol=4)
pairs(X)

您还可以拥有不同组的数据并希望比较它们的分布。可以将两个这样的分布与 qq-plot 进行比较。

set.seed(4415)
group1 <- rnorm(20)
group2 <- rnorm(20)
qqplot(group1, group2)
abline(c(0,1))

当您有多个组时，如果有一个显示 qq-plots 矩阵的对类型图会很方便。

colnames(X) <- c("group1", "group2", "group3", "group4")
qq.pairs(X)

有这样的功能吗？有没有一种直接的方法可以从头开始编码？

Answer 1

您可以通过*panelargs添加您喜欢的任何功能。例如：

set.seed(8092)
X <- matrix(rnorm(80), ncol=4)

panel.qq <- function(x, y, ...) {
  usr <- par("usr"); on.exit(par(usr))
  par(usr = c(0, 1, 0, 1), new = TRUE)
  qqplot(x, y, xlab = deparse(substitute(x)), ylab = deparse(substitute(y)))
  abline(c(0,1), ...)
}

pairs(X, lower.panel = panel.qq)

当然你可以panel=panel.qq用来制作所有的方块qq图