tests/testthat/testthat_post.R

In POSTm: Phylogeny-Guided OTU-Specific Association Test for Microbiome Data

test_that("`post()` properly handles intercept only", {
  data("POSTmData")
  
  y <- as.integer(x = metadata[,"GC"] == "BV")
  X <- metadata[,"mRace"]
  
  expect_error(post(y = y, 
               X = matrix(1, length(y), 1L), 
               OTU = otu[,1:20], 
               tree = otutree,
               cValues = seq(0,0.05,by=0.01)),
               "an intercept only model should be defined as `X = NULL`")
})