-## compute.brtime.R (2011-07-15)
+## compute.brtime.R (2011-07-26)
## Compute and Set Branching Times
## x: branching times (aka, node ages or heights)
if (identical(method, "coalescent")) { # the default
- x <- 2 * rexp(n - 1)/(as.double(n:2) * as.double((n - 1):1))
- if (is.null(force.positive))
- force.positive <- TRUE
+ x <- 2 * rexp(m)/(as.double((m + 1):2) * as.double(m:1))
+ ## x <- 2 * rexp(n - 1)/(as.double(n:2) * as.double((n - 1):1))
+ if (is.null(force.positive)) force.positive <- TRUE
} else if (is.numeric(method)) {
x <- as.vector(method)
if (length(x) != m)