R/PEVMEAN0.R

In STPGA: Selection of Training Populations by Genetic Algorithm

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