function output=assess_vec_b(input,whichgenes)


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

load /home/jesse/psdstuff/OMIMgenemat

load /home/jesse/psdstuff/genelist

  OMIMgenemat=OMIMgenemat(whichgenes,:);



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

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

output=ROC_using_sorted;