[AniMov] Calculating spatial overlap using kerneloverlap() with BRB UDs

Guilherme De Miranda Mourao guilherme.mourao at embrapa.br
Wed Sep 28 15:34:16 CEST 2016


Hi List:
I am trying to estimate overlaps using the kerneloverlaphr function. The input objects are from estUDm class resulted from the BRB function. However, the overlap values seem odd, as low-overlapped "home ranges" had estimates as higher as the high-overlapped ones.
Anyone can say what I am doing wrong?
Following is an example:

####
library(adehabitatHR)

x<- runif(100,1,2)
y<-runif(100,1,2)
x1<-runif(50,2,3)
y1<-runif(50,2,3)
z<- seq(1:150)
x<-c(x,x1)
y<-c(y,y1)
xy <- as.matrix(data.frame(cbind(x,y)))
id <- as.factor(c(rep(1,50),rep(2,50),rep(3,50)))
class(z) <- "POSIXct"	
traject <- as.ltraj( xy, z, id, typeII=TRUE )
plot(traject)
D <- BRB.D(traject, Tmax=100, Lmin=0.05)
UD <- BRB(traject, D, tau=1, Tmax=100, Lmin=0.05, hmin=0.2) ## Biased random bridge

contr1 <- getverticeshr(UD[[1]], lev =50)
contr2 <- getverticeshr(UD[[2]], lev =50)
contr3 <- getverticeshr(UD[[3]], lev =50)

plot(contr1, col=rgb(0,0,0,0.2),xlim=c(0,3),ylim=c(0,3))  ## plot
plot(contr2, add=TRUE, col=rgb(0,0,0,0.2))
plot(contr3, add=TRUE, col=rgb(0,0,0,0.2))

kerneloverlaphr(UD, method = "PHR")   ## values seem wrong
###

Thank you in advance.
Guilherme

-- 
Guilherme de Miranda Mourão
Pesquisador A
Laboratório de Vida Selvagem

Embrapa Pantanal
Corumbá MS

guilherme.mourao at embrapa.br
www.cpap.embrapa.br

Fone:  +55 67   3234-5929
Fax:   +55 67   3234-5815

____________________________________________________________________________
Aviso de confidencialidade 

Esta mensagem da Empresa  Brasileira de Pesquisa  Agropecuaria  (Embrapa),
empresa publica federal  regida pelo disposto  na Lei Federal no. 5.851,  de
7 de dezembro de 1972,  e  enviada exclusivamente  a seu destinatario e pode
conter informacoes  confidenciais, protegidas  por sigilo profissional.  Sua
utilizacao desautorizada  e ilegal e  sujeita o infrator as penas da lei. Se
voce  a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente,
esclarecendo o equivoco.

Confidentiality note

This message from Empresa  Brasileira de Pesquisa  Agropecuaria  (Embrapa), a
government company  established under  Brazilian law  (5.851/72), is directed
exclusively to  its addressee  and may contain  confidential data,  protected
under  professional secrecy  rules. Its unauthorized  use is illegal and  may
subject the transgressor to the law's penalties. If you are not the addressee,
please send it back, elucidating the failure.



More information about the AniMov mailing list