a collection of bug fixes

[ape.git] / R / rtree.R

diff --git a/R/rtree.R b/R/rtree.R
index f331ca3101f0472575ef35b6e92f3abf84021764..6829b02d294be46753894cea6f074f5e09bfbd3d 100644 (file)
--- a/R/rtree.R
+++ b/R/rtree.R
@@ -106,7 +106,7 @@ rcoal <- function(n, tip.label = NULL, br = "coalescent", ...)
     nbr <- 2*n - 2
     edge <- matrix(NA, nbr, 2)
     ## coalescence times by default:
-    x <- if (is.character(br)) 2*rexp(n - 1)/(n:2 * (n - 1):1)
+    x <- if (is.character(br)) 2*rexp(n - 1)/(as.double(n:2) * as.double((n - 1):1))
     else br(n - 1, ...)
     if (n == 2) {
         edge[] <- c(3L, 3L, 1:2)