getTipDistribution {RPANDA}

Distribution of tip trait values.

Description

Computes the mean and variance of the tip trait distribution under a specified model of trait evolution.

Usage

getTipDistribution(object, params, v)

Arguments

object

an object of class 'PhenotypicModel'

params

vector of parameters, given in the same order as in the 'model' object.

v

boolean specifying the verbose mode. Default value : FALSE.

Value

mean

Expectation vector of the tip trait distribution.

Sigma

Variance-covariance matrix of the tip trait distribution.

Author(s)

M Manceau

References

Manceau M., Lambert A., Morlon H. (2017) A unifying comparative phylogenetic framework including traits coevolving across interacting lineages Systematic Biology

Examples

#Loading an example tree
newick <- "((((A:1,B:0.5):2,(C:3,D:2.5):1):6,E:10.25):2,(F:6.5,G:8.25):3):1;"
tree <- read.tree(text=newick)

#Creating a BM model
modelBM <- createModel(tree, 'BM')

#Tip trait distribution under the model :
getTipDistribution(modelBM, c(0,0,0,1))