# SEIR Connected model
seirconn <- ModelSEIRCONN(
name = "Disease",
n = 10000,
prevalence = 0.1,
contact_rate = 2.0,
transmission_rate = 0.8,
incubation_days = 7.0,
recovery_rate = 0.3
)
set.seed(937)
run(seirconn, 50)
# Retrieving the transition probability
get_transition_probability(seirconn)
# Get active cases
head(get_active_cases(seirconn))
# Get outbreak size
head(get_outbreak_size(seirconn))
Run the code above in your browser using DataLab