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

.. only:: html

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

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

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

.. _sphx_glr_auto_examples_swissmetro_plot_b10nested_bottom_old.py:


Nested logit model normalized from bottom
=========================================

Example of a nested logit model where the normalization is done at the
 bottom level. The specification is using the original syntax for
 nests.  Since biogeme 3.13, a new syntax, more explicit, has been
 adopted.

:author: Michel Bierlaire, EPFL
:date: Sun Apr  9 18:05:04 2023

.. GENERATED FROM PYTHON SOURCE LINES 15-21

.. code-block:: default


    import biogeme.biogeme_logging as blog
    import biogeme.biogeme as bio
    from biogeme import models
    from biogeme.expressions import Beta








.. GENERATED FROM PYTHON SOURCE LINES 22-23

See the data processing script: :ref:`swissmetro_data`.

.. GENERATED FROM PYTHON SOURCE LINES 23-40

.. code-block:: default

    from swissmetro_data import (
        database,
        CHOICE,
        SM_AV,
        CAR_AV_SP,
        TRAIN_AV_SP,
        TRAIN_TT_SCALED,
        TRAIN_COST_SCALED,
        SM_TT_SCALED,
        SM_COST_SCALED,
        CAR_TT_SCALED,
        CAR_CO_SCALED,
    )

    logger = blog.get_screen_logger(level=blog.INFO)
    logger.info('Example b10nested_bottom.py')





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Example b10nested_bottom.py 




.. GENERATED FROM PYTHON SOURCE LINES 41-42

Parameters to be estimated

.. GENERATED FROM PYTHON SOURCE LINES 42-48

.. code-block:: default

    ASC_CAR = Beta('ASC_CAR', 0, None, None, 0)
    ASC_TRAIN = Beta('ASC_TRAIN', 0, None, None, 0)
    ASC_SM = Beta('ASC_SM', 0, None, None, 1)
    B_TIME = Beta('B_TIME', 0, None, None, 0)
    B_COST = Beta('B_COST', 0, None, None, 0)








.. GENERATED FROM PYTHON SOURCE LINES 49-55

This is the scale parameter of the choice model. It is usually
normalized to one. In this example, we normalize the nest parameter
instead, and estimate the scale parameter for the model.  If the
lower bound is set to zero, the model cannot be evaluated.
Therefore, we set the lower bound to a small number, strictly larger
than zero.

.. GENERATED FROM PYTHON SOURCE LINES 55-57

.. code-block:: default

    MU = Beta('MU', 0.5, 0.000001, 1.0, 0)








.. GENERATED FROM PYTHON SOURCE LINES 58-59

Definition of the utility functions.

.. GENERATED FROM PYTHON SOURCE LINES 59-63

.. code-block:: default

    V1 = ASC_TRAIN + B_TIME * TRAIN_TT_SCALED + B_COST * TRAIN_COST_SCALED
    V2 = ASC_SM + B_TIME * SM_TT_SCALED + B_COST * SM_COST_SCALED
    V3 = ASC_CAR + B_TIME * CAR_TT_SCALED + B_COST * CAR_CO_SCALED








.. GENERATED FROM PYTHON SOURCE LINES 64-65

Associate utility functions with the numbering of alternatives.

.. GENERATED FROM PYTHON SOURCE LINES 65-67

.. code-block:: default

    V = {1: V1, 2: V2, 3: V3}








.. GENERATED FROM PYTHON SOURCE LINES 68-69

Associate the availability conditions with the alternatives.

.. GENERATED FROM PYTHON SOURCE LINES 69-71

.. code-block:: default

    av = {1: TRAIN_AV_SP, 2: SM_AV, 3: CAR_AV_SP}








.. GENERATED FROM PYTHON SOURCE LINES 72-76

Definition of nests. In this example, we create a nest for the
existing modes, that is train (1) and car (3).  Each nest is
associated with a tuple containing (i) the nest parameter and (ii)
the list of alternatives.

.. GENERATED FROM PYTHON SOURCE LINES 76-80

.. code-block:: default

    existing = 1.0, [1, 3]
    future = 1.0, [2]
    nests = existing, future








.. GENERATED FROM PYTHON SOURCE LINES 81-85

Definition of the model. This is the contribution of each
observation to the log likelihood function.
The choice model is a nested logit, with availability conditions,
where the scale parameter mu is explicitly involved.

.. GENERATED FROM PYTHON SOURCE LINES 85-87

.. code-block:: default

    logprob = models.lognestedMevMu(V, av, nests, CHOICE, MU)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    It is recommended to define the nests of the nested logit model using the objects OneNestForNestedLogit and NestsForNestedLogit defined in biogeme.nests. 




.. GENERATED FROM PYTHON SOURCE LINES 88-89

Create the Biogeme object.

.. GENERATED FROM PYTHON SOURCE LINES 89-92

.. code-block:: default

    the_biogeme = bio.BIOGEME(database, logprob)
    the_biogeme.modelName = 'b10nested_bottom'





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    File biogeme.toml has been parsed. 




.. GENERATED FROM PYTHON SOURCE LINES 93-94

Estimate the parameters.

.. GENERATED FROM PYTHON SOURCE LINES 94-96

.. code-block:: default

    results = the_biogeme.estimate()





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    *** Initial values of the parameters are obtained from the file __b10nested_bottom.iter 
    Parameter values restored from __b10nested_bottom.iter 
    Optimization algorithm: hybrid Newton/BFGS with simple bounds [simple_bounds] 
    ** Optimization: Newton with trust region for simple bounds 
    Results saved in file b10nested_bottom~00.html 
    Results saved in file b10nested_bottom~00.pickle 




.. GENERATED FROM PYTHON SOURCE LINES 97-99

.. code-block:: default

    print(results.short_summary())





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Results for model b10nested_bottom
    Nbr of parameters:              5
    Sample size:                    6768
    Excluded data:                  3960
    Final log likelihood:           -5236.9
    Akaike Information Criterion:   10483.8
    Bayesian Information Criterion: 10517.9





.. GENERATED FROM PYTHON SOURCE LINES 100-102

.. code-block:: default

    pandas_results = results.getEstimatedParameters()
    pandas_results





.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <div>
    <style scoped>
        .dataframe tbody tr th:only-of-type {
            vertical-align: middle;
        }

        .dataframe tbody tr th {
            vertical-align: top;
        }

        .dataframe thead th {
            text-align: right;
        }
    </style>
    <table border="1" class="dataframe">
      <thead>
        <tr style="text-align: right;">
          <th></th>
          <th>Value</th>
          <th>Rob. Std err</th>
          <th>Rob. t-test</th>
          <th>Rob. p-value</th>
        </tr>
      </thead>
      <tbody>
        <tr>
          <th>ASC_CAR</th>
          <td>-0.343348</td>
          <td>0.118827</td>
          <td>-2.889482</td>
          <td>3.858768e-03</td>
        </tr>
        <tr>
          <th>ASC_TRAIN</th>
          <td>-1.051575</td>
          <td>0.164976</td>
          <td>-6.374109</td>
          <td>1.840299e-10</td>
        </tr>
        <tr>
          <th>B_COST</th>
          <td>-1.759647</td>
          <td>0.149310</td>
          <td>-11.785219</td>
          <td>0.000000e+00</td>
        </tr>
        <tr>
          <th>B_TIME</th>
          <td>-1.845915</td>
          <td>0.225653</td>
          <td>-8.180332</td>
          <td>2.220446e-16</td>
        </tr>
        <tr>
          <th>MU</th>
          <td>0.486839</td>
          <td>0.038918</td>
          <td>12.509269</td>
          <td>0.000000e+00</td>
        </tr>
      </tbody>
    </table>
    </div>
    </div>
    <br />
    <br />


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

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


.. _sphx_glr_download_auto_examples_swissmetro_plot_b10nested_bottom_old.py:

.. only:: html

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




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

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

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

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


.. only:: html

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

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