This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

A Study of an Atomic Mobility Game With Uncertainty Under Prospect Theory

Ioannis Vasileios Chremos1, Heeseung Bang1,2, Aditya Dave2, Viet-Anh Le1,3, Student Members, IEEE
Andreas A. Malikopoulos2, Senior Member, IEEE
This work was supported in part by NSF under Grants CNS-2149520 and CMMI-2219761, the Delaware Department of Transportation, and by the University of Delaware Graduate College through the University Doctoral Fellowship Award.1Department of Mechanical Engineering, University of Delaware, Newark, DE 19716 USA. ichremos@udel.edu2School of Civil and Environmental Engineering, Cornell University, Ithaca, NY 14850 USA. {hb489,a.dave,amaliko}@cornell.edu3Systems Engineering, Cornell University, Ithaca, NY 14850 USA. vl299@cornell.edu
Abstract

In this paper, we present a study of a mobility game with uncertainty in the decision-making of travelers and incorporate prospect theory to model travel behavior. We formulate a mobility game that models how travelers distribute their traffic flows in a transportation network with splittable traffic, utilizing the Bureau of Public Roads function to establish the relationship between traffic flow and travel time cost. Given the inherent non-linearities and complexity introduced by the uncertainties, we propose a smooth approximation function to estimate the prospect-theoretic cost functions. As part of our analysis, we characterize the best-fit parameters and derive an upper bound for the error. We then show the existence of an equilibrium and its its best-possible approximation.

I INTRODUCTION

Emerging mobility systems (e.g., connected and automated vehicles (CAVs), shared mobility) provide the most intriguing opportunity for enabling users to monitor transportation network conditions better and make efficient decisions for improving safety and transportation efficiency. The data and shared information of emerging mobility systems are associated with a new level of complexity in modeling and control [1]. The impact of selfish or irrational social behavior in routing networks of cars has been studied in recent years [2, 3, 4]. Other efforts have addressed how people learn and make routing decisions with behavioral dynamics [5]. The problem of how travelers often have to make decisions under the uncertainty of experiencing delays, especially when uncertainties directly affect travel time in a transportation network, has not been adequately approached yet. Hence, our problem of interest is to study in a game-theoretic setting these interactions and analyze the equilibrium of the travelers’ decisions under uncertainties [6]. We study the interactions of a finite group of players that seek to travel in a transportation network (with a unique origin-destination pair) comprised of roads with splittable traffic. A key characteristic of our approach is that we incorporate prospect theory, a behavioral model that captures the perceptions of utility under uncertainty (how likely and how much).

Some of the existing game-theoretical literature in control and transportation theory assumes that the players’ behavior follows the rational choice theory, i.e., each player is a risk-neutral and utility maximizer. This makes transportation models unrealistic, as unexpected travel delays can lead to uncertainty in a traveler’s utility. There is strong evidence from empirical experiments that show how humans’ choices and preferences systematically may deviate from the choice and preferences of a game-theoretic player under the rational choice theory [7]. For example, humans compare the outcomes of their choices to a known expected amount of utility (called reference) and make their final decision, using that reference to assess their losses or gains asymmetrically. Prospect theory has laid down the theoretical foundations to study such biases and the subjective perception of risk in the utility of humans [7, 8]. This theory has been recognized as a closer-to-reality behavioral model for the decision-making of humans in different engineering problems [9, 10, 11].

In general, one of the standard approaches to alleviate congestion in a transportation system has been managing the travel demand and supply while also taking into consideration the scarce resources. Such approaches focus primarily on traffic routing, which aims to optimize the routing decisions in a transportation network [12]. Another approach is game theory that allows us to investigate the impact of selfish routing on efficiency and congestion [13] and assign travelers routes to minimize travel time under a Nash Equilibrium (NE) [14, 15, 16, 17, 18]. A fundamental theoretical approach in alleviating congestion is routing/congestion games [19, 20], which are a generalization of the standard resource-sharing game of an arbitrary number of resources in a network.

In this paper, we use Prelec’s probability weighting function and an S-shaped value function to model how travelers perceive traffic uncertainties and their travel gains/losses. So, our first contribution is incorporating prospect theory into an atomic routing game with splittable traffic to capture a realistic version of the travelers’ decision-making regarding travel time costs. The S-shaped value function is adopted to represent the curvature of the travel time cost function and account for the travelers’ perception of gains/losses in travel time according to a reference point (defined using the US Bureau of Public Roads function). To address prospect theory’s mathematical intractabilities, our second contribution proposes a smooth approximation function that estimates the non-linear piecewise prospect-theoretic cost functions. Thus, we can estimate how travelers perceive gains/losses and probabilities in travel time costs. This work is focused on establishing the fitness of the approximation function, proving the existence of at least one NE in pure strategies.

The remainder of the paper is structured as follows. In Section II, we present the mathematical formulation of the proposed game-theoretic framework. In Section III, we derive the theoretical properties of the proposed framework, and finally, draw conclusions in Section IV.

II Modeling Framework

We consider a routing game with a finite non-empty set of players \mathcal{I}, ||=n|\mathcal{I}|=n\in\mathbb{N}. Each player ii may represents a class of travelers that could use connected and automated vehicles (CAVs) and who control a significant amount of traffic, say xi0x_{i}\in\mathbb{R}_{\geq 0}. Thus, we interpret xix_{i} as the representation of the flow of traffic that player ii contributes to a transportation network. We define traffic flow in this setting as the number of vehicles passing through each point in the network over time. This decision variable is non-negative as players (or the travelers) make trips using their vehicle over time in the transportation network. This is in contrast to non-atomic routing games, where players only control an infinitesimal amount of traffic. We also assume that traffic is splittable. Travelers seek to travel in a transportation network represented by a directed multigraph 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}), where each node in 𝒱\mathcal{V} may represent different city areas or neighborhoods (e.g., Braess’ paradox network). Each edge ee\in\mathcal{E} may represent a road. For our purposes, we think of 𝒢=(𝒱,)\mathcal{G}=(\mathcal{V},\mathcal{E}) as a representation of a smart city network with a road infrastructure. Any player ii\in\mathcal{I} seeks to travel from an origin o𝒱o\in\mathcal{V} to a destination d𝒱d\in\mathcal{V}. So, all players are associated with the same unique origin-destination pair (o,d)𝒱×𝒱(o,d)\in\mathcal{V}\times\mathcal{V}. Next, each player may use a sequence of edges that connects the OD pair (o,d)(o,d). We define 2\mathcal{R}\subset 2^{\mathcal{E}} as the set of routes available to any player ii\in\mathcal{I}, where their route rir_{i} consists of a sequence of edges connecting the origin-destination pair (o,d)(o,d). We are interested in how such players may compete over the routes in the network for routing their traffic flows (this is a multiple-route traffic flow decision-making problem).

Since each player ii\in\mathcal{I} seeks to route their traffic represented by flow xix_{i} in the network 𝒢\mathcal{G}, we define, for each ii\in\mathcal{I}, the set of actions as

𝒳i={xi0||:rixiri=x¯i},\mathcal{X}_{i}=\left\{x_{i}\in\mathbb{R}_{\geq 0}^{|\mathcal{R}|}:\sum_{r_{i}\in\mathcal{R}}x_{i}^{r_{i}}=\bar{x}_{i}\right\}, (1)

where xi=(xiri1,xiri2,,xiri||)x_{i}=(x_{i}^{r_{i}^{1}},x_{i}^{r_{i}^{2}},\dots,x_{i}^{r_{i}^{|\mathcal{R}|}}), x¯i0\bar{x}_{i}\in\mathbb{R}_{\geq 0} is the total flow of player ii, and rikr_{i}^{k} denotes the kk-th route in the network.

Remark 1.

Note here that each player ii controls their traffic flow xix_{i}, which we represent as a vector since player ii may choose to use different routes in the transportation network, thus sending traffic xirikx_{i}^{r_{i}^{k}} for some kk. The total traffic flow controlled by player ii is finite, though. And so, we represent this by introducing x¯i\bar{x}_{i}.

We write 𝒳=𝒳1×𝒳2××𝒳n\mathcal{X}=\mathcal{X}_{1}\times\mathcal{X}_{2}\times\dots\times\mathcal{X}_{n} for the Cartesian product of all the players’ action sets. We also write xi=(x1,x2,,xi1,xi+1,,xn)x_{-i}=(x_{1},x_{2},\dots,x_{i-1},x_{i+1},\dots,x_{n}) for the action profile that excludes player ii\in\mathcal{I}. Next, for the aggregate action profile, we write x=(xi,xi)x=(x_{i},x_{-i}), x𝒳x\in\mathcal{X}.

Definition 1.

The flow on edge ee\in\mathcal{E} is the sum of relevant components of all players’ traffic flows that have chosen a route that includes edge ee, i.e., fe(x)=iriexirif_{e}(x)=\sum_{i\in\mathcal{I}}\sum_{r_{i}\ni e}x_{i}^{r_{i}}.

In our routing game where each player ii\in\mathcal{I} chooses their traffic flow vector xix_{i} over a common set of routes \mathcal{R}, if player ii chooses to send traffic xirix_{i}^{r_{i}} along route rir_{i}, then this traffic will be distributed along all the edges in this route rir_{i}. This is because a traveler’s traffic on some route rir_{i} is a single quantity among all the route’s edges.

Next, we introduce a travel time latency function to capture the cost that players may experience. Intuitively, we capture the players’ preferences for different outcomes using a “cost function,” in which players are expected to act as cost minimizers. For each ee\in\mathcal{E}, we consider non-negative cost functions ce:00c_{e}:\mathbb{R}_{\geq 0}\to\mathbb{R}_{\geq 0}. We assume that the cost functions at each edge ee are convex, continuous, and differentiable with respect to fef_{e}. One standard way to define in an exact form cec_{e} is by the US Bureau of Public Roads (BPR) function, as it is a commonly used model for the relationship between flow and travel time. Mathematically, we have, for any edge ee\in\mathcal{E},

ce(fe)=ce0(1+320(fefeCRT)4),c_{e}(f_{e})=c_{e}^{0}\left(1+\frac{3}{20}\left(\frac{f_{e}}{f_{e}^{CRT}}\right)^{4}\right), (2)

where ce0c_{e}^{0} is the free-flow travel time and feCRTf_{e}^{CRT} is the critical capacity of traffic flow on road ee. Note that the BPR function is non-linear, continuous, differentiable, strictly increasing, and strictly convex for fe0f_{e}\geq 0.

Definition 2.

If the maximum flow on edge ee\in\mathcal{E} is femax>0f_{e}^{\max}\in\mathbb{R}_{>0}, then for the critical flow, feCRTf_{e}^{CRT}, on edge ee\in\mathcal{E} we have feCRT<femaxf_{e}^{CRT}<f_{e}^{\max}.

Next, for some route rir_{i} of any player ii, its cost is the sum of the costs on the edges that constitute route rir_{i}, i.e., cri(x)=erice(fe(x))c_{r_{i}}(x)=\sum_{e\in r_{i}}c_{e}(f_{e}(x)). The total cost for player ii is

ci(x)=ricri(x)=ri[erice(fe)],c_{i}(x)=\sum_{r_{i}\in\mathcal{R}}c_{r_{i}}(x)=\sum_{r_{i}\in\mathcal{R}}\left[\sum_{e\in r_{i}}c_{e}(f_{e})\right], (3)

where fe(x)=iriexirif_{e}(x)=\sum_{i\in\mathcal{I}}\sum_{r_{i}\ni e}x_{i}^{r_{i}}.

The game is fully characterized by the tuple =,(𝒳i)i,(ci)i\mathcal{M}=\langle\mathcal{I},(\mathcal{X}_{i})_{i\in\mathcal{I}},(c_{i})_{i\in\mathcal{I}}\rangle. This non-cooperative routing game is a simultaneous-move game where players make decisions simultaneously and commute in (o,d)(o,d) of network 𝒢\mathcal{G}. Players behave selfishly and aim to minimize their costs (e.g., travel time latencies). Naturally, players compete with each other over the available yet limited routes and how to utilize them in the transportation network. Indirectly, players make route choices that satisfy their travel needs (modeled through traffic flow). Next, we clarify “who knows what?” in \mathcal{M}. All players have complete knowledge of the game and the network. Each player knows their own information (action and cost) as well as the information of other players. At equilibrium, we want to ensure that no player has an incentive to unilaterally deviate from their chosen decisions and change how they distribute their traffic flows over the available routes in the network. So, for our purposes, we observe that an NE in terms of the players’ traffic flows in pure strategies is the most appropriate a solution concept for our game.

Definition 3.

A feasible flow profile x=(xiri)iri𝒳x^{*}=(x_{i}^{r_{i}})_{i\in\mathcal{I}}^{r_{i}\in\mathcal{R}}\in\mathcal{X} constitutes a NE if for each player ii\in\mathcal{I}, ci(xi,xi)ci(xi,xi)c_{i}(x_{i}^{*},x_{-i}^{*})\leq c_{i}(x_{i},x_{-i}^{*}), for all xi𝒳ix_{i}\in\mathcal{X}_{i}.

In other words, a flow profile xx^{*} is an NE if no player can reduce their total cost by unilaterally changing how they distribute their total traffic flow over the available routes in the network. In an NE, each player’s specific xix_{i} has the lowest possible cost among all possible distributions over the routes, given the choices made by other players.

II-A Prospect Theory Analysis

In this subsection, we briefly introduce prospect theory and its main concepts [21]. Prospect theory attempts to answer one of the main questions of how a decision-maker may evaluate different possible actions/outcomes under uncertain and risky circumstances. Thus, prospect theory is a descriptive behavioral model and focuses on three main behavioral factors: (i) Reference dependence: decision-makers make decisions based on their utility, which is measured from the “gains” or “losses.” However, the utility is a gain or loss relative to a reference point that may be unique to each decision-maker. (ii) Diminishing sensitivity: changes in value have a greater impact near the reference point than away from the reference point. (iii) Loss aversion: decision-makers are more conservative in gains and riskier in losses. One way to mathematize the above behavioral factors (1) - (3) is to consider an action by a decision-maker as a “gamble” with objective utility value zz\in\mathbb{R}. We say that this decision-maker perceives zz subjectively using a value function [7, 22]

v(z)={(zz0)β1,if zz0,λ(z0z)β2,if z<z0,v(z)=\begin{cases}(z-z_{0})^{\beta_{1}},&\text{if }z\geq z_{0},\\ -\lambda(z_{0}-z)^{\beta_{2}},&\text{if }z<z_{0},\end{cases} (4)

where z0z_{0} represents a reference point, β1,β2(0,1)\beta_{1},\beta_{2}\in(0,1) are parameters that represent the diminishing sensitivity. Both β1,β2\beta_{1},\beta_{2} shape (4) in a way that the changes in value have a greater impact near the reference point than away from the reference point. We observe that (4) is concave in the domain of gains and convex in the domain of losses. Moreover, λ1\lambda\geq 1 reflects the level of loss aversion of decision-makers. To the best of our knowledge, a widely agreed theory does not exist that determines and defines the reference dependence [7]. In engineering [23, 11], it is assumed that z0=0z_{0}=0 captures a decision-maker’s expected status-quo level of the resources.

Prospect theory models the subjective behavior of decision-makers under uncertainty and risk. Each objective utility zz\in\mathbb{R} is associated with a probabilistic occurrence, say p[0,1]p\in[0,1]. Decision-makers are subjective and perceive pp differently depending on its value. To capture this behavior, we introduce a strictly increasing function w:[0,1][0,1]w:[0,1]\to[0,1] with w(0)=0w(0)=0 and w(1)=1w(1)=1 called the probability weighting function. This function allows us to model how decision-makers may overestimate small probabilities of objective utilities, i.e., w(p)>pw(p)>p if pp is close to 0, or underestimate high probabilities, i.e., w(p)<pw(p)<p if pp is close to 11. We use Prelec’s probability weighting function first introduced in [24], w(p)=exp((log(p))β3),p[0,1]w(p)=\exp\left(-(-\log(p))^{\beta_{3}}\right),\quad p\in[0,1], where β3(0,1)\beta_{3}\in(0,1) represents a rational index, i.e., the distortion of a decision-maker’s probability perceptions. Mathematically, β3\beta_{3} controls the curvature of the weighting function.

Definition 4.

Suppose that there are KK\in\mathbb{N} possible outcomes available to a decision-maker and zkz_{k}\in\mathbb{R} is the kk-th gain/loss of objective utility. Then a prospect k\ell_{k} is a tuple of the utilities and their respective probabilities, i.e., k=(z0,z1,z2,,zK;p0,p1,p2,,pK)\ell_{k}=(z_{0},z_{1},z_{2},\dots,z_{K};p_{0},p_{1},p_{2},\dots,p_{K}), where k=0,1,2,,Kk=0,1,2,\dots,K. We denote the kk-th prospect more compactly as k=(zk,pk)\ell_{k}=(z_{k},p_{k}). We have that k=0Kpk=1\sum_{k=0}^{K}p_{k}=1 and k\ell_{k} is well-ordered, i.e., z0z1zKz_{0}\leq z_{1}\leq\cdots\leq z_{K}. Under prospect theory, the decision-maker evaluates their “subjective utility” as u()=0kKv(zk)w(pk)u(\ell)=\sum_{0\leq k\leq K}v(z_{k})w(p_{k}), where =(k)k=1K\ell=(\ell_{k})_{k=1}^{K} is the profile of prospects of KK outcomes.

In the remainder of this subsection, we apply the prospect theory to our modeling framework, clearly define the mobility outcomes (objective and subjective utilities), and then show that the prospect-theoretic game \mathcal{M} admits a NE.

Players may be uncertain about the value of the traffic disturbances as it is affected by unexpected factors, and so we use Prelec’s probability weighting function w:[0,1][0,1]w:[0,1]\to[0,1] to capture how different traveler populations “perceive” probabilities. In addition, we are interested in capturing how players may perceive their gains or losses regarding their travel time costs with respect to the costs at critical density. Hence, we define the mobility prospect as whether fef_{e} will reach its critical or jammed point. Formally, πe\pi_{e} is the probability that fe(0,feCRT)f_{e}\in(0,f_{e}^{CRT}), and 1πe1-\pi_{e} is the probability for fe(feCRT,femax]f_{e}\in(f_{e}^{CRT},f_{e}^{\max}]. We then use the prospect-theoretic S-shaped value function v(ce(fe)):0v(c_{e}(f_{e})):\mathbb{R}_{\geq 0}\to\mathbb{R} to capture how players may perceive such costs. Hence, we have

v(ce(fe))={λ(ce0ce(fe))β,if ce(fe)ce0,(ce(fe)ce0)β,if ce(fe)>ce0,v(c_{e}(f_{e}))=\begin{cases}\lambda(c_{e}^{0}-c_{e}(f_{e}))^{\beta},&\text{if }c_{e}(f_{e})\leq c_{e}^{0},\\ -(c_{e}(f_{e})-c_{e}^{0})^{\beta},&\text{if }c_{e}(f_{e})>c_{e}^{0},\end{cases} (5)

where the reference dependence is represented by ce0=ce(feCRT)c_{e}^{0}=c_{e}(f_{e}^{CRT}), β1=β2=β(0,1)\beta_{1}=\beta_{2}=\beta\in(0,1), and for each ee\in\mathcal{E}, we have πe[0,1]\pi_{e}\in[0,1]. We justify β1=β2\beta_{1}=\beta_{2} in above as it has been verified to produce extremely good results, and the outcomes are consistent with the original data [8]. We define

c~e(fe)={ce0ce(fe),if ce(fe)<ce0,ce(fe)ce0,if ce(fe)>ce0.\tilde{c}_{e}(f_{e})=\begin{cases}c_{e}^{0}-c_{e}(f_{e}),&\text{if }c_{e}(f_{e})<c_{e}^{0},\\ c_{e}(f_{e})-c_{e}^{0},&\text{if }c_{e}(f_{e})>c_{e}^{0}.\end{cases} (6)
Remark 2.

It is important to note that our prospect-theoretic value function is “reversed,” capturing the way a traveler will perceive the gains in travel time through a cost function. Using as a reference point the critical traffic flow on edge ee, we can pinpoint the exact point that any more delays become socially unacceptable, i.e., a higher flow causes a higher travel time that a traveler will not tolerate.

The new cost function is

cePT(fe)=w(πe)λ[c~e(fe|ce(fe)<ce0)]βw(1πe)[c~e(fe|ce(fe)>ce0)]β.c_{e}^{\text{PT}}(f_{e})=w(\pi_{e})\cdot\lambda\cdot[\tilde{c}_{e}(f_{e}\,|\,c_{e}(f_{e})<c_{e}^{0})]^{\beta}\\ -w(1-\pi_{e})\cdot[\tilde{c}_{e}(f_{e}\,|\,c_{e}(f_{e})>c_{e}^{0})]^{\beta}. (7)

Observe that this particular formulation allows only two main outcomes for any player. One outcome may represent an easy commute (no traffic), and the other may represent traffic. For our purposes, we naturally expect two probabilities for these two outcomes. Future work will ensure to allow a larger distribution of probabilities for many different outcomes for any player (in such cases, cumulative prospect theory would be a more appropriate model [8]).

The total cost on some route rir_{i} for player ii under prospect theory is criPT(x)=ericePT(fe)c_{r_{i}}^{\text{PT}}(x)=\sum_{e\in r_{i}}c_{e}^{\text{PT}}(f_{e}). Now, the total cost of some player ii is given by

ciPT(x)=riericePT(fe).c_{i}^{\text{PT}}(x)=\sum_{r_{i}\in\mathcal{R}}\sum_{e\in r_{i}}c_{e}^{\text{PT}}(f_{e}). (8)

Note, however, that in this case, the prospect-theoretic cost is capturing the gains and losses of travel. Thus, the aim is to maximize this function to maximize the gains (by minimizing the actual cost of travel latencies).

What we observe in (8) is that it is rather cumbersome to analyze it analytically as issues in its smoothness arise quickly. The problem in analyzing such a function is that the exponent takes values in (0,1)(0,1). To address this theoretical obstacle, we propose a new function that approximates the prospect-theoretic function and, most importantly, can be shown to have useful properties. Hence, we define the following function

σ(fe)=δ11+exp(δ2feδ3)+δ4,\sigma(f_{e})=\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-f_{e}}{\delta_{3}}\right)}}+\delta_{4}, (9)

where δ1,δ2,δ3,δ4\delta_{1},\delta_{2},\delta_{3},\delta_{4}\in\mathbb{R}, and fe[0,femax]f_{e}\in[0,f_{e}^{\max}]. Hence, we can approximately evaluate (8) with the following:

ciPT(x)=rieriσ(fe).c_{i}^{\text{PT}}(x)=\sum_{r_{i}\in\mathcal{R}}\sum_{e\in r_{i}}\sigma(f_{e}). (10)

III Analysis and Properties of the Game

In this section, we provide a formal analysis of the properties of our proposed modeling framework, characterize the coefficients of σ\sigma function, and show that our game admits an NE in pure strategies.

Lemma 1.

The strategy space of the game \mathcal{M} is non-empty, compact, and convex.

Proof.

The proof has been omitted here due to space-constraints. ∎

Lemma 2.

The approximation function given by (9) in the interval [0,κ][0,\kappa], κ<femax\kappa<f_{e}^{\max}, is strictly concave with respect to fef_{e} when δ3>0\delta_{3}>0, δ4\delta_{4}\in\mathbb{R}, and (i) δ1>0\delta_{1}>0, δ2>fe\delta_{2}>f_{e}, or alternatively (ii) δ1<0\delta_{1}<0, δ2<fe\delta_{2}<f_{e}.

Proof.

Given that fe0f_{e}\geq 0, we analyze the second-order derivative of the function σ(fe)=δ11+exp(δ2feδ3)+δ4\sigma(f_{e})=\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-f_{e}}{\delta_{3}}\right)}}+\delta_{4} to determine the conditions for strict concavity. First, let us find the first and second-order derivatives of σ\sigma with respect to fef_{e}, i.e.,

σ(fe)\displaystyle\sigma^{\prime}(f_{e}) =δ1exp(δ2feδ3)δ3(1+exp(δ2feδ3))2,\displaystyle=\frac{-\delta_{1}\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})}}{\delta_{3}(1+\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})})^{2}}, (11)
σ′′(fe)\displaystyle\sigma^{\prime\prime}(f_{e}) =2δ1exp(δ2feδ3)(1exp(δ2feδ3))δ32(1+exp(δ2feδ3))3.\displaystyle=\frac{2\delta_{1}\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})}(1-\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})})}{\delta_{3}^{2}(1+\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})})^{3}}. (12)

Now, we examine the conditions for σ′′(fe)<0\sigma^{\prime\prime}(f_{e})<0. First, δ1\delta_{1} controls the sign of the second-order derivative as follows: if δ1<0\delta_{1}<0 and δ3>0\delta_{3}>0, σ′′(fe)\sigma^{\prime\prime}(f_{e}) will be negative when 1exp(δ2feδ3)<01-\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})}<0, which simplifies to δ2>fe\delta_{2}>f_{e}. If δ3<0\delta_{3}<0 in either of the cases, then the signs are reversed. We do require though that δ32\delta_{3}^{2} is well-defined, so δ30\delta_{3}\neq 0. On greater detail, 1exp(δ2feδ3)1-\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})} determines the conditions for σ′′(fe)\sigma^{\prime\prime}(f_{e}) to be negative. If δ1<0\delta_{1}<0, we need 1exp(δ2feδ3)>01-\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})}>0, which implies that fe>δ2δ3log(1)f_{e}>\delta_{2}-\delta_{3}\log{(1)} (since fe0f_{e}\geq 0). If δ1<0\delta_{1}<0, we need 1exp(δ2feδ3)<01-\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})}<0, which implies that fe<δ2δ3log(1)f_{e}<\delta_{2}-\delta_{3}\log{(1)}.

Combining these insights, we can conclude that the function σ(fe)=δ11+exp(δ2feδ3)+δ4\sigma(f_{e})=\frac{\delta_{1}}{1+\exp{(\frac{\delta_{2}-f_{e}}{\delta_{3}})}}+\delta_{4} becomes strictly concave in the entire interval. So, it is strictly concave for fe0f_{e}\geq 0 if: (i) δ1>0,δ3>0\delta_{1}>0,\delta_{3}>0 and fe<δ2f_{e}<\delta_{2}; (ii) δ1<0\delta_{1}<0, δ3>0\delta_{3}>0, and fe>δ2f_{e}>\delta_{2}. If δ3<0\delta_{3}<0, then the relation between fef_{e} and δ2\delta_{2} is naturally reversed. Note that the parameter δ4\delta_{4} does not affect the convexity of the function, as it only shifts the function vertically. Therefore, we have derived the necessary conditions that ensure σ′′(fe)\sigma^{\prime\prime}(f_{e}) is negative for all fef_{e}, making σ(fe)\sigma(f_{e}) strictly concave. ∎

It follows that it is strictly decreasing, continuous, and (continuously) differentiable with respect to the traffic flow fe[0,femax]f_{e}\in[0,f_{e}^{\max}] for any edge ee\in\mathcal{E}.

Now, we discuss the error characterization of our approximation function. Let us define the error function Φ\Phi as the squared difference between cePT(fe)c_{e}^{\text{PT}}(f_{e}) and σ(fe)\sigma(f_{e}), integrated over the interval [0,κ][0,\kappa]:

Φ(δ1,δ2,δ3,δ4)=0κ(cePT(fe)σ(fe))2dx.\Phi(\delta_{1},\delta_{2},\delta_{3},\delta_{4})=\int_{0}^{\kappa}\left(c_{e}^{\text{PT}}(f_{e})-\sigma(f_{e})\right)^{2}\mathrm{d}x. (13)

The goal is to minimize Φ\Phi with respect to the parameters δ1,δ2,δ3\delta_{1},\delta_{2},\delta_{3}, and δ4\delta_{4}. First, we find the critical points of Φ\Phi by setting its gradient to zero and solving the resulting system of equations: Φ(δ1,δ2,δ3,δ4)=0\nabla\Phi(\delta_{1},\delta_{2},\delta_{3},\delta_{4})=0. This results in a system of equations involving the partial derivatives of Φ\Phi with respect to each of the parameters, i.e., Φδ1=0,Φδ2=0,Φδ3=0\frac{\partial\Phi}{\partial\delta_{1}}=0,\quad\frac{\partial\Phi}{\partial\delta_{2}}=0,\frac{\partial\Phi}{\partial\delta_{3}}=0, and Φδ4=0\frac{\partial\Phi}{\partial\delta_{4}}=0. To compute these partial derivatives, we need to differentiate the integrand with respect to each parameter and then integrate it again, for example, Φδ1=0κδ1(cePT(fe)σ(fe))2dx\frac{\partial\Phi}{\partial\delta_{1}}=\int_{0}^{\kappa}\frac{\partial}{\partial\delta_{1}}\left(c_{e}^{\text{PT}}(f_{e})-\sigma(f_{e})\right)^{2}\mathrm{d}x. This process needs to be repeated for all parameters. However, due to the complexity of the function cePT(fe)c_{e}^{\text{PT}}(f_{e}) (being a non-linear piecewise function), it is not possible to obtain an explicit analytical expression for these partial derivatives. For our purposes, we rely on numerical optimization techniques to find the exact best-fit parameters that minimize the error function, as these methods can easily handle complex and non-linear optimization.

Theorem 1.

The error ϕ()=(cePT(fe)σ(fe))\phi(\cdot)=\left(c_{e}^{\text{PT}}(f_{e})-\sigma(f_{e})\right) is upper bounded by γ+ε\gamma+\varepsilon, where γ\gamma is some real number and ε>0\varepsilon>0.

Proof.

For the purposes of this proof we assume that β=0.5\beta=0.5, λ=2\lambda=2 and ce0=13c_{e}^{0}=13 and feCRT=1f_{e}^{\text{CRT}}=1 and ce0=14.95c_{e}^{0}=14.95. We substitute now the known equations to get ϕ(δ1,δ2,δ3,δ4)=w(1πe)[c~e(fe|ce(fe)>ce0)]βδ11+exp(δ2feδ3)δ4\phi(\delta_{1},\delta_{2},\delta_{3},\delta_{4})=-w(1-\pi_{e})\cdot[\tilde{c}_{e}(f_{e}\,|\,c_{e}(f_{e})>c_{e}^{0})]^{\beta}-\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-f_{e}}{\delta_{3}}\right)}}-\delta_{4}. Using a straightforward computation of the second-order derivative, we can get the inflection point of ϕ\phi, which will lie in (1,1+ε)(1,1+\varepsilon). This means that it is sufficient for us to compute ϕ\phi at fe=1f_{e}=1 and focus on ϕ\phi for fe>1f_{e}>1. Since σ\sigma is smooth and strictly concave in that interval, it approximates the worst cePTc_{e}^{\text{PT}} around the inflection point. So, we have the following ϕ(δ1,δ2,δ3,δ4)=w(1πe)(ce(fe)ce0)βδ11+exp(δ2feδ3)δ4\phi(\delta_{1},\delta_{2},\delta_{3},\delta_{4})=-w(1-\pi_{e})(c_{e}(f_{e})-c_{e}^{0})^{\beta}-\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-f_{e}}{\delta_{3}}\right)}}-\delta_{4}. This expression simplifies to

ϕ=w(1πe)[13(1+320(fe)4)29920]βδ11+exp(δ2feδ3)δ4,\phi=-w(1-\pi_{e})\left[13\Bigg{(}1+\frac{3}{20}(f_{e})^{4}\Bigg{)}-\frac{299}{20}\right]^{\beta}\\ -\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-f_{e}}{\delta_{3}}\right)}}-\delta_{4}, (14)

where we have δ1<0\delta_{1}<0 and δ2,δ3,δ4>0\delta_{2},\delta_{3},\delta_{4}>0, and δ2>fe\delta_{2}>f_{e}. Since w(1πe)w(1-\pi_{e}) is only a positive parameter constant, it is negligible, so we drop it from our analysis. The first component simplifies to [3920(1(fe)4)]β\left[\frac{39}{20}\left(1-(f_{e})^{4}\right)\right]^{\beta}, which is negative when we evaluate near the inflection point. Next, it follows that the second component is positive for small values of δ2\delta_{2} and δ3\delta_{3}. We use the Taylor series expansion evaluated at fe=1+εf_{e}=1+\varepsilon, where ε\varepsilon is a small positive number to get

[3920((fe)41)]β=395ε34395ε32732395ε52+O(ε72),-\left[\frac{39}{20}\left((f_{e})^{4}-1\right)\right]^{\beta}=\\ -\sqrt{\frac{39}{5}}\sqrt{\varepsilon}-\frac{3}{4}\frac{\sqrt{39}}{5}\varepsilon^{\frac{3}{2}}-\frac{7}{32}\frac{\sqrt{39}}{5}\varepsilon^{\frac{5}{2}}+O(\varepsilon^{\frac{7}{2}}), (15)

which is clearly negative. For the second component, we use the Taylor series expansion at fe=1+εf_{e}=1+\varepsilon to get

δ11+exp(δ2(1+ε)δ3)=δ11+exp(δ21δ3)+δ1exp(δ21δ3)δ3(1+exp(δ21δ3))2ε12δ1exp(1δ3+δ2δ3)(exp(1δ3)exp(δ2δ3))δ32(exp(1δ3)+exp(δ2δ3))3ε2+O(ε3).\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-(1+\varepsilon)}{\delta_{3}}\right)}}=\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}\\ +\frac{\delta_{1}\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}{\delta_{3}\left(1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}\right)^{2}}\cdot\varepsilon\\ \frac{1}{2}\frac{\delta_{1}\exp{\left(\frac{1}{\delta_{3}}+\frac{\delta_{2}}{\delta_{3}}\right)}\left(\exp{\left(\frac{1}{\delta_{3}}\right)}-\exp{\left(\frac{\delta_{2}}{\delta_{3}}\right)}\right)}{\delta_{3}^{2}\left(\exp{\left(\frac{1}{\delta_{3}}\right)}+\exp{\left(\frac{\delta_{2}}{\delta_{3}}\right)}\right)^{3}}\cdot\varepsilon^{2}+O(\varepsilon^{3}). (16)

We combine the expressions for the first and second components. Next, we have

ϕ=395ε34395ε32732395ε52+O(ε72)+δ11+exp(δ21δ3)+δ1exp(δ21δ3)δ3(1+exp(δ21δ3))2ε+12δ1exp(1δ3+δ2δ3)(exp(1δ3)exp(δ2δ3))δ32(exp(1δ3)+exp(δ2δ3))3ε2+O(ε3).\phi=-\sqrt{\frac{39}{5}}\sqrt{\varepsilon}-\frac{3}{4}\frac{\sqrt{39}}{5}\varepsilon^{\frac{3}{2}}-\frac{7}{32}\frac{\sqrt{39}}{5}\varepsilon^{\frac{5}{2}}+O(\varepsilon^{\frac{7}{2}})\\ +\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}+\frac{\delta_{1}\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}{\delta_{3}\left(1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}\right)^{2}}\cdot\varepsilon\\ +\frac{1}{2}\frac{\delta_{1}\exp{\left(\frac{1}{\delta_{3}}+\frac{\delta_{2}}{\delta_{3}}\right)}\left(\exp{\left(\frac{1}{\delta_{3}}\right)}-\exp{\left(\frac{\delta_{2}}{\delta_{3}}\right)}\right)}{\delta_{3}^{2}\left(\exp{\left(\frac{1}{\delta_{3}}\right)}+\exp{\left(\frac{\delta_{2}}{\delta_{3}}\right)}\right)^{3}}\cdot\varepsilon^{2}+O(\varepsilon^{3}). (17)

We want to find an upper bound for the error, which means we need to show that (17) is less than or equal to γ+ε\gamma+\varepsilon for some γ\gamma\in\mathbb{R}. Note that for any a,ba,b\in\mathbb{R} with a<0a<0 and b>0b>0, it is always true that a+bmax{a,b}a+b\leq\max\{a,b\}. So,

ϕmax{395ε34395ε32732395ε52+O(ε72),δ11+exp(δ21δ3)+δ1exp(δ21δ3)δ3(1+exp(δ21δ3))2ε+12δ1exp(1δ3+δ2δ3)(exp(1δ3)exp(δ2δ3))δ32(exp(1δ3)+exp(δ2δ3))3ε2+O(ε3)}.\phi\leq\max\bigg{\{}-\sqrt{\frac{39}{5}}\sqrt{\varepsilon}-\frac{3}{4}\frac{\sqrt{39}}{5}\varepsilon^{\frac{3}{2}}-\frac{7}{32}\frac{\sqrt{39}}{5}\varepsilon^{\frac{5}{2}}+O(\varepsilon^{\frac{7}{2}}),\\ \frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}+\frac{\delta_{1}\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}{\delta_{3}\left(1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}\right)^{2}}\cdot\varepsilon\\ +\frac{1}{2}\frac{\delta_{1}\exp{\left(\frac{1}{\delta_{3}}+\frac{\delta_{2}}{\delta_{3}}\right)}\left(\exp{\left(\frac{1}{\delta_{3}}\right)}-\exp{\left(\frac{\delta_{2}}{\delta_{3}}\right)}\right)}{\delta_{3}^{2}\left(\exp{\left(\frac{1}{\delta_{3}}\right)}+\exp{\left(\frac{\delta_{2}}{\delta_{3}}\right)}\right)^{3}}\cdot\varepsilon^{2}+O(\varepsilon^{3})\bigg{\}}. (18)

As ε\varepsilon is positively small, we take the limit as ε0\varepsilon\to 0. We note that the term 395ε-\sqrt{\frac{39}{5}}\sqrt{\varepsilon} dominates as ε0\varepsilon\to 0, and so the first component approaches -\infty as ε0\varepsilon\to 0. For the second component, the term δ11+exp(δ21δ3)\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}} dominates as ε0\varepsilon\to 0, and since δ1<0\delta_{1}<0 and δ2,δ3>0\delta_{2},\delta_{3}>0, the second component is positive. Hence, we can write

limε0ϕlimε0max{395ε,δ11+exp(δ21δ3)}.\lim_{\varepsilon\to 0}\phi\leq\lim_{\varepsilon\to 0}\max\left\{-\sqrt{\frac{39}{5}}\sqrt{\varepsilon},\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}\right\}. (19)

As ε0\varepsilon\to 0, we have 395ε-\sqrt{\frac{39}{5}}\sqrt{\varepsilon}\to-\infty, hence

limε0ϕ\displaystyle\lim_{\varepsilon\to 0}\phi limε0max{395ε,δ11+exp(δ21δ3)}\displaystyle\leq\lim_{\varepsilon\to 0}\max\left\{-\sqrt{\frac{39}{5}}\sqrt{\varepsilon},\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}\right\}
=δ11+exp(δ21δ3).\displaystyle=\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}. (20)

Now, let γ=δ11+exp(δ21δ3)\gamma=\frac{\delta_{1}}{1+\exp{\left(\frac{\delta_{2}-1}{\delta_{3}}\right)}}. Since the second component is positive, we have γ>0\gamma>0, thus ϕγ+ε\phi\leq\gamma+\varepsilon. Therefore, the error ϕ\phi is upper bounded by γ+ε\gamma+\varepsilon, where γ\gamma\in\mathbb{R} and ε>0\varepsilon>0. ∎

Theorem 2.

The game \mathcal{M} admits at least one NE.

Proof.

We formally prove the existence of an NE in the prospect-theoretic routing game using Brouwer’s fixed point theorem. Recall that for any player ii, ciPT(x)=rieriσ(fe(x))c_{i}^{\text{PT}}(x)=\sum_{r_{i}\in\mathcal{R}}\sum_{e\in r_{i}}\sigma(f_{e}(x)), where σ\sigma is our smooth and monotonic approximation function. We define the best-response correspondence for each player ii as: bi(xi)=argmaxxiciPT(x)b_{i}(x_{-i})=\arg\max_{x_{i}}c_{i}^{\text{PT}}(x). Smoothness in the approximation function σ\sigma implies that it is continuous and has continuous derivatives. This implies that we can estimate the utility function ciPT(x)c_{i}^{\text{PT}}(x) continuously with respect to the traffic vector xx. To show that the best-response correspondence bi(xi)b_{i}(x_{-i}) is continuous, we need the argmax\arg\max operator to be continuous. Since the set of maximizers is compact, which actually follows from the compactness of the strategy space by Lemma 1. By Lemma 2, we have that σ\sigma is concave on a specific interval [0,κ][0,\kappa]. This implies that we can estimate the utility function ciPT(x)c_{i}^{\text{PT}}(x) within the interval [0,κ][0,\kappa] pointwise in a strictly decreasing and strictly concave curve with respect to xix_{i} for any player ii\in\mathcal{I}. However, a strictly concave function has at most one unique maximum, which ensures the single-validness of the best-response correspondence bi(xi)b_{i}(x_{-i}). We now define the combined best-response correspondence B(x)=(b1(x1),b2(x2),,bn(xn))B(x)=(b_{1}(x_{-1}),b_{2}(x_{-2}),\dots,b_{n}(x_{-n})). Since each bi(xi)b_{i}(x_{-i}) is continuous, B(x)B(x) is also continuous, and thus it maps the strategy space to itself. Hence, now we can apply Brouwer’s fixed point theorem, which guarantees that there exists a fixed point x=B(x)x^{*}=B(x^{*}); the result then follows. ∎

IV CONCLUSIONS AND FUTURE WORK

In this paper, we presented our mobility game incorporating an atomic splittable routing game with prospect theory to study travel behavior in mobility systems. We modeled the overestimation/underestimation of probabilities using Prelec’s probability weighting function, and we considered the traffic uncertainties and travelers’ perception of gains/losses in travel time using a prospect-theoretic S-shaped value function. We proposed an approximation function to address the non-linear and piecewise nature of the prospect-theoretic cost functions and showed that at least one NE exists. Finally, we also derived an upper bound for the error.

In future research, we can explore how to analyze a convex-concave piecewise non-linear optimization problem using optimization techniques, such as sequential convex programming or cutting plane methods. Developing such an optimization framework can enhance our ability to predict travel decisions in mobility systems under prospect theory. Another direction is to incorporate prospect theory and a taxation mechanism and study using artificial intelligence how we can incentivize prospect-theoretic travelers and the trade-offs of efficiency in the mobility systems [25, 26, 27].

References

  • [1] A. A. Malikopoulos, “A duality framework for stochastic optimal control of complex systems,” IEEE Transactions on Automatic Control, vol. 61, no. 10, pp. 2756–2765, 2016.
  • [2] N. Mehr and R. Horowitz, “How will the presence of autonomous vehicles affect the equilibrium state of traffic networks?” IEEE Transactions on Control of Network Systems, vol. 7(1), pp. 96–105, 2019.
  • [3] D. A. Lazar, S. Coogan, and R. Pedarsani, “Routing for traffic networks with mixed autonomy,” IEEE Transactions on Automatic Control, vol. 66(6), pp. 2664–2676, 2021.
  • [4] E. Bıyık, D. Lazar, R. Pedarsani, and D. Sadigh, “Incentivizing efficient equilibria in traffic networks with mixed autonomy,” IEEE Transactions on Control of Network Systems, vol. 8(4), pp. 1717–1729, 2021.
  • [5] W. Krichene, B. Drighès, and A. M. Bayen, “Online learning of nash equilibria in congestion games,” SIAM Journal on Control and Optimization, vol. 53(2), pp. 1056–1081, 2015.
  • [6] I. V. Chremos and A. A. Malikopoulos, “Mechanism design theory in control engineering: A tutorial and overview of applications in communication, power grid, transportation, and security systems,” arXiv preprint arXiv:2212.00756, 2022.
  • [7] D. Kahneman and A. Tversky, “Prospect theory: An analysis of decision under risk,” Econometrica, vol. 47(2), pp. 363–391, 1979.
  • [8] A. Tversky and D. Kahneman, “Advances in prospect theory: Cumulative representation of uncertainty,” Journal of Risk and Uncertainty, vol. 5(4), pp. 297–323, 1992.
  • [9] K. Nar, L. J. Ratliff, and S. Sastry, “Learning prospect theory value function and reference point of a sequential decision maker,” in IEEE 56th Annual Conference on Decision and Control (CDC), 2017, pp. 5770–5775.
  • [10] I. V. Chremos and A. A. Malikopoulos, “A traveler-centric mobility game: Efficiency and stability under rationality and prospect theory,” PLoS ONE, vol. 18(5):e0285322, 2023.
  • [11] S. R. Etesami, W. Saad, N. B. Mandayam, and H. V. Poor, “Smart routing of electric vehicles for load balancing in smart grids,” Automatica, vol. 120, no. 109148, 2020.
  • [12] A. Silva, H. Tembine, E. Altman, and M. Debbah, “Optimum and equilibrium in assignment problems with congestion: Mobile terminals association to base stations,” IEEE Transactions on Automatic Control, vol. 58(8), pp. 2018–2031, 2013.
  • [13] J. R. Marden and J. S. Shamma, “Game theory and distributed control,” in Handbook of Game Theory With Economic Applications, H. Peyton Young and S. Zamir, Eds.   Elsevier, 2015, vol. 4, ch. 16, pp. 861–899.
  • [14] P. N. Brown and J. R. Marden, “Optimal mechanisms for robust coordination in congestion games,” IEEE Transactions on Automatic Control, vol. 63(8), pp. 2437–2448, 2017.
  • [15] I. V. Chremos and A. A. Malikopoulos, “Socioeconomic impact of emerging mobility markets and implementation strategies,” in AI-enabled Technologies for Autonomous and Connected Vehicles, I. Kolmanovsky, Y. Murphey, and P. Watta, Eds.   Springer, 2023.
  • [16] ——, “An analytical study of a two-sided mobility game,” in 2022 American Control Conference (ACC), 2022, pp. 1254–1259.
  • [17] I. V. Chremos, L. E. Beaver, and A. A. Malikopoulos, “A game-theoretic analysis of the social impact of connected and automated vehicles,” in 2020 23rd International Conference on Intelligent Transportation Systems (ITSC).   IEEE, 2020, pp. 2214–2219.
  • [18] I. V. Chremos and A. A. Malikopoulos, “Design and stability analysis of a shared mobility market,” in 2021 European Control Conference (ECC), 2021, pp. 375–380.
  • [19] R. Rosenthal, “A class of games possessing pure-strategy nash equilibria,” International Journal of Game Theory, vol. 2, pp. 65–67, 1973.
  • [20] J. B. Rosen, “Existence and uniqueness of equilibrium points for concave n-person games,” Econometrica: Journal of the Econometric Society, vol. ), pp. 520–534, 1965.
  • [21] P. P. Wakker, Prospect Theory: For Risk and Ambiguity.   Cambridge University Press, 2010.
  • [22] A. Al-Nowaihi, I. Bradley, and S. Dhami, “A note on the utility function under prospect theory,” Economics Letters, vol. 99(2), pp. 337–339, 2008.
  • [23] A. R. Hota, S. Garg, and S. Sundaram, “Fragility of the commons under prospect-theoretic risk attitudes,” Games and Economic Behavior, vol. 98, pp. 135–164, 2016.
  • [24] D. Prelec, “The probability weighting function,” Econometrica, vol. 66(3), pp. 497–527, 1998.
  • [25] I. V. Chremos and A. A. Malikopoulos, “Social resource allocation in a mobility system with connected and automated vehicles: A mechanism design problem,” in Proceedings of the 59th IEEE Conference on Decision and Control (CDC), 2020, 2020, pp. 2642–2647.
  • [26] ——, “Mobility equity and economic sustainability using game theory,” in 2023 American Control Conference (ACC), 2023, pp. 1698–1703.
  • [27] J. Bharadiya, “Artificial intelligence in transportation systems a critical review,” American Journal of Computing and Engineering, vol. 6, no. 1, pp. 34–45, 2023.