Skip to content

Commit 60092df

Browse files
committed
[DOC] Trying to fix theory of the second prob tutorial
1 parent 417bc8a commit 60092df

File tree

7 files changed

+5957
-31
lines changed

7 files changed

+5957
-31
lines changed

docs/source/_static/Thickness_example.svg

Lines changed: 5919 additions & 0 deletions
Loading
-44.2 KB
Binary file not shown.
134 KB
Loading

examples/2-basic_geology/2-thickness_problem_gempy.py

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,13 @@
1+
"""
2+
2.2 - Including GemPy
3+
=====================
4+
5+
Complex probabilistic model
6+
---------------------------
7+
8+
"""
19
from pyro.infer.autoguide import init_to_mean
10+
211
from pyro.infer.inspect import get_dependencies
312

413
import gempy as gp

examples/2-basic_geology/README.rst

Lines changed: 1 addition & 31 deletions
Original file line numberDiff line numberDiff line change
@@ -13,48 +13,18 @@ Lucky for us, after we perform the first inference on the thickness, :math:`\til
1313

1414
.. math::
1515
:nowrap:
16-
1716
\pi(\theta_{thickness}) = \pi(\theta_{top}) - \pi(\theta_{bottom})
1817
1918
or,
2019

2120
.. math::
2221
:nowrap:
23-
2422
\pi(\theta_{bottom}) = \pi(\theta_{top}) - \pi(\theta_{thickness})
2523
2624
now the question is which probabilistic model design is more suitable. In the end this relates directly to the question the model is trying to answer---and possible limitations on the algorithms used---since joint probability follows the commutative and associative properties.
2725

28-
.. figure:: Thickness_example.eps
26+
.. figure:: ../_static/Thickness_example.svg
2927
:align: center
3028
:width: 80%
3129

3230
Prior (in blue) and posterior distributions (in red) of all the parameters of the probabilistic model :math:`\theta` (the full list of values of the simulation can be found in Appendix X). Case A) correspond to the subtraction equation while Case B) has been generated with the summation equation. Both cases yield the exact same posterior for the random variables---within the Monte Carlo error---as expected.
33-
34-
**Correlations**
35-
36-
Furthermore, keep in mind that the inference happens in a multidimensional manifold---although for obvious visualization reasons we display each parameter independently. This means that although in the previous figures the posterior distributions seem independent, a closer examination---e.g. joint plots---uncovers correlation in the parameters.
37-
38-
.. figure:: correlation.eps
39-
:align: center
40-
:width: 80%
41-
42-
Correlations between parameters can be uncovered with joint plots.
43-
44-
.. list-table::
45-
:widths: 10 30
46-
:header-rows: 0
47-
48-
* -
49-
- Here we are showing how to combine geometric data but with the same pattern we could encapsulate geological knowledge in form of geometrical rules
50-
51-
License
52-
-------
53-
The code in this case study is copyrighted by Miguel de la Varga and licensed under the new BSD (3-clause) license:
54-
55-
https://opensource.org/licenses/BSD-3-Clause
56-
57-
The text and figures in this case study are copyrighted by Miguel de la Varga and licensed under the CC BY-NC 4.0 license:
58-
59-
https://creativecommons.org/licenses/by-nc/4.0/
60-
Make sure to replace the links with actual hyperlinks if you're using a platform that supports it (e.g., Markdown or HTML). Otherwise, the plain URLs work fine for plain text.
Lines changed: 28 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,28 @@
1+
Correlations
2+
````````````
3+
4+
Furthermore, keep in mind that the inference happens in a multidimensional manifold---although for obvious visualization reasons we display each parameter independently. This means that although in the previous figures the posterior distributions seem independent, a closer examination---e.g. joint plots---uncovers correlation in the parameters.
5+
6+
.. figure::../_static/correlation.png
7+
:align: center
8+
:width: 80%
9+
10+
Correlations between parameters can be uncovered with joint plots.
11+
12+
.. list-table::
13+
:widths: 10 30
14+
:header-rows: 0
15+
16+
* -
17+
- Here we are showing how to combine geometric data but with the same pattern we could encapsulate geological knowledge in form of geometrical rules
18+
19+
License
20+
-------
21+
The code in this case study is copyrighted by Miguel de la Varga and licensed under the new BSD (3-clause) license:
22+
23+
https://opensource.org/licenses/BSD-3-Clause
24+
25+
The text and figures in this case study are copyrighted by Miguel de la Varga and licensed under the CC BY-NC 4.0 license:
26+
27+
https://creativecommons.org/licenses/by-nc/4.0/
28+
Make sure to replace the links with actual hyperlinks if you're using a platform that supports it (e.g., Markdown or HTML). Otherwise, the plain URLs work fine for plain text.

examples/2-basic_geology/next/__init__.py

Whitespace-only changes.

0 commit comments

Comments
 (0)