R/immig_event.R

In joshwlambert/DAISIEmainland: Simulates Data for Phylogenetic Data on Islands for DAISIE with the Addition of Mainland Dynamics

#' Updates the island_tbl when an immigration event occurs
#'
#' @inheritParams default_params_doc
#'
#' @return A data frame
#' @author Joshua W. Lambert
immig_event <- function(timeval,
                        island_tbl,
                        mainland_spec) {
  testit::assert(is.numeric(timeval))
  testit::assert(is.data.frame(island_tbl))
  testit::assert(ncol(island_tbl) == 7)
  testit::assert(is.numeric(mainland_spec))

  colonist <- DDD::sample2(mainland_spec, 1)
  if (length(island_tbl[, "spec_id"]) != 0) {
    isitthere <- which(island_tbl[, "spec_id"] == colonist)
  } else {
    isitthere <- c()
  }
  if (length(isitthere) == 0) {
    island_tbl <- rbind(
      island_tbl,
      data.frame(
        spec_id = colonist,
        main_anc_id = colonist,
        col_t = timeval,
        spec_type = "I",
        branch_code = as.character(NA),
        branch_t = NaN,
        ana_origin = as.character(NA)
      )
    )
  } else {
    island_tbl[isitthere, ] <- data.frame(
      spec_id = colonist,
      main_anc_id = colonist,
      col_t = timeval,
      spec_type = "I",
      branch_code = as.character(NA),
      branch_t = NaN,
      ana_origin = as.character(NA)
    )
  }
  return(island_tbl)
}