tests/testthat/test-aggregate.R

In EmilHvitfeldt/MethCon5: Identify and Rank CpG DNA Methylation Conservation Along the Human Genome

data <- data.frame(x = c(1, 1, 1, 2, 2, 3, 3, 3, 3, 4),
                   y = c(1, 2, 3, 4, 5, 6, 7, 8, NA, 0),
                   z = c(0, 9, 8, 7, 6, 5, 4, 3, 2, 1))

test_that("Base function work", {
  data_sum <- meth_aggregate(data, x, y)

  expect_equal(dim(data_sum), c(length(unique(data$x)), 3))
  expect_equal(data_sum$n, tabulate(data$x))
  expect_equal(data_sum$y, vapply(split(data$y, data$x), mean,
                                  FUN.VALUE = numeric(1), USE.NAMES = FALSE))

  expect_s3_class(data_sum, "methcon")
})

test_that("doens't work with multiple selections", {
  expect_error(meth_aggregate(data, x, y:z))
})

test_that("works with custom functions", {
  data_sum <- meth_aggregate(data, x, y, fun = min)

  expect_equal(dim(data_sum), c(length(unique(data$x)), 3))
  expect_equal(data_sum$y, vapply(split(data$y, data$x), min,
                                  FUN.VALUE = numeric(1), USE.NAMES = FALSE))

  data_sum <- meth_aggregate(data, x, y, fun = min, na.rm = TRUE)

  expect_equal(dim(data_sum), c(length(unique(data$x)), 3))
  expect_equal(data_sum$y, vapply(split(data$y, data$x), min,
                                  FUN.VALUE = numeric(1),
                                  na.rm = TRUE, USE.NAMES = FALSE))
})