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

  # Boundary occurrences
  bo <- length(which(occdf[, "lat"] %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"))

})