function [roc]=GMer2(index,kernel,folds) %[roc]=GMer2(index,{kernel},folds) rand('twister',sum(100*clock)); for i=1:length(kernel) kernel{i}(isnan(kernel{i}))=0; end scores1=sparse(length(index),1); scores2=sparse(length(index),1); path(path,'/home/jesse/psdstuff/GeneMANIACode'); path(path,'/home/jesse'); index=(real(index)); whichers=find(index); sizer=length(whichers); whichers=whichers(randperm(sizer)); output=[]; output2=[]; for i=1:folds whicherout{i}=whichers(floor((i-1)/folds*(sizer)+1):floor(sizer*i/folds)); end for i=1:folds newindex=index; newindex(whicherout{i})=0; newindex2=newindex; newindex2(newindex2==0)=-1; [p, k, wts]=predictClassesCG(newindex2,kernel); p(logical(newindex))=nan; p=tiedrank(p); output=[output;p(whicherout{i})]; end roc=sum(output(:))/folds; np=sum(index)/folds; nn=length(index)-sum(index); roc=(roc/np-(np+1)/2)/nn;