demo/bern-SIRC.R

In kypraios/epiABC: ABC methods for infectious disease epidemics

library(epiABC)
sin=proc.time()

N=10000
OUT=matrix(0,ncol=89,nrow=N)
PARA=matrix(0,ncol=5,nrow=N)

for(i in 1:N)
{
p=runif(1)
beta=rexp(1,2)
gamma=rgamma(1,2,1)
#print(PARA[i,])
Xout=bernSIRC(89,beta,gamma,p,40)
OUT[i,]=Xout$output
PARA[i,]=c(p,beta,gamma,Xout$count,OUT[i,Xout$count])
}


sout=proc.time()
sout-sin