具有第二个 x 轴和第二个 y 轴的格子图?

Lattice plot with both a second x- and a second y-axis?

我想在下面的(格子)水平图中添加第二个 y 轴(右)和第二个 x 轴(上)。这些轴应该只指示某些行和列(无标签),因此模仿基本图形的地毯功能。如何才能做到这一点?

library(lattice)
library(latticeExtra)

## Generate a correlation matrix
d <- 50
L <- diag(1:d)
set.seed(271)
L[lower.tri(L)] <- runif(choose(d,2))
Sigma <- L %*% t(L)
P <- cor(Sigma)

## Panel function
my_panel <- function(...) {
    panel.levelplot(...)
    panel.abline(h = (1:4)*10, v = (1:4)*10, lty = 2)
    panel.axis(side = "top", at = (1:50)-0.5, draw.labels = FALSE) # maybe a panel axis could do it? why not centered?
}

## Plot
obj1 <- levelplot(P, xlab = "Column", ylab = "Row",
                  col.regions = grey(c(seq(1, 0, length.out = 600))),
                  panel = my_panel)
obj2 <- xyplot(NA~NA, ylim = c(0, 50),
               scales = list(x = list(at = (1:50)-0.5, labels = rep("", 50)),
                             y = list(at = (1:50)-0.5, labels = rep("", 50))))
doubleYScale(obj1, obj2, use.style = FALSE) # idea based on latticeExtra; only gives a 2nd y-axis, though

您对 panel.rug() 有了一个好主意,但由于 lattice 将其绘图默认裁剪到面板内部而受阻。要解决这个问题,您可以通过 par.settings= 参数关闭裁剪。如果你想禁止在右侧和顶部面板边框上绘制默认轴刻度线,你可以使用 tck= 参数来实现,如下所示。

my_panel <- function(...) {
    panel.levelplot(...)
    panel.abline(h = (1:4)*10, v = (1:4)*10, lty = 2)
    ## Have panel.rug print tick marks starting at 1 npc (edge of panel)
    ## and extending to 1.02 npc (slightly outside of panel). (See ?unit)
    panel.rug(x = (1:51)-0.5, y = (1:51)-0.5, 
              start = 1, end = 1.02,
              col="black")
}

levelplot(P, xlab = "Column", ylab = "Row",
          col.regions = grey(c(seq(1, 0, length.out = 600))),
          ## Suppress default scales on right and top sides, by setting their
          ## tick lengths to zero 
          scales = list(tck=c(1,0)),
          ## Turn off clipping, so that panel.rug can plot outside of the panel
          par.settings = list(clip = list(panel = "off")),          
          panel = my_panel)