function [Y]=snormalise(X,sX); 
% median quantile normalisation of a to-be-predicted sample X against
% training sample sX

 m=median(sX,2); [m,i]=sort(m); Y=X; [N,n]=size(Y); 
 for j=1:n
   [a,i]=sort(X(:,j));
   Y(i,j)=m;
 end;