alpha=0.3)
ax1.text(
dateCache.parse("2020-04-15"),
45500,
"Geen smoesjes,\nje weet het best.\nAls je niet ziek wordt,\nhoef je ook niet getest.",
color="gray",
bbox=dict(facecolor='white', alpha=1.0, edgecolor='white'),
zorder=10)
ax1.plot(totaaltests['x'],
totaaltests['y'],
color='lightblue',
linestyle='-',
label='Uitgevoerde tests (nu: ' +
decimalstring(totaaltests['y'][-1]) + ')')
# Plot cases per dag
totaal_percentage = round(100 * totaal_positief / totaaltests['total'], 1)
ax1.plot(positief['x'][:-positief_voorspeld['avgsize']],
positief['y'][:-positief_voorspeld['avgsize']],
color='fuchsia',
label="Tests positief (nu: %s)" % decimalstring(positief['y'][-1]))
#ax1.plot(positief['x'][-11:], positief['y'][-11:], color='steelblue', linestyle='--', alpha=0.3, label='onvolledig')
ax1.plot(positief_voorspeld['x'][-positief_voorspeld['avgsize'] - 7:],
positief_voorspeld['y'][-positief_voorspeld['avgsize'] - 7:],
color='fuchsia',
linestyle=':')
huidigpercentage = decimalstring(round(positief_percentage['y'][-1], 1))
ax2.plot(positief_percentage['x'],
# ax1.grid(which='both', axis='both', linestyle='-.',
# color='gray', linewidth=1, alpha=0.3)
ax1.set_xlabel("Datum")
ax1.set_ylabel("Prikken per dag")
ax1.stackplot(
vaccins_delta['x'],
vaccins_delta['astra_zeneca'],
vaccins_delta['pfizer'],
vaccins_delta['cure_vac'],
vaccins_delta['janssen'],
vaccins_delta['moderna'],
vaccins_delta['sanofi'],
labels=('COVID-19 Vaccine AstraZeneca ® (%s)' %
decimalstring(vaccins_delta['astra_zeneca'][-1]),
'Comirnaty® (BioNTech/Pfizer) (%s)' %
decimalstring(vaccins_delta['pfizer'][-1]),
'CVnCoV (CureVac) (%s)' %
decimalstring(vaccins_delta['cure_vac'][-1]),
'COVID-19 Vaccine Janssen (%s)' %
decimalstring(vaccins_delta['janssen'][-1]),
'COVID-19 Vaccine Moderna ® (%s)' %
decimalstring(vaccins_delta['moderna'][-1]),
'Sanofi/GSK (%s)' % decimalstring(vaccins_delta['sanofi'][-1])),
colors=(
'mediumblue',
'deepskyblue',
'darkorange',
'tomato',
'yellow',
color='gray',
linewidth=1,
alpha=0.3)
ax2.grid(which='both',
axis='both',
linestyle='-.',
color='gray',
linewidth=1,
alpha=0.3)
nu_opgenomen = opgenomen['y'][-1]
ax1.plot(opgenomen['x'],
opgenomen['y'],
color='orange',
label='aantal op verpleegafdeling (nu: ' +
decimalstring(nu_opgenomen) + ')')
ax1.plot(opgenomen_voorspeld['x'],
opgenomen_voorspeld['y'],
color='orange',
linestyle=':')
nu_op_ic = ic['y'][-1]
ax1.plot(ic['x'],
ic['y'],
color='red',
label='aantal op IC (nu: ' + decimalstring(nu_op_ic) + ')')
ax1.plot(ic_voorspeld['x'], ic_voorspeld['y'], color='red', linestyle=':')
# Plot ziek based on RNA
ax2.plot(geschat_ziek_rna['x'],
geschat_ziek_rna['y'],
y=0,
s=datetime.datetime.now().strftime("%d"),
color="white",
fontsize=8,
ha="center",
va="center",
bbox=dict(boxstyle='round,pad=0.1', facecolor='red', alpha=1, edgecolor='red'),
zorder=10
)
plt.axvline(dateCache.today(), color='red', linewidth=0.5)
# Laat afkappunt R zien (is twee weken geleden)
plt.axvline(dateCache.today() - datetime.timedelta(days=14), color='blue', linestyle='--', linewidth=0.5)
plt.annotate(
decimalstring(Rt_avg['y'][-1]),
xy=(Rt_avg['x'][-1], Rt_avg['y'][-1]),
xytext=(Rt_avg['x'][-1], Rt_avg['y'][-1]+0.5),
fontsize=8,
bbox=dict(boxstyle='round,pad=0.4', fc='ivory', alpha=1),
arrowprops=dict(arrowstyle='->', connectionstyle='arc3,rad=0.1')
)
axes = plt.gca()
axes.set_ylim([0,3])
axes.set_xlim([date_range[0],date_range[-1]])
axes.set_xlabel("Datum")
axes.set_ylabel("Reproductiegetal")
gegenereerd_op=datetime.datetime.now().strftime("%Y-%m-%d %H:%M")
linestyle='-.',
color='gray',
linewidth=1,
alpha=0.3)
nu_opgenomen = opgenomen['y'][-1]
nu_op_ic = ic['y'][-1]
nu_totaal = totaal[-1]
plt.stackplot(
opgenomen['x'],
ic['y'],
opgenomen['y'],
colors=['red', 'orange'],
labels=[
'aantal op IC (nu: ' + decimalstring(nu_op_ic) + ')',
'aantal op verpleegafdeling (nu: ' + decimalstring(nu_opgenomen) + ')'
])
plt.plot(opgenomen['x'],
totaal,
color='black',
label='Totaal (%s)' % decimalstring(nu_totaal))
# laat huidige datum zien met vertikale lijn
plt.figtext(0.885,
0.125,
datetime.datetime.now().strftime("%d"),
color="red",
fontsize=8,
bbox=dict(facecolor='white', alpha=0.9, pad=0, edgecolor='white'),
# Plot cases per dag
ax1.plot(alpha['x'], alpha['y'], color='red', label='Schatting zieken RIVM')
# ax1.plot(beta['x'], beta['y'], color='steelblue', label='beta')
# ax2.plot(alpha['x'], alpha['y'], c olor='steelblue', label='alpha')
# ax1.plot(beta['x'], beta['y'], color='blue', label='Schatting zieken @rolfje obv RNA in RWZI')
ax1.plot(beta['x'], beta['y'], color='blue', label='Schatting zieken @rolfje obv tests en RNA in RWZI')
ax1.fill_between(
beta['x'],
beta['min'],
beta['max'],
facecolor='blue', alpha=0.1, interpolate=True)
ax1.plot(teta['x'], teta['y'], color='cyan', label='Infectieradar COVID-19 klachten (percentage x'+decimalstring(tetascale)+')')
ax2.plot(gamma['x'], gamma['y'], color='lightgreen', label='Afwijking @rolfje t.o.v. RIVM schatting in %')
ax1.set_ylim([0, 400000])
ax2.set_ylim([-200, 200])
ax1.set_xlabel("Datum")
gegenereerd_op=datetime.datetime.now().strftime("%Y-%m-%d %H:%M")
plt.title('Afwijking schatting @rolfje t.o.v. RIVM schatting, '+gegenereerd_op)
ax1.legend(loc="upper left")
ax2.legend(loc="upper right")
plt.savefig("../docs/graphs/schattingafwijking.svg", format="svg")
#plt.show()
with open('prikafspraak/geboortejaar.txt') as f:
vaccingeboortejaar = f.read().rstrip()
filename = '../cache/stats.csv'
with open(filename, 'r') as csv_file:
csv_reader = csv.reader(csv_file, delimiter=';')
for row in csv_reader:
if row[0] == 'python_lines':
python_lines = int(row[1])
if row[0] == 'cache_size':
cache_size = int(row[1]) * 1000
substitutes = {
'totaal_positief':
decimalstring(totaal_positief),
'totaal_positief_num':
totaal_positief,
'totaal_positief_datum':
totaal_positief_datum,
'geschat_ziek_rivm':
decimalstring(geschat_ziek_nu),
'geschat_ziek_nu_datum':
geschat_ziek_nu_datum,
'geschat_ziek_rna':
decimalstring(round(geschat_ziek_nu_rna)),
'geschat_ziek_nu_rna_datum':
geschat_ziek_nu_rna_datum,
'geschat_ziek_rolf':
decimalstring(round(geschat_ziek_nu_rolf)),
'geschat_ziek_nu_rolf_datum':
# Variants in order of today's highest
today = {}
for code in varianten_totaal:
if code != 'x':
today[code] = varianten_totaal[code][-1]
today = dict(sorted(today.items(), key=lambda x: x[1]))
yrray = []
ylabels = []
# yrray.append(varianten_totaal['totaal'])
# ylabels.append('Totaal')
yrray.append(varianten_totaal['overig'])
ylabels.append("Overig (nu: %s)" %
(decimalstring(round(varianten_totaal['overig'][-1]))))
for code in today:
if code in varianten_totaal and code in variantcodes:
yrray.append(varianten_totaal[code])
ylabels.append("%s (nu: %s)" %
(variantcodes[code],
decimalstring(round(varianten_totaal[code][-1]))))
ax1.stackplot(
varianten_totaal['x'],
*yrray,
labels=ylabels,
colors=
( # note these are in reverse order because we revert the labels later
'gray',