Hi,
I am attempting to estimate the ancestral states in a continuous character
under a single-optimum OU model. Following your latest post, I used
geiger::fitContinuous(model="OU") and transformed the branch lengths using
ouTree(
)
Unfortunately on putting these into ace I get an error message:
Error in nlm(function(p) dev.BM(p), p = c(1, rep(mean(x), nb.node)), hessian =
TRUE) :
non-finite value supplied by 'nlm'
And when I try using anc.ML I get this error message:
Error in solve.default(C) :
system is computationally singular: reciprocal condition number = 1.75196e-148
It works with my original data but for some reason the transformation of the
branch lengths is giving odd values. The alpha that fitContinuous gives me is
3.831327.
Here are the first few of the original branch lengths:
1] 10.82 26.39 18.01 18.01 14.40
And here are those values transformed using ouTree:
[1] 2.468450e-149 1.638125e-61 1.408962e-01 1.408962e-01 2.057990e-101
Thanks in advance,
Zac
[[alternative HTML version deleted]]
_______________________________________________
R-sig-phylo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
