tests/testthat/test-lcScore.R

In HajkD/myTAI: Evolutionary Transcriptomics

context("Test: lcScore() ")


data(PhyloExpressionSetExample)
data(DivergenceExpressionSetExample)

tai.vals <- TAI(PhyloExpressionSetExample)
tdi.vals <- TDI(DivergenceExpressionSetExample)

test_that("lcScore computes correct early conservation scores...",{
  expect_equal(lcScore(tai.vals,
                       early = 1:2, mid = 3:5, late = 6:7),min(c(mean(tai.vals[1:2] - mean(tai.vals[6:7])), mean(tai.vals[3:5] - mean(tai.vals[6:7])))))
  
  expect_equal(lcScore(tdi.vals,
                       early = 1:2, mid = 3:5, late = 6:7),min(c(mean(tdi.vals[1:2] - mean(tdi.vals[6:7])), mean(tdi.vals[3:5] - mean(tdi.vals[6:7])))))
})