function output=assess_vec_b(input,whichgenes)
%output=assess_vec_b(input,whichgenes)
%human

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

load /home/jesse/psdstuff/GOgenemat

load /home/jesse/psdstuff/genelist

  GOgenemat=GOgenemat(whichgenes,:);



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

for i=1:10127
  TF=logical(GOgenemat(:,i));
if sum(TF)>=1
ROC_using_sorted(i)=(sum(ranks(TF))/sum(TF)-(sum(TF)+1)/2)/sum(~TF);
end
end

output=ROC_using_sorted;