-## rtree.R (2008-01-08)
+## rtree.R (2008-01-13)
## Generates Random Trees
phy
}
-rcoal <- function(n, tip.label = NULL, br = rexp, ...)
+rcoal <- function(n, tip.label = NULL, br = "coalescent", ...)
{
nbr <- 2*n - 2
edge <- matrix(NA, nbr, 2)
- x <- br(n - 1, ...) # coalescence times
+ ## coalescence times by default:
+ x <- if (is.character(br)) 2*rexp(n - 1)/(n:2 * (n - 1):1)
+ else br(n - 1, ...)
if (n == 2) {
edge[] <- c(3, 3, 1:2)
edge.length <- rep(x, 2)