如何使用 2 个不同的 y 轴进行绘图？

我想在 R 中叠加两个散点图，以便每组点都有自己的（不同的）y 轴（即在图上的位置 2 和 4 中），但这些点看起来叠加在同一个图上。

是否可以这样做plot?

Edit显示问题的示例代码

# example code for SO question
y1 <- rnorm(10, 100, 20)
y2 <- rnorm(10, 1, 1)
x <- 1:10
# in this plot y2 is plotted on what is clearly an inappropriate scale
plot(y1 ~ x, ylim = c(-1, 150))
points(y2 ~ x, pch = 2)

update：复制 R wiki 上的材料http://rwiki.sciviews.org/doku.php?id=tips:graphics-base:2yaxes http://rwiki.sciviews.org/doku.php?id=tips:graphics-base:2yaxes，链接现已损坏：也可从回程机器 https://web.archive.org/web/20140118053202/http://rwiki.sciviews.org/doku.php?id=tips:graphics-base:2yaxes

同一图上有两个不同的 y 轴

（部分材料最初由 Daniel Rajdl 2006/03/31 15:26 提供）

请注意，很少有情况适合在同一绘图上使用两种不同的比例。很容易误导图形的观看者。检查以下两个示例以及对此问题的评论（example1 http://junkcharts.typepad.com/junk_charts/2006/06/illusion_of_suc.html, example2 http://junkcharts.typepad.com/junk_charts/2006/05/the_crossover_l.html from 垃圾图表 http://junkcharts.typepad.com/junk_charts/）， 也斯蒂芬·费斯的这篇文章 http://www.perceptualedge.com/articles/visual_business_intelligence/dual-scaled_axes.pdf（其结论是“我当然不能一劳永逸地得出结论，具有双尺度轴的图表永远没有用；只是我无法考虑到其他更好的解决方案来保证它们的情况。”）另见第#点4进这个卡通 http://www.smbc-comics.com/?id=3167#comic ...

如果您有决心，基本方法是创建您的第一个图，设置par(new=TRUE)为了防止 R 清除图形设备，创建第二个图axes=FALSE（并设置xlab and ylab留空——ann=FALSE也应该有效），然后使用axis(side=4)在右侧添加一个新轴，以及mtext(...,side=4)在右侧添加轴标签。这是一个使用一点虚构数据的示例：

set.seed(101)
x <- 1:10
y <- rnorm(10)
## second data set on a very different scale
z <- runif(10, min=1000, max=10000) 
par(mar = c(5, 4, 4, 4) + 0.3)  # Leave space for z axis
plot(x, y) # first plot
par(new = TRUE)
plot(x, z, type = "l", axes = FALSE, bty = "n", xlab = "", ylab = "")
axis(side=4, at = pretty(range(z)))
mtext("z", side=4, line=3)

twoord.plot() in the plotrix包自动化了这个过程，就像doubleYScale() in the latticeExtra包裹。

另一个例子（改编自 Robert W. Baer 的 R 邮件列表帖子）：

## set up some fake test data
time <- seq(0,72,12)
betagal.abs <- c(0.05,0.18,0.25,0.31,0.32,0.34,0.35)
cell.density <- c(0,1000,2000,3000,4000,5000,6000)

## add extra space to right margin of plot within frame
par(mar=c(5, 4, 4, 6) + 0.1)

## Plot first set of data and draw its axis
plot(time, betagal.abs, pch=16, axes=FALSE, ylim=c(0,1), xlab="", ylab="", 
   type="b",col="black", main="Mike's test data")
axis(2, ylim=c(0,1),col="black",las=1)  ## las=1 makes horizontal labels
mtext("Beta Gal Absorbance",side=2,line=2.5)
box()

## Allow a second plot on the same graph
par(new=TRUE)

## Plot the second plot and put axis scale on right
plot(time, cell.density, pch=15,  xlab="", ylab="", ylim=c(0,7000), 
    axes=FALSE, type="b", col="red")
## a little farther out (line=4) to make room for labels
mtext("Cell Density",side=4,col="red",line=4) 
axis(4, ylim=c(0,7000), col="red",col.axis="red",las=1)

## Draw the time axis
axis(1,pretty(range(time),10))
mtext("Time (Hours)",side=1,col="black",line=2.5)  

## Add Legend
legend("topleft",legend=c("Beta Gal","Cell Density"),
  text.col=c("black","red"),pch=c(16,15),col=c("black","red"))

类似的方法可用于叠加不同类型的图——条形图、直方图等。