"sagatNumRepsAUCEval" <-
function(data.mat, U.mat, reps.vec=c(1,2,5,15), nDE.sca=1122, nEG.sca=NULL, truth.vec=as.numeric(rank(-1*abs(prostate.data.EG.eklund.PLM$coefficients))<=nDE.sca))
{
# 2008.12.8:  Test dilution of compendium to determine minimal number of arrays needed for SAGAT performance

  m <- ncol(data.mat)
  if(is.null(nEG.sca)) nEG.sca <- round(ncol(U.mat)/2)
  results.lis <- list()

  for(k in reps.vec) {
    ind.vec <- seq(1, m-(m%%k)-k+1, k)
    auc.vec <- c()

    for(i in 1:length(ind.vec)) {
      scores.vec <- sagat(data.mat[,ind.vec[i]:(ind.vec[i]+k-1), drop=FALSE],U.mat[,1:nEG.sca],design=NULL,contrasts=1)
      auc.vec <- c(auc.vec, performance(prediction(predictions=abs(scores.vec),labels=truth.vec),measure="auc",fpr.stop=1)@y.values[[1]])
      cat(paste(ind.vec[i]," ",sep=""))
    }
    cat("\n")

    results.lis <- c(results.lis, list(auc.vec))
  }

  return(results.lis)
}