makeStandardizeF <- function(X) {
  if (missing(X)) {
    cat("Usage:
         standardize <- makeStandardizeF(X)  ## X is nSamples x nDimensions
         Xs <- standardize(X)
         X2s <- standardize(X2)\n")
    return(invisible())
  }
  ## X is nSamples x nDimensions
  mu <- colMeans(X)
  sigma <- sd(X) ##sd should be named colSds

  function(newX) {
    nr <- nrow(newX)
    nc <- ncol(newX)
    (newX - matrix(mu,nr,nc,byrow=TRUE)) / matrix(sigma,nr,nc,byrow=TRUE)
  }
}

combine <- function(results,parameterColumns,valueColumns,func) {
  uniqueCombos <- unique(results[,parameterColumns,drop=FALSE])
  comboResult <- c()
  for (uci in 1:nrow(uniqueCombos)) {
    parameters <- uniqueCombos[uci,]
    mask <- apply(results[,parameterColumns,drop=FALSE], 1,
                  function(ps) all(ps==parameters))
    comboResult <- rbind(comboResult,
                         c(parameters,
                           func(results[mask,valueColumns,drop=FALSE])))
  }
  comboResult
}