matrixProfile (version 0.5.0)

stamp: Scalable Time series Anytime Matrix Profile (stamp)

Description

Calculates a matrix profile of given data using STAMP algorithm.

Usage

stamp(q, t, by = 10, isPlot = FALSE)

Arguments

q

A query data for dot product.

t

A timeseries data for analysis.

by

A parameter that indicates the progress of the process in the process of calculating the matrix profile. For example, if by is 10, the text is displayed in the console window every 10 percent.

isPlot

A parameter that determines whether or not to draw a plot in the middle of calculating a matrix profile. The default value is FALSE.

Value

An object of class stamp.models.

MP

A matrix profile computed by given data.

MPI

A matrix profile index computed by given data.

MTI

A motif index of matrix profile. Unlike in the original paper, it denotes the pair of motif index with the smallest value of matrix profile.

Details

The matrix profile is calculated by the self join method using the STAMP algorithm. One of the key features of the STAMP algorithm is the attribute anytime. In other words, because the matrix profile is computed rather than randomly, the computation speed is the same, but it is quickly optimized.

References

Yeh, C. C. M., Zhu, Y., Ulanova, L., Begum, N., Ding, Y., Dau, H. A., ... & Keogh, E. (2016) <DOI:10.1109/ICDM.2016.0179>. 2016 IEEE 16th International Conference on Data Mining (ICDM), Barcelona, 2016, pp. 1317-1322.

http://www.cs.ucr.edu/~eamonn/MatrixProfile.html

See Also

mass (in package matrixProfile)

Examples

# NOT RUN {
# data input
dt = AirPassengers
dt = as.vector(dt)

# generates matrix profile
stamp <- stamp(q = dt[1:12], t = dt[-c(1:12)])

# plotting
par(mfrow = c(2,1))
plot(dt, type = "l", main = "Original Timeseries")
plot(stamp$MP, type = "l", main = "Matrix Profile", xlim = c(0, length(dt)))
# }