From 2a1e80f8ebca7ee936e730602dff0da2a968eeac Mon Sep 17 00:00:00 2001 From: markus Date: Mon, 25 Jan 2021 13:22:30 +0100 Subject: [PATCH] use new distribution for viral load --- cara/montecarlo.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/cara/montecarlo.py b/cara/montecarlo.py index 21e1361b..e245cd1f 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -77,7 +77,7 @@ class MCInfectedPopulation(models.Population): emissions = sct.weibull_min.isf(sct.norm.sf(np.random.normal(size=self.samples)), e_k, loc=0, scale=e_lambda) diameters = sct.weibull_min.isf(sct.norm.sf(np.random.normal(size=self.samples)), d_k, loc=0, scale=d_lambda) if self.viral_load is None: - viral_loads = np.random.normal(loc=7.8, scale=1.7, size=self.samples) + viral_loads = self._generate_viral_loads() else: viral_loads = np.full(self.samples, self.viral_load)