Nothing
test_that("sampled stt, 1 type, no geodynamics, oceanic island (same arguments
as geodynamics, 5 pars)", {
pars <- c(0.5, 0.1, 10, 1, 0.5)
time <- 1
mainland_n <- 1
verbose <- FALSE
sample_freq <- 1
area_pars <- create_area_pars(
max_area = 1,
current_area = 1,
proportional_peak_t = 0,
total_island_age = 0,
sea_level_amplitude = 0,
sea_level_frequency = 0,
island_gradient_angle = 0)
hyper_pars <- create_hyper_pars(d = 0, x = 0)
nonoceanic_pars <- c(0, 0)
set.seed(1)
island_replicates <- list()
out <- list()
out[[1]] <- DAISIE_sim_core_cr(
time = time,
pars = pars,
mainland_n = mainland_n,
area_pars = area_pars,
hyper_pars = hyper_pars,
nonoceanic_pars = nonoceanic_pars
)
island_replicates[[1]] <- out
testthat::expect_silent(
formatted_CS_sim <- DAISIE_format_CS_sampled_stt(
island_replicates = island_replicates,
time = time,
M = mainland_n,
sample_freq = sample_freq,
verbose = verbose
)
)
expected_CS_format <- list()
expected_CS_format[[1]] <- list()
stt_all <- matrix(ncol = 5, nrow = 2)
colnames(stt_all) <- c("Time", "nI", "nA", "nC", "present")
stt_all[1, ] <- c(1, 0, 0, 0, 0)
stt_all[2, ] <- c(0, 0, 0, 3, 1)
expected_CS_format[[1]][[1]] <- list(island_age = 1,
not_present = 0,
stt_all = stt_all)
expected_CS_format[[1]][[2]] <- list(branching_times = c(1.000000000000000,
0.244818166871655,
0.173128288990374,
0.029668240213840),
stac = 2,
missing_species = 0)
testthat::expect_equal(formatted_CS_sim, expected_CS_format)
})
test_that("sampled stt, 1 type, geodynamics, oceanic island (same arguments as
no geodynamics, 5 pars)", {
total_time <- 5
mainland_n <- 1
verbose <- FALSE
sample_freq <- 25
set.seed(3)
island_replicates <- list()
out <- list()
ext_pars <- c(1, 100)
island_ontogeny <- 1
sea_level <- 0
pars <- c(0.0001, 2.2, 0.005, 1, 1)
area_pars <- create_area_pars(
max_area = 5000,
current_area = 2500,
proportional_peak_t = 0.5,
total_island_age = 15,
sea_level_amplitude = 0,
sea_level_frequency = 0,
island_gradient_angle = 0
)
hyper_pars <- create_hyper_pars(d = 0.2, x = 0.1)
peak <- calc_peak(total_time = total_time,
area_pars = area_pars)
Amax <- get_global_max_area(total_time = total_time,
area_pars = area_pars,
peak = peak,
island_ontogeny = island_ontogeny,
sea_level = sea_level)
Amin <- get_global_min_area(total_time = total_time,
area_pars = area_pars,
peak = peak,
island_ontogeny = island_ontogeny,
sea_level = sea_level)
nonoceanic_pars <- c(0, 0)
out[[1]] <- DAISIE_sim_core_time_dep(
time = total_time,
pars = pars,
mainland_n = mainland_n,
island_ontogeny = island_ontogeny,
sea_level = sea_level,
hyper_pars = hyper_pars,
area_pars = area_pars,
peak = peak,
Amax = Amax,
Amin = Amin,
nonoceanic_pars = nonoceanic_pars
)
island_replicates[[1]] <- out
testthat::expect_silent(
formatted_CS_sim <- DAISIE_format_CS_sampled_stt(
island_replicates = island_replicates,
time = total_time,
M = mainland_n,
sample_freq = sample_freq,
verbose = verbose
)
)
testthat::expect_equal(
formatted_CS_sim[[1]][[1]]$island_age,
5
)
testthat::expect_equal(
formatted_CS_sim[[1]][[1]]$not_present,
0
)
testthat::expect_equal(
formatted_CS_sim[[1]][[1]]$stt_all[5, ],
c(Time = 4.2, nI = 0.0, nA = 0.0, nC = 0.0, present = 0.0)
)
testthat::expect_equal(
formatted_CS_sim[[1]][[1]]$stt_all[12, ],
c(Time = 2.8, nI = 1.0, nA = 0.0, nC = 0.0, present = 1.0)
)
testthat::expect_equal(
formatted_CS_sim[[1]][[1]]$stt_all[25, ],
c(Time = 0.2, nI = 0.0, nA = 2.0, nC = 0.0, present = 1.0)
)
testthat::expect_equal(
formatted_CS_sim[[1]][[2]]$branching_times,
c(5.00000000000000, 0.57929405393928002)
)
testthat::expect_equal(
formatted_CS_sim[[1]][[2]]$stac,
2
)
testthat::expect_equal(
formatted_CS_sim[[1]][[2]]$missing_species,
0
)
})
test_that("sampled stt, 2 type, no geodynamics, oceanic island (same arguments
as geodynamics, 10 pars)", {
pars <- c(0.4, 0.1, 10, 1, 0.5, 0.4, 0.1, 10, 1, 0.5)
total_time <- 1
M <- 10
mainland_n <- M
verbose <- FALSE
replicates <- 1
sample_freq <- 25
area_pars <- create_area_pars(
max_area = 1,
current_area = 1,
proportional_peak_t = 0,
total_island_age = 0,
sea_level_amplitude = 0,
sea_level_frequency = 0,
island_gradient_angle = 0)
hyper_pars <- create_hyper_pars(d = 0, x = 0)
set.seed(1)
island_replicates <- list()
prop_type2_pool <- 0.4
island_replicates <- DAISIE_sim_min_type2(
time = total_time,
M = M,
pars = pars,
replicates = replicates,
prop_type2_pool = prop_type2_pool,
area_pars = area_pars,
hyper_pars = hyper_pars,
verbose = FALSE
)
testthat::expect_silent(
formatted_CS_sim <- DAISIE_format_CS_sampled_stt(
island_replicates = island_replicates,
time = total_time,
M = mainland_n,
sample_freq = sample_freq,
verbose = verbose
)
)
testthat::expect_equal(
names(formatted_CS_sim[[1]][[1]]),
c(
"island_age",
"not_present_type1",
"not_present_type2",
"stt_all",
"stt_type1",
"stt_type2"
)
)
testthat::expect_equal(
names(formatted_CS_sim[[1]][[2]]),
c(
"branching_times",
"stac",
"missing_species",
"all_colonisations",
"type1or2"
)
)
testthat::expect_equal(
length(formatted_CS_sim[[1]]), 8
)
# Sampled STT has the correct size
testthat::expect_equal(
dim(formatted_CS_sim[[1]][[1]]$stt_all), c(26, 5)
)
testthat::expect_equal(
dim(formatted_CS_sim[[1]][[1]]$stt_type1), c(26, 5)
)
testthat::expect_equal(
dim(formatted_CS_sim[[1]][[1]]$stt_type2), c(26, 5)
)
})
test_that("sampled stt, 1 type, no geodynamics, nonoceanic (same arguments as
geodynamics, 5 pars)", {
total_time <- 5
island_age <- 0.4
clado_rate <- 2.550687345 # cladogenesis rate
ext_rate <- 2.683454548 # extinction rate
clade_carr_cap <- 10.0 # clade-level carrying capacity
imm_rate <- 0.00933207 # immigration rate
ana_rate <- 1.010073119 # anagenesis rate
pars <- c(clado_rate, ext_rate, clade_carr_cap, imm_rate, ana_rate)
area_pars <- create_area_pars(
max_area = 1,
current_area = 1,
proportional_peak_t = 0,
total_island_age = 0,
sea_level_amplitude = 0,
sea_level_frequency = 0,
island_gradient_angle = 0)
hyper_pars <- create_hyper_pars(d = 0, x = 0)
nonoceanic_pars <- c(0.1, 0.9)
sample_freq <- 25
mainland_n <- 1
verbose <- FALSE
set.seed(1)
island_replicates <- list()
out <- list()
out[[1]] <- DAISIE_sim_core_cr(
time = total_time,
pars = pars,
mainland_n = mainland_n,
area_pars = area_pars,
hyper_pars = hyper_pars,
nonoceanic_pars = nonoceanic_pars
)
island_replicates[[1]] <- out
testthat::expect_silent(
formatted_CS_sim <- DAISIE_format_CS_sampled_stt(
island_replicates = island_replicates,
time = total_time,
M = mainland_n,
sample_freq = sample_freq,
verbose = verbose
)
)
})
test_that("sampled stt, 1 type, no geodynamics, oceanic (same arguments as
geodynamics, 5 pars) - 3 replicates", {
pars <- c(0.4, 0.2, 10, 2, 0.5)
total_time <- 5
mainland_n <- 2
verbose <- FALSE
sample_freq <- 25
area_pars <- create_area_pars(
max_area = 1,
current_area = 1,
proportional_peak_t = 0,
total_island_age = 0,
sea_level_amplitude = 0,
sea_level_frequency = 0,
island_gradient_angle = 0)
hyper_pars <- create_hyper_pars(d = 0, x = 0)
nonoceanic_pars <- c(0, 0)
set.seed(1)
replicates <- 3
island_replicates <- list()
for (rep in 1:replicates) {
island_replicates[[rep]] <- list()
full_list <- list()
out <- list()
for (m_spec in 1:mainland_n) {
out$branching_times <- c(10)
while (length(out$branching_times) == 1) {
out <- DAISIE_sim_core_cr(
time = total_time,
mainland_n = 1,
pars = pars,
area_pars = area_pars,
hyper_pars = hyper_pars,
nonoceanic_pars = nonoceanic_pars
)
}
full_list[[m_spec]] <- out
}
island_replicates[[rep]] <- full_list
}
testthat::expect_silent(
formatted_CS_sim <- DAISIE_format_CS_sampled_stt(
island_replicates = island_replicates,
time = total_time,
M = mainland_n,
sample_freq = sample_freq,
verbose = verbose
)
)
})
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.