out <- epiDCM(type='SIR', s.num=1000, i.num=1, r.num=0,
beta=0.2, cont=2:8, nu=1/3,
b=0.011, ms=0.01, mi=0.03, mr=0.01,
dt=500, verbose=TRUE)
statsDCM(out, time=25, run=1)
statsDCM(out, time=25, run=7)
statsDCM(out, time=26, run=7)
Run the code above in your browser using DataLab