tests/testthat/test-SEIR_ode.R

context("SEIR_ode")



##Model Input
S_0 <- 989
E_0 <- 10
I_0 <- 1
R_0 <- 0
beta <- 3
gamma <-  1/2
tau <- 2
dt <- 1

parameters <- c(beta = beta, gamma = gamma, tau = tau)
inits <- c(S = S_0, E = E_0, I = I_0, R_0 = R_0)

test <- SEIR_ode(1, inits, parameters)


test_that("SEIR_ode can correctly return rates of change over one timestep", {
  skip_on_cran()
  expect_known_output(test, file = "../../tests/test-files/SEIR_ode/test-01.rds")
})

Try the idmodelr package in your browser

Any scripts or data that you put into this service are public.

idmodelr documentation built on Sept. 2, 2022, 5:06 p.m.