[R-br] Coordenadas e distâncias

[Roberto Guevara]

Olá pessoal do R-br.

Meus dados consistem basicamente em coordenadas e o horário em que
cada ponto foi obtido, como exemplificado a seguir:

obs  dia  min    x    y
  1    1    400  425 475
  2    1    410  425 425
  3    1    420  475 325
  4    1    430  525 325
  5    1    435  525 275
  6    1    443  525 275

obs se refere a cada unidade de amostra, dia a cada dia de coleta (um
total de 73), min é o horário em unidades de minutos e x e y cada
elemento do ponto cartesiano.

O que preciso obter desses dados é a distribuição de distâncias
euclidianas entre um ponto escolhido aleatoriamente e todos os pontos
seguintes de um mesmo dia. Ou seja, preciso escolher aleatoriamente um
ponto por dia, calcular, para todos os dias, a euclidiana entre este
ponto e todos os pontos posteriores do mesmo dia e relacionar a
distância espacial com a diferença de tempo entre os pontos e o dia em
que cada distância está associada.

Eu já elaborei um código pra isso. Caso tenham sugestões para melhorar
o código, ficarei agradecido. Porém, não consigo passar para uma
próxima etapa. Como poderia fazer para repetir esse procedimento o
número de vezes que fosse necessário? Por exemplo, existem milhares de
combinações de pontos escolhidos aleatoriamente ao longo dos 73 dias.
Eu gostaria de repetir o procedimento 1000 vezes para ter uma ideia da
variação de distribuição de distâncias que posso obter em decorrência
da escolha aleatória de um ponto de origem por dia.

Código já elaborado:

arquivo<-read.table("arquivo.txt",head=T)

fixo<-numeric()	#vetor com origens de cada dia
ult<-numeric()	#vetor com últimos pontos de cada dia
xy<-list()	#lista para armazenar a "parte dos pontos de cada dia"
delimitadas por fixo e ult
hor<-list() #lista para armazenar a "parte dos horarios de cada dia"
delimitadas por fixo e ult
distancias<-list()#distancias euclidianas
tempo<-list()#diferença de tempo entre os pontos
dia<-list()#dia de coleta
for (i in 1:max(unique(arquivo$dia))){	#repetindo até o total de dias
	fixo[i]<-sample(arquivo$obs[arquivo$dia==i],1)#sorteio ponto de origem
	ult[i]<-max(arquivo$obs[arquivo$dia==i])#armazenando último ponto de cada dia
	xy[[i]]<-arquivo[fixo[i]:ult[i],4:5]#armazenando parte dos pontos de cada dia
	hor[[i]]<-arquivo[fixo[i]:ult[i],6]#armazenando parte dos horários de cada dia
	distancias[[i]]<-dist(xy[[i]])[1:nrow(xy[[i]])-1]#somente as
distâncias em relação ao primeiro ponto interessam
	tempo[[i]]<-dist(horarios[[i]])[1:length(horarios[[i]])-1]#somente
diferenças ao primeiro ponto
	dia[[i]]<-rep(i,length(arquivo$dia[fixo[i]:ult[i]])-1)$cada dia de coleta
	}

distancias<-unlist(distancias)
tempo<-unlist(tempo)
dia<-unlist(dia)

dtd<-data.frame(distancias,tempo,dia)

Este data.frame seria como uma unidade de amostra das milhares de
distribuições que são possíveis de obter. Como otimizar para que
qualquer número desejado não seja dispendioso de obter?

Obrigado pela atenção

-- 
Roberto Guevara Ferreira Lima

Biólogo - CRBio 6ª Região: 73146/06-D
Mestre em Zoologia
Laboratório de Ecologia e Zoologia de Vertebrados
Instituto de Ciências Biológicas
Universidade Federal do Pará
Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110
Currículo Lattes: http://lattes.cnpq.br/7363382159007600