R/Min.Max.Multivar.PCA.R

Defines functions Min.Max.Multivar.PCA

Documented in Min.Max.Multivar.PCA

Min.Max.Multivar.PCA <- function(gamma, Sigma_SS, Sigma_T0T0, Sigma_T1T1){ 

  Max.PCA <- min(t(gamma) %*% Sigma_SS %*% gamma / ( (sqrt(Sigma_T0T0) - sqrt(Sigma_T1T1)   )**2), 1) 
  Min.PCA <- max(t(gamma) %*% Sigma_SS %*% gamma / ( (sqrt(Sigma_T0T0) + sqrt(Sigma_T1T1)   )**2), 0)
  
  cat("\n \nMin PCA: ", Min.PCA)
  cat("\n \nMax PCA: ", Max.PCA, "\n \n")  
  
fit <- 
    list(Call=match.call())
  
  class(fit) <- "Min.Max.Multivar.PCA"
  fit
}

Try the EffectTreat package in your browser

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

EffectTreat documentation built on July 8, 2020, 7:17 p.m.