#olá pessoal, #Estou querendo montar uma função que me mostre pontos que estão á uma distância menor que uma pré-estabelecida. #o script e explicações seguem abaixo: ############################ ##Conjunto de dados para testar a<-rep(c(1,4),each=4) b<-c(1,2,4,5,1,2,4,5) cc<-c("p1","p2","p3","p4","p5","p6","p7","p8") mat<-data.frame(cc,a,b) rownames(mat)<-cc coord<-mat[,c(2,3)] plot(coord,pch="",xlim=c(0,6),ylim=c(0,6));text(coord[,"a"],coord[,"b"],rownames(coord))#esse é o meu conjunto de pontos. #note que se eu for juntar os pontos que estão a uma distância menor ou igual á dois, formarei dois grupos. #um formado pelos pontos: p1, p2, p3 e p4 e, outro grupo com os pontos p5, p6, p7 e p8 #################### ##iniciando os cálculos de distância library(vegan) (mat.dist<-vegdist(coord,"euclidean"))#calculo de distância entre os pontos resu.col1<-numeric(0) resu.col2<-numeric(0) (distancias<-as.numeric(mat.dist)) locais<-rownames(coord) n.linhas<-nrow(coord) ##Agora vou criar uma tabela para comparar as distâncias par-a-par for(i in 2:n.linhas){ linhas<-c(i:n.linhas) resu.col2<-c(resu.col2,locais[linhas]) resu.col1<-c(resu.col1,rep(locais[i-1],length(linhas))) } (mat.resu<-data.frame(resu.col1,resu.col2,distancias))#tabela com as distâncias par-a-par dist.corte<-2##argumentos usado para estabelecer distância para juntar pontos (min.dists<-mat.resu[which(mat.resu[,"distancias"]<=dist.corte),])##tabela com todos os pontos que estão a uma distância menor, entre si, que á pré estabelecida. #note que se for juntando os pontos forma-se dois grupos # resu.col1 resu.col2 distancias *#1 p1 p2 1#p1 com p2* *#8 p2 p3 2#o p2 que estava juntado com o p1 junta-se ao p3* *#14 p3 p4 1**#o p3 que estava juntado com o p2 junta-se ao p4* *#23 p5 p6 1#idem* *#26 p6 p7 2* *#28 p7 p8 1* ######################################## *#Alguém, por favor, sabe me dizer como de agora em diante eu faço para saber quais que são os pontos que formam cada um dos grupos?* Atenciosamente, Luciano -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732