R/dbRDA.R

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

`dbRDA` <-
    function(D, X = NULL, Z = NULL)
{
    SOL <- PCoA(D)
    pos <- SOL$values > sqrt(.Machine$double.eps)
    Lambda <- SOL$values[pos]
    Y <- SOL$vectors[, pos] %*% diag(sqrt(Lambda))
    pRDA(Y, X, Z)
}