%SCATTERD Display scatterplot% %   H = SCATTERD(A)%   H = SCATTERD(A,DIM,S,CMAP,FONTSIZE,'label','both','legend','gridded')%% INPUT%   A     Dataset or matrix%   DIM   Number of dimensions: 1,2 or 3 (optional; default: 2)%   S     String specifying the colors and markers (optional)%   CMAP  Matrix with a color map (optional)% % OUTPUT%   H     Vector of handles%% DESCRIPTION% SCATTERD(A) displays a 2D scatterplot of the first two features of the % dataset A. If the number of dimensions DIM is specified (1..3), it plots % the first D features in a D-dimensional plot (D<4). If the plot string S % is provided, e.g. S = 'w+', all points are plotted accordingly. If given,% different plot strings are used for different classes. See PLOT for the% specification of plot strings.%   % If CMAP is specified, the color of the object symbols is determined by % CMAP indexed by the object labels. A colormap has the size [C x 3], where % C is the number of classes. The three components of CMAP(I,:) determine % the red, green and blue components of the color. For instance: %%   MAP = HSV; [M,K] = SIZE(A); LABELS = CEIL(64*[1:M]'/M); %   A = DATASET(A,LABELS); SCATTERD(A,'.','COLORMAP',MAP); %% This may be used for tracking ordered objects.%  % FONTSIZE may be a vector with three elements: fontsize, markersize and% size of the label font in case of a label plot.%  % Various other options are:%   'LABEL'  : plot labels instead of symbols%   'BOTH'   : plot labels next to each sample%   'LEGEND' : place a legend in the figure%   'GRIDDED': make a grid of 2D scatterplots of each pair of features%  % All the parameters, except for the dataset A can be specified in any % order or can be left out.%  % Classifiers can be plot in the scatterplot by PLOTC.% Note that PLOTC does not work for 1D and 3D scatterplots.%% EXAMPLES% See PREX_CONFMAT, PREX_DENSITY, PREX_PLOTC, PREX_MCPLOT.%% SEE ALSO% DATASETS, COLORMAP, PLOT, PLOTC% Copyright: D. de Ridder, R.P.W. Duin, duin@ph.tn.tudelft.nl% Faculty of Applied Sciences, Delft University of Technology% P.O. Box 5046, 2600 GA Delft, The Netherlands% $Id: scatterd.m,v 1.14 2004/10/27 09:16:31 duin Exp $% REVISIONS% DR1 - Dick, 05-10-2004% Added plotting of unlabeled data as 'k.'.function handle = scatterd(a,p1,p2,p3,p4,p5,p6,p7,p8)	prtrace(mfilename);	% Defaults	d = min(size(a,2),2);		% Dimensionality of plot	s = []; 					% Plot symbol(s)	map = []; 					% Color map	plotlab = 0;				% Put text labels instead of or next to samples	plotsym = 1;				% Plot symbols for samples	plotlegend = 0; 		    % Plot legend	gridded = 0;				% Make a gridded plot	gridrun = 0;				% Inner loop in a gridded plot?	font_size = [];	mark_size = [];	lab_size  = [];	hold_axis = ishold; % A flag to check if 'hold on' is set for the current axis	a = dataset(a); 		% Allow for a non-dataset data	if (nargin < 9), par{8} = []; else, par{8} = p8; end	if (nargin < 8), par{7} = []; else, par{7} = p7; end	if (nargin < 7), par{6} = []; else, par{6} = p6; end	if (nargin < 6), par{5} = []; else, par{5} = p5; end	if (nargin < 5), par{4} = []; else, par{4} = p4; end	if (nargin < 4), par{3} = []; else, par{3} = p3; end	if (nargin < 3), par{2} = []; else, par{2} = p2; end	if (nargin < 2), par{1} = []; else, par{1} = p1; end	% Set up default values.	for i = 1:5		if (~isempty(par{i})) 			if (length(par{i}) == 1) & par{i} < 5 & (~ischar(par{i})) 	% Dimensionality				d   = par{i}; par{i} = 2; 		% Necessary for gridded: D needs to be 2.			elseif ((size(par{i},1) > 1) & (size(par{i},2)==3) & (~ischar(par{i}))) 	% Color map				map = par{i};			elseif (strcmp(par{i},'label'))				plotlab = 1; plotsym = 0;			elseif (strcmp(par{i},'both'))				plotlab = 1; plotsym = 1;			elseif (strcmp(par{i},'legend'))				plotlegend = 1;			elseif (strcmp(par{i},'gridded'))				gridded = 1; 				par{i} = 'gridrun'; 					% Necessary for gridded: otherwise an infinite recursion :)			elseif (strcmp(par{i},'gridrun'))				gridrun = 1;			elseif ~isstr(par{i}) & length(par{i}) <= 3 & par{i}(1) >= 5				font_size = par{i}(1); 				if length(par{i}) >= 2					mark_size = par{i}(2);				end				if length(par{i}) == 3					lab_size = par{i}(3);				end			else				s	= par{i};			end		end	end	if (gridrun)		if isempty(font_size), font_size = 10; end		if isempty(mark_size), mark_size =  5; end		if isempty(lab_size),   lab_size =  8; end	elseif (~hold_axis)		%clf; 		cla;		if isempty(font_size), font_size = 16; end		if isempty(mark_size), mark_size =  7; end		if isempty(lab_size),   lab_size = 14; end	else		if isempty(font_size), font_size = get(gca,'fontsize'); end		if isempty(mark_size), mark_size = font_size/2; end		if isempty(lab_size),   lab_size = font_size-2; end	end	feats = getfeatlab(a,'string');	if size(feats,2) == 1		feats = [repmat('Feature ',size(feats,1),1) feats];	end	if (gridded)		clf;		gs = size(a,2);		for i = 1:gs			for j = 1:gs				subplot(gs,gs,(i-1)*gs+j);				gridrun = 1;				h = feval(mfilename,a(:,[j i]),par{1},par{2},par{3},par{4},par{5},par{6},par{7});				gridrun = 0;				if (i~=gs), xlabel(''); end				if (j~=1),  ylabel(''); end			end			subplot(111); 		% Takes care of clf for the next scatterplot.		end		return;	end	if (isa(a,'dataset')) & (~isempty(getlablist(a)))		[m,k,c] = getsize(a);		lab = getnlab(a); 		lablist = getlablist(a,'string');		dataset_name = getname(a);		a = double(a);	else		[m,k] = size(a);		lab = ones(m,1);		dataset_name = [];		c = 1;		a = double(a);	end	% Character string defining the plotting setup in terms of symbols and colors.	if (isempty(s))		vers = version;		if (str2num(vers(1)) < 5)			col = 'gbrm';			sym = ['+*xo.']';			i   = [1:20];			ss  = [col(i-floor((i-1)/4)*4)' sym(i-floor((i-1)/5)*5)];		else			col = 'gbrm';			sym = ['+*oxsdv^<>p']';			i   = [1:44];			ss  = [col(i-floor((i-1)/4)*4)' sym(i-floor((i-1)/11)*11)];		end    ss  = ['k.'; ss];  								% DR1 - Add symbol for "unlabeled" data.		[ms,ns] = size(ss);		if ms == 1, ss = setstr(ones(m,1)*ss); end	else		if size(s,1) == 1			ss = repmat(s,c,1); s = [];		else			ss = s; s = [];		end    ss  = ['k.'; ss];  								% DR1 - Add symbol for "unlabeled" data.	end	% Define some 'space' OY to be added around the data plotted in symbols.	oy = zeros(1,3);			if (plotsym)		oy = 0.03*(max(a)-min(a));	else		s = 'w.'; 				% Plot white spot instead of symbols.	end	oy(2) = 0;	% Make a plot.	lhandle = []; thandle = [];  % Also plot label "0" (unlabeled).	for i = 0:c		J = find(lab==i);		if (isempty(s)), symbol = ss(i+1,:); else, symbol = s; end		if ((d == 1) & ~isempty(J))			h = plot(a(J,1),zeros(length(J),1),symbol);			hold on;			set(h,'markersize',mark_size);			lhandle = [lhandle h];			if (plotlab)				for j = 1:length(J)					h = text(a(J(j),1)+oy(1),oy(2),lablist(lab(J(j)),:));					set(h,'fontsize',lab_size);					thandle = [thandle h]; if (~isempty(map)), set (h, 'color', map(i+1,:)); end				end			end		elseif ((d == 2) & ~isempty(J))			h = plot(a(J,1),a(J,2),symbol);			hold on;			set(h,'markersize',mark_size);			lhandle = [lhandle h];			if (plotlab)				for j = 1:length(J)					h = text(a(J(j),1)+oy(1),a(J(j),2)+oy(2),lablist(lab(J(j)),:));					set(h,'fontsize',lab_size);					thandle = [thandle h]; if (~isempty(map)), set (h, 'color', map(i+1,:)); end				end			end		elseif (~isempty(J))			h = plot3(a(J,1),a(J,2),a(J,3),symbol);			hold on;			set(h,'markersize',mark_size);			lhandle = [lhandle h];			if (plotlab)				for j = 1:length(J)					h = text(a(J(j),1)+oy(1),a(J(j),2)+oy(2),a(J(j),3)+oy(3),lablist(lab(J(j)),:));					set(h,'fontsize',lab_size);					thandle = [thandle h]; if (~isempty(map)), set (h, 'color', map(i+1,:)); end							end			end		end	end	if (plotsym)		if (~isempty(map))			for i = 0:c, set (lhandle(i+1), 'color', map(i+1,:)); end		end		if (plotlegend), 			[ht, hl] = legend(lhandle(:),lablist);			hl = hl(:)';			%set(hl(1:2*c),'markersize',mark_size);			lhandle = [lhandle hl];			thandle = [thandle ht(:)'];		end	end	% !%_%*!_% Matlab	set(gca,'fontsize',font_size);	if (~hold_axis)		dd = (max(a) - min(a))*0.05;		% offset, avoiding points on plot box.		J = find(dd==0);		dd(J) = dd(J)+1;		if (d == 1),     			axis ([min(a(:,1))-dd(1) max(a(:,1))+dd(1) -0.5 0.5]); 			hx = xlabel(feats(1,:));			thandle = [thandle hx];		elseif (d == 2),			axis ([min(a(:,1))-dd(1) max(a(:,1))+dd(1) min(a(:,2))-dd(2) max(a(:,2))+dd(2)]); 			hx = xlabel(feats(1,:)); hy = ylabel(feats(2,:));			thandle = [thandle(:)' hx hy];			view(2); 		else		% D = 3            			axis ([min(a(:,1))-dd(1) max(a(:,1))+dd(1) min(a(:,2))-dd(2) max(a(:,2))+dd(2) min(a(:,3))-dd(3) max(a(:,3))+dd(3)]); 			hx = xlabel(feats(1,:)); hy = ylabel(feats(2,:)); hz = zlabel(feats(3,:));			thandle = [thandle hx hy hz];			view(3); 		end	end	if (~gridrun) & (length(get(gcf,'children')) == 1 & any(get(gca,'position')>0.80))		set(gca,'position',[0.10 0.13 0.82 0.78]); % make axis labels visible	end			if (~gridrun) & (~isempty(dataset_name))		title(dataset_name);	end	hold off;	if (nargout > 0)		handle = [lhandle thandle];	endreturn;