Hi all,
I'm having problems computing the inverse of the expected variance-
covariance evolutionary matrix (C) on simulated trees generated with
birthdeath.tree().
solve(vcv.phylo(birthdeath.tree(1,0, taxa.stop=10)))
These matrices seem to be singular!
Matrices generated using rtree() and
Hi Santiago,
The problem is that this function, birth.death(,,taxa.stop=X), simulates
speciation and extinction under a Yule process until it gets to taxa=X
and then stops, without adding any length to the branches leading to the
last two taxa. This means that the distances from the root of