This function transforms a variance-covariance matrix into a phylogenetic tree.
vcv2phylo(mat, tolerance = 1e-7)
an object of class "phylo"
.
a square symmetric (positive-definite) matrix.
the numeric tolerance used to compare the branch lengths.
Simon Blomberg
The function tests if the matrix is symmetric and positive-definite (i.e., all its eigenvalues positive within the specified tolerance).
vcv
, corPhyl
tr <- rtree(10)
V <- vcv(tr) # VCV matrix assuming Brownian motion
z <- vcv2phylo(V)
identical(tr, z) # FALSE
all.equal(tr, z) # TRUE
Run the code above in your browser using DataLab