R/T2.R

In chengvt/cheng: cheng

#' @export
T2 <- function(prcomp, k){
    scores <- prcomp$x
    if (k == 1) { lambda <- prcomp$sdev[1:k]^2
    } else lambda <- diag(prcomp$sdev[1:k]^2)
    
    T2 <- vector() #initialize
    for (i in 1:nrow(scores)){
        Ti <- scores[i,1:k, drop = FALSE]
        Ti2 <- Ti %*% lambda %*% t(Ti)
        T2[i] <- Ti2 
    }
    return(T2)
}