added mira distribution loading

ciemss · Jul 3, 2023 · ba31252 · ba31252
1 parent 7fc261a
commit ba31252
Show file tree

Hide file tree

Showing 8 changed files with 1,651 additions and 1,597 deletions.
diff --git a/notebook/integration_demo/demo.ipynb b/notebook/integration_demo/demo.ipynb
@@ -2,40 +2,20 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 1,
    "metadata": {},
    "outputs": [],
    "source": [
     "import os\n",
     "from pyciemss.PetriNetODE.interfaces import (\n",
     "    load_and_sample_petri_model,\n",
-    "    load_and_calibrate_and_sample_petri_model,\n",
-    "    load_petri_model,\n",
-    "    setup_petri_model,\n",
-    "    sample\n",
-    ")\n",
-    "import numpy as np\n",
-    "from typing import Iterable\n",
-    "from pyciemss.utils.interface_utils import (\n",
-    "    assign_interventions_to_timepoints,\n",
-    "    interventions_and_sampled_params_to_interval,\n",
-    "    convert_to_output_format\n",
-    ")\n",
-    "from torch import tensor"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "cd ../.."
+    "    load_and_calibrate_and_sample_petri_model\n",
+    ")"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 2,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -53,7 +33,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 3,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -69,7 +49,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 4,
    "metadata": {
     "scrolled": true
    },
@@ -86,8 +66,7 @@
     "intervened_samples.to_csv(\n",
     "    os.path.join(DEMO_PATH, \"results_petri/sample_results_w_interventions.csv\"),\n",
     "    index=False,\n",
-    ")\n",
-    "intervened_samples"
+    ")"
    ]
   },
   {
@@ -100,9 +79,56 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 5,
    "metadata": {},
-   "outputs": [],
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "iteration 0: loss = 36.72989881038666\n",
+      "iteration 25: loss = 35.35893213748932\n",
+      "iteration 50: loss = 32.896029233932495\n",
+      "iteration 75: loss = 33.40498435497284\n",
+      "iteration 100: loss = 33.192803621292114\n",
+      "iteration 125: loss = 31.314371824264526\n",
+      "iteration 150: loss = 32.697070837020874\n",
+      "iteration 175: loss = 33.36095118522644\n",
+      "iteration 200: loss = 32.870508432388306\n",
+      "iteration 225: loss = 32.567596197128296\n",
+      "iteration 250: loss = 32.83535122871399\n",
+      "iteration 275: loss = 32.43224883079529\n",
+      "iteration 300: loss = 32.42908525466919\n",
+      "iteration 325: loss = 32.09095215797424\n",
+      "iteration 350: loss = 32.357651472091675\n",
+      "iteration 375: loss = 32.677024126052856\n",
+      "iteration 400: loss = 32.180010080337524\n",
+      "iteration 425: loss = 32.30828094482422\n",
+      "iteration 450: loss = 32.41364550590515\n",
+      "iteration 475: loss = 32.47804093360901\n",
+      "iteration 500: loss = 32.483015298843384\n",
+      "iteration 525: loss = 32.19251108169556\n",
+      "iteration 550: loss = 32.23183989524841\n",
+      "iteration 575: loss = 32.15587258338928\n",
+      "iteration 600: loss = 32.164281606674194\n",
+      "iteration 625: loss = 32.375654220581055\n",
+      "iteration 650: loss = 32.25739574432373\n",
+      "iteration 675: loss = 32.28496503829956\n",
+      "iteration 700: loss = 32.22272062301636\n",
+      "iteration 725: loss = 32.60892415046692\n",
+      "iteration 750: loss = 32.329415798187256\n",
+      "iteration 775: loss = 32.2548828125\n",
+      "iteration 800: loss = 32.61113715171814\n",
+      "iteration 825: loss = 32.33934807777405\n",
+      "iteration 850: loss = 32.35854434967041\n",
+      "iteration 875: loss = 32.08857345581055\n",
+      "iteration 900: loss = 32.40586733818054\n",
+      "iteration 925: loss = 32.314961194992065\n",
+      "iteration 950: loss = 32.23428010940552\n",
+      "iteration 975: loss = 32.46005034446716\n"
+     ]
+    }
+   ],
    "source": [
     "data_path = os.path.join(DEMO_PATH, \"data.csv\")\n",
     "num_samples = 100\n",
@@ -122,20 +148,13 @@
     "    os.path.join(DEMO_PATH, \"results_petri/calibrated_sample_results.csv\"), index=False\n",
     ")"
    ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": []
   }
  ],
  "metadata": {
   "kernelspec": {
-   "display_name": "pyciemss-main",
+   "display_name": "pyciemss-deriv",
    "language": "python",
-   "name": "pyciemss-main"
+   "name": "python3"
   },
   "language_info": {
    "codemirror_mode": {
@@ -147,7 +166,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.10.9"
+   "version": "3.10.11"
   }
  },
  "nbformat": 4,