如何使用R制作以下图？

Question

我想用R中的4个轴绘制一个图，使其看起来与此图相似：

我查看了Quick-R website的建议并修改了他们的一个示例（称为A Silly Axis Example）：

    # specify the data
x <- c(1:5); y <- x/2; 
w <- c(2:4)
z <- c(1:5)

# create extra margin room on the right for an axis
par(mar=c(5, 4, 4, 8) + 0.1)

# plot x vs. y
plot(x, y,type="b", pch=21, col="red",
   yaxt="n", lty=3, xlab="", ylab="")

# add x vs. 1/x
lines(x, z, type="b", pch=22, col="blue", lty=2)

# draw an axis on the left
axis(2, at=x,labels=x, col.axis="red", las=2)

# draw an axis on the right, with smaller text and ticks
axis(4, at=w,labels=round(w,digits=2),
  col.axis="blue", las=2, cex.axis=0.7, tck=-.01)

# draw an axis on the top
axis(3, at=z,labels=round(z,digits=2),
  col.axis="blue", las=2, cex.axis=0.7, tck=-.01)

# add a title for the right axis
mtext("L", side=3, line=3, cex.lab=1,las=2, col="blue")

# add a title for the right axis
mtext("OSR", side=4, line=3, cex.lab=1,las=2, col="red")

# add a main title and bottom and left axis labels
title("", xlab="GSI",   ylab="FSI")

此代码生成以下图表：

我很难弄清楚不同的轴如何有不同的比例。例如，我希望顶轴L从5到13，但如果我设置z <-c(5:13)，它将不会将轴设置为这些值。但是，我可以覆盖标签的内容：

axis(3, at=z,labels=round(c(9:13),digits=2), col.axis="blue", 
las=2, cex.axis=0.7, tck=-.01)

但是如果我想用这四个参数绘制一个点，该点将不会显示在正确的位置。我该怎么做？

Answer 1

一个（也许是繁琐的）选项是编写转换函数，在两个尺度之间转换值。假设您提前知道了顶轴和底轴的数据范围，您可以编写如下函数：

convertScaleToBottom <- function(x,botRange,topRange){
    temp <- (x - topRange[1]) / (topRange[2] - topRange[1])
    return(botRange[1] + (temp * (botRange[2] - botRange[1])))
}

在顶轴刻度中获取一组值x，并将它们转换为底轴刻度。然后，您可以绘制转换后的值并将原件保留为标签：

z1 <- 5:13
z1Adj <- convertScaleToBottom(z1,range(x),range(z1))

# draw an axis on the top
axis(3, at=z1Adj,labels=round(z1,digits=2),
     col.axis="blue", las=2, cex.axis=0.7, tck=-.01)

可以轻松修改此方法以反转顶轴的顺序：

convertScaleToBottomRev <- function(x,botRange,topRange){
    temp <- (x - topRange[1]) / (topRange[2] - topRange[1])
    return(botRange[2] - (temp * (botRange[2] - botRange[1])))
}