tests/testthat/test-sdraw.SpatialPolygons.R

In tmcd82070/SDraw: Spatially Balanced Samples of Spatial Objects

# test-sdraw.SpatialPolygons.R
context("Testing sdraw.SpatialPolygons()")


# create spatial polygons object
# make "squares" with integer values rounded in a list
square <- rbind( c(2, 4, 3, 4, 3, 5,
                   2, 5, 2, 4, 2, 4),
                 c(6, 9, 7, 9, 7, 8,
                   6, 8, 6, 9, 6, 9))

# give these "squares" an identification
ID <- c("shape1", "shape2")

# create SpatialPolygon object from these squares 
spatPoly <- SpatialPolygons(list(
  Polygons(list(Polygon(matrix(square[1, ], ncol = 2, byrow = TRUE))), ID[1]),
  Polygons(list(Polygon(matrix(square[2, ], ncol = 2, byrow = TRUE))), ID[2])
))


# check if the function stops with message
test_that("invalid SpatialPolygons sample type", {
  capture.output(expect_error(sdraw.SpatialPolygons(spatPoly, 15, type="HOP"),"Invalid SpatialPolygons sample type = HOP",fixed=TRUE))
  capture.output(expect_error(sdraw.SpatialPolygons(spatPoly, 20, type="BTS"),"Invalid SpatialPolygons sample type = BTS",fixed=TRUE))
  capture.output(expect_error(sdraw.SpatialPolygons(spatPoly, 25, type="SAS"),"Invalid SpatialPolygons sample type = SAS",fixed=TRUE))
  capture.output(expect_error(sdraw.SpatialPolygons(spatPoly, 30, type="STS"),"Invalid SpatialPolygons sample type = STS",fixed=TRUE))
  capture.output(expect_error(sdraw.SpatialPolygons(spatPoly, 35, type="RTS"),"Invalid SpatialPolygons sample type = RTS",fixed=TRUE))
})

# check the output and length
test_that("length sdraw.SpatialPolygons(spatPoly,10, type = GRTS is 10", {
  capture.output(expect_length(sdraw.SpatialPolygons(spatPoly,10, type ="GRTS"),length(primes(10))))
})