plot3dproj.m

来自「一个关于数据聚类和模式识别的程序,在生物化学,化学中因该都可以用到.希望对大家有」· M 代码 · 共 41 行

function plot2dProj(DS)
%plotClassVsFeature: Plot of class vs. feature

%	Roger Jang, 20041209

if nargin == 0, selfdemo; return, end

[featureNum, dataNum]=size(DS.input);
distinctClass = elementCount(DS.output);
classNum = length(distinctClass);
sequence = combine(1:featureNum, 3);
plotNum = size(sequence,1);
side = ceil(sqrt(plotNum));
k = 1;
for i = 1:side,
	for j = 1:side,
		if k <= plotNum,
			subplot(side, side, k);
			for p=1:classNum
				index=find(DS.output==distinctClass(p));
				line(DS.input(sequence(k,1),index), DS.input(sequence(k,2),index), DS.input(sequence(k,3),index),'marker', '.', 'lineStyle', 'none', 'color', getColor(p));
			end
			box on; axis([-inf inf -inf inf]); axis image; view(-37.5, 30);
			if isfield(DS, 'inputName')
				xlabel(['x', num2str(sequence(k,1)), ': ', DS.inputName{sequence(k,1)}]);
				ylabel(['x', num2str(sequence(k,2)), ': ', DS.inputName{sequence(k,2)}]);
				zlabel(['x', num2str(sequence(k,3)), ': ', DS.inputName{sequence(k,3)}]);
			else
				xlabel(['x', num2str(sequence(k,1))]);
				ylabel(['x', num2str(sequence(k,2))]);
				zlabel(['x', num2str(sequence(k,3))]);
			end
			k = k+1;
		end
	end
end

% ====== Self demo
function selfdemo
DS=prData('iris');
feval(mfilename, DS);