我想用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)
但是如果我想用这四个参数绘制一个点,该点将不会显示在正确的位置。我该怎么做?
答案 0 :(得分:3)
一个(也许是繁琐的)选项是编写转换函数,在两个尺度之间转换值。假设您提前知道了顶轴和底轴的数据范围,您可以编写如下函数:
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])))
}