Nothing
source("helper-diversitree.R")
## Will be useful to have a function for testing tolerance to within
## 1e-7, as that works out to be how accurate most things actually
## are.
equals7 <- function(...)
equals(..., tolerance=1e-7)
context("GeoSSE (time-dependent)")
test_that("geosse-t", {
set.seed(4)
pars <- c(0.1, 0.2, 0.001, 0.03, 0.03, 0.01, 0.01)
phy <- tree.geosse(pars, max.t=30, x0=0)
## Plain GeoSSE function for comparison:
lik.base <- make.geosse(phy, phy$tip.state)
control.0 <- list()
control.d <- list(backend="deSolve")
control.g <- list(backend="gslode", compiled=FALSE)
control.G <- list(backend="gslode", compiled=TRUE)
functions <- rep(c("sigmoid.t", "constant.t"), c(3, 4))
lik.0 <- make.geosse.t(phy, phy$tip.state, functions, control=control.0)
lik.d <- make.geosse.t(phy, phy$tip.state, functions, control=control.d)
lik.g <- make.geosse.t(phy, phy$tip.state, functions, control=control.g)
lik.G <- make.geosse.t(phy, phy$tip.state, functions, control=control.G)
## First, evaluate the functions with no time effect and check that they
## are the same as the base GeoSSE model
t <- max(branching.times(phy))/2
p.t <- c(pars[1], pars[1], t, .5,
pars[2], pars[2], t, .5,
pars[3], pars[3], t, .5,
pars[4:7])
ll <- lik.base(pars)
expect_that(ll, equals(-293.082650678509))
expect_that(lik.0(p.t), equals(ll))
expect_that(lik.d(p.t), equals7(ll))
expect_that(lik.g(p.t), equals(ll))
expect_that(lik.G(p.t), equals(ll))
## With some random parameters, still no time dependence:
set.seed(1)
pars2 <- pars + runif(7, 0, .2)
p2.t <- c(pars2[1], pars2[1], t, .5,
pars2[2], pars2[2], t, .5,
pars2[3], pars2[3], t, .5,
pars2[4:7])
ll <- lik.base(pars2)
expect_that(ll, equals(-344.946604073163))
## And expect that the time dependent versions agree again:
expect_that(lik.0(p2.t), equals(ll))
expect_that(lik.d(p2.t), equals7(ll))
expect_that(lik.g(p2.t), equals(ll))
expect_that(lik.G(p2.t), equals(ll))
## Add some time dependence:
p3.t <- p2.t
p3.t[2] <- p3.t[2] + .2
ll <- lik.0(p3.t)
expect_that(ll, equals(-346.858778360156))
expect_that(lik.d(p3.t), equals7(ll))
expect_that(lik.g(p3.t), equals(ll))
expect_that(lik.G(p3.t), equals(ll))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.