2次元PLOT(散布図)に等高線を重ねる

plotした後に、contour関数をadd=TRUEとして使えばOK。contour関数に食わせるデータはouter関数で作るヨロシ。

x <- seq(-4, 4, 0.1)
y <- x
z <- outer(x, y, function(x,y){exp(-0.5*(x^2+y^2))})
plot(rnorm(200), rnorm(200), xlim=c(-4, 4), ylim=c(-4, 4))
contour(x, y, z, add=TRUE)