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 rcoal() are fine.
Any insights on this?
Cheers,
Santiago
Santiago Claramunt
Museum of Natural Science,
119 Foster Hall,
Louisiana State University,
Baton Rouge, LA70803
scla...@tigers.lsu.edu
http://www.museum.lsu.edu/Claramunt/Home.html
_______________________________________________
R-sig-phylo mailing list
R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
