tests/testthat/test-jstransientmodel.R

In r-glennie/openpopscr: Open Population Spatial Capture-Recapture

context("JsTransientModel test")

# set truth 
true_par <- list(D = 1000, lambda0 = 2, sigma = 20, phi = 0.5, beta = 0.25, sd = 10)
# make detectors array 
detectors <- make.grid(nx = 7, ny = 7, spacing = 20, detector = "count")
# make mesh 
mesh <- make.mask(detectors, buffer = 100, nx = 64, ny = 64, type = "trapbuffer")
# set number of occasions to simulate
n_occasions <- 5 
# simulate ScrData 
scrdat <- simulate_js_openscr(true_par, n_occasions, detectors, mesh, move = TRUE, print = FALSE, seed = 225622)
# formula 
par <- list(lambda0 ~ 1, 
            sigma ~ 1, 
            beta ~ 1, 
            phi ~ 1, 
            sd ~ 1, 
            D ~ 1)
# get automatic start values 
start <- get_start_values(scrdat, model = "JsTransientModel")

test_that("Basic model computes correct llk", {
  obj <- JsTransientModel$new(par, scrdat, start, print = FALSE)
  expect_equal(signif(obj$calc_llk(), 6), -806.707)
}  
)