-## CDF.birth.death.R (2010-09-27)
+## CDF.birth.death.R (2012-09-14)
-## Functions to simulate and fit
+## Functions to Simulate and Fit
## Time-Dependent Birth-Death Models
-## Copyright 2010 Emmanuel Paradis
+## Copyright 2010-2012 Emmanuel Paradis
## This file is part of the R-package `ape'.
## See the file ../COPYING for licensing issues.
if (!case %in% c(1, 3, 6)) Pi <- Vectorize(Pi)
- denom <- if (fast) integrateTrapeze(Pi, 0, Tmax) else integrate(Pi, 0, Tmax)$value
+ denom <-
+ if (fast) integrateTrapeze(Pi, 0, Tmax)
+ else integrate(Pi, 0, Tmax)$value
n <- length(x)
p <- numeric(n)
if (fast) {
phy <- list(edge = edge, edge.length = edge.length,
tip.label = paste("t", 1:(i + 1), sep = ""), Nnode = i)
class(phy) <- "phylo"
+ attr(phy, "order") <- "cladewise"
phy
}
.makePhylo(edge[1:j, ], edge.length[1:j], i)
}
-drop.fossil <- function(phy, tol = 0)
+drop.fossil <- function(phy, tol = 1e-8)
{
n <- Ntip(phy)
x <- dist.nodes(phy)[n + 1, ][1:n]