.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "auto_data_analysis/statistical_tests/plot_fitted_distribution_ranking.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note Click :ref:`here ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_auto_data_analysis_statistical_tests_plot_fitted_distribution_ranking.py: Select fitted distributions =========================== .. GENERATED FROM PYTHON SOURCE LINES 7-14 In this example help to make a choice between several distributions fitted to a sample. Several methods can be used: - the ranking by the Kolmogorov p-values (for continuous distributions), - the ranking by the ChiSquared p-values (for discrete distributions), - the ranking by BIC values. .. GENERATED FROM PYTHON SOURCE LINES 16-22 .. code-block:: default from __future__ import print_function import openturns as ot import openturns.viewer as viewer from matplotlib import pylab as plt ot.Log.Show(ot.Log.NONE) .. GENERATED FROM PYTHON SOURCE LINES 23-24 Create a sample from a continuous distribution .. GENERATED FROM PYTHON SOURCE LINES 24-29 .. code-block:: default distribution = ot.Beta(2.0, 2.0, 0.0, 1.) sample = distribution.getSample(1000) graph = ot.UserDefined(sample).drawCDF() view = viewer.View(graph) .. image-sg:: /auto_data_analysis/statistical_tests/images/sphx_glr_plot_fitted_distribution_ranking_001.png :alt: X0 CDF :srcset: /auto_data_analysis/statistical_tests/images/sphx_glr_plot_fitted_distribution_ranking_001.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 30-31 **1. Specify the model only** .. GENERATED FROM PYTHON SOURCE LINES 33-34 Create the list of distribution estimators .. GENERATED FROM PYTHON SOURCE LINES 34-36 .. code-block:: default factories = [ot.BetaFactory(), ot.TriangularFactory()] .. GENERATED FROM PYTHON SOURCE LINES 37-38 Rank the continuous models by the Lilliefors p-values: .. GENERATED FROM PYTHON SOURCE LINES 38-42 .. code-block:: default estimated_distribution, test_result = ot.FittingTest.BestModelLilliefors( sample, factories) test_result .. raw:: html

class=TestResult name=Unnamed type=Lilliefors Beta binaryQualityMeasure=false p-value threshold=0.5 p-value=0.006 statistic=0.0327766 description=[Beta(alpha = 1.72649, beta = 1.66568, a = 0.00526109, b = 0.970313) vs sample Beta]



.. GENERATED FROM PYTHON SOURCE LINES 43-44 Rank the continuous models wrt the BIC criteria (no test result): .. GENERATED FROM PYTHON SOURCE LINES 44-46 .. code-block:: default ot.FittingTest.BestModelBIC(sample, factories) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none [class=Beta name=Beta dimension=1 alpha=1.72649 beta=1.66568 a=0.00526109 b=0.970313, -0.19254944819710879] .. GENERATED FROM PYTHON SOURCE LINES 47-48 Rank the continuous models wrt the AIC criteria (no test result) .. GENERATED FROM PYTHON SOURCE LINES 48-50 .. code-block:: default ot.FittingTest.BestModelAIC(sample, factories) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none [class=Beta name=Beta dimension=1 alpha=1.72649 beta=1.66568 a=0.00526109 b=0.970313, -0.21218046931303733] .. GENERATED FROM PYTHON SOURCE LINES 51-52 Rank the continuous models wrt the AICc criteria (no test result): .. GENERATED FROM PYTHON SOURCE LINES 52-54 .. code-block:: default ot.FittingTest.BestModelAICC(sample, factories) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none [class=Beta name=Beta dimension=1 alpha=1.72649 beta=1.66568 a=0.00526109 b=0.970313, -0.2121402683080122] .. GENERATED FROM PYTHON SOURCE LINES 55-56 **2. Specify the model and its parameters** .. GENERATED FROM PYTHON SOURCE LINES 58-59 Create a collection of the distributions to be tested .. GENERATED FROM PYTHON SOURCE LINES 59-61 .. code-block:: default distributions = [ot.Beta(2.0, 2.0, 0.0, 1.0), ot.Triangular(0.0, 0.5, 1.0)] .. GENERATED FROM PYTHON SOURCE LINES 62-63 Rank the continuous models by the Kolmogorov p-values: .. GENERATED FROM PYTHON SOURCE LINES 63-67 .. code-block:: default estimated_distribution, test_result = ot.FittingTest.BestModelKolmogorov( sample, distributions) test_result .. raw:: html

class=TestResult name=Unnamed type=Kolmogorov Beta binaryQualityMeasure=true p-value threshold=0.05 p-value=0.127302 statistic=0.0369407 description=[Beta(alpha = 2, beta = 2, a = 0, b = 1) vs sample Beta]



.. GENERATED FROM PYTHON SOURCE LINES 68-69 Rank the continuous models wrt the BIC criteria: .. GENERATED FROM PYTHON SOURCE LINES 69-71 .. code-block:: default ot.FittingTest.BestModelBIC(sample, distributions) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none [class=Beta name=Beta dimension=1 alpha=2 beta=2 a=0 b=1, -0.21804827501286062] .. GENERATED FROM PYTHON SOURCE LINES 72-73 Rank the continuous models wrt the AIC criteria: .. GENERATED FROM PYTHON SOURCE LINES 73-75 .. code-block:: default ot.FittingTest.BestModelAIC(sample, distributions) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none [class=Beta name=Beta dimension=1 alpha=2 beta=2 a=0 b=1, -0.21804827501286062] .. GENERATED FROM PYTHON SOURCE LINES 76-77 Rank the continuous models wrt the AICc criteria: .. GENERATED FROM PYTHON SOURCE LINES 77-79 .. code-block:: default ot.FittingTest.BestModelAICC(sample, distributions) .. rst-class:: sphx-glr-script-out Out: .. code-block:: none [class=Beta name=Beta dimension=1 alpha=2 beta=2 a=0 b=1, -0.21804827501286062] .. GENERATED FROM PYTHON SOURCE LINES 80-81 **Discrete distributions** .. GENERATED FROM PYTHON SOURCE LINES 83-84 Create a sample from a discrete distribution .. GENERATED FROM PYTHON SOURCE LINES 84-89 .. code-block:: default distribution = ot.Poisson(2.0) sample = distribution.getSample(1000) graph = ot.UserDefined(sample).drawCDF() view = viewer.View(graph) .. image-sg:: /auto_data_analysis/statistical_tests/images/sphx_glr_plot_fitted_distribution_ranking_002.png :alt: X0 CDF :srcset: /auto_data_analysis/statistical_tests/images/sphx_glr_plot_fitted_distribution_ranking_002.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 90-91 Create the list of distribution estimators .. GENERATED FROM PYTHON SOURCE LINES 91-93 .. code-block:: default distributions = [ot.Poisson(2.0), ot.Geometric(0.1)] .. GENERATED FROM PYTHON SOURCE LINES 94-95 Rank the discrete models wrt the ChiSquared p-values: .. GENERATED FROM PYTHON SOURCE LINES 95-99 .. code-block:: default estimated_distribution, test_result = ot.FittingTest.BestModelChiSquared( sample, distributions) test_result .. raw:: html

class=TestResult name=Unnamed type=ChiSquared Poisson binaryQualityMeasure=true p-value threshold=0.05 p-value=0.184085 statistic=8.81784 description=[Poisson(lambda = 2) vs sample Poisson]



.. GENERATED FROM PYTHON SOURCE LINES 100-101 Rank the discrete models wrt the BIC criteria: .. GENERATED FROM PYTHON SOURCE LINES 101-103 .. code-block:: default ot.FittingTest.BestModelBIC(sample, distributions) plt.show() .. rst-class:: sphx-glr-timing **Total running time of the script:** ( 0 minutes 0.797 seconds) .. _sphx_glr_download_auto_data_analysis_statistical_tests_plot_fitted_distribution_ranking.py: .. only :: html .. container:: sphx-glr-footer :class: sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_fitted_distribution_ranking.py ` .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_fitted_distribution_ranking.ipynb ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_