markov chain notebooke changed

3bbcc4d7 · Mirko Birbaumer · cdb6ea7c · 3bbcc4d7
Commit 3bbcc4d7 authored 1 year ago by Mirko Birbaumer
--- a/notebooks/MCMC/markov_chain_weather.ipynb
+++ b/notebooks/MCMC/markov_chain_weather.ipynb
@@ -47,14 +47,6 @@
    "w"
   ]
  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "b9b378fb-493c-4dee-a007-2c0c18760ab0",
-   "metadata": {},
-   "outputs": [],
-   "source": []
-  },
  {
   "cell_type": "code",
   "execution_count": 10,

 %% Cell type:markdown id:d2be01c2-bb2a-490b-bade-56c0defc4a31 tags:

 ### Weather Simulation with Markov Chain

 %% Cell type:code id:0938a32d-af05-440e-a207-b1af52425e98 tags:

 ``` python
 import random
 # List with sunny and rainy
 d = ['S', 'R']
 # Define empty list
 w = []
 # Initialize first entry: Chooses randomly a value of list # d and adds (append) it to the empty list w
 w.append(random.choice(d))
 # Number of days
 n = 1000
 ```

 %% Cell type:code id:fd5f38eb-b522-4859-805b-d9cd4b0bd3b2 tags:

 ``` python
 w
 ```

 %% Output

    ['S']

-%% Cell type:code id:b9b378fb-493c-4dee-a007-2c0c18760ab0 tags:
-
-``` python
-```
-
 %% Cell type:code id:e98e5870-9f9e-446c-9714-e25d8585ff87 tags:

 ``` python
 # Markov chain
 for i in range(n):
 # If it's sunny on day i add randomly 'S' from d with
 # absweight .9 or 'R' with weight 0.1
    if w[i] == 'S':
        w.append(random.choices(d, weights=[0.9, 0.1])[0])
    # If it's rainy on day i add randomly 'S' with weight .5
    # or 'R' with weight 0.5
    else:
        w.append(random.choices(d, weights=[0.5, 0.5])[0])
        # Removes quotes ' ' in list w
 w_1 = ' '.join(w)
 w_1
 ```

 %% Output

    'S S S S S S S R S S S S S S S S S S S S S S S S S S S R S S S S S S S R R S S S S S S S S S S S S R R R R S S S S S S S S R S S S R R R S R R S S R R S S S S S S S S R S S S S S S S S S S S S S S S S S R S R R S S S S S S S R R S R R R R R S S S S S S S S R S S S S S S S S S S S S S S R S S S R R R R R R S S S S S S S S S S S S S S S R S S S S R R S S S S S S S S S S S S S S S S S R S S S S S S S S S S S R R S S S S R S S S S S S S S S S S S S S S S S S S S S S R R R S S R R S S S S S S S S S S S S S S S S S S R R R R S S S S S S S S S S S S S R R R R R R R S S S S S S S S S R R S S S S S S S S S S S S S S S S S S S S S R S S S S S S S S S S S S S S S S S S S S S R R S S R S S S S S S S S S S S S S S S S S S S S S S S S R R R R S S S S S R R S R S S S R R S S S S S S S S S S R S S S S S S S S S S S S S S S S R R R R S S S S S R S S R R S S S S S S R S S R S S S S S S S S S S S S R R S S S S S S S S S S S S S S S S S R S S S S S S S S S S S S S S R S S S S S S S S S S S S S S R S S R S S S S S S S S S S S S S S S R R S S S S S S S S S R S S S S S S S S S S S S S S S S S S S S S R S S S S S S S S S S S S S S S S R R R R R S S S S S S R S S R R R R S S S S S S S S S R R R S S S S S R R R S S S S S S S S S S S S S S S S S S S S S S S S S S S R S S S S S S S S S R S S S R R R R R S S S R S S S S S R R R R S S S S S S S S S S S S S S S S S S S S S S S S R S S S S S S R R R R S R S S R R R R R R R S S S S S S S S S S R R R R R S S S S S S S S S S S S S R S S S S S S S S S S S S S S S S S R S S S S S S R R S S S S S S S S S S S S R S S S R S S S S S S S S S S S S S S S S S S S S S S S R S S S R R S S S S S S S S S S S S S R R R R S S R S S S S S S S S S S S S S S S S S S S S S S S S R S S S S S R S S R S S S S S S S S S R S S S S S S S S S S S S S S S S S S S R R S S S R S S S S S S S S S S S S S R R R R R S S S S R S S S S S S S S S S S S S S S S S S S S R R S S S S S R R S S S S S S S R R S S S S S R S S S S S S S S S R R S S S S S S S S S S S S S S S S S R R S S S S S S S'

 %% Cell type:code id:6c8e96be-29f4-4ac7-a230-38fbccdcaa71 tags:

 ``` python
 # Probability of sunny days
 w.count('S') / len(w)
 ```

 %% Output

    0.8231768231768232