Nothing
### * Setup
n_cores <- min(2, parallel::detectCores())
n_chains <- max(n_cores, 2)
run_mcmc <- function(...) {
isotracer:::run_mcmc(..., cores = n_cores, chains = n_chains)
}
new_networkModel <- function() {
isotracer::new_networkModel(quiet = TRUE)
}
### * Ops.topology()
test_that("Ops.topology() works", {
expect_error({topo1 <- make_topology(links = "NH4, NO3 -> epi -> pseph, tricor")},
NA)
# A larger foodweb
links <- c("NH4, NO3 -> seston, epi, CBOM, FBOM",
"seston -> lepto", "epi -> petro, pseph",
"CBOM, FBOM -> eudan", "CBOM -> phyllo",
"FBOM -> tricor -> arg, euthy")
expect_error({topo2 <- make_topology(links = links, split = "epi")},
NA)
# Tests
expect_false(topo1 == topo2)
expect_true(topo1 != topo2)
# Using a data frame to specify the links
links <- data.frame(source = c("NH4", "NO3", "epi"),
consumer = c("epi", "epi", "petro"))
expect_error({topo3 <- make_topology(links, from = "source", to = "consumer")},
NA)
expect_false(topo1 == topo3)
expect_true(topo3 == topo3)
expect_length(unique(list(topo1, topo2, topo1, topo3, topo1)), 3)
topo4 <- topo1
attr(topo4, "split") <- "epi"
expect_false(topo4 == topo1)
expect_length(unique(list(topo1, topo2, topo1, topo3, topo1, topo4)), 4)
})
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.