Phylogenetic reduced major axis regression (RMA).

Input

result<-phyl.RMA(x,y,tree,method="BM")

"tree" is a phylogenetic tree in "phylo" format;
"x" and "y" are data vectors with names;
"method" is the method to obtain the correlation structure for the model: can be "BM" or "lambda".

Output

phyl.pca() returns a list with the following elements:
"RMA.beta" a 2x1 vector of regression coefficients (intercept and slope);
"V" a VCV matrix for the traits;
"lamba" fitted version of lambda (or 1.0, if method="BM");
"logL" log-likelihood;
"resid" a vector of residuals for y given x.

Citation

Function written by Liam J. Revell 2010 (modified 2011).