R/AOPT.R

In STPGA: Selection of Training Populations by Genetic Algorithm

AOPT <-
  function(Train,Test=NULL, P, lambda=1e-5, C=NULL){
    PTrain<-P[rownames(P)%in%Train,]
    if (!is.null(C)){
    D<-mean(diag(C%*%solve(crossprod(PTrain)+lambda*diag(ncol(PTrain)), t(C))))
    }else {
      D<-mean(diag(solve(crossprod(PTrain)+lambda*diag(ncol(PTrain)))))
    }
    
    return(D)
  }