tests/testthat/test_value.R

In fedeago/NewWave: Negative binomial model for scRNA-seq

context("Test numerical correctness of functions.")
set.seed(1234)


test_that("Estimates are reasonable when data is Poisson", {
  
    counts <- matrix(rpois(10000, lambda=50), nrow=100, ncol=100)
    m1 <- newFit(counts)
    expect_true(all(newPhi(m1) < 1e-4))
    
    m2 <- newFit(counts, commondispersion = FALSE)
    expect_true(all(newPhi(m2) < 1e-4))
    
  
    expect_true(abs(mean(newMu(m1)) - 50) < 1)
    expect_true(abs(mean(newMu(m2)) - 50) < 1)
})

test_that("Estimates are reasonable when data is Negative Binomial", {
  
    counts <- matrix(rnbinom(10000, mu=50, size = 10), nrow=100, ncol=100)
  
    m1 <- newFit(counts, commondispersion = TRUE)
  
    expect_true(abs(mean(newMu(m1)) - 50) < 1)
    expect_true(abs(mean(newTheta(m1)) - 10) < 2)
})