piecewise.py

Go to the documentation of this file.

from biogeme import *



# \param V A <a
# href="http://docs.python.org/py3k/tutorial/datastructures.html#dictionaries"
# target="_blank">dictionary</a> mapping each alternative id with the
# expression of the utility function.
# @param Gi A <a
# href="http://docs.python.org/py3k/tutorial/datastructures.html#dictionaries"
# target="_blank">dictionary</a> mapping each alternative id with the function
# \f[
#   \frac{\partial G}{\partial y_i}(e^{V_1},\ldots,e^{V_J})
#\f]
# where \f$G\f$ is the MEV generating function. If an alternative \f$i\f$ is not available, then \f$G_i = 0\f$.
# @param av A <a
# href="http://docs.python.org/py3k/tutorial/datastructures.html#dictionaries"
# target="_blank">dictionary</a> mapping each alternative id with its
# availability condition.
# @param choice Expression producing the id of the chosen alternative.
# @return Choice probability of the MEV model, given by
#  \f[
#    \frac{e^{V_i + \ln G_i(e^{V_1},\ldots,e^{V_J})}}{\sum_j e^{V_j + \ln G_j(e^{V_1},\ldots,e^{V_J})}}
#  \f]
#
# \code
# def mev(V,Gi,av,choice) :
#     H = {}
#     for i,v in V.items() :
#        H[i] =  Elem({0:0, 1: v + log(Gi[i])},Gi[i]!=0)  
#     P = bioLogit(H,av,choice)
#     return P
# \endcode

def piecewise(x,beta,slopes) :

    H = {}
    for i,v in V.items() :
        H[i] =  Elem({0:0, 1: v + log(Gi[i])},av[i]!=0)  
    P = bioLogit(H,av,choice)
    return P


def logmev(V,Gi,av,choice) :
    H = {}
    for i,v in V.items() :
        H[i] =  Elem({0:0, 1: v + log(Gi[i])},av[i]!=0)  
    logP = bioLogLogit(H,av,choice)
    return logP




def mev_selectionBias(V,Gi,av,correction,choice) :
    H = {}
    for i,v in V.items() :
        H[i] = v + log(Gi[i]) + correction[i]

    P = bioLogit(H,av,choice)
            
    return P




def logmev_selectionBias(V,Gi,av,correction,choice) :
    H = {}
    for i,v in V.items() :
        H[i] = v + log(Gi[i]) + correction[i]

    P = bioLogLogit(H,av,choice)
            
    return P