renaming MCMC notebooks

101f5b5e · Mirko Birbaumer · 0db1b502 · 101f5b5e · 101f5b5e · 101f5b5e
Commit 101f5b5e authored 6 months ago by Mirko Birbaumer
--- a/notebooks/MCMC/markov_chain_weather.ipynb
+++ b/notebooks/MCMC/markov_chain_weather.ipynb
--- a/notebooks/MCMC/metro_walk.ipynb
+++ b/notebooks/MCMC/metro_walk.ipynb
--- a/notebooks/MCMC/MCMC_3_5_2.ipynb
+++ b/notebooks/MCMC/MCMC_3_5_2.ipynb
--- a/notebooks/MCMC/metro_hist.ipynb
+++ b/notebooks/MCMC/metro_hist.ipynb
--- a/notebooks/MCMC/Untitled.ipynb
+++ b/notebooks/MCMC/Untitled.ipynb
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "id": "9b27cc58-6e9d-4e3a-9c5b-a9d1f2080022",
+   "metadata": {
+    "tags": []
+   },
+   "outputs": [],
+   "source": [
+    "from scipy.stats import uniform\n",
+    "n = 1000\n",
+    "x = uniform.rvs(size=1000, loc=0, scale=1)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 2,
+   "id": "98bb71fa-bc54-4e6b-a6c5-5639ffc3b16d",
+   "metadata": {
+    "tags": []
+   },
+   "outputs": [],
+   "source": [
+    "import numpy as np\n",
+    "from scipy.stats import uniform\n",
+    "n = 1000\n",
+    "x = uniform.rvs(size=1000, loc=0, scale=1)\n",
+    "integrand = np.exp(-np.square(x)/2)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 7,
+   "id": "e9d73b09-de8f-472b-b448-c20b35462143",
+   "metadata": {
+    "tags": []
+   },
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "0.34320423614573237\n"
+     ]
+    }
+   ],
+   "source": [
+    "from math import pi\n",
+    "import numpy as np\n",
+    "from scipy.stats import uniform\n",
+    "n = 1000\n",
+    "x = uniform.rvs(size=1000, loc=0, scale=1)\n",
+    "integrand = np.exp(-np.square(x)/2)\n",
+    "integral = 1/(np.sqrt(2*pi))*np.sum(integrand)/n\n",
+    "print(integral)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 4,
+   "id": "7eee71cd-16c8-484e-8171-ee506ae2b254",
+   "metadata": {
+    "tags": []
+   },
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "0.3413447460685429\n"
+     ]
+    }
+   ],
+   "source": [
+    "from scipy.stats import norm\n",
+    "integral = norm.cdf(1)-norm.cdf(0)\n",
+    "print(integral)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "d7083f3a-e905-417a-9afb-7ac22debd3a1",
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.10.11"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}
+%% Cell type:code id:9b27cc58-6e9d-4e3a-9c5b-a9d1f2080022 tags:
+
+``` python
+from scipy.stats import uniform
+n = 1000
+x = uniform.rvs(size=1000, loc=0, scale=1)
+```
+
+%% Cell type:code id:98bb71fa-bc54-4e6b-a6c5-5639ffc3b16d tags:
+
+``` python
+import numpy as np
+from scipy.stats import uniform
+n = 1000
+x = uniform.rvs(size=1000, loc=0, scale=1)
+integrand = np.exp(-np.square(x)/2)
+```
+
+%% Cell type:code id:e9d73b09-de8f-472b-b448-c20b35462143 tags:
+
+``` python
+from math import pi
+import numpy as np
+from scipy.stats import uniform
+n = 1000
+x = uniform.rvs(size=1000, loc=0, scale=1)
+integrand = np.exp(-np.square(x)/2)
+integral = 1/(np.sqrt(2*pi))*np.sum(integrand)/n
+print(integral)
+```
+
+%% Output
+
+    0.34320423614573237
+
+%% Cell type:code id:7eee71cd-16c8-484e-8171-ee506ae2b254 tags:
+
+``` python
+from scipy.stats import norm
+integral = norm.cdf(1)-norm.cdf(0)
+print(integral)
+```
+
+%% Output
+
+    0.3413447460685429
+
+%% Cell type:code id:d7083f3a-e905-417a-9afb-7ac22debd3a1 tags:
+
+``` python
+```
--- a/notebooks/MCMC/metro_beta_walk_hist.ipynb
+++ b/notebooks/MCMC/metro_beta_walk_hist.ipynb