% 05-22-2008
% Script file: clusterbio.m
%              Drawing 2d dendrogram 
d=16; N=55;    % 55 patients, each has 16 genes
fin=fopen('biodata.txt');
fgetl(fin); fgetl(fin); fgetl(fin);  % skip 3 headlines
A=fscanf(fin,'%f',[1+N d]);
%
X=A(2:56,1:d);   % 38 patients of NVI,  17 patients of HCC VI
R={'Pat101','Pat102','Pat103','Pat104','Pat105','Pat106',...
   'Pat107','Pat108','Pat109','Pat110','Pat111','Pat112',...
   'Pat113','Pat114','Pat115','Pat116','Pat117','Pat118',...
   'Pat119','Pat120','Pat121','Pat122','Pat123','Pat124',...
   'Pat125','Pat126','Pat127','Pat128','Pat129','Pat130',...
   'Pat131','Pat132','Pat133','Pat134','Pat135','Pat136',...
   'Pat137','Pat138',...  
'Pat22201','Pat22202','Pat22203','Pat22204','Pat22205','Pat22206',...
'Pat22207','Pat22208','Pat22209','Pat22210','Pat22211','Pat22212',...
'Pat22213','Pat22214','Pat22215','Pat22216','Pat22217'};
C={'Gene01','Gene02','Gene03','Gene04','Gene05','Gene06','Gene07','Gene08',...
   'Gene09','Gene10','Gene11','Gene12','Gene13','Gene14','Gene15','Gene16'};
clustergram(X,'Pdist','euclidean','Linkage','complete','Dimension',2,...
'ROWLABELS',R,'COLUMNLABELS',C,'Dendrogram',{'color',5})