Merge branch 'main' gf github.com:Priesemann-Group/testing-artefacts-pharmaco-multipath

forgot to push

Author: laura-mueller

notebooks/bigHIVmodel.ipynb

@@ -99,8 +99,11 @@ def duration2rate(x):
     [14.0,  24.0,  167.0,  203.0],
     [15.0,  40.0,  120.0,  203.0],
     [10.0,  38.0,  141.3, 203.0],
+    [2.0, 80.0, 53.3, 200.0],
+    [2.0, 63.5, 100.0, 200.0],
+    [75, 39, 18.5, 2],
 ])
-c = sets_of_c[2]
+c = sets_of_c[2] # 2, 5, 7
 
 #H = 5.0 # HIV hazard
 #P = 50.0 # PrEP fraction

notebooks/min_PrEP_VS_lambda_P.ipynb

@@ -19,6 +19,22 @@ def save_results(results, filename):
     np.save(f"../results/{filename}.npy", results)
     logger.info(f"Results saved to ../results/{filename}.npy")
 
+# Load the results from a file
+def load_results(filename):
+    """
+    Load the results from a file.
+
+    Parameters:
+    filename (str): Name of the file to load the results from.
+
+    Returns:
+    dict: Loaded results.
+    """
+    # Load the results from a file
+    results = np.load(f"../results/{filename}.npy", allow_pickle=True).item()
+    logger.info(f"Results loaded from ../results/{filename}.npy")
+    return results
+
 
 # Define a function to read and evaluate the file
 # PD: using np and jnp in calculation of args in HIVandSTI.py (in fraction2rate() etc.)