Algoritmo das Projeções Sucessivas Para Seleção de ... - PPGQ

Anexos
end
save(filename,'L')
else % Carregar
filename = input('Entre com o nome do arquivo (entre apóstrofos e sem extensão): ');
load(filename)
end
% Determinação do custo de classificação p/ todas as cadeias de variaveis
disp('Determinando o custo de classificação para determinar a melhor cadeia de variaveis ...')
disp(' ')
R = zeros(1,N2);
Lopt = zeros(N2,N2);
custo = zeros(N2,Nlambdas);
N2
for N = N1:N2 % Para a cadeia de comprimento N
N
for i = 1:Nlambdas % partindo da variável i
lambdas = L(1:N,i); % Variáveis da cadeia
% Respostas instrumentais nas variáveis da cadeia
S = cov(Xpooled(:,lambdas),1);
invS = inv(S);
Val2 = Val(:,lambdas);
custoaux = 0;
for j=1:size(Val,1) % Para cada objeto de Val
x = Val2(j,:);
grupox = Group_Val(j); % Classificação correta do objeto
for k=1:C % Para cada classe
mu = media{k};
mu = mu(lambdas);
r(k) = (x - mu)*invS*(x - mu)';
end
num = r(grupox); % Distancia de Mahalanobis a classe correta
remaining = setdiff([1:C],grupox); % Demais classes
den = min(r(remaining)); % Menor Distancia de Mahalanobis as classes restantes
custoaux = custoaux + num/den;
end
% Calculo do custo associado a cadeia
custo(N,i)= custoaux/size(Val,1); % Custo medio
end
[R(N) imin] = min(custo(N,:));
Lopt(1:N,N)=L(1:N,imin);
end
[Rbest,Nbest] = min(R(N1:N2));
Nbest = Nbest+N1-1;
custoopt = custo(Nbest,:);
l = (Lopt(1:Nbest,Nbest))';
disp(['Menor custo obtido no conjunto de validação: ' num2str(Rbest)])
disp(['Número ideal de variáveis: ' num2str(Nbest)])
figure,plot([N1:N2],R(N1:N2)),grid
axis([N1 N2 -inf inf])
119