Nothing
context("Internals")
spatstat.options(npixel = 2^3)
#creating a test that involves NA values in an image
xifull <- as.mask(heather$coarse, eps = 5 * heather$coarse$xstep)
w <- setminus.owin(Frame(xifull), square(5))
xiowin <- intersect.owin(xifull, w)
phat <- coverageprob(xiowin, obswin = w)
#creating logically valued image to pass
xiim.l <- as.im(xiowin, W = w, value = TRUE, na.replace = FALSE)
xiim.l[ setminus.owin(Frame(xiim.l), w) ] <- NA
#numerically valued image of 0, 1 and NA.
xiim.n <- xiim.l * 1
#numerically valued image of 0, 3 and NA.
xiim.n3 <- xiim.l * 3
test_that("isbinarymap operates correctly", {
expect_true(isbinarymap(xiim.n))
expect_true(isbinarymap(xiim.l))
expect_false(isbinarymap(xiim.n3))
expect_error(isbinarymap(xiim.n3, requiretrue = TRUE))
})
test_that("cppicka operates when passed im (logical and numeric) and fails correctly when im values not binary", {
cpest <- cppicka(xiowin, obswin = w)
expect_equal(max(abs(cpest - cppicka(xiim.l))), 0)
expect_equal(max(abs(cpest - cppicka(xiim.n))), 0)
expect_error(cppicka(xiim.n3), regexp = "has values other than")
})
test_that("plugincvc operates when passed im (logical and numeric) and fails correctly when im values not binary", {
cvcest <- plugincvc(xiowin, obswin = w)
expect_equal(max(abs(cvcest - plugincvc(xiim.l))), 0)
expect_equal(max(abs(cvcest - plugincvc(xiim.n))), 0)
expect_error(plugincvc(xiim.n3), regexp = "has values other than")
})
test_that("balancedracscovariance operates when passed im (logical and numeric) and fails correctly when im values not binary", {
cvcest.o <- racscovariance(xiowin, obswin = w)[["pickaH"]]
cvcest.l <- racscovariance(xiim.l)[["pickaH"]]
cvcest.n <- racscovariance(xiim.n)[["pickaH"]]
expect_equal(max(abs(cvcest.o - cvcest.l)), 0)
expect_equal(max(abs(cvcest.o - cvcest.n)), 0)
expect_error(racscovariance(xiim.n3), regexp = "has values other than")
})
test_that("paircorr operates when passed im (logical and numeric) and fails correctly when im values not binary", {
pclnest.o <- paircorr(xiowin, obswin = w)[["pickaint"]]
pclnest.l <- paircorr(xiim.l)[["pickaint"]]
pclnest.n <- paircorr(xiim.n)[["pickaint"]]
expect_equal(max(abs(pclnest.o - pclnest.l)), 0)
expect_equal(max(abs(pclnest.o - pclnest.n)), 0)
expect_error(paircorr(xiim.n3), regexp = "has values other than")
})
test_that("gblemp operates as expected when pass im object", {
boxwidths <- seq(1, 5, by = 1)
gblemp.l <- gblemp(boxwidths, xiim.l)
gblemp.n <- gblemp(boxwidths, xiim.n)
expect_equal(max(abs(gblemp.l - gblemp.n)), 0)
expect_error(gblemp(boxwidths, xiim.n3), regexp = "has values other than")
})
test_that("gbl() operates well when passed im or owin boxes", {
boxwidths <- seq(1, 5, by = 1)
expect_warning(gblest.l <- gbl(xiim.l, boxwidths,
c("GBLg.mattfeldt", "GBLg.pickaint", "GBLg.pickaH", "GBLcc.mattfeldt",
"GBLcc.pickaint", "GBLc", "GBLemp")),
regexp = "harmon")
expect_warning(gblest.n <- gbl(xiim.n, boxwidths,
c("GBLg.mattfeldt", "GBLg.pickaint", "GBLg.pickaH", "GBLcc.mattfeldt",
"GBLcc.pickaint", "GBLc", "GBLemp")),
regexp = "harmon")
expect_equal(gblest.l, gblest.n)
expect_error(gbl(xiim.n3, boxwidths), regexp = "isbinarymap")
})
reset.spatstat.options()
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.