pingu_98/cara
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

use new distribution for viral load

Browse source
This commit is contained in:
markus 2021-01-25 13:22:30 +01:00
parent 4fe50ea22a
commit 2a1e80f8eb
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
cara/montecarlo.py
View file

@ -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)