From: Milton Cezar Ribeiro <milton_ruser_at_yahoo.com.br>

Date: Fri, 8 Feb 2008 20:01:09 -0800 (PST)

df.0<-cbind(x0,y0)

df.1<-cbind(x1,y1)

plot(df.0,xlim=range(c(x0,x1)),ylim=range(c(y0,y1))) points(df.1,col=2)

rm(x0,x1,y0,y1)

R-help_at_r-project.org mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Sat 09 Feb 2008 - 04:05:16 GMT

Date: Fri, 8 Feb 2008 20:01:09 -0800 (PST)

x0<-rnorm(n=500,mean=1,sd=runif(1)) y0<-rnorm(n=500,mean=3,sd=runif(1)) x1<-rnorm(n=700,mean=8,sd=runif(1)) y1<-rnorm(n=700,mean=5,sd=runif(1))

df.0<-cbind(x0,y0)

df.1<-cbind(x1,y1)

plot(df.0,xlim=range(c(x0,x1)),ylim=range(c(y0,y1))) points(df.1,col=2)

rm(x0,x1,y0,y1)

#merge two data.frames of points

#### IT SPEND many time

df.merge<-merge(df.0,df.1,all=T)

#compute distances between each pair of points

attach(df.merge)

df.merge$distance<-((x0-x1)^2+(y0-y1)^2)^0.5
detach(df.merge)

#identify the minimum distance

df.merge.distance.min<-min(df.merge$distance)

#select the pair of points (x0,y0,x1,y1) with shortest distance

df.merge.distance.min.subset<-subset(df.merge,df.merge$distance<=df.merge.distance.min)

#trace a arrow between the points with shortest distance

arrows(df.merge.distance.min.subset[1,1],df.merge.distance.min.subset[1,2],df.merge.distance.min.subset[1,3],df.merge.distance.min.subset[1,4],code=0,col=3,lwd=2,lty=1)

Any help are welcome

Miltinho

Brazil.

para armazenamento!

[[alternative HTML version deleted]]

R-help_at_r-project.org mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Sat 09 Feb 2008 - 04:05:16 GMT

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.2.0, at Sat 09 Feb 2008 - 17:30:12 GMT.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*