evolvcv.lite | R Documentation |
Fits a multi-regime model for differences in the evolutionary variance-covariance structure or correlation between two continuous traits, following Revell & Collar (2009) and Revell et al. (2022).
evolvcv.lite(tree, X, maxit=2000, tol=1e-10, ...)
tree |
an object of class |
X |
an |
maxit |
an optional integer value indicating the maximum number of iterations for optimization - may need to be increased for large trees. |
tol |
tolerance value for |
... |
other optional arguments. The most important optional argument is probably |
This function takes an object of class "simmap"
with a mapped binary or multi-state trait and data for two and only two continuously valued character. It then fits (by default) four different evolutionary models: common rates and correlation; different rates, common correlation; different correlations, common rates; no common structure.
In addition to the four default models specified above, evolvcv.lite
now fits an additional four additional models.
The set of models to be fit can be specified using the optional argument models
in multiple ways.
First, if left unspecified, then the four models listed above will be fit.
Second, if models
is set to "all models"
than eight models will be fit.
Lastly, one or more (up to all eight) models can be fit by encoding the models to be fit into a single vector containing a subset or all of the following elements: "1"
, "2"
, "2b"
, "3"
, "3b"
, "3c"
, and "4"
. These codes correspond to the following eight models: 1. common rates, common correlation; 2. different rates, common correlation; 2b. different rates for trait 1 only, common correlation; 2c. different rates for trait 2 only, common correlation; 3. common rates, different correlations; 3b. different rates for trait 1 only, different correlations; 3c. different rates for trait 2 only, different correlation; and 4. no common structure.
A list with the results summarized for each model.
Liam Revell liam.revell@umb.edu
Revell, L. J. (2024) phytools 2.0: an updated R ecosystem for phylogenetic comparative methods (and other things). PeerJ, 12, e16505.
Revell, L. J., and D. C. Collar (2009) Phylogenetic analysis of the evolutionary correlation using likelihood. Evolution, 63, 1090-1100.
Revell, L. J., K. S. Toyama, and D. L. Mahler (2022) A simple hierarchical model for heterogeneity in the evolutionary correlation on a phylogenetic tree. PeerJ, 10, e13910.
brownie.lite
, evol.vcv
## load data from Revell & Collar (2009)
data(sunfish.tree)
data(sunfish.data)
## fit heirarchical common-structure models
sunfish.fit<-evolvcv.lite(sunfish.tree,
sunfish.data[,2:3],models="all models")
## print fitted models
print(sunfish.fit)
## compare models
anova(sunfish.fit)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.