tests/testthat/test-depthgram.R

In astamm/roahd: Robust Analysis of High Dimensional Data

test_that("`depthgram()` works as expected for all supported input types", {
  # Arrange
  withr::local_seed(1234)
  N <- 100
  P <- 100
  grid <- seq(0, 1, length.out = P)
  Cov <- exp_cov_function(grid, alpha = 0.3, beta = 0.4)

  Data <- list()
  Data[[1]] <- generate_gauss_fdata(
    N,
    centerline = sin(2 * pi * grid),
    Cov = Cov
  )
  Data[[2]] <- generate_gauss_fdata(
    N,
    centerline = sin(2 * pi * grid),
    Cov = Cov
  )
  names <- paste0("id_", 1:nrow(Data[[1]]))
  fD <- fData(grid, Data[[1]])
  mfD <- mfData(grid, Data)

  # Act
  actual_list <- depthgram(Data, marginal_outliers = TRUE, ids = names)
  actual_fData <- depthgram(fD, marginal_outliers = TRUE, ids = names)
  actual_mfData <- depthgram(mfD, marginal_outliers = TRUE, ids = names)

  # Assert
  expect_snapshot_value(actual_list, style = "json2")
  expect_snapshot_value(actual_fData, style = "json2")
  expect_snapshot_value(actual_mfData, style = "json2")
})