tests/testthat/test-removeNA.R

In CenterForStatistics-UGent/RCM: Fit row-column association models with the negative binomial distribution for the microbiome

context("Correct for missingness")

n = 50; p = 100
X = matrix(rnbinom(n*p, size = 1, mu = 2), n,p)
X[sample(seq_len(n*p), n)] = NA
mu = outer(rowSums(X, na.rm = TRUE), colSums(X, na.rm = TRUE))/
  sum(X, na.rm = TRUE)

test_that("correctXMissingness() removes all NAs", {
  expect_true(anyNA(RCM:::correctXMissingness(X, mu, FALSE)))
  expect_false(anyNA(RCM:::correctXMissingness(X, mu, TRUE)))
})