function output=assess_vec(input) path(path,'/home/jesse') load /home/jesse/psdstuff/figures/forfig5_m fig5_GO2gene_m load /home/jesse/psdstuff/figures/forfig5_m list_fig5_m_n %load /home/jesse/psdstuff/figures/forfig5_m fig5_numingo_m %load /home/jesse/psdstuff/figures/forfig5_m fig5_gene2GO_m ranks=tiedrank(input); ROC_using_sorted=nan(length(fig5_GO2gene_m),1); for i=1:8907 [TF,loc]=ismember(list_fig5_m_n, fig5_GO2gene_m{i}); if sum(TF)>20 & sum(TF)<1000 ROC_using_sorted(i)=1-(sum(ranks(TF))/sum(TF)-(sum(TF)+1)/2)/sum(~TF); end end output=1-ROC_using_sorted;