tlda.m

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

%load iris.dat
%feature = normal(iris(:, 1:end-1));
%class = iris(:,end);

%load wine.dat;
%feature = normal(wine(:,2:end));
%class = wine(:,1);

load glass.dat;
feature = normal(glass(:,1:end-1));
class = glass(:,end);

dim = size(feature,2);

recog = ldainsel(feature, class, dim);
k2 = 1:5;
for i=1:length(k2),
	this = fldainsel(feature, class, size(feature,2), 1, k2(i));
	recog = [recog, this]; 
end

marker = {'+', 'o', 'x', '*', 'd', 's', '.', 'p', 'h', 'v', '^', '<', '>'};
color = {'m'; 'b'; 'r'; 'g'; 'c'; 'k'};
figure
%plot(1:dim, recog);
hold on
for i=1:(1+length(k2)),
	plot(1:dim, recog(:,i), [color{i}, '-', marker{i}]); 
end
hold off
xlabel('Dimensions used in LDA');
ylabel('Leave-one-out recognition rates (%)');
box on;
axis tight
legend('LDA', ...
	'Fuzzy LDA (k=1)', ...
	'Fuzzy LDA (k=2)', ...
	'Fuzzy LDA (k=3)', ...
	'Fuzzy LDA (k=4)', ...
	'Fuzzy LDA (k=5)');