njit-etd2003-081 - New Jersey Institute of Technology

297
The MATLAB program below performs the cluster analysis on the data set created using
the principal components obtained from the PCA on the 12 different HRV parameters
and actual heart rate, blood pressure and breathing (respiration) rate.
function [RITrue,RISelf,EstIndex]=DoClustering(Algorithm,Data,k,TrueID,Repetitions)
% DOCLUSTERING Performs cluster analysis of specified data using specified
% algorithm.
% [RITrue,RISelf,Index]=DOCLUSTERING(Algorithm,Data,k,TrueID,Repetitions)
% where Algorithm is the name of the clustering algorithm (string),
% Data is 2-d matrix, with variables as columns and records as rows,
% k is number of clusters to search for,
% TrueID is a vector of class labels, one for each row in "Data", ([]->not known)
% and Repetitions is the number of times to repeat the algorithm (default 1).
% RITrue is the (adjusted) Rand Index comparing the solution to the true labels
% (where known). RISelf is the (adjusted) Rand Index between pairs of solutions
% (when Repetitions > l).
% Index is the estimated labels of each record.
% AlgorithmList=DOCLUSTERING with no parameters returns list of all
% the aglorithms available
% Copyright (C) David Corney 2000 D.Corney@cs.ucl.ac.uk
%If no inputs, then return list of all aglorithms used (arbitrary order).
if nargin == 0
RITrue={'kmeans','EM_spherical','EM_elliptical',
'hier_centroid','hier_complete','hier_single',
'fuzzy','random'};
return
end
if nargin < 3
error(DOCLUSTERING requires at least 3 arguments');
help(mfilename)
return
end
if nargin < 4, TrueID=[]; end %no known class labels
if nargin < 5, Repetitions = 1;end %only do it once
if nargin < 6, options = 0;