xlogit 0.2.7
A Python package for GPU-accelerated estimation of mixed logit models.
License: MIT
Requires Python: >=3.5

===========================================================================================
xlogit: A Python Package for GPU-Accelerated Estimation of Mixed Logit Models.
===========================================================================================

.. image:: https://raw.githubusercontent.com/arteagac/xlogit/master/docs/xlogit_logo_1000px.png
  :width: 400

|Build| |Coverage| |Community| |Docs| |PyPi| |License|

.. _Mixed Logit: https://xlogit.readthedocs.io/en/latest/api/mixed_logit.html
.. _Multinomial Logit: https://xlogit.readthedocs.io/en/latest/api/multinomial_logit.html

`Examples <https://xlogit.readthedocs.io/en/latest/examples.html>`__ | `Docs <https://xlogit.readthedocs.io/en/latest/index.html>`__ | `Installation <https://xlogit.readthedocs.io/en/latest/install.html>`__ | `API Reference <https://xlogit.readthedocs.io/en/latest/api/index.html>`__ | `Contributing <https://xlogit.readthedocs.io/en/latest/contributing.html>`__ | `Contact <https://xlogit.readthedocs.io/en/latest/index.html#contact>`__ 

Quick start
===========
The following example uses ``xlogit`` to estimate a mixed logit model for choices of electricity supplier (`See the data here <https://github.com/arteagac/xlogit/blob/master/examples/data/electricity_long.csv>`__). The parameters are:

* ``X``: 2-D array of input data (in long format) with choice situations as rows, and variables as columns
* ``y``: 1-D array of choices (in long format)
* ``varnames``: List of variable names that matches the number and order of the columns in ``X``
* ``alts``:  1-D array of alternative indexes or an alternatives list
* ``ids``:  1-D array of the ids of the choice situations
* ``panels``: 1-D array of ids for panel formation
* ``randvars``: dictionary of variables and their mixing distributions (``"n"`` normal, ``"ln"`` lognormal, ``"t"`` triangular, ``"u"`` uniform, ``"tn"`` truncated normal)

The current version of `xlogit` only supports input data in long format.

.. code-block:: python

    # Read data from CSV file
    import pandas as pd
    df = pd.read_csv("examples/data/electricity_long.csv")
    
    # Fit the model with xlogit
    from xlogit import MixedLogit
    
    varnames = ['pf', 'cl', 'loc', 'wk', 'tod', 'seas']
    model = MixedLogit()
    model.fit(X=df[varnames],
              y=df['choice'],
              varnames=varnames,
              ids=df['chid'],
              panels=df['id'],
              alts=df['alt'],
              n_draws=600,
              randvars={'pf': 'n', 'cl': 'n', 'loc': 'n',
                        'wk': 'n', 'tod': 'n', 'seas': 'n'})
    model.summary()


::

    GPU processing enabled.
    Optimization terminated successfully.
             Current function value: 3888.413414
             Iterations: 46
             Function evaluations: 51
             Gradient evaluations: 51
    Estimation time= 2.6 seconds
    ----------------------------------------------------------------------
    Coefficient         Estimate      Std.Err.         z-val         P>|z|
    ----------------------------------------------------------------------
    pf                -0.9996286     0.0331488   -30.1557541     9.98e-100 ***
    cl                -0.2355334     0.0220401   -10.6865870      1.97e-22 ***
    loc                2.2307891     0.1164263    19.1605300      5.64e-56 ***
    wk                 1.6251657     0.0918755    17.6887855      6.85e-50 ***
    tod               -9.6067367     0.3112721   -30.8628296     2.36e-102 ***
    seas              -9.7892800     0.2913063   -33.6047603     2.81e-112 ***
    sd.pf              0.2357813     0.0181892    12.9627201      7.25e-31 ***
    sd.cl              0.4025377     0.0220183    18.2819903      2.43e-52 ***
    sd.loc             1.9262893     0.1187850    16.2166103      7.67e-44 ***
    sd.wk             -1.2192931     0.0944581   -12.9083017      1.17e-30 ***
    sd.tod             2.3354462     0.1741859    13.4077786      1.37e-32 ***
    sd.seas           -1.4200913     0.2095869    -6.7756668       3.1e-10 ***
    ----------------------------------------------------------------------
    Significance:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
    
    Log-Likelihood= -3888.413
    AIC= 7800.827
    BIC= 7847.493


For more examples of ``xlogit`` see `this Jupyter Notebook in Google Colab <https://colab.research.google.com/github/arteagac/xlogit/blob/master/examples/mixed_logit_model.ipynb>`__. **Google Colab provides GPU resources for free**, which will significantly speed up your model estimation using ``xlogit``.

Quick install
=============
Install ``xlogit`` using ``pip`` as follows:

.. code-block:: bash

   pip install xlogit


.. hint::

   To enable GPU processing, you must install the `CuPy Python library <https://docs.cupy.dev/en/stable/install.html>`__.  When ``xlogit`` detects that CuPy is properly installed, it switches to GPU processing without any additional setup. If you use Google Colab, CuPy is usually installed by default.


For additional installation details check xlogit installation instructions at: https://xlogit.readthedocs.io/en/latest/install.html


No GPU? No problem
==================
``xlogit`` can also be used without a GPU. However, if you need to speed up your model estimation, there are several low cost and even free options to access cloud GPU resources. For instance:

- `Google Colab <https://colab.research.google.com>`_ offers free GPU resources with no setup required, as the service can be accessed using a web browser. Using xlogit in Google Colab is very easy as it runs out of the box without having to to install CUDA or CuPy, which are installed by default. For examples of xlogit running in Google Colab `see this link <https://colab.research.google.com/github/arteagac/xlogit/blob/master/examples/mixed_logit_model.ipynb>`_.
- `Amazon Sagemaker Studio Lab <https://studiolab.sagemaker.aws/>`_ offers Python runtime environments with free GPUs.
- `Google Cloud platform <https://cloud.google.com/compute/gpus-pricing>`_ offers GPU processing at less than $1 USD per hour for NVIDIA Tesla K80 GPU with 4,992 CUDA cores.
- `Amazon Sagemaker <https://aws.amazon.com/ec2/instance-types/p2/>`_ offers virtual machine instances with the same TESLA K80 GPU at a similar price range of less than $1 USD per hour.

Benchmark
=========
As shown in the plots below, ``xlogit`` is significantly faster than existing estimation packages. Also, ``xlogit`` provides convenient scaling when the number of random draws increases. These results were obtained using a modest and low-cost NVIDIA GTX 1060 graphics card. More sophisticated graphics cards are expected to provide even faster estimation times. For additional details about this benchmark and for replication instructions check https://xlogit.readthedocs.io/en/latest/benchmark.html.

.. image:: https://raw.githubusercontent.com/arteagac/xlogit/master/examples/benchmark/results/time_benchmark_artificial.png
  :width: 300

.. image:: https://raw.githubusercontent.com/arteagac/xlogit/master/examples/benchmark/results/time_benchmark_apollo_biogeme.png
  :width: 300

Notes
=====
The current version allows estimation of:

- `Mixed Logit`_ with several types of mixing distributions (normal, lognormal, triangular, uniform, and truncated normal)
- `Mixed Logit`_ with panel data
- `Mixed Logit`_ with unbalanced panel data
- `Mixed Logit`_ with Halton draws
- `Multinomial Logit`_ models
- `Conditional logit <https://xlogit.readthedocs.io/en/latest/api/multinomial_logit.html>`_ models
- Handling of unbalanced availability of choice alternatives for all of the supported models 
- Post-estimation tools for prediction and specification testing
- Inclusion of sample weights for all of the supported models

Contributors
============
The following contributors have tremendously helped in the enhancement and expansion of `xlogit`'s features.  

- `@crforsythe <https://github.com/crforsythe>`__
- John Helveston (`@jhelvy  <https://github.com/jhelvy>`__)

Contact
=======

If you have any questions, ideas to improve ``xlogit``, or want to report a bug, `chat with us on gitter <https://gitter.im/xlogit/community>`__ or open a `new issue in xlogit's GitHub repository <https://github.com/arteagac/xlogit/issues>`__.

Citing ``xlogit``
=================
Please cite ``xlogit`` as follows:

    Arteaga, C., Park, J., Beeramoole, P. B., & Paz, A. (2022). xlogit: An open-source Python package for GPU-accelerated estimation of Mixed Logit models. Journal of Choice Modelling, 42, 100339. https://doi.org/10.1016/j.jocm.2021.100339
    
Or using BibTex as follows::

    @article{xlogit,
        title = {xlogit: An open-source Python package for GPU-accelerated estimation of Mixed Logit models},
        author = {Cristian Arteaga and JeeWoong Park and Prithvi Bhat Beeramoole and Alexander Paz},
        journal = {Journal of Choice Modelling},
        volume = {42},
        pages = {100339},
        year = {2022},
        issn = {1755-5345},
        doi = {https://doi.org/10.1016/j.jocm.2021.100339},
    }


.. |Build| image:: https://github.com/arteagac/xlogit/actions/workflows/python-tests.yml/badge.svg
   :target: https://github.com/arteagac/xlogit/actions/workflows/python-tests.yml

.. |Docs| image:: https://readthedocs.org/projects/xlogit/badge/?version=latest
   :target: https://xlogit.readthedocs.io/en/latest/?badge=latest
   :alt: Documentation Status

.. |Community| image:: https://badges.gitter.im/xlogit/community.svg
   :target: https://gitter.im/xlogit/community?utm_source=badge&utm_medium=badge&utm_campaign=pr-badge
   :alt: Community

.. |Coverage| image:: https://coveralls.io/repos/github/arteagac/xlogit/badge.svg?branch=master
   :target: https://coveralls.io/github/arteagac/xlogit?branch=master

.. |PyPi| image:: https://badge.fury.io/py/xlogit.svg
   :target: https://badge.fury.io/py/xlogit

.. |License| image:: https://img.shields.io/github/license/arteagac/xlogit
   :target: https://github.com/arteagac/xlogit/blob/master/LICENSE