wolffd@0: function som_projections_plot(pmode,varargin)
wolffd@0: 
wolffd@0: % SOM_PROJECTIONS_PLOT Projection visualizations.
wolffd@0: %
wolffd@0: % som_projections_plot(pmode,varargin)
wolffd@0: %
wolffd@0: %   [cPCAarg, Pdata, Pproto] = som_projections(D,sM);
wolffd@0: %   som_projections_plot('scatter',Pdata(:,1:3))
wolffd@0: %   som_projections_plot('scatter',Pproto(:,1:3),Pproto(:,5:7),5,sM)
wolffd@0: %   som_projections_plot('residuals',Pdata(:,1:4))
wolffd@0: %   som_projections_plot('scree',cPCAarg{3})
wolffd@0: %
wolffd@0: % The other arguments depend on the pmode:
wolffd@0: %
wolffd@0: %  pmode = 'scatter'
wolffd@0: %
wolffd@0: %     arg1: Co     (matrix) coordinates
wolffd@0: %     arg2: color  (matrix) color vectors
wolffd@0: %                  (string) colorstring ('k' by default)
wolffd@0: %     arg3: psize  (scalar) point size
wolffd@0: %     arg4: sT     (struct) topology struct, if map grid is drawn
wolffd@0: %  
wolffd@0: %  pmode = 'residuals'
wolffd@0: %
wolffd@0: %     arg1: Co     (matrix) coordinates (2 first columns) + residuals (the rest)
wolffd@0: %     arg2: color  (string) colorstring ('k' by default)
wolffd@0: %
wolffd@0: %  pmode = 'scree'
wolffd@0: %
wolffd@0: %     arg1: eigval (vector) vector of eigenvalues
wolffd@0: %      
wolffd@0: % See also  SOM_PROJECTIONS.
wolffd@0: 
wolffd@0: switch pmode, 
wolffd@0: case 'scatter', 
wolffd@0:     Co = varargin{1}; 
wolffd@0:     if length(varargin)>1, color = varargin{2}; else color = 'k'; end
wolffd@0:     if length(varargin)>2, psize = varargin{3}; else psize = 5; end
wolffd@0:     if length(varargin)>3, sT = varargin{4}; else sT = []; end
wolffd@0:     if isstruct(sT) & strcmp(sT.type,'som_map'), sT = sT.topol; end
wolffd@0:     
wolffd@0:     if isempty(sT),
wolffd@0:         som_grid({'rect',[size(Co,1) 1]},'Coord',Co,'Markercolor',color,'line','none','markersize',psize);
wolffd@0:     else
wolffd@0:         if ischar(color), lcolor = color; else lcolor = 'k'; end
wolffd@0:         som_grid(sT,'Coord',Co,'Markercolor',color,'markersize',psize,'linecolor',lcolor); 
wolffd@0:     end 
wolffd@0:     
wolffd@0: case 'residuals',
wolffd@0:     CoRes = varargin{1}; 
wolffd@0:     n = size(CoRes,1); 
wolffd@0:     if length(varargin)>1, color = varargin{2}; else color = 'k'; end
wolffd@0:     Co = CoRes(:,1:2); Co(end+1,:) = NaN; 
wolffd@0:     res = sqrt(sum(CoRes(:,3:end).*CoRes(:,3:end),2)); 
wolffd@0:     h=plot(Co(:,1),Co(:,2),'k.'); set(h,'color',color);
wolffd@0:     Co(end+1,:) = NaN;
wolffd@0:     res = [res; 0; NaN];  
wolffd@0:     i = [1:n; 1:n; (n+1)+zeros(1,n)]; i = i(:); 
wolffd@0:     j = [(n+1)+zeros(1,n); 1:n; (n+2)+zeros(1,n)]; j = j(:);  
wolffd@0:     h = line(Co(i,1),Co(i,2),res(j)); set(h,'color',color);
wolffd@0:     axis tight, axis equal, view(3)
wolffd@0:     
wolffd@0: case 'scree', 
wolffd@0:     eigval = varargin{1};
wolffd@0:     if size(eigval,1)>1, eigval = eigval'; end
wolffd@0:     eigval = eigval / sum(eigval); 
wolffd@0:     cumeig = cumsum(eigval);
wolffd@0:     bar(cumeig,0.1)
wolffd@0:     i = find(cumeig>=0.75); i75 = i(1); 
wolffd@0:     hold on
wolffd@0:     plot([0 2],cumeig([2 2]),'r-',[0 i75],cumeig([i75 i75]),'r-');
wolffd@0:     set(gca,'YTick',unique([0:0.1:1 cumeig(2) cumeig(i75)]));  
wolffd@0:     axis tight, grid on
wolffd@0: 
wolffd@0: end 
wolffd@0: