Switch to lognormal for defaults

pchlenski · Feb 17, 2022 · dae28e9 · dae28e9
1 parent 39e7710
commit dae28e9
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/src/OmicsGenerator.py b/src/OmicsGenerator.py
@@ -961,7 +961,8 @@ def generate_multiple(
         # Initialize first sample distribution
         if initial_distribution is None:
             def initial_distribution(size):
-                return np.random.exponential(size=size) * np.random.binomial(1, 1-extinct_fraction, size=size)
+                # return np.random.exponential(size=size) * np.random.binomial(1, 1-extinct_fraction, size=size)
+                return np.random.lognormal(size=size) * np.random.binomial(1, 1-extinct_fraction, size=size)
 
         # Generation loop
         for i in range(n):
@@ -1115,7 +1116,8 @@ def _random(size) -> np.ndarray:
 
     def _init_full(
         self,
-        initial_distribution : callable = np.random.exponential,
+        # initial_distribution : callable = np.random.exponential,
+        initial_distribution : callable = np.random.lognormal,
         growth_rate_distribution : callable = _random,
         **kwargs) -> None:
         """