set.seed(234)
mydata <- data.frame(infected = rchainbinom(n = 15, s0 = 5, sar = 0.2,
i0 = 1, generations = Inf),
s0 = 5, i0 = 1, generations = Inf)
xmat <- model.matrix(~ 1, data = mydata)
res <- cbmod(y = mydata$infected, s0 = mydata$s0, x = xmat, i0 = mydata$i0,
generations = mydata$generations, link = 'identity')
summary(res)
predict(res, x = xmat, type = 'sar')
Run the code above in your browser using DataLab