tests/testthat/test_momentuHMM.R

context("momentuHMM")

test_that("Exceptions are thrown",{
  data <- data.frame(1,2,3)
  states <- rep(1,10)
  mle <- list(1,2,3)
  stepDist <- "gamma"
  angleDist <- "vm"
  mod <- list(1,2,3)
  conditions <- list(dist=list(step=stepDist,angle=angleDist),zeroInflation=list(step=FALSE),estAngleMean=list(angle=TRUE),stationary=FALSE,formula=~c1)
  rawCovs <- data.frame(c1=rnorm(10))
  stateNames <- c("s1","s2")

  m <- list(data=data,mle=mle,mod=mod,
            conditions=conditions,rawCovs=rawCovs,stateNames=stateNames)
  expect_error(momentuHMM(m),NA)

  m <- list(data=data,states=states,mle=mle,stepDist=stepDist)
  expect_error(momentuHMM(m))
})

test_that("The output has the right class attribute",{
  data <- data.frame(1,2,3)
  states <- rep(1,10)
  mle <- list(1,2,3)
  stepDist <- "gamma"
  angleDist <- "vm"
  mod <- list(1,2,3)
  conditions <- list(dist=list(step=stepDist,angle=angleDist),zeroInflation=list(step=FALSE),estAngleMean=list(angle=TRUE),stationary=FALSE,formula=~c1)
  rawCovs <- data.frame(c1=rnorm(10))
  stateNames <- c("s1","s2")

  m <- list(data=data,mle=mle,mod=mod,
            conditions=conditions,rawCovs=rawCovs,stateNames=stateNames)
  m <- momentuHMM(m)

  expect_equal(length(which(class(m)=="momentuHMM")),1)
})
bmcclintock/momentuHMM documentation built on Oct. 26, 2022, 1 a.m.