Dear Members,
My objective is to establish weekly space/time kernels (using kernelkc),
separately for 20 animals (saved as ltraj class). Having done that I would
like to extract the 95% home range for each specified time period and
animal (using getverticeshrk) and save each as as a separate spatialpolygon
(sp class) - maintaining the id of the animal and timestamp indicating the
specific date of that kernelkc output. I am planning to use the
spatialpolygon to extract environmental variables at that specific
time/animal.
This example/attempt is based from an example found in the adehabitat
package:
#data
library(adehabitat)
puech <-data(puechcirc)
#time vector for kernel estimates to be calculated
vv <- seq(min(puechcirc[[2]]$date), max(puechcirc[[2]]$date),
length=50)
## The estimate the UD
re <- lapply(1:length(vv),
function(i) {
uu <- kernelkc(impkoball, h=c(3000,3000,30*24*3600),
tcalc = vv[i], same4all=TRUE,
extent=0.1)
})
#The above example does not include output for a home range
contours/polygon, i have tried various ways of including the following
functions inside the above function after the kernelkc function, which did
not work...:
hr <-getverticeshrk(uu, lev=95)
#and
lapply(1:length(uu), function(i) {
getverticeshrk(uu[[i]], clev=95)
})
Any suggestions how I can make this work will be much appreciated. Thank
you for your time.
Regards,
Don
[[alternative HTML version deleted]]
_______________________________________________
R-sig-ecology mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
