Skip to content
Snippets Groups Projects
Commit 0647ae89 authored by freba's avatar freba :call_me:
Browse files

Updated documentation.org (improved clarity).

parent 470ed0dc
No related branches found
No related tags found
No related merge requests found
......@@ -483,11 +483,18 @@ filepath = plot_one_step_prediction_map(mp, xs, us, cluster=[0])
#+end_src
* TODO Markov Chain Monte Carlo
* Markov Chain Monte Carlo
The code blocks in this section associate with Figures 4, 5 and 6. The code can be found in the listed scripts, respectively:
- =c0_stan.py= :: this code block is associated with subfigure (a) in Figure 4.
- =c0_stan_ppc.py= :: this code block is associated with subfigure (a) in Figure 5.
- =c1_stan.py= :: this code block is associated with subfigure (b) in Figure 4.
- =c1_stan_ppc.py= :: this code block is associated with subfigure (b) in Figure 5.
- =stan_cosimulate_ID_c0= :: this code block is associated with Figure 6.
This code block is associated with subfigure (a) in Figure 4.
#+NAME: c0_stan
#+begin_src python :noweb strip-export :exports both :results file :tangle ./Scripts/c0_stan.py :return filepath
#+begin_src python :noweb strip-export :exports none :results file :tangle ./Scripts/c0_stan.py :return filepath
cluster=0
<<sindyc_header>>
<<code_stan_fitModel>>
......@@ -537,9 +544,8 @@ filepath = Plotter.plot_density(
#+end_src
This code block is associated with subfigure (a) in Figure 5.
#+NAME: c0_stan_ppc
#+begin_src python :exports both :results file :tangle ./Scripts/c0_stan_ppc.py :noweb yes :return filepath
#+begin_src python :exports none :results file :tangle ./Scripts/c0_stan_ppc.py :noweb yes :return filepath
import numpy as np
<<paper_plotting_header>>
import arviz as az
......@@ -588,9 +594,8 @@ filepath = Plotter.plot_ppc(
#+end_src
This code block is associated with subfigure (b) in Figure 4.
#+NAME: c1_stan
#+begin_src python :noweb strip-export :exports both :results file :tangle ./Scripts/c1_stan.py :return filepath
#+begin_src python :noweb strip-export :exports none :results file :tangle ./Scripts/c1_stan.py :return filepath
cluster=1
<<simo_IDCL_sindy_header>>
<<code_stan_fitModel>>
......@@ -637,8 +642,6 @@ filepath = Plotter.plot_density(stan_results,
#+END_SRC
This code block is associated with subfigure (b) in Figure 5.
#+NAME: c1_stan_ppc
#+begin_src python :exports code :results file :tangle ./Scripts/c1_stan_ppc.py :noweb yes :return filepath
import numpy as np
......@@ -688,10 +691,8 @@ filepath = Plotter.plot_ppc(
#+end_src
This code block is associated with Figure 6.
#+NAME: stan_cosimulate_ID_c0
#+begin_src python :results file :exports both :return filename :tangle ./Scripts/stan_cosimulate_ID_c0.py :noweb yes
#+begin_src python :results file :exports none :return filename :tangle ./Scripts/stan_cosimulate_ID_c0.py :noweb yes
from collections import OrderedDict
import numpy as np
......@@ -743,7 +744,16 @@ filename = Plotter.plot_sectors(
#+end_src
* WAIT Stan submodels
* Stan submodels
The code blocks in this section represent the Stan models with increasing complexity. The code can be found in the listed scripts, respectively:
- =sindyc_stan_code_12.stan= :: Stan model with active terms 1 and 2.
- =sindyc_stan_code_23.stan= :: Stan model with active terms 2 and 3.
- =sindyc_stan_code_34.stan= :: Stan model with active terms 3 and 4.
- =sindyc_stan_code_345.stan= :: Stan model with active terms 3, 4 and 5.
- =sindyc_stan_code_124.stan= :: Stan model with active terms 1, 2 and 4.
- =sindyc_stan_code_1234.stan= :: Stan model with active terms 1, 2, 3 and 4.
#+NAME: sindyc_stan_code_12
#+begin_src stan :exports code :results none :tangle ./Scripts/sindyc_stan_code_12.stan
......@@ -1383,7 +1393,7 @@ def read_from_disc(filename="save.pkl"):
Common prosumer response dynamics.
#+NAME: paper_prosumerDynamics
#+begin_src python :results none :exports code
#+begin_src python :results none :exports none
class Prosumer:
"""Common prosumer response model.
"""
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment