Note to self: sim_island always returns an uninteresting result?

[richelbilderbeek]

Will continue on this tomorrow...

test_that("use", {
  set.seed(
    1,
    kind = "Mersenne-Twister",
    normal.kind = "Inversion",
    sample.kind = "Rejection"
  )
  mainland <- DAISIEmainland:::sim_mainland(
    total_time = 1,
    m = 10,
    mainland_ex = 2)
  plot_mainland(mainland = mainland)

  count_extant_mainland_species <- function(mainland) {
    purrr::map_dbl(mainland, function(x) { return(sum(x$spec_ex_t == 1.0)) } )
  }
  clade_with_most_species <- which(
    count_extant_mainland_species(mainland) ==
      max(count_extant_mainland_species(mainland))
  )

  mainland_clade <- mainland[[clade_with_most_species]]
  cladogenesis_rate <- 100.0
  extinction_rate <- 100.0
  carrying_capacity <- 100.0
  immigration_rate <- 100.0
  anagenesis_rate <- 100.0

  island_pars <- c(
    cladogenesis_rate,
    extinction_rate,
    carrying_capacity,
    immigration_rate,
    anagenesis_rate
  )
  mainland_sample_type <- "unsampled"
  mainland_sample_type <- "undiscovered"
  mainland_sample_type <- "complete"
  testit::assert(mainland_sample_type == "unsampled" ||
                   mainland_sample_type == "undiscovered" ||
                   mainland_sample_type == "complete")

  island <- sim_island(
    total_time = 100,
    island_pars = island_pars,
    mainland_clade = mainland_clade,
    mainland_sample_prob = 1,
    mainland_sample_type = mainland_sample_type)
  island
  plot_island(island)
})

[richelbilderbeek]

Nope, see test-plot_ideal_island.R for plenty of interesting results.