d <- read.table("out_MDS_map.csv",skip=15,sep=",") dat <- d[,2:ncol(d)] plot(dat[,2],dat[,1],pch=19,col="red",xlim=c(-3000,3000),ylim=c(-3000,3000)) par(xpd=NA) text(dat[,2],dat[,1],labels=d[,1],pos=2,cex=1.2) segments(dat[1,2],dat[1,1],dat[2,2],dat[2,1],col="red") segments(dat[1,2],dat[1,1],dat[3,2],dat[3,1],col="red") segments(dat[1,2],dat[1,1],dat[4,2],dat[4,1],col="red") segments(dat[1,2],dat[1,1],dat[5,2],dat[5,1],col="red") segments(dat[1,2],dat[1,1],dat[6,2],dat[6,1],col="red") segments(dat[1,2],dat[1,1],dat[7,2],dat[7,1],col="red") segments(dat[1,2],dat[1,1],dat[8,2],dat[8,1],col="red") segments(dat[1,2],dat[1,1],dat[9,2],dat[9,1],col="red")