swh:1:snp:a4c99a50dc49f82b591f268001b320f8c3ca0041
Tip revision: ce60f670aac0a708d3a9af5f0cf46a752d46ba8a authored by John M Chambers on 28 October 2020, 07:59:48 UTC
version 1.0-6.1
version 1.0-6.1
Tip revision: ce60f67
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
}