swh:1:snp:a4c99a50dc49f82b591f268001b320f8c3ca0041
Tip revision: e3fb24ebd97f7686caee78bcef4e558b908a465e authored by jmc on 01 December 2007, 00:00:00 UTC
version 0.97-1
version 0.97-1
Tip revision: e3fb24e
trackSpeed.R
trackSpeed <- function(coords, time) {
dt <- diff(as.numeric(time))
d2 <- diff(as.matrix(coords))^2
speed <- d2[,1]
ndim <- ncol(d2)
if(ndim > 1)
speed <- d2 %*% rep(1,ndim)
else
speed <- d2
sqrt(as.vector(speed))/dt
}