r - 将矩阵与 qq 图配对

Question

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

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

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

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

当你有几个组时，如果有一个显示 qq-plots 矩阵的 pairs-type plot 会很方便。

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

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

Answer 1

您可以通过 *panel 参数添加您喜欢的任何功能。例如:

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来制作所有方 block 的qq图