someData <- data.frame(
  region = c(
    "transverse temporal",
    "insula",
    "precentral",
    "superior parietal",
    "transverse temporal",
    "insula",
    "precentral",
    "superior parietal"
  ),
  p = sample(seq(0, .5, .001), 8),
  Group = c(rep("G1", 4), rep("G2", 4)),
  stringsAsFactors = FALSE
) %>%
  group_by(Group)

test_that("Check that merging with grouped data works", {
  dk2 <- as_ggseg_atlas(dk)
  testData <- brain_join(someData, unnest(dk2, cols = ggseg))

  expect_equal(names(testData)[1], "Group")
  expect_equal(unique(testData$Group), c("G1", "G2"))
})

test_that("Check that plotting with grouped data works", {
  pp <- ggseg(.data = someData, mapping = aes(fill = p)) +
    facet_wrap(~Group)

  expect_true("Group" %in% names(pp$data))
})

test_that("Check that simple brain_join works", {
  someData <- someData %>%
    #tidyr::unnest() %>%
    dplyr::filter(Group == "G1")
  dk2 <- as_ggseg_atlas(dk)

  testData <- brain_join(someData, unnest(dk2, ggseg))

  expect_equal(names(testData)[1], "Group")
  expect_equal(unique(testData$Group), "G1")
})

test_that("Check that simple-features brain_join works", {
  someData <- someData %>%
    dplyr::filter(Group == "G1")

  testData <- brain_join(someData, dk)
  expect_true(inherits(testData, "sf"))
  expect_equal(names(testData)[1], "Group")
  expect_equal(unique(testData$Group), "G1")
})