mrmr-ftest.zip_feature selection_image registration_mrmr_mrmr ma

% Function: Minimum redundancy-Maximum relevance (T-test) % for continuous variables\ % X should include more than 2 attributes % Written by Yi Zhang % 07/24/2007 function Wnum=mRMRC(X,y,Topn,TypeD,TypeM) % Wnum is the attribute index of the top Topn % if TypeD is 1, TypeM can be chosen as 1 or 2 % if TypeD is 2, TypeM can be chosen as 3 or 4 [m,n] = size(X); Wnum = zeros(1,Topn); % Calculate F-test value F = Ftest(X,y,n); F1 = zeros(1,n); W = zeros(1,n); X1 = zeros(1,n); X2 = zeros(m,n); % Reduce redundance between two attributes if TypeD == 1 % WD = abs(PearsonCorrD(X))*(10^(-15)); WD = abs(PearsonCorrD(X)); else % WD = EuclideanD(X); % scale X into [0,1]; X1 = min(X); X3 = max(X)-min(X); for i = 1:n for j = 1:m X2(j,i) = (X(j,i)-X1(i))/X3(i); end end WD = EuclideanD(X2); end Wnum(1) = 1; Ft = F(1); for i = 1:n-1 if F(i+1) >= Ft Wnum(1) = i+1; Ft = F(i+1); end end F1(Wnum(1)) = 1; Ftemp1 = F(Wnum(1)); Ftemp2 = 0; for i = 2:Topn for j = 1:n if F1(j) == 0 if TypeM == 1 % FCD: max(F-Wd); Wd: Pearson correlation for s = 1:i-1 W(j) = Ftemp2-WD(j,Wnum(s)); end Ftempa = Ftemp1+F(j); Ftempb = Ftemp2+W(j); W(j) = Ftempa/i+Ftempb/(i^2); elseif TypeM == 2 %FCQ: max(F/Wd); Wd: Pearson correlation for s = 1:i-1 W(j) = Ftemp2+abs(WD(j,Wnum(s))); end Ftempa = Ftemp1+F(j); Ftempb = Ftemp2+W(j); W(j) = Ftempa/Ftempb*i; elseif TypeM == 3 %FDM: max(F*Wd); Wd: EuclideanD for s = 1:i-1 W(j) = Ftemp2+WD(j,Wnum(s)); end Ftempa = Ftemp1+F(j); Ftempb = Ftemp2+W(j); W(j) = Ftempa/i*Ftempb/(i^2); else %FSQ: max(D/Wd); Wd=Wd^(-1) for s = 1:i-1 W(j) = Ftemp2+WD(j,Wnum(s))^(-1); end Ftempa = Ftemp1+F(j); Ftempb = Ftemp2+W(j); W(j) = Ftempa/Ftempb*i; end end end Ft2 = 0; for j = 1:n if F1(j) == 0 Ft2 = W(j); break end end for j = 1:n if F1(j) == 0 if W(j) >= Ft2 Wnum(i) = j; Ft2 = W(j); end end end F1(Wnum(i)) = 1; Ftemp1 = Ftemp1+F(Wnum(i)); if TypeM == 1 for s = 1:i-1 Ftemp2 = Ftemp2-WD(Wnum(i),Wnum(s)); end elseif TypeM == 2 for s = 1:i-1 Ftemp2 = Ftemp2+WD(Wnum(i),Wnum(s)); end elseif TypeM == 3 for s = 1:i-1 Ftemp2 = Ftemp2+WD(Wnum(i),Wnum(s)); end else for s = 1:i-1 Ftemp2 = Ftemp2+WD(Wnum(i),Wnum(s))^(-1); end end end