R/fitted.rda.R

`fitted.rda` <-
    function (object, model = c("CCA", "CA", "pCCA"), type = c("response", "working"), ...) 
{
    type <- match.arg(type)
    model <- match.arg(model)
    if (is.null(object[[model]]))
        stop("component ", model, " does not exist")
    if (model == "pCCA")
        Xbar <- object$pCCA$Fit
    else
        Xbar <- object[[model]]$Xbar
    if (model == "CCA") 
        Xbar <- qr.fitted(object$CCA$QR, Xbar)
    if (type == "response") {
        cent <- attr(Xbar, "scaled:center")
        scal <- attr(Xbar, "scaled:scale")
        if (!is.null(scal)) {
            Xbar <- sweep(Xbar, 2, scal, "*")
            attr(Xbar, "scaled:scale") <- NULL
        }
        Xbar <- sweep(Xbar, 2, cent, "+")
        attr(Xbar, "scaled:center") <- NULL
    } else {
        Xbar <- Xbar/sqrt(nrow(Xbar)-1)
    }
    Xbar
}

Try the vegan package in your browser

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

vegan documentation built on May 2, 2019, 5:51 p.m.