18b. Ordinal probit model

Example of an ordinal probit model. This is just to illustrate the syntax, as the data are not ordered. But the example assume, for the sake of it, that the alternatives are ordered as 1->2->3

Michel Bierlaire, EPFL Thu Jun 26 2025, 15:54:37

from IPython.core.display_functions import display

import biogeme.biogeme_logging as blog
from biogeme.biogeme import BIOGEME
from biogeme.expressions import Beta, OrderedLogProbit
from biogeme.results_processing import (
    EstimationResults,
    get_pandas_estimated_parameters,
)

See the data processing script: Data preparation for Swissmetro.

from swissmetro_data import CHOICE, TRAIN_COST_SCALED, TRAIN_TT_SCALED, database

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

Example b18b_ordinal_probit.py

Parameters to be estimated

b_time = Beta('b_time', 0, None, None, 0)
b_cost = Beta('b_cost', 0, None, None, 0)

Threshold parameters for the ordered probit.

\(\tau_1 \leq 0\).

tau1 = Beta('tau1', -1, None, 0, 0)

\(\delta_2 \geq 0\).

delta2 = Beta('delta2', 2, 0, None, 0)

\(\tau_2 = \tau_1 + \delta_2\)

tau2 = tau1 + delta2

Utility

utility = b_time * TRAIN_TT_SCALED + b_cost * TRAIN_COST_SCALED

Associate each discrete indicator with an interval.

log_probability = OrderedLogProbit(
    eta=utility,
    cutpoints=[tau1, tau2],
    y=CHOICE,
    categories=[1, 2, 3],
    neutral_labels=[],
)

Create the Biogeme object.

the_biogeme = BIOGEME(database, log_probability)
the_biogeme.model_name = 'b18b_ordinal_probit'

Biogeme parameters read from biogeme.toml.

Estimate the parameters.

try:
    results = EstimationResults.from_yaml_file(
        filename=f'saved_results/{the_biogeme.model_name}.yaml'
    )
except FileNotFoundError:
    results = the_biogeme.estimate()

print(results.short_summary())

Results for model b18b_ordinal_probit
Nbr of parameters:              4
Sample size:                    6768
Excluded data:                  3960
Final log likelihood:           -5789.055
Akaike Information Criterion:   11586.11
Bayesian Information Criterion: 11613.39

pandas_results = get_pandas_estimated_parameters(estimation_results=results)
display(pandas_results)

{'Estimated parameters':      Name     Value  Robust std err.  Robust t-stat.  Robust p-value
0  b_time  0.018053         0.023389        0.771831        0.440215
1  b_cost  0.687183         0.036818       18.664467        0.000000
2    tau1 -0.604797         0.038571      -15.680045        0.000000
3  delta2  1.754346         0.029600       59.269266        0.000000}