% Example 10.12
% Shows how to get Andrews curves for a 
% the data used previously.

load example104
% This loads up the reduced BPM features using ISOMAP.
% Use the 3-D data.
X = Yiso.coords{3}';
% Find the observations for the two topics: 6 and 9.
ind6 = find(classlab==6);
ind9 = find(classlab==9);
% Use the csandrews function from the Computational
% Statistics Toolbox.
% Construct the plot for topic 6.
csandrews(X(ind6,:),'-','r')
hold on
% Construct the plot for topic 9.
csandrews(X(ind9,:),':','g')
hold off
