powered by
This function transforms a variance-covariance matrix into a phylogenetic tree.
vcv2phylo(mat, tolerance = 1e-7)
a square symmetric (positive-definite) matrix.
the numeric tolerance used to compare the branch lengths.
an object of class "phylo".
"phylo"
The function tests if the matrix is symmetric and positive-definite (i.e., all its eigenvalues positive within the specified tolerance).
vcv, corPhyl
vcv
corPhyl
# NOT RUN { 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