some bug fixes and '...' in rTrait*()

[ape.git] / man / rlineage.Rd

diff --git a/man/rlineage.Rd b/man/rlineage.Rd
index d4412454cf78c7a056d5cd965d0eef01e1ecfae4..b2b1aed354c2fb780689216142543d5ef444dc22 100644 (file)
--- a/man/rlineage.Rd
+++ b/man/rlineage.Rd
@@ -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}