Encontrei uma função feita por J.Ahumada e achei super interessante, tudo a ver com meu trabalho. A função é para separar registros fotográficos de uma espécie em uma determinada unidade amostral(ua) a partir de um intervalos de independência desejado.
Criei um objeto com as informações de uma determinada espécie em uma única unidade amostral. O objeto chama-se "paca".
paca <- filter(meus.dados, ua==" GF", data, hora, especie=="Cuniculus paca")
paca
ua data hora especie
(chr) (time) (chr) (chr)
1 GF 2012-06-02 01:12:00 Cuniculus paca
2 GF 2012-06-11 23:50:00 Cuniculus paca
3 GF 2012-06-12 00:06:00 Cuniculus paca
4 GF 2012-06-12 01:16:00 Cuniculus paca
5 GF 2012-07-11 20:35:00 Cuniculus paca
6 GF 2012-07-24 23:52:00 Cuniculus paca
7 GF 2012-08-01 21:39:00 Cuniculus paca
8 GF 2012-08-09 02:37:00 Cuniculus paca
9 GF 2012-08-11 00:24:00 Cuniculus paca
10 GF 2012-08-13 00:55:00 Cuniculus paca
11 GF 2012-08-13 19:47:00 Cuniculus paca
12 GF 2012-08-15 19:16:00 Cuniculus paca
13 GF 2012-08-18 02:35:00 Cuniculus paca
14 GF 2012-08-18 22:28:00 Cuniculus paca
15 GF 2012-08-24 02:27:00 Cuniculus paca
Ao rodar a função, retorna a sequência de números correspondente ao número da linha criada pelo R (1 a 15). E quando o registro não respeita o intervalo de 60min ele repete o número da linha onde está o registro.
reg.independentes<-function(dados,independencia){
l<-length(dados$data)
intervalo<-diff(dados$data)
intervalo<-intervalo/60 #informar intervalo de independência em minutos
intervalo<-as.numeric(intervalo)
ev<-1;res<-numeric()
cond<-intervalo> independencia
for(i in 1:(l-1)){
if(!cond[i]) ev<-ev
else ev<-ev+1
res<-c(res,ev)
}
c(1,res)
}
reg.independentes(paca, 60)
[1] 1 2 3 3 4 5 6 7 8 9 9 10 11 11 12 12 13 14 15
A função não considera o fato do registro ter sido em datas diferentes, está considerando apenas o horário. Gerando duas situações:
Primeira: repetindo linhas onde o registro foi na mesma data,porém com intervalos maior que 60min. Por exemplo, ele repete a linha 3, ao conferir os registros, eles possuem intervalo maior que 60min, conforme o desejado (data é igual e horário diferente - 00:06 e 01:16). Não entendi o por que, não era para este registro estar sendo sinalizado!!
Segunda: repetindo linhas onde o registro é em datas diferentes apenas do horário ser semelhante. A função não está considerando a questão da data ser diferente,por exemplo, sinalizou a linha (9, 11 e 12) porém os registros são em datas diferentes, tornando-se independente.
O registro é considerado Não Independente se ocorrer no mesmo dia e em intervalo menor que 60min. Se o registro for em horários semelhantes, mas em datas diferentes são considerados Independentes (isso que eu preciso que a função faça)
Tentei alterar a fórmula da função mas não estou conseguindo. Gostaria de acrescentar que a função retornasse uma tabela com apenas os registros independentes.... Alguém pode me ajudar?