phytools (version 0.2-30)

phyl.resid: Phylogenetic size-correction via GLS regression

Description

This function fits one or multiple phylogenetic regressions (depending on the number of columns in Y) and computes the residuals. Designed for phylogenetic size correction using GLS regression (e.g., Revell 2009; Evolution).

Usage

phyl.resid(tree, x, Y, method="BM")

Arguments

tree

a phylogenetic tree in "phylo" format.

x

vector containing the single independent variable (e.g., size), or matrix with multiple independent variables in columns.

Y

vector or matrix with one or multiple dependent variables in columns.

method

method to obtain the correlation structure: can be "BM" or "lambda".

Value

• a list with the following components:
• beta
  a vector or matrix of regression coefficients.
• resid
  a vector or matrix of residuals for species.
• lambda
  a vector of lambda values (method="lambda" only).
• logL
  a vector of log-likelihoods (method="lambda" only).

Details

Optionally fits $\lambda$ for each regression model. Likelihood optimization of $\lambda$ is performed for method="lambda" using optimize on the interval (0,1). This function is redundant with residuals(gls(...,correlation=corPagel(...))), but may be easier to use to use for some users.

References

Revell, L. J. 2009. Size-correction and principal components for interspecific comparative studies. Evolution, 63,3258--3268. Revell, L. J. 2010. Phylogenetic signal and linear regression on species data. Methods in Ecology and Evolution, 1, 319--329.

See Also

phyl.pca, gls