Skip to content
Prev 139999 / 398502 Next

minimum distances

Dave Depew <ddepew <at> sciborg.uwaterloo.ca> writes:
Sewers<-data.frame(x=c(10,20,35,50),y=c(100,150,200,300))
transect <- data.frame(x=0:50,y=100:150)


dist <- sqrt((outer(Sewers$x,transect$x,"-"))^2+
             (outer(Sewers$y,transect$y,"-"))^2)
mindist = apply(dist,2,min)
closest = apply(dist,2,which.min)

library(MASS)
eqscplot(Sewers,pch=16,col=1:4)
points(transect,col=closest)