already found a bug with bursts containing 2 data points... here is a code corrected: mindistkeep2 <- function(x, threshold) { if (!inherits(x,"ltraj")) stop("x should be of class 'ltraj'") foo <- function(y) { ul <- 1
if (nrow(y) > 2){ for (i in 2:nrow(y)) { if(y$dist[i-1]<threshold){ ul[i]= ul[i-1] y$dist[i] = sqrt( (y$x[i+1] -y$x [ul[i]])*(y$x[i+1] -y$x [ul[i]])+(y$y[i+1] -y$y [ul[i]])*(y$y[i+1] -y$y [ul[i]]) ) } else ul[i]={i} } } else print(nrow (y)) if (!is.null(attr(y, "infolocs"))) { infol <- attr(y, "infolocs") } z <- y[ul,c("x","y")] da <- y$date if (!is.null(attr(y, "infolocs"))) { infol <- attr(y, "infolocs") lt <- as.ltraj(z, da, id=attr(y,"id"), burst=attr(y,"burst"), typeII=attr(x,"typeII"), infolocs=infol) } else { lt <- as.ltraj(z, da, id=attr(y,"id"), burst=attr(y,"burst"), typeII=attr(x,"typeII")) } return(lt) } res <- do.call("c.ltraj", lapply(x, foo)) class(res) <- c("ltraj","list") return(res) } _______________________________________________ AniMov mailing list AniMov@faunalia.it http://lists.faunalia.it/cgi-bin/mailman/listinfo/animov