Я пытаюсь смоделировать базовую c модель SIR, а затем рассчитать (и построить график) коэффициент рождаемости на душу населения в хостах ... это правильно или мне нужно что-то более сложное? Что, если, например, я хочу изменить параметр (например, бета-версию), а затем рассчитать / изучить, как изменения рождаются при увеличении / уменьшении беты? В таком случае мне понадобится a для l oop (или что-то подобное), или приведенный ниже код даст правильный ответ?
Для ясности, этот код работает, но построение дБ показывает экспоненциальный рост рождаемости во времени ... это не кажется правильным.
Большое спасибо заранее.
dS = b*S + b*I - d*S - beta*S*p # susceptible hosts
dI = beta*S*p - (d + v)*I # infected hosts
db <- (b*S + b*I)/(S + I) # attempt to capture birth rate, weighted by total number of hosts
res <- c(dS, dI, db)