Produces a data.frame of speed and distance estimates to analyze, as well as a plot highlighting potential speed and distance outliers in telemetry data.
outlie(data,plot=TRUE,by='d',...)# S3 method for outlie
plot(x,level=0.95,units=TRUE,axes=c('d','v'),xlim=NULL,ylim=NULL,...)
Returns an outlie object, which is a data.frame of distance and speed information. Can also produce a plot as a side effect.
telemetry object.
Output a plot highlighting high speeds (blue) and distant locations (red).
Color and size side-effect plot points by 'd', 'v', 'dz', 'vz', for distance from center, minimum speed, vertical distance from center, and minimum vertical speed.
Arguments passed to plot.
outlie object to plot.
Confidence level for error bars.
Convert axes to natural units.
\(x\)-\(y\) axes to plot. Can be any of 'd', 'v', 'dz', 'vz', for time, distance from center, minimum speed, vertical distance from center, and minimum vertical speed.
\(x\)-axis plotting range in SI units.
\(y\)-axis plotting range in SI units.
C. H. Fleming.
If plot=TRUE in outlie(), intervals of high speed are highlighted with blue segments, while distant locations are highlighted with red points.
When plotting the outlie object itself, `median deviation' denotes distances from the geometric median, while `minimum speed' denotes the minimum speed required to explain the location estimate's displacement as straight-line motion. Both estimates account for telemetry error and condition on as few data points as possible. The speed estimates furthermore account for timestamp truncation and assign each timestep's speed to the most likely offending time, based on its other adjacent speed estimate.
The output outlie object contains the above noted speed and distance estimates in a data.frame, with rows corresponding to those of the input telemetry object.
C. H. Fleming et al, ``A comprehensive framework for handling location error in animal tracking data'', bioRxiv 2020.06.12.130195 (2020) tools:::Rd_expr_doi("10.1101/2020.06.12.130195").
as.telemetry.
# Load package and data
library(ctmm)
data(turtle)
# look for outliers in a turtle
OUT <- outlie(turtle[[3]])
# look at the distribution of estimates
plot(OUT)
Run the code above in your browser using DataLab