SI_model <- odemodel(
name = "SI",
model = list(
S ~ - beta*S*I/N,
I ~ beta*S*I/N - gamma*I
),
observation = list(
susceptible ~ dnorm(mean=S, sd=sigma1),
infected ~ dnorm(mean=I, sd=sigma2)
),
initial = list(
S ~ N * (1 - i0),
I ~ N * i0
),
par = c("beta", "gamma", "N", "i0", "sigma1", "sigma2"),
link = c(i0="logit")
)
Run the code above in your browser using DataLab