Low dimensional representation
The software also provides a basic means of visualizing low-dimensional representations of the data, using PCA as
TS_PlotLowDim.This can be done for a time-series dataset as follows:
1
TS_PlotLowDim('norm','pca');
Copied!
This uses the normalized data (specifying
'norm'), plotting time series in the reduced, two-dimensional principal components space of operations (the leading two principal components of the data matrix).By default, the user will be prompted to select 10 points on the plot to annotate with their corresponding time series, which are annotated as the first 300 points of that time series (and their names by default).
After selecting 10 points, we have the following:
pca_image
The proportion of variance explained by each principal component is provided in parentheses in the axis label.
Customizing annotations
Annotation properties can be altered with some detail by specifying properties as the
annotateParams input variable, for example:1
% Set whether to plot marginal distributions:
2
showDistributions = false;
3
% Set up the annotateParams structure:
4
annotateParams = struct;
5
annotateParams.n = 8; % annotate 8 time series
6
annotateParams.maxL = 150; % annotates the first 150 samples of time series
7
annotateParams.userInput = false; % points not selected by user but allocated randomly
8
annotateParams.textAnnotation = false; % don't display names of annotated time series
9
​
10
% Generate a plot using these settings:
11
TS_PlotLowDim('norm','pca',showDistributions,'',annotateParams);
Copied!
which yields:
annotated plot
Plotting grouped time series
If groups of time series have been specified (using
TS_LabelGroups), then these are automatically recognized by TS_PlotLowDim, which will then distinguish the labeled groups in the resulting 2-dimensional annotated time-series plot.Consider the sample dataset containing 20 periodic signals with additive noise (given the keyword noisy in the database), and 20 purely periodic signals (given the keyword periodic in the database). After retrieving and normalizing the data, we store the two groups in the metadata for the normalized dataset HCTSA_N.mat:
1
TS_LabelGroups('norm',{'noisy','periodic'},'ts');
Copied!
1
We found:
2
noisy -- 20 matches
3
periodic -- 20 matches
4
Saving group labels and information back to HCTSA_N.mat... Saved.
Copied!
Now when we plot the dataset in
TS_PlotLowDim, it will automatically distinguish the groups in the plot and attempt to classify the difference in the reduced principal components space.Running the following:
1
annotateParams = struct('n',6); % annotate 6 time series
2
showDistributions = true; % plot marginal distributions
3
TS_PlotLowDim('norm','pca',showDistributions,'',annotateParams);
Copied!
The function then directs you to select 6 points to annotate time series to, producing the following:
Notice how the two labeled groups have been distinguished as red and blue points, and a linear classification boundary has been added (with in-sample misclassification rate annotated to the title and to each individual principal component). If marginal distributions are plotted (setting
showDistribution = true above), they are labeled according to the same colors.Last modified 3yr ago
Copy link