Mark number of medically relevant cases in graph

plot.py

@@ -19,6 +19,7 @@ data_source               = config['most_recent_values']['data_source']
 
 mortality_rate_percent    = float(config['assumptions']['mortality_rate_percent'])
 severe_cases_rate_percent = float(config['assumptions']['severe_cases_rate_percent'])
+medically_relevant_cases_percent = float(config['assumptions']['medically_relevant_cases_percent'])
 available_icus            = int(config['assumptions']['available_icus'])
 available_beds            = int(config['assumptions']['available_beds'])
 recovery_period           = int(config['assumptions']['recovery_period'])
@@ -50,10 +51,12 @@ linewidth=0.5
 # 0.6% of treated cases die
 # 5% of untreated cases die
 infected_to_severe = dead / infected * severe_cases_rate_percent / mortality_rate_percent
+infected_to_medically_relevant = dead / infected * medically_relevant_cases_percent / mortality_rate_percent
 
 ax1.plot(data['actual'], color='g', label="Nachgewiesene Infizierte -> " + data_source, marker='2', markersize=4, linestyle="None", linewidth=linewidth)
 ax1.plot(data['gompertz'], color='g', label="Schätzung Infizierte -> Gompertz-Kurve\n" + str(recovery_period) + " Tage bis Genesung, Sättigung bei " + str(saturation_percentage) + "%", linewidth=linewidth)
 ax1.plot(data['gompertz'] * infected_to_severe, color='b', label="Schätzung schwere Fälle (" + str(severe_cases_rate_percent) + "%)", linestyle="dotted", linewidth=linewidth)
+ax1.plot(data['gompertz'] * infected_to_medically_relevant, color='b', label="Schätzung auf Behandlung angewiesene Fälle (" + str(medically_relevant_cases_percent) + "%)", linestyle="dashdot", linewidth=linewidth)
 #
 # http://www.gbe-bund.de/oowa921-install/servlet/oowa/aw92/dboowasys921.xwdevkit/xwd_init?gbe.isgbetol/xs_start_neu/&p_aid=3&p_aid=58665717&nummer=838&p_sprache=D&p_indsp=-&p_aid=83716383
 #

src/main.rs

@@ -130,6 +130,8 @@ fn main() -> std::io::Result<()> {
     let dead = days[last_day_index].dead;
     let mortality_rate_percent = 0.6;
     let severe_cases_rate_percent = 5.0;
+    // https://www.medicalnewstoday.com/articles/coronavirus-81-of-cases-are-mild-study-says
+    let medically_relevant_cases_percent = 19.1;
     let available_icus = 28000 + 8000;
     let available_beds = 497000;
     conf.with_section(Some("most_recent_values"))
@@ -140,6 +142,7 @@ fn main() -> std::io::Result<()> {
     conf.with_section(Some("assumptions"))
         .set("mortality_rate_percent", mortality_rate_percent.to_string())
         .set("severe_cases_rate_percent", severe_cases_rate_percent.to_string())
+        .set("medically_relevant_cases_percent", medically_relevant_cases_percent.to_string())
         .set("available_icus", available_icus.to_string())
         .set("available_beds", available_beds.to_string())
         .set("recovery_period", recovery_period.to_string());