test_that("validation", {
  
  env <- new.env()
  
  
  
  
  # assert_integer_counts() ====
  
  env$counts <- counts_p
  expect_error(assert_integer_counts(env))
  
  
  
  # validate_pairs() ====
  
  env$counts <- matrix(nrow = 5, ncol = 0)
  expect_error(validate_pairs(env))
  
  env$counts <- counts
  
  env$pairs <- function (i,j) list()
  expect_error(validate_pairs(env))
  
  env$pairs <- TRUE; expect_error(validate_pairs(env))
  env$pairs <- 1.5;  expect_error(validate_pairs(env))
  env$pairs <- -1;   expect_error(validate_pairs(env))
  env$pairs <- 100;  expect_error(validate_pairs(env))
  
  env$pairs <- 2:5
  expect_silent(validate_pairs(env))
  
  env$pairs <- c(F,F,T,T,T,F,F,F,F,F)
  expect_silent(validate_pairs(env))
  
  
  
  
  # validate_alpha() ====
  
  env$alpha <- 1L
  expect_silent(validate_alpha(env))
  
  
  
  
  # validate_cpus() ====
  
  env$cpus <- 1
  expect_silent(validate_cpus(env))
  
  
  
  
  # validate_power() ====
  
  env$power <- 1L
  expect_silent(validate_power(env))
  
  
  
  
  # validate_pseudocount() ====
  
  env$pseudocount <- NULL
  expect_silent(validate_pseudocount(env))
  
  env$pseudocount <- 1L
  expect_silent(validate_pseudocount(env))
  
  
  
  
  # validate_rescale() ====
  
  env$rescale <- NULL
  expect_error(validate_rescale(env))
  
  env$rescale <- 'a'
  expect_error(validate_rescale(env))
  
  
  
  
  # validate_tree() ====
  
  env$counts <- counts[1:3,,drop=FALSE]
  tree2      <- tree
  tree2$edge.length <- as.integer(tree$edge.length * 100)
  tree2$edge <- matrix(
    data     = as.numeric(tree$edge), 
    nrow     = nrow(tree$edge), 
    ncol     = ncol(tree$edge),
    dimnames = dimnames(tree$edge) )
  env$tree <- tree2
  expect_silent(validate_tree(env))
  
  env$counts <- counts[,1:4,drop=FALSE]
  env$tree   <- tree
  expect_silent(validate_tree(env))
  
  
  
  
  # validate_counts() ====
  
  env$tree   <- NULL
  env$counts <- counts
  attr(env$counts, 'tree') <- tree
  expect_silent(validate_counts(env))
  
  skip_on_cran()
  skip_if_not_installed('rbiom')
  
  hmp50               <- do.call(`::`, list('rbiom', 'hmp50'))
  convert_to_phyloseq <- do.call(`::`, list('rbiom', 'convert_to_phyloseq'))
  convert_to_TSE      <- do.call(`::`, list('rbiom', 'convert_to_TSE'))
  convert_to_SE       <- do.call(`::`, list('rbiom', 'convert_to_SE'))
  
  env$tree   <- NULL
  env$counts <- 1:10
  expect_silent(validate_counts(env))
  
  env$tree   <- NULL
  env$counts <- t(hmp50$counts)
  expect_silent(validate_counts(env))
  
  env$tree   <- NULL
  env$counts <- hmp50
  expect_silent(validate_counts(env))
  
  env$tree   <- NULL
  env$counts <- convert_to_phyloseq(hmp50)
  expect_silent(validate_counts(env))
  
  env$tree   <- NULL
  env$counts <- convert_to_TSE(hmp50)
  expect_silent(validate_counts(env))
  
  env$tree   <- NULL
  env$counts <- convert_to_SE(hmp50)
  expect_silent(validate_counts(env))
  
  
})