new image.DNAbin()

[ape.git] / man / rlineage.Rd

diff --git a/man/rlineage.Rd b/man/rlineage.Rd
index d4412454cf78c7a056d5cd965d0eef01e1ecfae4..5b2056808ebfd31e8e1070fb82e15da3004c4cd2 100644 (file)
--- a/man/rlineage.Rd
+++ b/man/rlineage.Rd
@@ -15,7 +15,7 @@ rlineage(birth, death, Tmax = 50, BIRTH = NULL,
          DEATH = NULL, eps = 1e-6)
 rbdtree(birth, death, Tmax = 50, BIRTH = NULL,
         DEATH = NULL, eps = 1e-6)
-drop.fossil(phy, tol = 0)
+drop.fossil(phy, tol = 1e-8)
 }
 \arguments{
   \item{birth, death}{a numeric value or a (vectorized) function
@@ -49,8 +49,8 @@ drop.fossil(phy, tol = 0)
   process. \emph{Annals of Mathematical Statistics}, \bold{19}, 1--15.
 
   Paradis, E. (2010) Time-dependent speciation and extinction from
-  phylogenies: a least squares approach. (under revision)
-  %\emph{Evolution}, \bold{59}, 1--12.
+  phylogenies: a least squares approach. \emph{Evolution} (in press)
+  %, \bold{59}, 1--12.
 }
 \author{Emmanuel Paradis}
 \seealso{
@@ -60,9 +60,16 @@ drop.fossil(phy, tol = 0)
 \examples{
 plot(rlineage(0.1, 0)) # Yule process with lambda = 0.1
 plot(rlineage(0.1, 0.05)) # simple birth-death process
-b <- function(t) 1/(1 + exp(0.1*t - 2)) # logistic
-layout(matrix(1:2, 1))
-plot(rlineage(b, 0.01))
-plot(rbdtree(b, 0.01))
+b <- function(t) 1/(1 + exp(0.2*t - 1)) # logistic
+layout(matrix(0:3, 2, byrow = TRUE))
+curve(b, 0, 50, xlab = "Time", ylab = "")
+mu <- 0.07
+segments(0, mu, 50, mu, lty = 2)
+legend("topright", c(expression(lambda), expression(mu)),
+       lty = 1:2, bty = "n")
+plot(rlineage(b, mu), show.tip.label = FALSE)
+title("Simulated with 'rlineage'")
+plot(rbdtree(b, mu), show.tip.label = FALSE)
+title("Simulated with 'rbdtree'")
 }
 \keyword{datagen}