tests/testthat/test_glasso_NA_methods.R

In MalteLond/rfcd: Random Forest Change Point Detection

context("NA methods")

test_that("complete", {
  
  x_del <- readRDS('../data/reference_glasso_x_del_blockwise_002')
  tree_complete <- hdcd(x_del, method = 'glasso', optimizer = 'section_search', delta = 0.2, lambda = 0.1, control = hdcd_control(glasso_NA_method = 'complete_observations'))
  expect_equal(get_change_points_from_tree(tree_complete, 'max_gain'), c(70, 130))
  
})

test_that("average", {
  
  x_del <- readRDS('../data/reference_glasso_x_del_blockwise_03')
  tree_average <- hdcd(x_del, method = 'glasso', optimizer = 'section_search', delta = 0.1, lambda = 0.1, control = hdcd_control(glasso_NA_method = 'average'))
  expect_equal(get_change_points_from_tree(tree_average, 'max_gain'), c(56, 79, 105, 139, 173))
  
})
  
test_that("LW", {
  
  x_del <- readRDS('../data/reference_glasso_x_del_blockwise_03')
  tree_LW <- hdcd(x_del, method = 'glasso', optimizer = 'section_search', delta = 0.1, lambda = 0.01, control = hdcd_control(glasso_NA_method = 'loh_wainwright_bias_correction'))
  expect_equal(get_change_points_from_tree(tree_LW, 'max_gain'), c(35, 70, 95, 120, 157, 178))
  
})

test_that("pairwise", {
  
  x_del <- readRDS('../data/reference_glasso_x_del_blockwise_03')
  tree_pairwise <- hdcd(x_del, method = 'glasso', optimizer = 'section_search', delta = 0.1, lambda = 0.1, control = hdcd_control(glasso_NA_method = 'pairwise'))
  expect_equal(get_change_points_from_tree(tree_pairwise, 'max_gain'),  c(71, 95, 120, 157, 178))
  
})