tests/testthat/test_spacing.r

In valr: Genome Interval Arithmetic

test_that("start intervals are NA", {
  x <- tibble::tribble(
    ~chrom, ~start, ~end,
    "chr1", 1, 100,
    "chr1", 150, 200,
    "chr2", 200, 300
  )

  res <- interval_spacing(x)

  first <- res |>
    group_by(chrom) |>
    slice(1) |>
    select(chrom:end)

  nas <- filter(res, is.na(.spacing)) |>
    select(chrom:end)

  expect_true(all(first == nas))
})

# from bedtools2
test_that("bt test succeeds", {
  x <- tibble::tribble(
    ~chrom, ~start, ~end,
    "chr1", 20, 30,
    "chr1", 25, 40,
    "chr1", 40, 50,
    "chr1", 60, 80,
    "chr1", 75, 100,
    "chr1", 105, 110,
    "chr2", 115, 130,
    "chr2", 120, 160,
    "chr2", 170, 180
  )

  res <- interval_spacing(x)

  exp <- c(NA, -5, 0, 10, -5, 5, NA, -10, 10)
  expect_equal(res$.spacing, exp)
})