.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "gallery/miscellaneous/dirichlet_pdf.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_gallery_miscellaneous_dirichlet_pdf.py>`
        to download the full example code

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_gallery_miscellaneous_dirichlet_pdf.py:


======================
Dirichlet distribution
======================

.. GENERATED FROM PYTHON SOURCE LINES 6-43



.. image-sg:: /gallery/miscellaneous/images/sphx_glr_dirichlet_pdf_001.svg
   :alt: ${\mathbf{\alpha}}$ = (1.5, 1.5, 1.5), ${\mathbf{\alpha}}$ = (5.0, 5.0, 5.0), ${\mathbf{\alpha}}$ = (1.0, 2.0, 2.0), ${\mathbf{\alpha}}$ = (2.0, 4.0, 8.0)
   :srcset: /gallery/miscellaneous/images/sphx_glr_dirichlet_pdf_001.svg
   :class: sphx-glr-single-img





.. code-block:: Python

    import matplotlib.pyplot as plt
    from mpltern.datasets import get_dirichlet_pdfs

    fig = plt.figure(figsize=(10.8, 8.8))
    fig.subplots_adjust(
        left=0.1,
        right=0.9,
        bottom=0.1,
        top=0.9,
        wspace=0.5,
        hspace=0.5,
    )

    alphas = ((1.5, 1.5, 1.5), (5.0, 5.0, 5.0), (1.0, 2.0, 2.0), (2.0, 4.0, 8.0))
    for i, alpha in enumerate(alphas):
        ax = fig.add_subplot(2, 2, i + 1, projection="ternary")
        t, l, r, v = get_dirichlet_pdfs(n=61, alpha=alpha)
        cmap = "Blues"
        shading = "gouraud"
        cs = ax.tripcolor(t, l, r, v, cmap=cmap, shading=shading, rasterized=True)
        ax.tricontour(t, l, r, v, colors="k", linewidths=0.5)

        ax.set_tlabel("$x_1$")
        ax.set_llabel("$x_2$")
        ax.set_rlabel("$x_3$")

        ax.taxis.set_label_position("tick1")
        ax.laxis.set_label_position("tick1")
        ax.raxis.set_label_position("tick1")

        ax.set_title("${\\mathbf{\\alpha}}$ = " + str(alpha))

        cax = ax.inset_axes([1.05, 0.1, 0.05, 0.9], transform=ax.transAxes)
        colorbar = fig.colorbar(cs, cax=cax)
        colorbar.set_label("PDF", rotation=270, va="baseline")

    plt.show()


.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (0 minutes 4.117 seconds)


.. _sphx_glr_download_gallery_miscellaneous_dirichlet_pdf.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: dirichlet_pdf.ipynb <dirichlet_pdf.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: dirichlet_pdf.py <dirichlet_pdf.py>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_