Question

我正在通过以下方式将ID定义的几个图表绘制到不同的文件中。

xy <- data.frame(NAME=c("NAME1","NAME1","NAME1","NAME2","NAME2","NAME2","NAME3","NAME3"),ID=c(87,87,87,199,199,199,233,233), X_START_YEAR=c(1950,1988,1994,1899,1909,1924,1945,1948),Y_START_VALUE=c(75,25,-90,-8,-55,-10,-9,12),X_END_YEAR=c(1985,1994,1999,1904,1924,1987,1946,1949), Y_END_VALUE=c(20,50,-15,-70,-80,-100,24,59))

  NAME  ID X_START_YEAR Y_START_VALUE X_END_YEAR Y_END_VALUE
1 NAME1  87         1950            75       1985          20
2 NAME1  87         1988            25       1994          50
3 NAME1  87         1994           -90       1999         -15
4 NAME2 199         1899            -8       1904         -70
5 NAME2 199         1909           -55       1924         -80
6 NAME2 199         1924           -10       1987        -100
7 NAME3 233         1945            -9       1946          24
8 NAME3 233         1948            12       1949          59

# Split glacier names by ID
ind <- split(x = xy,f = xy[,'ID'])

# Plot
for (i in ind){
  xx = unlist(i[,grep('X_',colnames(i))])
  yy = unlist(i[,grep('Y_',colnames(i))])    
  filename <- paste0(i[1, 'ID'], '.png') 
  png(filename, width=1679, height=1165, res=150)
  if(any(xx < 1946)) {my_x_lim <- c(min(xx), 2014)} else {my_x_lim <- c(1946, 2014)} 
  par(mar=c(6,8,6,5))
  plot(xx, yy)
  i <- i[,-1]
  segments(i[,2],i[,3],i[,4],i[,5],lwd=2)
  points(xx, yy, pch=21,cex=1, bg='white') 
  dev.off()
}

我的问题：我正在寻找的是一个基本的R解决方案，将xy的整个数据包含在一个图中（而不是像上面代码那样的几个图）并且图例中的多个线和点（由ID定义的组）。我想通过ggplot找到了一种非常简单的方法，但我感兴趣的是如何在基数R中做到这一点。任何想法？

以下ggplot代码可能会让您了解我尝试做的事情：

library(ggplot2)
p = ggplot(xy)
p = p + geom_segment(aes(x = X_START_YEAR, xend = X_END_YEAR, y = Y_START_VALUE, yend = Y_END_VALUE,linetype=NAME))
print(p)

Answer 1

noms<-levels(xy$NAME)
nblev<-nlevels(xy$NAME)
linetype<-(1:9)[1:nblev]

par(mar=c(5,4,4,10),xpd=TRUE)
plot(c(xy[,3],xy[,5]),c(xy[,4],xy[,6]),type="p",pch=21,xlab="YEAR",ylab="VALUE",las=1)
segments(xy[,3],xy[,4],xy[,5],xy[,6],lty=linetype[as.numeric(xy$NAME)])
coord<-par("usr")
legend(coord[2]+10,sum(coord[3:4])/2,lty=linetype,legend=noms,bty="n")

R：在基础R图中联合多个段（）和点（）

1 个答案: