%  loads features from ccSchzSchz and ccSchzNorms
 %  male subjects only

 fprintf(fid,'Schizophrenics vs Normals - males\n');
 fprintf(fid,'Data: from ccSchzSchz and ccSchzNorms\n');
 fprintf(fid, '# Schizophrenics = 23; # Normals = 23\n\n');
 

 ClassA =  zeros(35,153);
 ClassB =  zeros(48,153);
 
 cd  ../../ccSchzSchz
 ClassA(:,1:3) = [load('rel_gn') load('rel_spl') load('rel_rbn')];
 ClassA(:,4:75) = load('axis_rho');

 cd  ../ccSchzNorms
 ClassB(:,1:3) = [load('rel_gn') load('rel_spl') load('rel_rbn')];
 ClassB(:,4:75) = load('axis_rho');

 ClassA = ClassA(13:35,:);
 ClassB = ClassB(22:44,:);

 cd ../classification/1FeatureAnalysis
  
 divno = [36 18 9 6 4 3 2];
 featureno = [75 111 129 138 144 148 151];
 for nn = 1:7
    N = divno(nn);
    Last = featureno(nn);
    I = round(72/N);
    for k=1:N
        kk = round(4+(k-1)*I);
        ClassA(:,Last+k)=sum([ClassA(:,kk:kk+I-1)],2); 
        ClassB(:,Last+k)=sum([ClassB(:,kk:kk+I-1)],2); 
     end
 end
 
 

 % first 3 features: gn,spl,rbn
 % features 4-75: samples of rho along the main axis
 % remaining: averages of the 72 rho samples,
 % features 76-111: averages over groups of 2
 % features 112-129: averages over groups of 4
 % features 130-138: averages over groups of 8
 % features 139-144: averages over groups of 12
 % features 145-148: averages over groups of 18
 % features 149-151: averages over groups of 24
 % features 152-153: averages over groups of 36