如何在R的下面板和上面板中使用不同的数据集？

Question

想在R中创建一个散点图矩阵，所以我使用函数pairs，如下所示：

pairs(X, lower.panel=panel.contour1, gap=0.5)

如何将一组不同的数据（我们称之为Y）引入panel.contour1函数？基本上我需要为矩阵图的上三角和下三角使用两组不同的数据。

Answer 1

也许有办法传递其他数据集，但如果失败，您可以rbind X和Y，然后参考lower.panel中的相关子集和upper.panel。

例如：

X <- matrix(runif(300), ncol=3)
Y <- matrix(c(sort(runif(100, 0, 10)), 
              sort(runif(100, 0, 10)),
              sort(runif(100, 0, 10))), ncol=3)
XY <- rbind(X, Y)

pairs(XY,
      lower.panel=function(x, y, ...) {
        Xx <- x[seq_len(nrow(X))] # corresponds to X subset
        Xy <- y[seq_len(nrow(X))] # corresponds to X subset
        usr <- par("usr"); on.exit(par(usr))
        par(usr = c(range(X[, -ncol(X)]), range(X[, -1]))) # set up limits
        points(Xx, Xy)
        if(par('mfg')[2] == 1) axis(2) # if left plot, add left axis
        if(par('mfg')[1] == ncol(X)) axis(1) # if bottom plot add bottom axis
      }, 
      upper.panel=function(x, y, ...) {
        Yx <- x[(nrow(X) + 1):length(x)] # Y subset
        Yy <- y[(nrow(X) + 1):length(y)] # Y subset
        cntr <- outer(Yx, Yx, FUN='*') # arbitrary function for contour
        usr <- par("usr"); on.exit(par(usr))
        par(usr = c(range(Y[, -1]), range(Y[, -ncol(Y)]))) # set up limits
        contour(Yx, Yy, cntr, add=TRUE)
        if(par('mfg')[2] == ncol(Y)) axis(4) # if right plot, add right axis
        if(par('mfg')[1] == 1) axis(3) # if top plot, add top axis
      }, 
      tick=FALSE, # suppress the default tick marks
      line=10) # move the default tick labels off the plot