function output=assess_vec_KEGG_b(input,whichgenes)


  whichgenes=logical(whichgenes);
path(path,'/home/jesse')

load /home/jesse/psdstuff/KEGGgenemat



  KEGGgenemat=KEGGgenemat(whichgenes,:);


ranks=tiedrank(input);
ROC_using_sorted=nan(size(KEGGgenemat,2),1);

for i=1:217
  TF=logical(KEGGgenemat(:,i));
if sum(TF)>20 & sum(TF)<1000
ROC_using_sorted(i)=(sum(ranks(TF))/sum(TF)-(sum(TF)+1)/2)/sum(~TF);
end
end

output=ROC_using_sorted;