tests/testthat/test-Mbound.R

In nickbrazeau/discent: Estimation of Deme Inbreeding Spatial Coefficients with Gradient Descent

test_that("M is properly bounded", {
  # sim data
  data("IBD_simulation_data", package = "discent")
  dat <- IBD_simulation_data
  # start params
  our_start_params <- rep(0.2, 3)
  names(our_start_params) <- 1:3
  our_start_params <- c(our_start_params, "m" = 1e3)
  # run model w/ EXTREMELY TIGHT bounds
  inputdisc <- dat %>%
    dplyr::filter(deme1 != deme2)
  mod1 <- disc(discdat = inputdisc,
               start_params = our_start_params,
               learningrate = 1e-3,
               m_lowerbound = 999,
               m_upperbound = 1001,
               b1 = 0.9,
               b2 = 0.999,
               e = 1e-8,
               steps = 1e2,
               report_progress = TRUE,
               normalize_geodist = FALSE,
               return_verbose = TRUE)
  testthat::expect_gte(min(mod1$m_run), 999)
  testthat::expect_lte(max(mod1$m_run), 1001)

})