getEVAF: getEVAF

In EliGurarie/smoove: Simulation and Estimation of Correlated Velocity Movement (CVM) Models

getEVAF

Description

Obtain empirical velocity autocorrelation function.

Usage

getEVAF(Z, XY = NULL, V = NULL, T, lagmax = max(T)/2)

Arguments

Z	complex vector of locations
XY	two-column matrix of x-y coordinates of locations
V	complex vector of velocities
T	times of observations
lagmax	maximum lag to compute the empirical VAF for

Details

Requires at least one of complex locations, XY coordinates, or velocties to compute.

Value

A data frame with columns: lag and vaf

Examples

nu <- 2; tau <- 5
ucvm <- simulateUCVM(nu=nu, tau = tau, T.max = 1000, dt = .1)
evaf <- with(ucvm, getEVAF(Z=Z,T=T, lagmax = 20))
plot(evaf)
eta <- 2/sqrt(pi)*nu
curve(eta^2 * exp(-x/tau), add=TRUE, col=2, lwd=2)

EliGurarie/smoove documentation built on Aug. 2, 2022, 10:26 p.m.