tests/testthat/test_carp_sparsity_pattern.R

In jjn13/clustRviz: Convex Clustering and Bi-Clustering via Algorithmic Regularization

context("Test CARP Sparsity Invariants")

test_that("CARP begins with no fusions", {
  clustRviz_options(keep_debug_info = TRUE)

  carp_fit <- CARP(presidential_speech, back_track = FALSE)
  expect_zeros(carp_fit$debug$row$v_zero_indices[, 1])

  clustRviz_reset_options()
})

test_that("CARP ends with all fusions", {
  clustRviz_options(keep_debug_info = TRUE)

  carp_fit <- CARP(presidential_speech, back_track = FALSE)
  expect_ones(carp_fit$debug$row$v_zero_indices[, NCOL(carp_fit$debug$row$v_zero_indices)])

  clustRviz_reset_options()
})

test_that("CARP-VIZ begins with no fusions", {
  clustRviz_options(keep_debug_info = TRUE)

  carp_fit <- CARP(presidential_speech, back_track = TRUE)
  expect_zeros(carp_fit$debug$row$v_zero_indices[, 1])

  clustRviz_reset_options()
})

test_that("CARP-VIZ ends with all fusions", {
  clustRviz_options(keep_debug_info = TRUE)

  carp_fit <- CARP(presidential_speech, back_track = TRUE)
  expect_ones(carp_fit$debug$row$v_zero_indices[, NCOL(carp_fit$debug$row$v_zero_indices)])

  clustRviz_reset_options()
})

test_that("CARP begins with no fusions (uniform weights)", {
  clustRviz_options(keep_debug_info = TRUE)

  weight_mat <- matrix(1, nrow=NROW(presidential_speech), ncol=NROW(presidential_speech))
  carp_fit   <- CARP(presidential_speech, weights = weight_mat, back_track = FALSE)
  expect_zeros(carp_fit$debug$row$v_zero_indices[, 1])

  clustRviz_reset_options()
})

test_that("CARP ends with all fusions (uniform weights)", {
  clustRviz_options(keep_debug_info = TRUE)

  weight_mat <- matrix(1, nrow=NROW(presidential_speech), ncol=NROW(presidential_speech))
  carp_fit   <- CARP(presidential_speech, weights = weight_mat, back_track = FALSE)
  expect_ones(carp_fit$debug$row$v_zero_indices[, NCOL(carp_fit$debug$row$v_zero_indices)])

  clustRviz_reset_options()
})

test_that("CARP-VIZ begins with no fusions (uniform weights)", {
  clustRviz_options(keep_debug_info = TRUE)

  weight_mat <- matrix(1, nrow=NROW(presidential_speech), ncol=NROW(presidential_speech))
  carp_fit   <- CARP(presidential_speech, weights = weight_mat, back_track = TRUE)
  expect_zeros(carp_fit$debug$row$v_zero_indices[, 1])

  clustRviz_reset_options()
})

test_that("CARP-VIZ ends with all fusions (uniform weights)", {
  clustRviz_options(keep_debug_info = TRUE)

  weight_mat <- matrix(1, nrow=NROW(presidential_speech), ncol=NROW(presidential_speech))
  carp_fit   <- CARP(presidential_speech, weights = weight_mat, back_track = TRUE)
  expect_ones(carp_fit$debug$row$v_zero_indices[, NCOL(carp_fit$debug$row$v_zero_indices)])

  clustRviz_reset_options()
})