tests/testthat/test-bin_lat.R

In palaeoverse: Prepare and Explore Data for Palaeobiological Analyses

test_that("bin_lat works", {
  # Load occurrence data
  occdf <- tetrapods
  # Generate latitudinal bins
  bins <- lat_bins_degrees(size = 10)
  # Bin data
  expect_equal(nrow(bin_lat(occdf = occdf, bins = bins, lat = "lat")),
               nrow(occdf))

  # Boundary occurrences
  bo <- length(which(occdf[, "lat", drop = TRUE] %in% c(bins$max, bins$min)))
  expect_equal(nrow(bin_lat(occdf = occdf, bins = bins, lat = "lat",
                            boundary = TRUE)),
               nrow(occdf) + bo)

  # Expect error
  expect_error(bin_lat(occdf = 2, bins = bins, lat = "lat"))
  expect_error(bin_lat(occdf = occdf, bins = 2, lat = "lat"))
  expect_error(bin_lat(occdf = occdf, bins = 2, lat = "plat"))
  bins <- bins[, -c(1)]
  expect_error(bin_lat(occdf = occdf, bins = bins, lat = "lat"))
  occdf$lat[1] <- NA
  expect_error(bin_lat(occdf = occdf, bins = bins, lat = "lat"))
  occdf$lat[1] <- 91
  expect_error(bin_lat(occdf = occdf, bins = bins, lat = "lat"))
  expect_error(bin_lat(occdf = occdf, bins = bins, lat = "latitude"))

})