Morning all,
Whenever I try to use ape's (v 3.0-2) birthdeath function e.g.
> birthdeath(rcoal(sample(10:100, 1)))
I get the following error:
Error in while (foo(up) < 0) up <- up + inc :
missing value where TRUE/FALSE needed
The warnings all appear to be this sort of thing:
1: In log(1 - a) : NaNs produced
2: In log(exp(r * x[2:N]) - a) : NaNs produced
3: In nlm(function(p) dev(p[1], p[2]), c(0.1, 0.2), hessian = TRUE) :
NA/Inf replaced by maximum positive value
and seem to be generated from this line:
out <- nlm(function(p) dev(p[1], p[2]), c(0.1, 0.2), hessian = TRUE)
Is there a fix/workaround for this?
Kind regards,
Simon
_______________________________________________
R-sig-phylo mailing list
R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
