tests/testthat/tests_cv2.R

In reumandc/tsvr: Timescale-Specific Variance Ratio for Use in Community Ecology

context("cv2")

test_that("test error catching",{
  X<-matrix(runif(10*20)+1,20,20)
  expect_error(cv2(X,"test"),"Error in cv2: type must be com, comip, or pop")
})

test_that("test cases where it actually provides results",{
  set.seed(401)
  X<-matrix(runif(10*20)+1,10,20)
  h<-cv2(X, "com")
  Xtot<-colSums(X)
  expect_equal(h,(sd(Xtot)/mean(Xtot))^2)
  h<-cv2(X, "comip")
  vars<-apply(FUN=var,X=X,MARGIN=1)
  expect_equal(sum(vars)/((mean(Xtot))^2),h)
  h<-cv2(X, "pop")
  expect_equal(h,(sum(sqrt(vars)))^2/(mean(Xtot)^2))
})