tests/egf_examples_day_of_week.R

In davidearn/epigrowthfit: Nonlinear Mixed Effects Models of Epidemic Growth

library(epigrowthfit)
options(warn = 2L, error = if (interactive()) recover)


## day_of_week #########################################################

r <- log(2) / 20
tinfl <- 100
K <- 25000
w <- c(1, 1, 1, 0.8, 0.8, 1.4)
disp <- 50

zz <- simulate(egf_model(curve = "logistic", family = "nbinom",
                         day_of_week = 3L),
               nsim = 1L,
               seed = 366465L,
               mu = log(c(r, tinfl, K, disp, w)),
               cstart = 10)
mm <- egf(zz)

stopifnot(all.equal(coef(zz), coef(mm), tolerance = 5e-02))