Skip to content

Commit

Permalink
Edit corner plots
Browse files Browse the repository at this point in the history
  • Loading branch information
@WWGolay
WWGolay committed Oct 18, 2022
1 parent e487f37 commit 25f63a1
Show file tree
Hide file tree
Showing 54 changed files with 66 additions and 166 deletions.
Binary file modified _build/.doctrees/environment.pickle
View file
Binary file not shown.
Binary file modified _build/.doctrees/notebooks/HD283572_corner-plot.doctree
View file
Binary file not shown.
Binary file modified _build/.doctrees/notebooks/V410Tau_corner-plot.doctree
View file
Binary file not shown.
Binary file modified _build/.doctrees/notebooks/algol.doctree
View file
Binary file not shown.
Binary file modified _build/.doctrees/notebooks/rs-cvn_corner-plots.doctree
View file
Binary file not shown.
Binary file modified _build/html/_images/HD283572_corner-plot_8_0.png
View file
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified _build/html/_images/V410Tau_corner-plot_8_0.png
View file
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified _build/html/_images/algol_11_0.png
View file
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified _build/html/_images/rs-cvn_corner-plots_14_0.png
View file
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified _build/html/_images/rs-cvn_corner-plots_9_0.png
View file
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
10 changes: 3 additions & 7 deletions _build/html/_sources/notebooks/HD283572_corner-plot.ipynb
View file

Large diffs are not rendered by default.

10 changes: 3 additions & 7 deletions _build/html/_sources/notebooks/V410Tau_corner-plot.ipynb
View file

Large diffs are not rendered by default.

10 changes: 3 additions & 7 deletions _build/html/_sources/notebooks/algol.ipynb
View file

Large diffs are not rendered by default.

20 changes: 6 additions & 14 deletions _build/html/_sources/notebooks/rs-cvn_corner-plots.ipynb
View file

Large diffs are not rendered by default.

8 changes: 2 additions & 6 deletions _build/html/notebooks/HD283572_corner-plot.html
View file
Original file line number Diff line number Diff line change
Expand Up @@ -504,8 +504,6 @@ <h2>8.4. Plot<a class="headerlink" href="#plot" title="Permalink to this headlin
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="n">thin</span> <span class="o">=</span> <span class="mi">10000</span>
<span class="n">plt_labels</span> <span class="o">=</span> <span class="p">[</span><span class="s1">&#39;$L_</span><span class="si">{pwr}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$\delta$&#39;</span><span class="p">,</span> <span class="s1">&#39;$n_{e,pwr}$&#39;</span><span class="p">,</span> <span class="s1">&#39;$B_</span><span class="si">{pwr}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$\phi_</span><span class="si">{pwr}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$L_</span><span class="si">{th}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$T_e$&#39;</span><span class="p">,</span> <span class="s1">&#39;$n_{e,th}$&#39;</span><span class="p">,</span> <span class="s1">&#39;$B_</span><span class="si">{th}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$\phi_</span><span class="si">{th}</span><span class="s1">$&#39;</span><span class="p">]</span>
<span class="n">median_values</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">median</span><span class="p">(</span><span class="n">samples</span><span class="p">[::</span><span class="n">thin</span><span class="p">,</span> <span class="p">:],</span> <span class="n">axis</span><span class="o">=</span><span class="mi">0</span><span class="p">)</span>
<span class="n">winnerWalker</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">argmax</span><span class="p">(</span><span class="n">sampler</span><span class="o">.</span><span class="n">get_log_prob</span><span class="p">(</span><span class="n">flat</span><span class="o">=</span><span class="kc">True</span><span class="p">))</span>
<span class="n">mostParams</span> <span class="o">=</span> <span class="n">samples</span><span class="p">[</span><span class="n">winnerWalker</span><span class="p">]</span>
<span class="n">cmap</span> <span class="o">=</span> <span class="n">mpl</span><span class="o">.</span><span class="n">cm</span><span class="o">.</span><span class="n">get_cmap</span><span class="p">(</span><span class="s1">&#39;Reds&#39;</span><span class="p">)</span>
<span class="n">cmap</span><span class="o">.</span><span class="n">set_under</span><span class="p">(</span><span class="n">color</span><span class="o">=</span><span class="s1">&#39;w&#39;</span><span class="p">)</span>

Expand All @@ -515,10 +513,8 @@ <h2>8.4. Plot<a class="headerlink" href="#plot" title="Permalink to this headlin
<span class="n">levels</span><span class="o">=</span><span class="p">(</span><span class="mf">0.393</span><span class="p">,</span> <span class="mf">0.865</span><span class="p">,</span> <span class="mf">0.989</span><span class="p">),</span> <span class="n">hist_kwargs</span><span class="o">=</span><span class="p">{</span><span class="s1">&#39;color&#39;</span><span class="p">:</span><span class="n">cmap</span><span class="p">(</span><span class="mf">0.75</span><span class="p">)},</span>
<span class="n">contourf_kwargs</span><span class="o">=</span><span class="p">{</span><span class="s1">&#39;colors&#39;</span><span class="p">:(</span><span class="n">cmap</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">),</span> <span class="n">cmap</span><span class="p">(</span><span class="mf">0.25</span><span class="p">),</span> <span class="n">cmap</span><span class="p">(</span><span class="mf">0.75</span><span class="p">),</span> <span class="n">cmap</span><span class="p">(</span><span class="mf">0.999</span><span class="p">)),</span> <span class="s1">&#39;alpha&#39;</span><span class="p">:</span><span class="mf">0.75</span><span class="p">})</span>

<span class="n">corner</span><span class="o">.</span><span class="n">overplot_lines</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">lmfit_params</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C1&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#39;Least squares&#39;</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_points</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="n">lmfit_params</span><span class="p">]),</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C1&#39;</span><span class="p">,</span> <span class="n">marker</span><span class="o">=</span><span class="s1">&#39;s&#39;</span><span class="p">,</span> <span class="n">ms</span><span class="o">=</span><span class="mi">3</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_lines</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">mostParams</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C0&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#39;Most probable&#39;</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_points</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="n">mostParams</span><span class="p">]),</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C0&#39;</span><span class="p">,</span> <span class="n">marker</span><span class="o">=</span><span class="s1">&#39;s&#39;</span><span class="p">,</span> <span class="n">ms</span><span class="o">=</span><span class="mi">5</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_lines</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">lmfit_params</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C0&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#39;LMFIT&#39;</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_points</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="n">lmfit_params</span><span class="p">]),</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C0&#39;</span><span class="p">,</span> <span class="n">marker</span><span class="o">=</span><span class="s1">&#39;s&#39;</span><span class="p">,</span> <span class="n">ms</span><span class="o">=</span><span class="mi">5</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_lines</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">median_values</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;black&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#39;Median&#39;</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_points</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="n">median_values</span><span class="p">]),</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;black&#39;</span><span class="p">,</span> <span class="n">marker</span><span class="o">=</span><span class="s1">&#39;s&#39;</span><span class="p">,</span> <span class="n">ms</span><span class="o">=</span><span class="mi">5</span><span class="p">)</span>

Expand Down
8 changes: 2 additions & 6 deletions _build/html/notebooks/V410Tau_corner-plot.html
View file
Original file line number Diff line number Diff line change
Expand Up @@ -504,8 +504,6 @@ <h2>7.4. Plot<a class="headerlink" href="#plot" title="Permalink to this headlin
<div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="n">thin</span> <span class="o">=</span> <span class="mi">10000</span>
<span class="n">plt_labels</span> <span class="o">=</span> <span class="p">[</span><span class="s1">&#39;$L_</span><span class="si">{pwr}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$\delta$&#39;</span><span class="p">,</span> <span class="s1">&#39;$n_{e,pwr}$&#39;</span><span class="p">,</span> <span class="s1">&#39;$B_</span><span class="si">{pwr}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$\phi_</span><span class="si">{pwr}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$L_</span><span class="si">{th}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$T_e$&#39;</span><span class="p">,</span> <span class="s1">&#39;$n_{e,th}$&#39;</span><span class="p">,</span> <span class="s1">&#39;$B_</span><span class="si">{th}</span><span class="s1">$&#39;</span><span class="p">,</span> <span class="s1">&#39;$\phi_</span><span class="si">{th}</span><span class="s1">$&#39;</span><span class="p">]</span>
<span class="n">median_values</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">median</span><span class="p">(</span><span class="n">samples</span><span class="p">[::</span><span class="n">thin</span><span class="p">,</span> <span class="p">:],</span> <span class="n">axis</span><span class="o">=</span><span class="mi">0</span><span class="p">)</span>
<span class="n">winnerWalker</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">argmax</span><span class="p">(</span><span class="n">sampler</span><span class="o">.</span><span class="n">get_log_prob</span><span class="p">(</span><span class="n">flat</span><span class="o">=</span><span class="kc">True</span><span class="p">))</span>
<span class="n">mostParams</span> <span class="o">=</span> <span class="n">samples</span><span class="p">[</span><span class="n">winnerWalker</span><span class="p">]</span>
<span class="n">cmap</span> <span class="o">=</span> <span class="n">mpl</span><span class="o">.</span><span class="n">cm</span><span class="o">.</span><span class="n">get_cmap</span><span class="p">(</span><span class="s1">&#39;Purples&#39;</span><span class="p">)</span>
<span class="n">cmap</span><span class="o">.</span><span class="n">set_under</span><span class="p">(</span><span class="n">color</span><span class="o">=</span><span class="s1">&#39;w&#39;</span><span class="p">)</span>

Expand All @@ -515,10 +513,8 @@ <h2>7.4. Plot<a class="headerlink" href="#plot" title="Permalink to this headlin
<span class="n">levels</span><span class="o">=</span><span class="p">(</span><span class="mf">0.393</span><span class="p">,</span> <span class="mf">0.865</span><span class="p">,</span> <span class="mf">0.989</span><span class="p">),</span> <span class="n">hist_kwargs</span><span class="o">=</span><span class="p">{</span><span class="s1">&#39;color&#39;</span><span class="p">:</span><span class="n">cmap</span><span class="p">(</span><span class="mf">0.75</span><span class="p">)},</span>
<span class="n">contourf_kwargs</span><span class="o">=</span><span class="p">{</span><span class="s1">&#39;colors&#39;</span><span class="p">:(</span><span class="n">cmap</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">),</span> <span class="n">cmap</span><span class="p">(</span><span class="mf">0.25</span><span class="p">),</span> <span class="n">cmap</span><span class="p">(</span><span class="mf">0.75</span><span class="p">),</span> <span class="n">cmap</span><span class="p">(</span><span class="mf">0.999</span><span class="p">)),</span> <span class="s1">&#39;alpha&#39;</span><span class="p">:</span><span class="mf">0.75</span><span class="p">})</span>

<span class="n">corner</span><span class="o">.</span><span class="n">overplot_lines</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">lmfit_params</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C1&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#39;Least squares&#39;</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_points</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="n">lmfit_params</span><span class="p">]),</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C1&#39;</span><span class="p">,</span> <span class="n">marker</span><span class="o">=</span><span class="s1">&#39;s&#39;</span><span class="p">,</span> <span class="n">ms</span><span class="o">=</span><span class="mi">3</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_lines</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">mostParams</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C0&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#39;Most probable&#39;</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_points</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="n">mostParams</span><span class="p">]),</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C0&#39;</span><span class="p">,</span> <span class="n">marker</span><span class="o">=</span><span class="s1">&#39;s&#39;</span><span class="p">,</span> <span class="n">ms</span><span class="o">=</span><span class="mi">5</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_lines</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">lmfit_params</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C0&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#39;LMFIT&#39;</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_points</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="n">lmfit_params</span><span class="p">]),</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;C0&#39;</span><span class="p">,</span> <span class="n">marker</span><span class="o">=</span><span class="s1">&#39;s&#39;</span><span class="p">,</span> <span class="n">ms</span><span class="o">=</span><span class="mi">5</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_lines</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">median_values</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;black&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#39;Median&#39;</span><span class="p">)</span>
<span class="n">corner</span><span class="o">.</span><span class="n">overplot_points</span><span class="p">(</span><span class="n">cornerFig</span><span class="p">,</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="n">median_values</span><span class="p">]),</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#39;black&#39;</span><span class="p">,</span> <span class="n">marker</span><span class="o">=</span><span class="s1">&#39;s&#39;</span><span class="p">,</span> <span class="n">ms</span><span class="o">=</span><span class="mi">5</span><span class="p">)</span>

Expand Down
Loading

0 comments on commit 25f63a1

Please sign in to comment.