From 4d5337377dae5c825b14730b741bbcf87d5a87b1 Mon Sep 17 00:00:00 2001 From: markus Date: Thu, 4 Feb 2021 12:31:13 +0100 Subject: [PATCH] add mask_leak_out --- cara/montecarlo.py | 13 +++++++++++++ 1 file changed, 13 insertions(+) diff --git a/cara/montecarlo.py b/cara/montecarlo.py index 62e83391..b86b633b 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -68,6 +68,19 @@ concentration_vs_diameter = ( ) +def mask_leak_out(d: float) -> float: + if d < 0.5: + return 1 + + if d < 1: + return 1 - (0.5893 * d + 0.1546) + + if d < 3: + return 1 - (0.0509 * d + 0.0664) + + return 1 - 0.8167 + + log_viral_load_frequencies = scoeh_vl_frequencies if USE_SCOEH else symptomatic_vl_frequencies