tests/testthat/test-pa.R

In ElliotDovers/scampr: Spatially Correlated, Approximate Modelling of Point patterns in R

context("Presence/Absence Model Function: pa()")
# usethis::use_test() when particular R function file is open to create these

test_that("ipp pa model creates scampr model", {
  dat_pa <- flora$pa
  mod <- pa(sp1 ~ MNT, dat_pa, model.type = "ipp")
  testthat::expect_equal(class(mod), "scampr")
})

test_that("lgcp (va) pa model doesn't work", {
  dat_pa <- flora$pa
  bfs <- simple_basis(nodes.on.long.edge = 5, data = dat_pa)
  testthat::expect_error(pa(sp1 ~ MNT, dat_pa, model.type = "variational", simple.basis = bfs))
})

test_that("lgcp (laplace) pa model creates scampr model", {
  dat_pa <- flora$pa
  bfs <- simple_basis(nodes.on.long.edge = 5, data = dat_pa)
  mod <- pa(sp1 ~ MNT, dat_pa, model.type = "laplace", simple.basis = bfs)
  testthat::expect_equal(class(mod), "scampr")
})

test_that("ipp pa model creates scampr model w. subset", {
  dat_pa <- flora$pa
  mod <- pa(sp1 ~ MNT, dat_pa, model.type = "ipp", subset = 1:nrow(dat_pa) < (0.5 * nrow(dat_pa)))
  testthat::expect_equal(class(mod), "scampr")
})

test_that("ipp pa model creates scampr model works with interaction terms in formula", {
  dat_pa <- flora$pa
  mod <- pa(sp1 ~ MNT*MXT, dat_pa, model.type = "ipp")
  testthat::expect_equal(class(mod), "scampr")
})

test_that("lgcp (laplace) pa model works with interaction terms in formula", {
  dat_pa <- flora$pa
  bfs <- simple_basis(nodes.on.long.edge = 5, data = dat_pa)
  mod <- pa(sp1 ~ MNT*MXT, dat_pa, model.type = "laplace", simple.basis = bfs)
  testthat::expect_equal(class(mod), "scampr")
})

test_that("lgcp (laplace) pa model works w. subset", {
  dat_pa <- flora$pa
  bfs <- simple_basis(nodes.on.long.edge = 5, data = dat_pa)
  mod <- pa(sp1 ~ MNT, dat_pa, model.type = "laplace", simple.basis = bfs, subset = 1:nrow(dat_pa) < (0.5 * nrow(dat_pa)))
  testthat::expect_equal(class(mod), "scampr")
})