# Deterministic SIR model with varying act.rate
mod <- epiDCM(type="SIR", s.num=1000, i.num=1, r.num=0,
trans.rate=0.2, act.rate=2:8, rec.rate=1/3,
b.rate=0.011, ds.rate=0.01, di.rate=0.03,
dr.rate=0.01, nsteps=500, verbose=TRUE)
summary(mod, time=25, run=1, comp.plot=TRUE)
summary(mod, time=25, run=7, comp.plot=TRUE)
summary(mod, time=26, run=7, comp.plot=TRUE)
Run the code above in your browser using DataLab