tests/testthat/test-gREM.R

In timcdlucas/encounteR: Estimate animal density and abundance using random encounter models

context("Test grem functions")


# calcProfileWidth() tests

test_that('Profile width is within 0 and 2*r', {
        alpha <- seq(0, 2*pi, length.out=100)
        theta <- seq(0, 2*pi, length.out=100)
        paras <- expand.grid(alpha, theta)
        paras <- cbind(paras, 1)
        colnames(paras) <- c('theta','alpha','r')
        p <- sapply(1:nrow(paras), function(x) do.call(calcProfileWidth,as.list(paras[x,])))
        expect_that(sum(p>2*pi) == 0, is_true())
        expect_that(sum(p<0) == 0, is_true())
})


test_that('Profile width is 0 when alpha is 0', {
        theta <- seq(0, 2*pi, length.out=100)
        paras <- cbind(theta, 0, 1)
        colnames(paras) <- c('theta','alpha','r')
        p <- sapply(1:nrow(paras), function(x) do.call(calcProfileWidth,as.list(paras[x,])))
        expect_that(all(p==0), is_true())
})


test_that('Profile width is 0 when r is 0', {
        alpha <- seq(0, 2*pi, length.out=100)
        theta <- seq(0, 2*pi, length.out=100)
        paras <- expand.grid(alpha, theta)
        paras <- cbind(paras, 0)
        colnames(paras) <- c('theta','alpha','r')
        p <- sapply(1:nrow(paras), function(x) do.call(calcProfileWidth,as.list(paras[x,])))
        expect_that(all(p==0), is_true())
})


# gremDensity() tests

test_that('Density returns error if p is 0', {
#        expect_that(gremDensity(2, 0.1, 0.1, 0, 2, 2), throws_error())
#        expect_that(gremDensity(2, 0.1, 0.1, 10, 0, 2), throws_error())
#        expect_that(gremDensity(2, 0.1, 0.1, 10, 2, 0), throws_error())
#        expect_that(gremDensity(2, 0, 0.1, 10, 2, 2), throws_error())
})

test_that('z of zero gives error', {
        expect_that(gremDensity(0, 0.1, 0.1, 2, 2, 2), throws_error())
})


# gremAbundance() tests