a few changes....

[ape.git] / man / rlineage.Rd

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