Dear R Phylo List, I'm trying to test an Early Burst model of evolution against a tree using the 'fitContinuous' function in 'geiger'. However, I get the following error message, and am unsure as to what I can do to fix this:
> fitContinuous(tree,trait,model="EB") Fitting EB model: Error in solve.default(phyvcv) : Lapack routine dgesv: system is exactly singular: U[4,4] = 0 My tree is based on fossil taxa and my branch lengths are calculated using time (in millions of years). I seem to be able to test every other model, only EB is giving me troubles. I'm pretty sure it has something to do with the branch lengths, cause if I used other measures, such as all=1 or grafen, I don't get an error. Has anybody had this issue before? Any assistance would be greatly appreciated. Cheers, -Nic- ---- Nicolás E. Campione M.Sc., Ph.D. Candidate Dept. Ecology & Evolutionary Biology University of Toronto 25 Wilcocks St. Toronto, ON Canada M5S 3B2 Royal Ontario Museum 100 Queen's Park Toronto, ON Canada M5S 2C6 Office: 416-586-5591 Email: [email protected] [[alternative HTML version deleted]]
_______________________________________________ R-sig-phylo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
