.. 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-21 .. code-block:: default 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 22-23 Create a sample from a continuous distribution .. GENERATED FROM PYTHON SOURCE LINES 23-28 .. 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 29-30 **1. Specify the model only** .. GENERATED FROM PYTHON SOURCE LINES 32-33 Create the list of distribution estimators .. GENERATED FROM PYTHON SOURCE LINES 33-35 .. code-block:: default factories = [ot.BetaFactory(), ot.TriangularFactory()] .. GENERATED FROM PYTHON SOURCE LINES 36-37 Rank the continuous models by the Lilliefors p-values: .. GENERATED FROM PYTHON SOURCE LINES 37-41 .. 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 42-43 Rank the continuous models wrt the BIC criteria (no test result): .. GENERATED FROM PYTHON SOURCE LINES 43-45 .. 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 46-47 Rank the continuous models wrt the AIC criteria (no test result) .. GENERATED FROM PYTHON SOURCE LINES 47-49 .. 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 50-51 Rank the continuous models wrt the AICc criteria (no test result): .. GENERATED FROM PYTHON SOURCE LINES 51-53 .. 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 54-55 **2. Specify the model and its parameters** .. GENERATED FROM PYTHON SOURCE LINES 57-58 Create a collection of the distributions to be tested .. GENERATED FROM PYTHON SOURCE LINES 58-60 .. 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 61-62 Rank the continuous models by the Kolmogorov p-values: .. GENERATED FROM PYTHON SOURCE LINES 62-66 .. 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 67-68 Rank the continuous models wrt the BIC criteria: .. GENERATED FROM PYTHON SOURCE LINES 68-70 .. 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 71-72 Rank the continuous models wrt the AIC criteria: .. GENERATED FROM PYTHON SOURCE LINES 72-74 .. 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 75-76 Rank the continuous models wrt the AICc criteria: .. GENERATED FROM PYTHON SOURCE LINES 76-78 .. 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 79-80 **Discrete distributions** .. GENERATED FROM PYTHON SOURCE LINES 82-83 Create a sample from a discrete distribution .. GENERATED FROM PYTHON SOURCE LINES 83-88 .. 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 89-90 Create the list of distribution estimators .. GENERATED FROM PYTHON SOURCE LINES 90-92 .. code-block:: default distributions = [ot.Poisson(2.0), ot.Geometric(0.1)] .. GENERATED FROM PYTHON SOURCE LINES 93-94 Rank the discrete models wrt the ChiSquared p-values: .. GENERATED FROM PYTHON SOURCE LINES 94-98 .. 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 99-100 Rank the discrete models wrt the BIC criteria: .. GENERATED FROM PYTHON SOURCE LINES 100-102 .. code-block:: default ot.FittingTest.BestModelBIC(sample, distributions) plt.show() .. rst-class:: sphx-glr-timing **Total running time of the script:** ( 0 minutes 0.509 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 `_