1
function res=mrses_orig(A,B,k,Niter,Ncycle,distmod)
16
error('As minimum two matrixes needed for MRSES');
19
error('Too much parameters');
23
sa=size(A);sb=size(B);
27
error('Features dimension mismatch');
31
optki=zeros(Ncycle,k);
33
% SELECT k GENES {ki} FOR TEST AND EXCLUDE THEM FROM ALL GENES {ke}
40
dist1=bmc(A(:,ki),B(:,ki));
43
xke=ceil(rand(1)*(genes-k));
49
dist2=bmc(A(:,ki),B(:,ki));
51
if(dist2(distmod)<dist1(distmod)) % COMPARES BHATA DISTANCES 1-bhata, 2-mahal, 3-corr
60
optki=reshape(optki,1,[]);
61
[n,g]=hist(optki,1:genes);
63
res=flipud(sortrows(H'));
67
function [rbhata,rmahal,rcorr]=bmc(x1,x2)
76
rmahal=((m2-m1)/c)*(m2-m1)';
77
rcorr=2.*log(det(c)./sqrt(det(c1).*det(c2)));
78
rbhata=rmahal./8+rcorr./4;