tests/testthat/test-wrap_cluster_graph.R

In dynverse/dynwrap: Representing and Inferring Single-Cell Trajectories

context("Testing add_cluster_graph")

# cell data
id <- "a"
cell_ids <- letters
group_ids <- LETTERS[1:5]
grouping <- sample(group_ids, length(cell_ids), replace = T) %>% set_names(cell_ids)
extras <- "banana"

wr_without_grouping <- wrap_data(
  id = id,
  cell_ids = cell_ids
)

wr_orig <- wr_without_grouping %>% add_grouping(
  group_ids = group_ids,
  grouping = grouping,
  extras = extras
)

milestone_network <- tibble(
  from = group_ids[1:4],
  to = group_ids[2:5],
  length = c(1, 2, 4, 5),
  directed = TRUE
)

test_that("Testing add_cluster_graph", {
  wr <- wr_orig %>%
    add_cluster_graph(
      milestone_network = milestone_network
    )

  # testing is_ti_data_wrapper
  expect_false(is_wrapper_with_trajectory(wr_orig))
  expect_true(is_wrapper_with_trajectory(wr))

  expect_equivalent(wr$id, id)
  expect_equivalent(wr$group_ids, group_ids)
  expect_setequal(wr$milestone_ids, group_ids)
  expect_equivalent(wr$cell_ids, cell_ids)
  expect_equivalent(wr$extras, extras)
  expect_equivalent(wr$grouping, grouping)
  expect_equivalent(wr$milestone_network, milestone_network)

  # percentages are either 0 or 1
  expect_true(all(abs(abs(wr$milestone_percentages$percentage - .5) - .5) < 1e-6))

  # test with providing a grouping in cluster_graph
  wr <- wr_without_grouping %>%
    add_cluster_graph(
      milestone_network = milestone_network,
      grouping = grouping
    )

  expect_true(is_wrapper_with_trajectory(wr))
})



test_that("Testing cluster_graph_add_explicit_splits", {
  milestone_network <- tribble(
    ~from, ~to, ~length, ~directed,
    "A", "B", 1, TRUE,
    "A", "C", 1, TRUE,
    "C", "D", 2, TRUE,
    "E", "F", 3, TRUE,
    "B", "G", 1, TRUE,
    "C", "G", 10, TRUE
  )

  milestone_network <- cluster_graph_add_explicit_splits(milestone_network)

  expect_true(
    all(
      c(
        "split_A -> B",
        "split_A -> C",
        "B -> convergence_G",
        "C -> convergence_G",
        "A -> split_A",
        "convergence_G -> G"
      ) %in%
        paste0(milestone_network$from, " -> ", milestone_network$to)
    )
  )
})