def __init__(self): self.population = POP_DENVER self.dayspergen = 6.8 self.r0 = 2.65 # Infections caused while infected self.beta = calc_beta(self.r0, self.dayspergen) # Recoveries of infected per day self.gamma = 1.0 / self.dayspergen self.susceptible = self.population - 1 self.infected = 1 self.recovered = 0 self.total_days = 0 self.S_domain = [] self.E_domain = [] self.I_domain = [] self.R_domain = [] self.incubation = 0 self.exposed = 0 self.daystoisolation = 0 self.daystorecovery = 0 self.alpha = 0.2 self.time_domain = None
def run_r0_set(self, date_offsets, r0_values):
day_counter = 0
for itr in range(0, len(date_offsets)):
self.set_r0(r0_values[itr])
self.beta = calc_beta(self.r0, self.dayspergen)
print(f"{self.beta} = {calc_beta(self.r0, self.dayspergen)})")
while day_counter < date_offsets[itr]:
self.step_day()
day_counter += 1
def run_r0_set(self, date_offsets, r0_values):
self.reset()
prev_date = 0
for itr in range(0, len(date_offsets)):
self.set_r0(r0_values[itr])
self.beta = calc_beta(self.r0, self.dayspergen)
span = date_offsets[itr] - prev_date + 1
self.run_period(span)
prev_date = date_offsets[itr]
def run_r0_set(self, date_offsets, beta_values):
self.susceptible = self.population - self.infected - self.recovered
self.beta = calc_beta(self.r0, self.dayspergen)
prev_date = 0
for itr in range(0, len(date_offsets)):
self.set_r0(beta_values[itr])
span = date_offsets[itr] - prev_date + 1
self.run_period(span)
prev_date = date_offsets[itr]
def run_r0_set(self, date_offsets, r0_values):
self.scenario.hospital_door_aggregator = []
day_counter = 0
for itr in range(0, len(date_offsets)):
self.set_r0(r0_values[itr])
self.beta = calc_beta(r0_values[itr], self.infectious.period)
while day_counter < date_offsets[itr]:
self.step_day()
day_counter += 1
self.scenario.hospital_door_aggregator.append(
self.scenario.hospital_door_aggregator[-1])
def run_r0_set(self, date_offsets, r0_values):
self.reset()
prev_date = 0
for itr in range(0, len(date_offsets)):
self.set_r0(r0_values[itr])
self.beta = calc_beta(self.r0, self.dayspergen)
self.recalculate()
span = date_offsets[itr] - prev_date + 1
self.run_period2(span)
prev_date = date_offsets[itr]
if len(self.unhospitalized.domain) != len(self.dead.domain):
raise ValueError(
f"oops, {len(self.unhospitalized.domain)} != {len(self.dead.domain)}"
)
def __init__(self): self.population = POP_DENVER self.dayspergen = 6.8 self.r0 = 2.65 self.susceptible = self.population - 1 self.infected = 1 self.recovered = 0 # Infections caused while infected self.beta = calc_beta(self.r0, self.dayspergen) # Recoveries of infected per day self.gamma = 1.0 / self.dayspergen self.total_days = 0 self.S_domain = [] self.I_domain = [] self.R_domain = []
def recalculate(self): days_infectious = (self.dayspergen - self.incubation) * 2 self.alpha = 1.0 / self.incubation self.beta = calc_beta(self.r0, self.dayspergen) self.gamma = 1.0 / self.daystoisolation
def recalculate(self):
self.beta = calc_beta(self.r0, self.dayspergen)
def recalculate(self):
self.beta = calc_beta(self.r0, self.dayspergen)
self.gamma = 1.0 / self.dayspergen