dtwclust (version 1.2.0)

DBA: DTW Barycenter Averaging

Description

A global averaging method for time series under DTW (Petitjean, Ketterlin and Gancarski, 2011).

Usage

DBA(X, center = NULL, max.iter = 20, norm = "L1", window.size = NULL,
  delta = 1e-06, error.check = TRUE, trace = FALSE)

Arguments

X

A data matrix where each row is a time series. Optionally, a list where each element is a time series.

center

Optionally, a time series to use as reference. It must be a numeric vector. Defaults to a random series of X if NULL.

max.iter

Maximum number of iterations allowed.

norm

Norm for the local cost matrix of DTW. Either "L1" for Manhattan distance or "L2" for Euclidean distance.

window.size

Window constraint for the DTW calculations. NULL means no constraint.

delta

At iteration i, if all(abs(center_{i} - center_{i-1}) < delta), convergence is assumed.

error.check

Should inconsistencies in the data be checked?

trace

If TRUE, the current iteration is printed to screen.

Value

• The average time series.

Details

This function tries to find the optimum average series between a group of time series in DTW space. Refer to the cited article for specific details on the algorithm.

If a given series reference is provided in center, the algorithm should always converge to the same result provided the rows of X keep the same values, although their order may change.

References

Petitjean F, Ketterlin A and Gancarski P (2011). ``A global averaging method for dynamic time warping, with applications to clustering.'' Pattern Recognition, 44(3), pp. 678 - 693. ISSN 0031-3203, http://dx.doi.org/10.1016/j.patcog.2010.09.013, http://www.sciencedirect.com/science/article/pii/S003132031000453X.

Examples

# Sample data
data(uciCT)

# Obtain an average for the first 5 time series
dtw.avg <- DBA(CharTraj[1:5], CharTraj[[1]], trace = TRUE)
plot(dtw.avg, type="l")

# Change the provided order
dtw.avg2 <- DBA(CharTraj[5:1], CharTraj[[1]], trace = TRUE)
all(dtw.avg == dtw.avg2)