R 轮廓水平不匹配 filled.contour

Question

希望这是一个简单的问题，但我使用 filled.contour() 在 R 中做了一个简单的图形。它看起来不错，以及它应该喜欢给定的数据。但是，我想沿着 0（级别 = 0）的等高线添加一条参考线，并且绘制的线与 filled.contour 图上的颜色不匹配。这条线很近，但与图形不匹配（并最终从 filled.contour 图中穿过另一个轮廓）。知道为什么会这样吗？

aa <- c(0.05843150,  0.11300040,  0.15280030,  0.183524400,  0.20772430,  0.228121000)
bb <- c(0.01561055,  0.06520635,  0.10196237,  0.130127650,  0.15314544,  0.172292410)
cc <- c(-0.02166599,  0.02306650,  0.05619421,  0.082193680,  0.10334837,  0.121156780)
dd <- c(-0.05356592, -0.01432910,  0.01546647,  0.039156660,  0.05858709,  0.074953650)
ee <- c(-0.08071987, -0.04654243, -0.02011676,  0.000977798,  0.01855881,  0.033651089)
ff <- c(-0.10343798, -0.07416114, -0.05111547, -0.032481132, -0.01683215, -0.003636035)
gg <- c(-0.12237798, -0.09753544, -0.07785126, -0.061607548, -0.04788856, -0.036169540)

hh <-rbind(aa,bb,cc,dd,ee,ff,gg)

z <- as.matrix(hh)
y <- seq(0.5,1.75,0.25)
x <- seq(1,2.5,0.25)

filled.contour(x,y,z,
           key.title = title(main=expression("log"(lambda))),
           color.palette = topo.colors) #This works

contour(x,y,z, level=0,add=T,lwd=3) #This line doesn't match plot

Answer 1

这在 ?filled.contour 帮助页面中得到了完整的回答。在注释部分它指出

The output produced by filled.contour is actually a combination of two plots; one is the filled contour and one is the legend. Two separate coordinate systems are set up for these two plots, but they are only used internally – once the function has returned these coordinate systems are lost. If you want to annotate the main contour plot, for example to add points, you can specify graphics commands in the plot.axes argument. See the examples.

并且该帮助页面中给出的示例显示了如何在主要情节之上进行注释。在这种特殊情况下，正确的方法是

filled.contour(x,y,z,
    key.title = title(main=expression("log"(lambda))),
    color.palette = topo.colors,
    plot.axes = {
        axis(1)
        axis(2)
        contour(x,y,z, level=0,add=T,lwd=3)             
    }
)

产生

R 轮廓水平不匹配 filled.contour

R contour levels don't match filled.contour

r

contour