R/fitted.procrustes.R

`fitted.procrustes` <-
  function(object, truemean = TRUE, ...)
{
  fit <- object$Yrot
  if (truemean)
    fit <- sweep(fit, 2, object$xmean, "+")
  fit
}

## Like above, except when takes 'newata'

`predict.procrustes` <-
    function(object, newdata, truemean = TRUE, ...)
{
    if (missing(newdata))
        return(fitted(object, truemean = truemean))
    if (object$symmetric)
        stop(gettextf("'predict' not available for symmetric procrustes analysis with 'newdata'"))
    Y <- as.matrix(newdata)
    ## scaling and rotation
    Y <- object$scale * Y %*% object$rotation
    ## translation: always
    Y <- sweep(Y, 2, object$translation, "+")
    if (!truemean)
        Y <- sweep(Y, 2, object$xmean, "-")
    Y
}

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.