R/CA.R

In barebone: Multivariate Methods in Ecology with R: Basic Algorithms

`CA` <-
    function(Y)
{
    Y <- Y/sum(Y)
    r <- rowSums(Y)
    c <- colSums(Y)
    Ybar <- (Y - r %o% c)/sqrt(r %o% c)
    SOL <- svd(Ybar)
    SOL$u <- diag(1/sqrt(r)) %*% SOL$u
    SOL$v <- diag(1/sqrt(c)) %*% SOL$v
    SOL
}