% Example 6.7
% This shows how to use the function for the whole MBC procedure.

load iris
data = [setosa;versicolor;virginica];
% Call the model-based clustering procedure with a
% maximum of 6 clusters.
[bics,bestmodel,allmodel,Z,clabs] = ...
    mbclust(data,6);

% Display the BIC curves.
plotbic(bics)

% Display the agglomerative part in the 
% One can use the clabs information along with gplotmatrix
% to do a grouped scatterplot matrix of the data. Do the groups seem
% reasonable?