Package: matrixprofile-ts
matrixprofile-tsmatrixprofile-ts is a Python 2 and 3 library for evaluating time series data using the Matrix Profile algorithms developed by theKeoghandMueenresearch groups at UC-Riverside and the University of New Mexico. Current implementations include MASS, STMP, STAMP, STAMPI, STOMP, SCRIMP++, and FLUSS.Read the Target blog posthere.Further academic description can be foundhere.The PyPi page for matrixprofile-ts ishereContentsMatrix Profile FoundationQuestionsInstallationQuick startDetailed ExampleAlgorithm ComparisonMatrix Profile in Other LanguagesContactCitationsMatrix Profile Foundationmatrixprofile-ts is part of the Matrix Profile Foundation, an organization dedicated to facilitating community awareness & adoption of the Matrix Profile through educational resources and high-quality code. More information can be found athttps://www.matrixprofile.org/.An interactive Matrix Profile UI can be found athttps://ui.matrixprofile.org/.QuestionsA key goal of the Matrix Profile Foundation is to provide assistance with deploying the Matrix Profile in real-world settings. In addition to e-mailing foundation members, MPF provides a dedicated Discord channel where practitioners can discuss Matrix Profile applications and seek specific assistance. To join the channel,click here.InstallationMajor releases of matrixprofile-ts are available on the Python Package Index:pip install matrixprofile-tsDetails about each release can be foundhere.Quick start>>> from matrixprofile import *
>>> import numpy as np
>>> a = np.array([0.0,1.0,1.0,0.0,0.0,1.0,1.0,0.0,0.0,1.0,1.0,0.0])
>>> matrixProfile.stomp(a,4)
(array([0., 0., 0., 0., 0., 0., 0., 0., 0.]), array([4., 5., 6., 7., 0., 1., 2., 3., 0.]))Note that SCRIMP++ is highly recommended for calculating the Matrix Profile due to its speed and anytime ability.ExamplesJupyter notebooks containing various examples of how to use matrixprofile-ts can be found underdocs/examples.As a basic introduction, we can take a synthetic signal and use STOMP to calculate the corresponding Matrix Profile (this is the same synthetic signal as in theGolang Matrix Profile library). Code for this example can be foundhereThere are several items of note:The Matrix Profile value jumps at each phase change. High Matrix Profile values are associated with "discords": time series behavior that hasn't been observed before.Repeated patterns in the data (or "motifs") lead to low Matrix Profile values.We can introduce an anomaly to the end of the time series and use STAMPI to detect itThe Matrix Profile has spiked in value, highlighting the (potential) presence of a new behavior. Note that Matrix Profile anomaly detection capabilities will depend on the nature of the data, as well as the selected subquery length parameter. Like all good algorithms, it's important to try out different parameter values.Algorithm ComparisonThis section shows the matrix profile algorithms and the time it takes to compute them. It also discusses use cases on when to use one versus another. The timing comparison is based on the synthetic sample data set to show run time speed.For a more comprehensive runtime comparison, please review the notebookdocs/examples/Algorithm Comparison.ipynb.All time comparisons were ran on a 4 core 2.8 ghz processor with 16 GB of memory. The operating system used was Ubuntu 18.04LTS 64 bit.AlgorithmTime to CompleteDescriptionSTAMP310 ms ± 1.73 ms per loop (mean ± std. dev. of 7 runs, 1 loop each)STAMP is an anytime algorithm that lets you sample the data set to get an approximate solution. Our implementation provides you with the option to specify the sampling size in percent format.STOMP79.8 ms ± 473 µs per loop (mean ± std. dev. of 7 runs, 10 loops each)STOMP computes an exact solution in a very efficient manner. When you have a historic time series that you would like to examine, STOMP is typically the quickest at giving an exact solution.SCRIMP++59 ms ± 278 µs per loop (mean ± std. dev. of 7 runs, 10 loops each)SCRIMP++ merges the concepts of STAMP and STOMP together to provide an anytime algorithm that enables "interactive analysis speed". Essentially, it provides an exact or approximate solution in a very timely manner. Our implementation allows you to specify the max number of seconds you are willing to wait for a solution to obtain an approximate solution. If you are wanting the exact solution, it is able to provide that as well. The original authors of this algorithm suggest that SCRIMP++ can be used in all use cases.Matrix Profile in Other LanguagesR:https://github.com/franzbischoff/tsmpGolang:https://github.com/aouyang1/go-matrixprofileC++:https://github.com/shapelets/khivaContactAndrew Van Benschoten (avbs89@gmail.com)Tyler Marrs (tylerwmarrs@gmail.com)CitationsChin-Chia Michael Yeh, Yan Zhu, Liudmila Ulanova, Nurjahan Begum, Yifei Ding, Hoang Anh Dau, Diego Furtado Silva, Abdullah Mueen, Eamonn Keogh (2016). Matrix Profile I: All Pairs Similarity Joins for Time Series: A Unifying View that Includes Motifs, Discords and Shapelets. IEEE ICDM 2016Matrix Profile II: Exploiting a Novel Algorithm and GPUs to break the one Hundred Million Barrier for Time Series Motifs and Joins.  Yan Zhu, Zachary Zimmerman, Nader Shakibay Senobari, Chin-Chia Michael Yeh, Gareth Funning, Abdullah Mueen, Philip Berisk and Eamonn Keogh (2016). EEE ICDM 2016Matrix Profile V: A Generic Technique to Incorporate Domain Knowledge into Motif Discovery. Hoang Anh Dau and Eamonn Keogh. KDD'17, Halifax, Canada.Matrix Profile XI: SCRIMP++: Time Series Motif Discovery at Interactive Speed. Yan Zhu, Chin-Chia Michael Yeh, Zachary Zimmerman, Kaveh Kamgar and Eamonn Keogh, ICDM 2018.Matrix Profile VIII: Domain Agnostic Online Semantic Segmentation at Superhuman Performance Levels. Shaghayegh Gharghabi, Yifei Ding, Chin-Chia Michael Yeh, Kaveh Kamgar, Liudmila Ulanova, and Eamonn Keogh. ICDM 2017.