#Estou querendo montar uma função que me mostre pontos que estão á uma distância menor que uma pré-estabelecida.
##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. SgarbiMestrando em Ecologia e Evolução - UFG
Laboratório
de Ecologia de InsetosCel.
(62)8174-2262 Lab. (62)3521-1732