我有一个data.frame xy
,我正在绘制像下面的代码。R:固定的x轴缩放和转换图()xy参数到点()
有没有一种方法,我可以怎样我xx
和在xy走向(yy
)转换为points()
命令,这样我可以设置type='n'
并添加points()
段,以便命令来控制它更完善?
xy <- data.frame(NAME=c("NAME1","NAME1","NAME1","NAME2","NAME2","NAME2"),ID=c(87,87,87,199,199,199), X_START_YEAR=c(1984,1986,1984,1899,1909,1924),Y_START_VALUE=c(75,25,-90,-8,-55,-10),X_END_YEAR=c(1986,1994,1999,1909,1924,1927), Y_END_VALUE=c(20,50,-15,-70,-80,-100))
xy
NAME ID X_START_YEAR Y_START_VALUE X_END_YEAR Y_END_VALUE
1 NAME1 87 1984 75 1986 20
2 NAME1 87 1986 25 1994 50
3 NAME1 87 1984 -90 1999 -15
4 NAME2 199 1899 -8 1909 -70
5 NAME2 199 1909 -55 1924 -80
6 NAME2 199 1924 -10 1927 -100
ind <- split(xy,xy$ID)
for (x in ind){
xx = unlist(x[,grep('X_',colnames(x))])
yy = unlist(x[,grep('Y_',colnames(x))])
fname <- paste0(x[1, 'ID'], '.png')
png(fname, width=1679, height=1165, res=150)
par(mar=c(6,8,6,5))
plot(xx,
yy,
main=unique(x[,1]),
xlab="Time [Years]",
ylab="Value [m]")
axis(1, at = seq(1000, 2050, 5), cex.axis=1, labels=FALSE, tcl=-0.3)
axis(2, at = seq(-100000, 100000, 500), cex.axis=1, labels=FALSE, tcl=-0.3)
x <- x[,-1]
segments(x[,2],x[,3],x[,4],x[,5],lwd=2)
dev.off()
}
如果有可能这将是巨大的,如果x轴可以是在一个固定的范围内(例如2040至14年),如果1940年之前的值是存在于x轴应该是自动的。 y轴的范围始终不同。我怎样才能将其纳入我的代码?
感谢您的回答!这是做我想做的,完美!一个小细节我遇到了问题,也许你有一个解决方案:我想用'pch = 21'和'bg ='white''来绘制点,从而产生一个带有黑边的白色圆点。然而,现在绘制的点被称为第一个,然后段被绘制在点上。这看起来不如在具有黑色边缘的白色点绘制在线段上时更好。这就是为什么我最初正在寻找一个带有“()”的解决方案,以便我可以定义要绘制的特征的顺序... – kurdtc 2014-09-22 14:12:22
我添加了带点的行......我希望它适用于您。 – 2014-09-22 14:57:03
谢谢!没有看到它是如此容易纳入。你帮了我很多! – kurdtc 2014-09-22 15:04:52