R/PEVMAX.R

Defines functions PEVMAX

Documented in PEVMAX

PEVMAX <-
  function(Train,Test, P, lambda=1e-5,C=NULL){
    PTrain<-P[rownames(P)%in%Train,]
    PTest<-P[rownames(P)%in%Test,]
     if (length(Test)==1){PTest=matrix(PTest, nrow=1)}
    if (!is.null(C)){ PTest<-C%*%PTest}
    PEVmax<-max(diag(PTest%*%solve(crossprod(PTrain)+lambda*diag(ncol(PTrain)),t(PTest))))
    return(PEVmax)
  }

Try the STPGA package in your browser

Any scripts or data that you put into this service are public.

STPGA documentation built on May 29, 2017, 3:44 p.m.