function res=mrses_orig(A,B,k,Niter,Ncycle,distmod) if (nargin<6) distmod=1; end if (nargin<5) Ncycle=1000; end if (nargin<4) Niter=500; end if (nargin<3) k=5; end if (nargin<2) error('As minimum two matrixes needed for MRSES'); end if (nargin>6) error('Too much parameters'); end sa=size(A);sb=size(B); if (sa(2)==sb(2)) genes=sa(2); else error('Features dimension mismatch'); end optki=zeros(Ncycle,k); for icycle=1:Ncycle % SELECT k GENES {ki} FOR TEST AND EXCLUDE THEM FROM ALL GENES {ke} tt=randperm(genes); ki=tt(1:k); ke=tt(k+1:end); for iter=1:Niter dist1=bmc(A(:,ki),B(:,ki)); xki=ceil(rand(1)*k); xke=ceil(rand(1)*(genes-k)); t=ki(xki); ki(xki)=ke(xke); ke(xke)=t; dist2=bmc(A(:,ki),B(:,ki)); if(dist2(distmod)