phyl.RMA {phytools} | R Documentation |
This function performs phylogenetic RMA regression.
phyl.RMA(x, y, tree, method="BM", lambda=NULL, fixed=FALSE, h0=1.0)
x |
vector with names. |
y |
vector with names. |
tree |
a phylogenetic tree in |
method |
method to obtain the correlation structure: can be |
lambda |
value of |
fixed |
logical value indicating whether or not λ should be optimized using likelihood. |
h0 |
null hypothesis for |
Optionally jointly estimates lambda if method="lambda"
. Likelihood optimization of lambda is performed using optimize
on the interval (0,1). Note that some statistician think there is never a condition in which a reduced-major-axis regression should be used.
A list with the following components:
RMA.beta |
a vector of RMA regression coefficients. |
V |
a VCV matrix for the traits. |
lambda |
fitted value of lambda ( |
logL |
log-likelihood ( |
test |
a vector containing results for hypothesis tests on |
resid |
a vector of residuals for |
Liam Revell liam.revell@umb.edu
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.
phyl.cca
, phyl.pca
, phyl.resid