## DISCOVERING CAR-FOLLOWING DYNAMICS FROM TRAJECTORY DATA THROUGH DEEP LEARNING

## A PREPRINT

## Ohay Angah

Department of Civil and Environmental Engineering University of Washington Seattle, WA USA

## Xuegang (Jeff) Ban

Department of Civil and Environmental Engineering University of Washington Seattle, WA USA banx@uw.edu

August 2, 2024

## ABSTRACT

This study aims to discover the governing mathematical expressions of car-following dynamics from trajectory data directly using deep learning techniques. We propose an expression exploration framework based on deep symbolic regression (DSR) integrated with a variable intersection selection (VIS) method to find variable combinations that encourage interpretable and parsimonious mathematical expressions. In the exploration learning process, two penalty terms are added to improve the reward function: (i) a complexity penalty to regulate the complexity of the explored expressions to be parsimonious, and (ii) a variable interaction penalty to encourage the expression exploration to focus on variable combinations that can best describe the data. We show the performance of the proposed method to learn several car-following dynamics models and discuss its limitations and future research directions.

Keywords Car-following dynamics · Deep symbolic regression · Deep learning · Variable interaction selection

## 1 Introduction

There are two recent trends in transportation and the broader science/engineering fields, which make the headlines almost every day. The first one is the emergence of connected/automated vehicles (CAVs) that i) may introduce new, complex traffic dynamics and interactions in the current and future traffic streams, and ii) generate increasingly available and massive datasets from both vehicles and the infrastructure. The second trend is the rapid development and application of deep learning techniques that seem to revolutionize almost every aspect of technology, science, engineering, and the entire society. While there have been numerous studies and applications of deep learning in transportation, in the paper, we are interested in the question of whether deep learning can help discover traffic dynamics (car-following models in particular) from data directly with no or little human involvement. An affirmative answer to this question will not only help discover/develop traffic dynamics models in this era but also have important implications for other science/engineering fields where dynamical systems and their governing equations are widely used and studied.

Car-following depicts the driving behavior of how a vehicle (driver) follows and interacts with the vehicle in front of it. It is one of the basic traffic models in revealing traffic dynamics characteristics at the microscopic traffic flow level Brackstone and McDonald [1999]. Car-following studies can be traced back to the 1950s and 1960s when Pipes [1953], Chandler et al. [1958], Kometani and Sasaki [1958], Gazis et al. [1959, 1961], and Helly [1959] initiated an era of modeling car-following and traffic dynamics. Pipes [1953] modeled car-following dynamics on the assumption that

## James Enouen

Department of Computer Science University of Southern California Los Angeles, CA USA

## Yan Liu

Department of Computer Science University of Southern California Los Angeles, CA USA

a following vehicle ( the follower ) is required to maintain a certain distance from the preceding vehicle ( the leader ). Chandler et al. [1958] further analyzed vehicle movements by a linear relationship between the acceleration of the follower and the velocity gap of the leader-follower pair. The model is given in Eq. (1) below. Here a f ( t + δt ) , v l ( ) t , and v f ( ) t , respectively, represent the follower's acceleration at the next time step, and the velocity of the leader and the velocity of the follower at the current time step, with c the so-called sensitivity term that may depend on the system states.

$$a _ { f } ( t + \delta t ) = c [ v _ { l } ( t ) - v _ { f } ( t ) ].$$

The model listed in Eq.(1) is named the General Motors (GM) model and is widely regarded as the foundational model based on which many later car-following models were developed, e.g., by specifying various forms of the sensitivity term c . It states that the follower will accelerate or decelerate depending on the speed difference between the leader and the follower, while the action is with a time delay δt representing the reaction time of the follower. The model is simple and intuitively makes sense, and was also validated by real-world data Chakroborty and Kikuchi [1999]. It can thus be considered as the "governing equation" of the underlying dynamical system of car-following. The example of the above GMmodel also illustrates the unique feature of existing, physics-based car-following models: they are concise, explicit analytical models that have clear physical meanings and generalize well with different datasets/sites.

The development of existing car-following models can be generally summarized as a scientific process Treiber and Helbing [2002], Wen-xing and Lei [2007], Pourabdollah et al. [2017], Kurtc [2020], which includes key steps such as problem definition, background research, experimental design, data collection, data analysis, hypothesis, and testing/validation. This is often an involved, iterative process that is time and resource-consuming. The data analysis and hypothesis steps are arguably the most intellectually challenging and most creative in the entire process. This is because the data analysis step requires innovative ways with intuition/knowledge to interpret the data, properly visualize, reveal salient patterns, and discover the right variables and parameters to build suitable car-following relations, as well as to apply knowledge to compare/contrast the patterns with existing models in the same or similar fields. This data analysis step draws useful insights that serve as the input to the hypothesis step. This later step then proposes a suitable dynamical model to best describe the discovered patterns and relations in the data analysis step, while at the same time ensuring that the proposed model is tractable for analysis and computation. This data analysis hypothesis -process often involves many iterations for trial and error, with heavy human involvement.

With the advances in deep learning and emerging technologies in transportation, we wonder if the data analysis -hypothesis steps in the existing scientific process can be replaced by a deep learning approach with minimum human interaction, for discovering the governing equations of car-following dynamics. The motivation for this is multifaceted. First, while predictive models Chen et al. [2021], Kamienny et al. [2022], Shi et al. [2021], including those predicting the car-following behavior Mo et al. [2021], have been increasingly developed in transportation based on deep learning techniques and outperformed traditional physics-based models, they are often a black box that lacks interpretations. We believe that concise models with clear physical interpretations, as traditionally done in transportation and other fields, are still crucial as "human knowledge" and in understanding and modeling new phenomena with emerging technologies and systems. This leads to the second motivation, i.e., important practical needs to develop such learning-based methods. Existing car-following models were developed mainly for human-driven vehicles (HDVs). With the emergence and rapid deployment of CAVs, it is imperative to understand the interactions, including car-following, among CA Vs and between CAVs and HDVs Zhong et al. [2020], Xiao et al. [2018]. Distinct to CAVs, CAV and HDV data may be collected by CAV themselves or from the infrastructure Gloudemans et al. [2023], which may make the classical scientific process cumbersome to discover the governing equations for new car-following and other traffic dynamics. Deep learning models may be a promising alternative method. Third, it is an intriguing research topic to develop a fully data-driven approach integrated with deep learning, with no or little human involvement, to directly produce a car-following model that shares the same features with the existing models: concise, interpretable, and generalizable. In the meantime, a successful development of such learning-based methods can provide useful insights into the interpretability of deep learning techniques, an open question that is drawing much attention nowadays.

Developing the aforementioned learning techniques for car-following from data can be viewed as a data-driven discovery problem Rudy et al. [2017], Brunton et al. [2016], Antoniou et al. [2011]. Early research has been primarily investigated on how to draw physical laws from ideal scientific experiments using classical methods such as regression models. In the early 2000s, Bongard and Lipson [2007] and Schmidt and Lipson [2009] proposed an approach named "symbolic regression (SR)" for identifying the underlying governing equation of a dynamical system directly from data. The purpose of the SR approach is to determine a mathematical expression f : R d → R , X i ∈ R d , y i ∈ R , and a dataset ( X,y , ) such that f is a short mathematical expression and y ≈ f ( X ) . The SR model was solved by Genetic Programming (GP) [Koza, 1994] based approaches, which were later widely extended [Bongard and Lipson, 2007,

Schmidt and Lipson, 2009, Bäck et al., 2018, La Cava et al., 2016, 2018, Virgolin et al., 2019, 2021, Kommenda et al., 2020, de Franca and Aldeia, 2021, Petersen et al., 2019]. We refer the readers to Section 2 for more detailed reviews.

Despite the above progress of learning-based discovery methods, discovering the underlying governing expressions from real-world traffic data that often involves noise and uncertainties can be challenging. Further, integrating with GP to solve SR problems may struggle with the downside of GP, such as sub-optimal and complex solutions. This study applies deep symbolic regression (DSR) Petersen et al. [2019] and GP to uncover the underlying car-following mathematical models from trajectory data. In particular, we make two improvements to the DSR-GP framework. First, DSR-GP is given some prior knowledge before learning the dynamics. The prior knowledge is a set of feature (variable) combinations extracted from the trajectory data through the variable intersection selection (VIS) algorithm proposed by Enouen and Liu [2022]. Through this prior knowledge, the framework is guided to mathematical expressions that have reasonable physical meanings. Second, we propose a new reward function to encourage exploring low-complexity mathematical expressions. These two improvements can lead to expressions that are parsimonious and interpretable. We call the proposed learning framework the VIS-enhanced DSR.

In the next section, we briefly review the DSR and VIS methods. The VIS-enhanced DSR framework is presented in Section 3. In this framework, we first apply VIS to discover variable combinations from a trajectory dataset and then use DSR to learn the dynamics based on the discovered variable combinations. The expression reconstruction performance of the proposed framework is compared with the performance using only DSR and DSR with GP. Further, we investigate the impact of data noise and the performance of the method on different types of car-following dynamics. The results and comparisons are included in Section 4. We discuss the limitation of the proposed method in Section 5, with conclusions and future research directions presented in Section 6.

The main contributions of this study are summarized as follows.

- · Deep learning models are proposed for identifying car-following dynamics from trajectory data with minimum human involvement.
- · VIS-enhanced DSR framework is developed to find mathematical expressions that tend to be parsimonious and interpretable.
- · Reward function in DSR is improved by new penalty terms to regulate the complexity of the explored expressions and to search expressions with desirable feature combinations.

## 2 Related Works

In this section, we briefly review car-following models and the data-driven models in solving SR problems, including the GP methods and the DSR framework Petersen et al. [2019]. We also briefly review the VIS approach Enouen and Liu [2022].

## 2.1 Car-following Models

Most of the pioneering car-following models Kometani and Sasaki [1958], Pipes [1966], Herman et al. [1959], Chandler et al. [1958], Edie [1961], Helly [1959], Gazis et al. [1959, 1961] between fifties to seventies modeled car-following dynamics as an acceleration function with the differences of positions and speeds between the leader and the follower. Later on, car-following models were extended to consider multiple scenarios including the follower's emergency braking, following the leader, approaching the leader, free moving, and accelerating from rest, where notable works include Wiedemann and Reiter [1992], Michaels [1963], Burnham and Bekey [1976], Lee [1976], Kumamoto et al. [1995], Krauß [1998], Treiber et al. [2000]. In particular, the Krauss model Krauß [1998] considers that the follower selects its action from its desired speed, its comfortable acceleration, or the speed that is able to maintain a safe distance from the leader. The Intelligent Driver Model (IDM) [Treiber et al., 2000] further considers a "soft braking strategy" when the follower approaches a slow-moving leader. It has been shown that the Krauss model performs better over IDM when modeling unsteady traffic conditions, whereas IDM performs better under steady traffic conditions Kanagaraj et al. [2013], Naveen Kumar [2013]. In the numerical experiments of this paper, we use a traffic simulation model that integrates the Krauss model for car following. The model is defined by Eq. (2), where t react , v f ( ) t , v l ( ) t , s f ( ) t , and s des respectively represent follower's reaction time, the velocities of the follower and the leader at the current time step, the spacing at the current time step, and the desired spacing ( s des ) that is defined by measuring the distance that the leader is estimated to drive forward in the next time interval, i.e., s des = v l ( ) t δt . The equation describes that given the current velocities of the leader and the follower, and the spacing between the pair, the follower selects its velocity from the choices of its current velocity plus its maximum acceleration v f ( )+ t a max δt , safe velocity v s , and the maximum velocity that this vehicle is capable of driving at, i.e., v max that is a known parameter. Other parameters that are pre-defined and

assumed fixed in this model include the maximum acceleration a max, comfortable deceleration b , driver's reaction time t react , driving uncertainty ε , and time step δt .

$$v _ { f } \left ( t + \delta t \right ) & = \max ( 0, \, v _ { d } - \varepsilon a _ { \max } ) \\ v _ { d } & = \min ( v _ { f } ( t ) + a _ { \max } \delta t, v _ { s }, v _ { \max } ) \\ v _ { s } & = v _ { l } ( t ) + \frac { s _ { f } ( t ) - s _ { \text{des} } ( t ) } { \frac { v _ { f } ( t ) + v _ { l } ( t ) } { 2 b } + t _ { \text{react} } }. \\. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

As this paper concerns about discovering existing car-following models using deep learning methods, we do not aim to conduct a comprehensive review of car-following models; interested readers may refer to Brackstone and McDonald [1999], Treiber and Kesting [2013] for more recent reviews.

## 2.2 Data-driven Discovery Methods

GP Koza [1994] was used to solve SR problems Bongard and Lipson [2007], Schmidt and Lipson [2009], which was shown to be effective by leveraging the evolutionary features of GP to improve fitting functions [Koza, 1994, Schmidt and Lipson, 2009, Bäck et al., 2018]. Many later approaches for solving SR problems were based on the concept of GP, including the epsilon-lexicase selection (EPLEX) La Cava et al. [2016], the feature engineering automation tool (FEAT) La Cava et al. [2018], semantic backpropagation GP (SBP-GP) Virgolin et al. [2019], the GP-based gene-pool optimal mixing evolutionary algorithm (GP-GOMEA) Virgolin et al. [2021], Operon Kommenda et al. [2020], and the interaction-transformation evolutionary algorithm (ITEA) de Franca and Aldeia [2021]. Despite the powerful solving capabilities from the nature of evolution, GP methods suffer from generating appropriate expressions that are parsimonious and general [O'Neill et al., 2010], and fitting may take a long time due to the large size of search space.

Brunton et al. [2016] developed the sparse identification of nonlinear dynamics (SINDy) framework that describes the expression discovery problem by a sparse regression problem. The framework is based on an underlying assumption that the dynamics of a physical system can be represented by a linear combination of a limited number of mathematical terms. The SINDy approach was integrated with autoencoders by Champion et al. [2019] and further with Bayesian approaches by Gao and Kutz [2022]. However, SINDy users need to provide all the possible feature combinations for the framework to discover important terms. Jin et al. [2019] incorporated Bayesian approaches in solving SR problems. In the deep learning domain, the AIFeynman approach [Udrescu et al., 2020] is also worth mentioning. It is inspired by the divide-and-conquer method powered by neural networks. Cranmer et al. [2020] and Kamienny et al. [2022] applied graph neural networks (GNNs) and transformers to solve SR problems, respectively. Similar to the works by Champion et al. [2019] and Gao and Kutz [2022], which discovered the governing behaviors from the intrinsic coordinates, Chen et al. [2021] proposed a two-stage autoencoder framework to estimate the intrinsic variables and use the variables to predict the behaviors at the next time step. The intrinsic patterns are however hidden in the autoencoders and thus are not interpretable.

DSR Petersen et al. [2019] is a recent framework that uses a recurrent neural network (RNN) with the proposed risk-seeking policy gradient approach for solving SR problems. DSR leverages the conditional probability feature in RNNs to increase the parsimony and accelerate the exploration efficiency. It was developed on the basis of the fact that a mathematical expression can be represented by an expression tree [Meurer et al., 2017] composed of mathematical tokens, where the tokens are sampled by RNNs from a pre-defined token pool. The sampling is restricted by fixed constraints to restrict token combinations and lengths of expressions, i.e., expression complexities. Given a range of expression complexity as an exploration constraint, the RNN will stop sampling if the expression has enough complexity and the expression tree is complete. This mechanism may result in expressions with a high level of complexity that lose parsimony and interpretability, which is also the main issue of GP. Data noise is another challenge: results showed that 10% of data noise affected over 10% of expression reconstruction, although it outperformed benchmark solvers, including GP, Eureqa Schmidt and Lipson [2009], etc. Expression searching in the DSR framework relies on training RNNs or long-short-term-memory (LSTM) neural networks, a special structure of RNNs, whose training can sometimes be computationally expensive. Combining GP with DSR has also been proposed Larma et al. [2021], for which GP is an option to assist the search.

## 2.3 Variable Interaction Selection (VIS)

To encourage the exploration to search for expressions that have reasonable lengths and physical meanings, we add a VIS step prior to DSR to provide knowledge to the exploration. This VIS step produces potential variable combinations extracted from data. The VIS algorithm was introduced in Enouen and Liu [2022] in their proposed sparse interaction additive networks (SIANs). SIANs leverages the generalized additive function, which speculates that a dependent

variable is linearly related to an addition of a set of shape functions. A shape function can be any specific function, e.g., a polynomial function of a single variable or variable combinations that have a strong interaction, where the interaction strength is quantified by Archipelago Tsang et al. [2020]. The algorithm first fits a deep neural network to the data, where the deep neural network acts as the reference function. The algorithm detects all the possible variable combination candidates and the selection is based on several key criteria including a pre-defined threshold of the interaction strength. That is, a variable combination is determined as a strong interaction if its interaction strength is greater than the threshold. We next present the VIS-enhanced DSR method to discover car-following models from data.

## 3 VIS-Enhanced DSR Method

Fig. 1a depicts the flowchart of our proposed VIS-enhanced DSR method. On the left of the flowchart, we use a VIS process to detect potential variable combinations from the trajectory data to increase the interpretability of the explored expressions. On the right of the flowchart, we use an expression checker to detect whether the expression determined by the RNN sampling has reasonable physical meanings based on the recommended variable combinations. The RNN training is governed by our proposed reward function that encourages parsimonious expressions for better interpretation. Fig. 1b illustrates a pipeline example of the VIS process. Given a set of trajectory data, the process trains a deep neural network to match the data and then produces a set of variable combinations that can capture the salient patterns in the data. The identified variable combinations can guide DSR to search for mathematical expressions that contain such variables. To introduce the framework and the reward function in detail, we first define the problem and notations as follows.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00251v1-figure-47af0772-4066-49bf-8207-1313bd49fb3d.jpg)

(a) Expression Exploration Flowchart

(b) Variable Interaction Pipeline Enouen and Liu [2022]

Figure 1: Expression Exploration Flowchart And An Example

Define the same mathematical expression f ( X ) and dataset ( X,y ) introduced in Section 1, i.e.,

$$y = f ( X ) + \epsilon,$$

where ϵ represents the approximation error term. We aim to discover the exact form of f ( X ) from the dataset. Given a set of mathematical tokens T , which can include mathematical operations ( + , -, sin , constant placeholders, etc.) and variables, the DSR framework samples each mathematical token τ i from T according to the probabilities

determined by the RNN. The token in T with the highest probability will be sampled and stored in the latest node spot. Each mathematical term is sampled without replacement. We denote the resulting expression by f s ( τ ) , where τ = { τ i } n s i =1 , ∀ τ i ∈ T , and n s represents the number of sampling occasions, i.e., expression complexity. We apply the DSR approach to find f s ( τ ) through sampling tokens from the given T such that ϵ is small. As an example, f s ( τ ) can be αx + sinx if the sample token set is { + , × , α x , , sin, x }, where x is the only variable and α is a constant.

As shown in Fig. 1a, we first feed the collected trajectory dataset, ( X,y ) , into both the VIS and DSR blocks. VIS starts to find the variable interactions from the trajectory data, as shown in Fig. 1b, by first fitting a deep neural network to the data, calculating the interaction strength for variables or variable combinations, and determining whether a variable or combination is significant to the data relationship. In the meantime, DSR starts to explore potential expressions that best describe the relationship between X and y . The potential expressions are then fed into the expression-checking process, which will check each expression candidate term by term for reward evaluation. The checking process will also assess whether the expression has recommended variable combinations by VIS and tag those without any recommended combinations (so that a penalty may be imposed as shown in Section 3.2). The DSR process will then continue evaluating the reward of the retained expression candidates, training the RNN, and exploring potential expressions.

## 3.1 Variable Interaction Selection (VIS)

In the VIS process, we consider an extended generalized additive model Enouen and Liu [2022], in which the dependent variable linearly depends on some functions consisting of single variables and/or groups of variable combinations, as defined in Eq.(4). In the equation, f ∅ is a constant term. Each h i ( X i ) function is a shape function that represents a specific function, e.g., a polynomial function. X i is an element of the variable set, i.e., X i ∈ I ⊆ { X i } d i =1 . The last term considers the variable combinations, where Z j (e.g., Z j = ( X ,X 1 2 ) ) represents an element in the set of variable combinations. The variable combination set is defined as J , Z j ∈ J .

$$y = f _ { \varnothing } + \sum _ { X _ { i } \in \mathcal { I } } h _ { i } ( X _ { i } ) + \sum _ { Z _ { j } \in \mathcal { J } } h _ { j } ( Z _ { j } ).$$

In VIS, a deep neural network is first trained to fit the data. We use f dnn ( X ) to denote the deep neural network. The interaction strengths of variables (in set I ) or variable combinations (in set J ), denoted as ψ ( I ) and ψ ( J ) respectively, can be measured by the Archipelago Tsang et al. [2020], which uses the approximation of the Hessian of f dnn ( X ) as shown in Eq. (5). The VIS algorithm proposed in Enouen and Liu [2022] is then used to detect strong variable interactions performed in the data based on f dnn ( X ) and Eq. (5).

$$\psi ( \mathcal { I } ) \coloneqq \mathbb { E } \left [ \frac { \partial ^ { | \mathcal { I } | } f ^ { \text{dnn} } ( X ) } { \partial X _ { i _ { 1 } } \partial X _ { i _ { 2 } } \cdots \partial X _ { i _ { | \mathcal { I } | } } } \right ] ^ { 2 } > 0. \\ \psi ( \mathcal { J } ) \coloneqq \mathbb { E } \left [ \frac { \partial ^ { | \mathcal { J } | } f ^ { \text{dnn} } ( X ) } { \partial Z _ { j _ { 1 } } \partial Z _ { j _ { 2 } } \cdots \partial Z _ { j _ { | \mathcal { J } | } } } \right ] ^ { 2 } > 0.$$

## 3.2 Improved DSR for Discovering Car-following Dynamics

In the DSR process, we use an LSTM neural network to sample tokens. The LSTM is trained through the risk-seeking policy gradient algorithm proposed by Petersen et al. [2019], which updates the LSTM weights through the topperforming cases (instead of the average-performing cases as normally done). We redesign the original reward function in Petersen et al. [2019] to encourage the LSTM to (i) explore expressions consisting of the recommended variable combinations generated by VIS, and (ii) explore expressions with reasonable expression complexity. The new reward function is shown in Eq. (6), where R c represents the combined reward function, L e represents an error-based loss function, p is the complexity of an expression f s , and β is a penalty parameter that penalizes expressions that do not contain the recommended variable combinations from VIS. The penalty only steps in for the first EPOCH epochs to initially guide the exploring direction and prevent the exploration from being time-consuming, as checking each term of an expression needs extra time.

$$R & = \begin{cases} R _ { c } ( 1 - \beta ), & \text{if no potential interactions and epoch} \leq \text{EPOCH} \\ R _ { c }, & \text{otherwise}, \end{cases}, \\ R _ { c } & = \frac { 2 } { ( 1 + L _ { e } ) + ( 1 + \text{norm} ( p ( f _ { s } ) ) ) },$$

As shown in Eq. (6), the combined reward function R c considers the error-based loss and the complexity of the expression candidate. R c is squashed to a range between 0 to 1. R c → 1 0 . when the error loss and the loss of complexity are both close to zero, i.e., L e → 0 and norm ( p f ( s )) → 0 . R c goes down when the predicted error and/or the complexity of an expression increase. To encourage the LSTM to explore expressions with reasonable complexities, we design the loss of expression complexity as a normalization function. The normalization relies on human knowledge to provide a reasonable range, as defined in Eq. (7), where p min and p max represent the minimum expression complexity and maximum expression complexity, respectively, which are pre-defined based on human knowledge.

$$\text{norm} ( p ( f _ { s } ) ) = \begin{cases} \frac { p ( f _ { s } ) - p _ { \min } } { p _ { \max } - p _ { \min } }, & \text{if $p_{\min} \leq p(f_{s})$} \\ \frac { p _ { \max } - p ( f _ { s } ) } { p _ { \max } - p _ { \min } }, & \text{otherwise}. \end{cases}$$

The equation's expression complexity is scaled into a base range of 0 to 1. When the complexity is greater than or equal to the minimum p min , it is scaled by the top conversion of the equation. To prevent the same conversion from scaling the complexity to a negative number when the complexity is smaller than the minimum, we use the bottom conversion of the equation to scale the complexity. The bottom conversion comes from the idea of treating the smaller complexity (smaller than the desired range) as the greater complexity (greater than the desired range) to

(

)

prevent producing negative losses. That is, if p f

(

p

max

+(

p

min

-

p

)

-

p

min

s

)

&lt; p min

,

then norm

(

p f

(

s

)) =

p

max

-

p

min

,

and rewriting

$$\frac { \left ( _ { \max } + ( p _ { \min } - p ) \right ) - p _ { \min } } { p _ { \max } - p _ { \min } } = \frac { \left ( _ { \max } - p \right ) } { p _ { \max } - p _ { \min } } \, \text{obtains the bottom conversion of $Eq.$} (7).$$

L e here captures the "accuracy" of the expression in describing the data. It can be any predicted loss function, such as the root-mean-squared error (RMSE), mean-squared-error (MSE), normalized MSE (NMSE), NRMSE, log MSE, etc. This study uses NRMSE as L e . Given the expression formed by sampled tokens f s ( τ ) and dataset ( X,y ) , L e is defined as Eq. (8), where σ y is the standard deviation of y .

$$L _ { e } = \frac { \sqrt { \frac { 1 } { n } \Sigma _ { i = 1 } ^ { n } ( y _ { i } - f _ { s } ( \tau, X _ { i } ) ) ^ { 2 } } } { \sigma _ { y } }.$$

Before concluding this section, we summarize the human involvements/knowledge that are required by the proposed VIS-enhanced DRS method. A set of tokens is first required, including variables ( X ) and mathematical operations. The variables can be relatively easily obtained as the key elements in the data: e.g., for trajectory data, this can be speeds, accelerations, locations, etc. Operations are a bit trickier and require certain level of human knowledge, e.g., the operations used in existing car-following models if new car-following models are to be developed. The good news is that we do not need to provide the exact set of operations, e.g., multiplication and subtraction for the car-following model in Eq. (1). Rather a superset of the needed operations is sufficient, e.g., the set of operations that contain × ÷ , , + , -for Eq. (1). Other human inputs include the min and max complexities in Eq (7) and the β parameter in Eq (6). These may be treated as hyberparameters for machine leanring, which may be pre-defined based on knowledge or learned by trial and error. Compared with traditional methods of deriving a car-following model from observed data, we can clearly see that the human involvement in the proposed method is marginal since the data analysis and hypothesis steps are largely automated via the VIS and DSR processes.

## 4 Numerical Results

We now use the proposed VIS-enhanced DSR framework to help discover the potential governing mathematical expressions from car-following trajectory datasets. We first present how the datasets are prepared.

## 4.1 Dataset and Problem Statement

The trajectory data is collected through numerical simulation with car-following pairs using the Simulation of Urban MObility (SUMO). To enhance the variety of the car-following behaviors, we control the leading vehicle by randomly selecting the velocity between 0 to 30 m/s, whereas the decisions of the following vehicle are described by the Krauss model Krauß [1998], Krajzewicz et al. [2002], as defined in Eq. (2), with predefined vehicle parameters: a max = 2 6 . m/s , 2 v max = 55 55 . m/s, b = 4 5 . m/s , and 5-m vehicle length. The follower is first assumed to have zero 2 driving uncertainty, i.e., ε = 0 0 . , and their reaction time is a second, i.e., t react = 1 0 . second. Later, we also impose data noises and test how the proposed method can deal with such noises. The trajectory dataset has in total ∼ 3.6k car-following samples (pairs). Given the parameters, the Krauss model can be simplified as Eq. 9. The function will

be the target expression in the following experiments. From the simulation, we collect the velocities of the follower and the leader at the current time step, the spacing at the current time step, the following velocity at the next time step, and the difference between current spacing and the desired spacing, i.e., , v f ( ) t , v l ( ) t , s f ( ) t , v f ( t + δt ) , and ds t ( ) = s f ( ) t -s des ( ) t , where s des ( ) = t v l ( ) t δt .

The data is organized into independent variables X = { v f ( ) t , v l ( ) t , s f ( ) t , ds ( ) t } , whereas the dependent variable is the velocity of the follower at the next time step, i.e., y = v f ( t + δt ) . The problem is to develop a mathematical expression that describes the relationship between the collected features and the next velocity of the follower, i.e., v f ( t + δt ) = f ( v f ( ) t , v l ( ) t , s f ( ) t , ds ( )) t .

$$v _ { f } \left ( t + \delta t \right ) & = \min ( v _ { f } ( t ) + 2. 6, v _ { s } ) \\ v _ { s } & = v _ { l } ( t ) + \frac { 9. 0 ( s _ { f } ( t ) - s _ { d e s } ( t ) ) } { v _ { f } ( t ) + v _ { l } ( t ) + 9. 0 }.$$

Note that certain simplifications are performed to derive (9), including only vehicles from a single class (cars), setting the reaction time to be 1s, as well as fixing the maximum acceleration and speed. In Section 5, we discuss these issues and possible ways to relax them.

## 4.2 Variable Interactions

The interaction selection process starts by training a 3-layer (in-32-64-32-out) neural network. The trajectory dataset is split into an 80% training set and a 20% validation set. After training, the mean squared error of the deep neural network estimation drops to 7.5E-04, 7.1E-04, and 4.6E-04 after 50, 100, and 300 epochs of training, as shown in Fig. 2. The training takes 5.8, 16.9, and 60.6 seconds at 50, 100, and 300 epochs on an Intel Core i7-12700 CPU.

Figure 2: DNN Training Loss

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00251v1-figure-6f100a1b-568d-4cfd-8303-02cc9a9e9989.jpg)

Given variables v f ( ) t , v l ( ) t , s f ( ) t , and ds t ( ) , the potential variable combinations along with their interaction strengths after 300 epochs of training are listed in Table 1, where variable combinations are listed in ascending order by the interaction strength presented in Fig. 3. In this figure, a few elbow points can be observed, including strengths of v l ( ) t , ( ds t ( ) , v l ( ) t , v f ( ) t ) , ( v l ( ) t , v f ( ) t ) , ( ds t ( ) , v f ( ) t ) , and ( ds t ( ) , s f ( ) t ) . Based on these elbow points, five sets of variables or variable combinations are recommended as follows.

Set 0 { v f ( ) t , v l ( ) t },

Set 1 { ds t ( ) , ( ds t ( ) , v l ( ) t , v f ( ) t ) },

Set 2 { ( ds t ( ) , v l ( ) t ) , ( v l ( ) t , v f ( ) t ) },

Set 3 { ( ds t ( ) , v f ( ) t ) },

Set 4 { ( ds t ( ) , s f ( ) t ) }.

The two variables in Set 0 are unlikely to be the only features (variables) describing the car-following behavior since the spacing difference may also be an important feature. We therefore combine Set 0-1 as the first possible recommended set of variable combinations, denoted as Scenario 1 . The rest sets are respectively coded up as Scenario 2 to Scenario 4 . Table 1 presents all variable combinations, their interaction strengths, and the scenarios. For instance, Scenario #1 uses Set 0 and Set 1, and Scenario #2 uses Set 0 - Set 2.

Table 1: Variable Interaction Strengths

| Interaction                                  |   Interaction Strength (IS) |   log (IS) | Scenario   | Set   |
|----------------------------------------------|-----------------------------|------------|------------|-------|
| v f ( t )                                    |                      347.07 |       5.85 |            |       |
| v l ( t )                                    |                       15.12 |       2.72 | #1         | 0 &1  |
| ds ( t )                                     |                       11.12 |       2.41 |            |       |
| ds ( t ) , v l ( t ) , v f ( t )             |                        9.6  |       2.26 |            |       |
| ds ( t ) , v l ( t )                         |                        5.59 |       1.72 | #2         | 0 - 2 |
| v l ( t ) , v f ( t )                        |                        3.21 |       1.17 |            |       |
| ds ( t ) , v f ( t )                         |                        3.04 |       1.11 | #3         | 0 - 3 |
| ds ( t ) , s f ( t )                         |                        0.41 |      -0.89 | #4         | 0 - 4 |
| s f ( t )                                    |                        0.29 |      -1.22 | -          |       |
| v l ( t ) , s f ( t )                        |                        0.28 |      -1.27 | -          |       |
| ds ( t ) , v l ( t ) , s f ( t )             |                        0.13 |      -2.08 | -          |       |
| v l ( t ) , s f ( t ) , v f ( t )            |                        0.09 |      -2.37 | -          |       |
| s f ( t ) , v f ( t )                        |                        0.09 |      -2.38 | -          |       |
| ds ( t ) , s f ( t ) , v f ( t )             |                        0.07 |      -2.67 | -          |       |
| ds ( t ) , v l ( t ) , s f ( t ) , v f ( t ) |                        0.06 |      -2.87 | -          |       |

Figure 3: Interaction Strength Elbow Plot

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00251v1-figure-434d3c17-09bd-4109-a6d4-ddc6d03d7aa3.jpg)

## 4.3 Expression Exploration

Based on the literature and analysis in the past [Pipes, 1953, Chandler et al., 1958, Gazis et al., 1959, 1961, Helly, 1959, Pipes, 1966, Krauß, 1998, Treiber et al., 2013], car-following behaviors are related to the velocity and gap differences with + , - × ÷ , , , min and parameters a max, and b . We thus provide + , - × ÷ , , , min , v f ( ) t , v l ( ) t , ds t ( ) , s f ( ) t , a max, and b as the set of tokens to the DSR process. If we allow exploration to focus solely on fitting accuracy without constraints and complexity regulations, we may end up with a complex expression that lacks physical meaning, e.g.,

the expression below, which may have a high level of accuracy but also a high level of expression complexity. Such expression is certainly not desirable.

(

```
& \left ( s _ { f } ( t ) \left ( - v _ { l } ( t ) / \left ( ( v _ { f } ( t ) \, - \, a _ { \max } / \left ( v _ { f } ( t ) s _ { f } ( t ) \, \min \{ - 2 b, \, v _ { l } ( t ) \} \right ) \right ) \min \{ v _ { f } ( t ), \, v _ { l } ( t ) \} \right ) \, - \, 2 b \right ) \, + \\ & \quad s _ { f } ( t ) \, - \min \{ a _ { \max }, \, s _ { f } ( t ) \} \, - \min \{ v _ { f } ( t ), \, a _ { \max } \} \, - \min \{ s _ { f } ( t ), \, - 2 b \, v _ { l } ( t ) \} \right ) / \left ( v _ { f } ( t ) \, + v _ { l } ( t ) \right ).
```

In the following expression exploration, min is limited to only being sampled once and should be the first token sampled. We configure a size of 300 to be the batch size of exploration. The reward penalty to encourage recommended variable interactions will be turned off after 10 epochs of exploration, i.e., EPOCH = 10 in Eq. (6), to prevent too much running time while gaining initial guidance. The classic GP method is set to be allowed to assist expression exploration. In the following, we use DSR based on the variable combinations recommended by the VIS method with the proposed complexity penalty reward strategy. All of the cases in the remainder of this subsection are tested using a noise-free dataset. We compare the performances of the pure DSR framework Petersen et al. [2019], the DSR framework with GP assistance Larma et al. [2021], and our proposed VIS-enhanced DSR framework with GP assistance. We use DSR-only, DSR-GP, and VIS-DSR-GP in the following to refer to the three methods, respectively. The DSR-GP method is configured with the same constraints of limiting min to only being sampled once and first. Additionally, DSR-GP is configured with a constraint of limiting the length of explored expressions to a range of 10 to 40.

## 4.3.1 Expression Exploration with Penalty

We first compare explorations using different penalty parameters of the complexity penalty strategy, i.e., β in Eq. (6). In the experiment, β s are tested with values from 0.00 to 1.0, with a 0.05 increment from 0.00 to 0.20 and then with a 0.20 increment from 0.20 to 1.0. All the tests are based on the recommended variable combinations in Scenario #1 shown in Table 1. Note that β = 0 indicates R = R c , i.e., there is no penalty if DSR explores expressions with variable combinations not recommended by the VIS.

We compare exploration epochs until converging using the VIS-DSR-GP and DSR-GP methods with different β 's in Fig. 4, where each result is averaged over 10 random seeds with the standard errors shown by the error bars in red. All have discovered the target expression as Eq. (9). DSR-only is excluded from the comparison, as it does not recover the target expression. Fig. 4a and Fig. 4b compare the numbers of learning epochs and running times over different β s, where each value is also the average of the results from running the test over 10 random seeds. The comparisons show that 0.15 is a sweet spot of β (VIS-DSR-GP), as it takes the fewest number of epochs and the shortest running time to discover the target expression. Even adding the running time of VIS (60.6 seconds from Section 4.2), VIS-DSR-GP with β = 0 15 . is still faster than using DSR-GP, which takes more than 150 epochs and over 200 seconds to discover the target expression. Choosing β = 0 00 . takes more epochs to learn the dynamics, but the running time is relatively short. It is because in the test of β = 0 00 . , the exploration does not check each expression term, and the reward takes only accuracy and parsimony into account, i.e., R = R c in Eq. (6). Using DSR-GP is found to tend to look for better expressions through higher complexity, as shown in Fig. 4c that the high range of expression complexities (20 to 40 given the desired range from 10 to 40) dominates the exploration compared to the range of low complexities (10 to 20). As shown in Fig. 4d, running multiple tests with the DSR-GP method, we found that this method tends to try exploring better expressions through a high level of expression complexity. Our proposed method, on the other hand, encourages the exploration of parsimonious expressions, as presented in Fig. 4c, where the expression complexities explored are mostly in a range of 10 to 20.

(a) Reward Comparisons with Standard Errors

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00251v1-figure-7946c519-5fad-4793-85ad-7e1d7a3f7142.jpg)

Exploration with Penalty Factors

(b) Elapsed Time Comparisons with Standard Errors Exploration with Penalty Factors

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00251v1-figure-57a5f3aa-9966-425b-a13b-5757ffa8cee8.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00251v1-figure-ea4eeeb2-97c4-4297-8e1a-074b98207b1d.jpg)

(c) Density of Complexities Explored

(d) Density of Complexities Explored (DSR-GP)

Figure 4: Exploration by Penalty Parameters

## 4.3.2 Expression Exploration with Recommended Variable Combinations

The above experiments indicate that using β = 0.15 performs the best in terms of the number of exploration epochs and running times. We adopt β = 0 15 . and apply the same constraints to test explorations with different recommended variable combinations listed in Table 1. The following tests still use the noise-free dataset. Fig. 5 presents a performance comparison using different sets of recommended variable combinations, where blue bars show the average of running 10 random seeds, and error bars in red demonstrate the standard errors. The target expression is discovered by all the test cases. The results show that a good set of variable combinations can help reduce the number of exploration epochs and further improve the running time. With that said, providing too many conditions to regulate expression terms may result in long running times, as the elapsed times of Scenarios #2 - #4 shown in Fig. 5b. In particular, Scenarios #3 and #4 take a longer time to find the target expression than DSR-GP.

Figure 5: Exploration Performance with Variable Interactions

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00251v1-figure-494d903b-1178-488b-af06-11a4c9806006.jpg)

## 4.4 Noise Impact Investigation

The above experiments are based on perfectly simulated data described by only one governing equation. Here we test how data noises may impact the performance of the proposed method. The following tests in this subsection are given the set of variable combinations recommended in Scenario #1 (see Table 1) with β = 0 15 . and the same constraint configuration, and a recommended expression complexity ranging from 10 to 40.

In the investigation, the dataset is contaminated 100% by different noise levels, which increase from 0 to 10% with a 1% increment of the standard deviation of the clean data. Fig. 6a shows an example of the clean data, data with 5%-level noise, and data with 10%-level noise. The VIS algorithm yields the same recommended variable combination sets as listed in Table 1 at all noise levels. Thus, the following focuses on comparing the performances of the DSR process under different noise levels. The performances in terms of mean percent errors (MPE) are presented in Fig. 6b, where we compare the MPEs between using DSR-only, DSR-GP, and VIS-DSR-GP.

Figure 6: Data Noise Example and Exploration Performance by Interactions

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00251v1-figure-4b8d609b-c031-4b01-8f31-7b612e00ab60.jpg)

Table 2: Best Recovered Expressions with Noise Levels

(

| Ground truth   | min v f ( t ) +2 . 6 , v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0   | min v f ( t ) +2 . 6 , v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0   | min v f ( t ) +2 . 6 , v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0   | min v f ( t ) +2 . 6 , v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0   |
|----------------|---------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------|
| Noise level    | VIS-DSR-GP                                                                                        |                                                                                                   | DSR-GP                                                                                            | DSR                                                                                               |
| 0%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          | v l ( t ) + 9 . 0 v f ( t )+ v l ( t )+9 . 0                                                      |
| 1%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          | v l ( t ) + 1 9 . 0 ( s f ( t ) - v l ( t ) )                                                     |
| 2%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          | v l ( t ) + ( 2 s f ( t )+ v l ( t ) ) 9 . 0                                                      |
| 3%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          | v l ( t ) + v f ( t ) v f ( t )+2 . 6 - 2 . 6                                                     |
| 4%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          | v l ( t ) + s des - 9 . 0 v f ( t )+9 . 0                                                         |
| 5%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          | v l ( t ) + s des 9 . 0 +0 . 3                                                                    |
| 6%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          | v l ( t ) + s des 9 . 0 +0 . 3                                                                    |
| 7%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          | v l ( t ) + s des 9 . 0 - 1 . 0                                                                   |
| 8%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+6 . 4                          | v l ( t ) - 2 . 6                                                                                 |
| 9%             | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          | v l ( t )                                                                                      | - s f ( t )+9 . 0 s f ( t ) ( s f ( t )+9 . 0 ) + v f ( t )+ v l ( t )+9 . 0                   | v l ( t ) - 2 . 6                                                                                 |
| 10%            | v l ( t ) + 9 . 0 ( s f ( t ) - s des ( t ) ) v f ( t )+ v l ( t )+9 . 0                          |                                                                                                   | ( t ) ( 1 - 9 . 0 v f ( t )+ v l ( t )+9 . 0 s f ( t ) )                                          | v l ( t ) - 2 . 6                                                                                 |

This figure shows MPEs (with respect to the ground truth) over different noise levels, where each result is the average of running 10 random seeds with error bars (vertical lines) indicating standard errors. The results demonstrate that DSR-only fails to discover a suitable expression that can adequately capture the behaviors performed in the data for all cases, even on clean data (0% noise level). The average MPEs are at around 3% given the noise levels below 7%. DSR-GP can reconstruct the target expression for all random seeds on clean data. The averaged MPEs drop below 2% for noise levels from 0% to 9%, while MPEs range 2-5% when the dataset is imposed 10% of noise. The results clearly show that GP helps reduce the impact of data noise. On the other hand, with prior knowledge of potential variable combinations, VIS-DSR-GP can uncover the target expression for all tests when noise levels are at or below 3%, and produce expressions with the highest MPEs below 2% (at 10% of noise level). The best expressions (out of the 10 random seeds) recovered using these methods are listed in Table 2, where we only document the component of the second term in the min function, as all the methods can recover the first term ( v f ( ) + 2 6 t . ). Here the expression is chosen as the best if it has the highest accuracy and/or the most similar expression structure as the target expression. The results show that DSR-only is not able to find the correct expression at all noise levels. It falls into a local optimum, where the expression structure of each term is limited. With the help of GP, DSR-GP can recover the target expression at or below 6%. On the other hand, VIS-DSR-GP can recover the target expression for all tests at least once at all noise levels. The results also show that learning the exact structure of the expression (a division of two terms in this case) and the correct values of the constant (e.g., 9.0) are challenging for DSR and DSR-GP.

## 4.5 Testing Performance on Other Car-following Models

In this subsection, we test the expression reconstruction capability of the proposed VIS-enhanced DSR framework (VIS-DSR-GP) on other car-following models. We choose the GM model [Chandler et al., 1958] and the Gazis-HermanRothery (GHR) model [Gazis et al., 1961] as the target expressions in the following experiments. The two expressions are respectively listed in Eq. (10) and Eq. (11). Here we simplify the sensitivity term c in the GM model as a constant.

$$v _ { f } ( t + \delta t ) = v _ { f } ( t ) + c [ v _ { l } ( t ) - v _ { f } ( t ) ].$$

$$v _ { f } ( t + \delta t ) = v _ { f } ( t ) + k _ { 1 } v _ { f } ^ { k _ { 2 } } ( t ) \frac { v _ { l } ( t - \delta t ) - v _ { f } ( t - \delta t ) } { s _ { f } ^ { k _ { 3 } } ( t - \delta t ) }.$$

Each car-following dataset is generated through the SUMO numerical simulation with the corresponding car-following model sharing the same vehicle parameters, including δt , a max, v max, b , and t react , as the Krauss model mentioned in Section 4.1. For constants in the tested car-following models, we use the recommended c = 0 368 . [Chandler et al., 1958] for the GM model, and k 1 = 1 2 . , k 2 = 1 0 . , k 3 = 1 1 . for the GHR model [Al-Jameel et al., 2009, May and Keller, 1967]. While collecting the GHR car-following data, we simplify the velocity difference as a term, i.e.,

∆ ( v t -δt ) = v l ( t -δt ) -v f ( t -δt ) . The simplified car-following models with pre-determined parameters are as Eq. (12) and Eq. (13) for GM and GHR, respectively.

$$v _ { f } ( t + 1 ) = v _ { f } ( t ) + 0. 3 6 8 \left [ v _ { l } ( t ) - v _ { f } ( t ) \right ].$$

$$v _ { f } ( t + 1 ) = v _ { f } ( t ) + 1. 2 v _ { f } ( t ) \frac { \Delta v ( t - 1 ) } { s _ { f } ^ { 1. 1 } ( t - 1 ) }.$$

The collected GM dataset shares the same variables as the Krauss model dataset, i.e., features in the GM dataset include v f ( ) t , v l ( ) t , s f ( ) t , and ds t ( ) . The collected GHR dataset, on the other hand, comprises variables v f ( ) t , v l ( ) t , s f ( t -1) , and ∆ ( v t -1) , as it has more specific features. Feeding both datasets respectively in the VIS algorithm obtains interaction strengths as shown in Fig. 7. We apply the elbow method to determine recommended combinations and document them in Table 3.

Now, we use the VIS-enhanced DSR framework to reconstruct the target expressions. After being simplified, the GM and GHR expressions have one and two constants, respectively. We reconstruct the two models by giving the framework the following mathematical token pool: { + , - × ÷ , , , min , power, constant}, parameters { a max, v max, and b }, and the variables as shown in Table 3. The hard constraints include: i) min and constants cannot occur more than twice, and ii) min should be the first sampled token. Soft regulations, such as the reward penalty factor and recommended complexity, are configured the same as the values in the previous subsections.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00251v1-figure-1fc94089-fa55-4bbc-aa08-6aa183abf38a.jpg)

Figure 7: Interaction Strengths

| Table 3: Variable Interaction Strengths   | Table 3: Variable Interaction Strengths      | Table 3: Variable Interaction Strengths              |
|-------------------------------------------|----------------------------------------------|------------------------------------------------------|
|                                           | GM                                           | GHR                                                  |
| Interactions                              | v f ( t ) v l ( t ) ( v f ( t ) , v l ( t )) | v f ( t ) ( v f ( t ) , ∆ v ( t - 1) , s f ( t - 1)) |

By imposing noise from level 0% to 10%, we document the best-reconstructed expressions in Table 4, to which we document only the second term in each expression since the first term can be correctly learned by the method for all noise levels. From the table, VIS-DSR-GP is able to reconstruct the GM expression with the given token pool and

constraints for the zero-noise dataset. Although the decimal values of the explored constants between noise levels 1% to 7% do not exactly match the target value, they are quite close (less than 0.5%). However, the constant starts to differ significantly from the target value when the noise level exceeds 8%. For the GHR, on the other hand, the best expression the framework can explore is v f ( t +1) = v f ( ) + 1 01 t . v f ( ) t v l ( t -1) -v f ( t -1) s f ( t -1) -0 16 . . Although the resulting expression has the same structure as the GHR: v f ( ) t , plus the product of a constant and a fraction term, the constant and the power term (also a constant) in the denominator is the most challenging to learn and discover. This becomes even more challenging when the data has a higher level of noise. Starting at 5% noise level, an additional term appears beside the primary expression structure. When the noise level exceeds 9%, the structure of the explored expression begins to noticeably differ from the target expression of GHR. The results here further confirms that learning the exact structure of the model and the constants (especially those with decimal values) are challening even for the proposed VIS-DSR-GP method.

Table 4: Best Recovered Expressions Using VIS-DSR-GP

|             | GM                               | GHR                                                                             |
|-------------|----------------------------------|---------------------------------------------------------------------------------|
| Noise level | 0 . 368 ( v l ( t ) - v f ( t )) | 1 . 20 v f ( t ) ∆ v ( t - 1) s 1 . 1 f ( t - 1)                                |
| 0%          | 0 . 368 ( v l ( t ) - v f ( t )) | 1 . 01 v f ( t ) ∆ v ( t - 1) s f ( t - 1) - 0 . 16                             |
| 1%          | 0 . 367 ( v l ( t ) - v f ( t )) | 1 . 01 v f ( t ) ∆ v ( t - 1) s f ( t - 1) - 0 . 16                             |
| 2%          | 0 . 366 ( v l ( t ) - v f ( t )) | 1 . 01 v f ( t ) ∆ v ( t - 1) s f ( t - 1) - 0 . 16                             |
| 3%          | 0 . 366 ( v l ( t ) - v f ( t )) | 1 . 01 v f ( t ) ∆ v ( t - 1) s f ( t - 1) - 0 . 16                             |
| 4%          | 0 . 366 ( v l ( t ) - v f ( t )) | 1 . 00 v f ( t ) ∆ v ( t - 1) s f ( t - 1) - 0 . 16                             |
| 5%          | 0 . 366 ( v l ( t ) - v f ( t )) | 1 . 00 v f ( t ) ∆ v ( t - 1) s f ( t - 1) - 0 . 17 + 0 . 21 ∆ v ( t - 1)       |
| 6%          | 0 . 366 ( v l ( t ) - v f ( t )) | 1 . 00 ∆ v ( t - 1) s f ( t - 1) - 0 . 17 - 0 . 01 v f ( t )                    |
| 7%          | 0 . 366 ( v l ( t ) - v f ( t )) | 1 . 00 v f ( t ) ∆ v ( t - 1) s f ( t - 1) - 0 . 17 + 0 . 02 2∆ v ( t - 1)      |
| 8%          | 0 . 385 ( v l ( t ) - v f ( t )) | 1 . 00 v f ( t ) ∆ v ( t - 1) s f ( t - 1) - 0 . 17 - 0 . 78                    |
| 9%          | 0 . 385 ( v l ( t ) - v f ( t )) | v f ( t )∆ v ( t - 1)( - 0 . 06 + 1 . 27 s f ( t ) )                            |
| 10%         | 0 . 385 ( v l ( t ) - v f ( t )) | [ v f ( t )(∆ v ( t - 1)) ] 2 v ( t )∆ v ( t - 1) ( s ( t ) - 0 . 06 ) - 0 . 57 |

## 5 Discussions

The results of the above numerical experiments illustrate that under relatively simple scenarios - Krauss or GM model, single vehicle class, noise-free or low-noise level (less than 8% as shown in Table 2 and Table 3) - the proposed VIS-enhanced DSR method can discover the target car-following models exactly or almost exactly from data directly with minimum human involvement. We think this is quite impressive considering the various structures/terms, math operations, variables, and especially constants in the expressions of the car-following models. This is especially so when compared with traditional methods for developing car-following models, for which extensive, time-consuming, and iterative processes of data analysis, visualization, hypothesis, validation, and testing are essential steps. The results, however, also demonstrate that the proposed framework faces significant challenges in exploring expressions correctly for dynamics with complex expression structures, constants, and data noises. Target expressions with intricate structures, such as the power term in the GHR model, can make it more difficult for the proposed method to find the correct expression. When there are many constants to be fitted, finding the exact expressions can also be challenging. Another challenge is data noises, which can have a significant impact on the exploration process, as they can obscure the desired predictions and thus make it difficult to obtain the target expression terms correctly. Furthermore, real-world traffic streams often contain multiple vehicle classes, such as cars, buses, and trucks, which may have different key model parameters (maximum speed, acceleration, etc.) even if they follow the same car-following dynamics. Our initial testing shows that this can bring great challenge to the proposed learning method, due to the large number of constants and the multimodal nature of the problem. As a result, there is much room for improvements in VIS-enhanced DSR to address these challenges.

More advanced deep learning techniques, such as large language model (LLM) based methods, need to be explored to address these challenges; see Section 6 for more discussions on this aspect. On the other hand, clearly defining objectives and developing measures that can properly capture such objectives are also crucial. With the goal of learning the governing equations of dynamical systems, such as the car-following dynamics in this paper, the learned mathematical expressions should accurately capture the observed data, which should also have proper physical meanings and can be generalized to different sites and traffic conditions. Learning accuracy is not the only objective, which is starkly different from other DL-based applications. This prompts the definition of four key measures that help evaluate the quality of

the derived expressions: i) accuracy, ii) interpretability, iii) parsimony, and iv) generalizability, where interpretability, parsimony (compactness), and generalizability are also discussed in [Tenachi et al., 2023].

Accuracy measures the 'goodness of fit' of the learned expressions with respect to the observed data, which is the basic measure to evaluate the quality of the learned expression. A good expression should provide reasonable accuracy in fitting the observed data. Metrics such as MSE, NMSR, or MPE have been widely used to measure accuracy; for example, DSR uses NMSE as one measure for the performance of the resulting expression. In addition, other mathematical properties of the model, such as stability, are also studied and used as metrics to evaluate the performance of the models.

Interpretability means that a learned expression should have human-understandable justifications about the 'physical meaning' of the model. It is important as it implies 'trust' - people often have more confidence in an interpretable model that can then be more easily adopted and used [Murdoch et al., 2019, Arrieta et al., 2020, Madiega, 2021]. Interpretability is an important objective for learning governing equations, as the main goal is to add new knowledge to our understanding of the physical systems by developing such models. This can be particularly essential (and challenging) when we learn the potentially new dynamics for mixed traffic flow of CA Vs and HDVs, as discussed in the Conclusion Section below. As we expect that explicit analytical expressions are produced in this research, we focus on two aspects of interpretability: i) the variables in the model represent clear physical phenomena or characteristics of the traffic system; and ii) the relations of the variables described in a model have clear, concise (see parsimony below) explanations. In this paper, we applied VIS to find the key (strong) interactions among variables, with the help of constraining the model complexity, which helps produce expressions that have reasonable physical meanings. As interpretability and the above two aspects are more subjective, it is imperative to develop specific measures for interpretability that can be properly integrated into the DSR learning process (e.g., the design of the reward function).

Parsimony is 'a good principle to explain the phenomena by the simplest hypotheses possible' [Gauch Jr, 2002, Lycan, 1975] or the simplest possible solution to a problem, which has been widely used as a guiding principle for daily life and scientific investigations [Lycan, 1975, Bentley and McIlroy, 1993]. In our specific context, parsimony means that the learned expressions are concise and straightforward, given all others, e.g., the metrics for other objectives being equal or very similar. In this paper, we use the complexity of the expression, i.e., the number of variables and operations included in the expression, as the measure to define and evaluate its parsimony. However, parsimony is also a subjective measure requiring substantial human judgment of the learned expression. Specific procedures/measures need to be developed to help evaluate the parsimony of an expression so that it is not too complex or over-parsimonious (i.e. too simple), and also help reach a proper balance with other objectives.

Generalizability means that the learned governing equations can be applied to data collected from broader traffic systems at various sites with different traffic mixes and conditions, which indicates that the learned expression can capture the underlying dynamics of the traffic system rather than only fitting to a specific dataset. Generalizability can be measured by the model performance (e.g., accuracy) across different datasets and different sites. An expression that can describe all scenarios better is considered a more generalizable model than those that only work well for certain datasets but not on others. This may be measured by the metrics such as their averages and variances of other three objectives discussed above across different datasets and test sites. Furthermore, expressions with no or minimum calibration but can well describe the dynamics of different datasets and/or sites also have better generalizability properties.

## 6 Conclusion

In this study, we proposed an expression exploration framework based on the DSR approach while encouraging parsimonious mathematical expressions that include potential variable combinations detected through a VIS algorithm. The variable combinations are used to guide DSR in focusing on expressions that involve the identified combinations, which can also increase the interpretability of the explored mathematical expressions. For this, we designed two penalty terms for the reward function of training RNNs in DSR: (i) a complexity penalty to encourage the exploration to crawl toward the desired range of expression complexity, and (ii) a variable interaction penalty to encourage the exploration of expressions consisting of the recommended variable combinations. The following results were discovered. First, the recommended variable combinations were found useful for directing the exploration to identify suitable mathematical expressions. Secondly, the proposed complexity penalty strategy was shown to be better for regulating the complexity range searched than the traditional complexity regulation method in DSR. Thirdly, the exploration based on the recommended variable combinations with the complexity penalty strategy helped identify the governing equations more efficiently than the traditional DSR. From the data noise investigation experiments, we recommended turning on the exploration assistance of GP, as it enables a wider search space for exploration.

The methods and results in this paper present the first step of discovering governing equations of car-following dynamics using deep learning. Limitations exist which motivate future research. We need to improve the current learning

framework to address the limitations summarized in Section 5, and apply the improved framework to more diversified, real trajectory data. In particular, the proposed framework can be applied to multimodal traffic systems, e.g., to capture the car-following dynamics among CAVs and between HDVs and CAVs, and to help develop interpretable dynamical models for the interactions between vehicles and vulnerable road users. This may be done by integrating multiple dynamic exploration frameworks as a coordinated learning system in which each exploration framework learns an interaction type in a coordinated way from the collected multimodal data.

There are two main directions that can help improve the current learning framework and address the challenges. First, as discussed in the previous section, accurately fitting the observed data is NOT the only objective for deriving proper governing equations. One has to balance with other objectives, including interpretability, parsimony, and generalizability. Proper metrics for each of these objectives need to be defined and integrated into the learning process, e.g., to design a comprehensive reward function that includes metrics for measuring parsimony, interpretability, generalization, and accuracy. An exploration incorporated with a reward that takes these features into account may also help find appropriate expressions from data with noise and complicated underlying dynamics. However, developing such a reward function is challenging, as there may be conflicts between accuracy and other objectives. Despite the challenges, this is an important research topic that has the potential to advance the capability of using deep learning models to "automatically" discover the underlying expressions from collected data. Second, given the advance of deep learning especially recent LLMs, one should leverage the ability of interaction and memory of LLMs to enhance search efficiency and expression quality in discovering the governing equations. Research has started in this direction Kamienny et al. [2022], which needs to be further improved and tailored to car-following and other engineering applications.

## References

Mark Brackstone and Mike McDonald. Car-following: a historical review. Transportation Research Part F: Traffic Psychology and Behaviour , 2(4):181-196, 1999.

Louis A Pipes. An operational analysis of traffic dynamics. Journal of applied physics , 24(3):274-281, 1953.

Robert E Chandler, Robert Herman, and Elliott W Montroll. Traffic dynamics: studies in car following. Operations research , 6(2):165-184, 1958.

EIJI Kometani and TSUNA Sasaki. On the stability of traffic flow (report-i). Journal of the Operations Research Society of Japan , 2(1):11-26, 1958.

Denos C Gazis, Robert Herman, and Renfrey B Potts. Car-following theory of steady-state traffic flow. Operations research , 7(4):499-505, 1959.

Denos C Gazis, Robert Herman, and Richard W Rothery. Nonlinear follow-the-leader models of traffic flow. Operations research , 9(4):545-567, 1961.

Walter Helly. Simulation of bottlenecks in single-lane traffic flow. Theory of traffic flow , 1959.

Partha Chakroborty and Shinya Kikuchi. Evaluation of the general motors based car-following models and a proposed fuzzy inference model. Transportation Research Part C: Emerging Technologies , 7(4):209-235, 1999.

Martin Treiber and Dirk Helbing. Reconstructing the spatio-temporal traffic dynamics from stationary detector data. Cooper@ tive Tr@ nsport@ tion Dyn@ mics , 1(3):3-1, 2002.

Zhu Wen-xing and Jia Lei. Analysis of a new car-following model with a consideration of multi-interaction of preceding vehicles. In 2007 IEEE Intelligent Transportation Systems Conference , pages 320-324. IEEE, 2007.

Mitra Pourabdollah, Eric Bjärkvik, Florian Fürer, Björn Lindenberg, and Klaas Burgdorf. Calibration and evaluation of car following models using real-world driving data. In 2017 IEEE 20th International conference on intelligent transportation systems (ITSC) , pages 1-6. IEEE, 2017.

Valentina Kurtc. Studying car-following dynamics on the basis of the highd dataset. Transportation research record , 2674(8):813-822, 2020.

Boyuan Chen, Kuang Huang, Sunand Raghupathi, Ishaan Chandratreya, Qiang Du, and Hod Lipson. Discovering state variables hidden in experimental data. arXiv preprint arXiv:2112.10755 , 2021.

Pierre-Alexandre Kamienny, Stéphane d'Ascoli, Guillaume Lample, and François Charton. End-to-end symbolic regression with transformers. arXiv preprint arXiv:2204.10532 , 2022.

Rongye Shi, Zhaobin Mo, Kuang Huang, Xuan Di, and Qiang Du. Physics-informed deep learning for traffic state estimation. arXiv preprint arXiv:2101.06580 , 2021.

Zhaobin Mo, Rongye Shi, and Xuan Di. A physics-informed deep learning paradigm for car-following models. Transportation research part C: emerging technologies , 130:103240, 2021.

| Zijia Zhong, Earl E Lee, Mark Nejad, and Joyoung Lee. Influence of cav clustering strategies on mixed traffic flow characteristics: An analysis of vehicle trajectory data. Transportation Research Part C: Emerging Technologies , 115: 102611, 2020.                  |
|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Lin Xiao, Meng Wang, Wouter Schakel, and Bart van Arem. Unravelling effects of cooperative adaptive cruise control deactivation on traffic flow characteristics at merging bottlenecks. Transportation research part C: emerging technologies , 96:380-397, 2018.       |
| Derek Gloudemans, Yanbing Wang, Junyi Ji, Gergely Zachar, William Barbour, Eric Hall, Meredith Cebelak, Lee Smith, and Daniel B Work. I-24 motion: An instrument for freeway traffic science. Transportation Research Part C: Emerging Technologies , 155:104311, 2023. |
| Samuel H Rudy, Steven L Brunton, Joshua L Proctor, and J Nathan Kutz. Data-driven discovery of partial differential equations. Science advances , 3(4):e1602614, 2017.                                                                                                  |
| Steven L Brunton, Joshua L Proctor, and J Nathan Kutz. Discovering governing equations from data by sparse identification of nonlinear dynamical systems. Proceedings of the national academy of sciences , 113(15):3932-3937, 2016.                                    |
| Constantinos Antoniou, Ramachandran Balakrishna, and Haris NKoutsopoulos. Asynthesis of emerging data collection technologies and their impact on traffic management applications. European Transport Research Review , 3(3): 139-148, 2011.                            |
| Josh Bongard and Hod Lipson. Automated reverse engineering of nonlinear dynamical systems. Proceedings of the National Academy of Sciences , 104(24):9943-9948, 2007.                                                                                                   |
| Michael Schmidt and Hod Lipson. Distilling free-form natural laws from experimental data. science , 324(5923):81-85, 2009.                                                                                                                                              |
| John R Koza. Genetic programming as a means for programming computers by natural selection. Statistics and computing , 4(2):87-112, 1994.                                                                                                                               |
| Thomas Bäck, David BFogel, and Zbigniew Michalewicz. Evolutionary computation 1: Basic algorithms and operators CRC press, 2018.                                                                                                                                        |
| William La Cava, Lee Spector, and Kourosh Danai. Epsilon-lexicase selection for regression. In Proceedings of the Genetic and Evolutionary Computation Conference 2016 , pages 741-748, 2016.                                                                           |
| William La Cava, Tilak Raj Singh, James Taggart, Srinivas Suri, and Jason H Moore. Learning concise representations for regression by evolving networks of trees. arXiv preprint arXiv:1807.00981 , 2018.                                                               |
| Marco Virgolin, Tanja Alderliesten, and Peter AN Bosman. Linear scaling with and within semantic backpropagation- based genetic programming for symbolic regression. In Proceedings of the genetic and evolutionary computation conference , pages 1084-1092, 2019.     |
| Marco Virgolin, Tanja Alderliesten, Cees Witteveen, and Peter AN Bosman. Improving model-based genetic program- ming for symbolic regression of small expressions. Evolutionary computation , 29(2):211-237, 2021.                                                      |
| Michael Kommenda, Bogdan Burlacu, Gabriel Kronberger, and Michael Affenzeller. Parameter identification for symbolic regression using nonlinear least squares. Genetic Programming and Evolvable Machines , 21(3):471-501, 2020.                                        |
| Fabricio Olivetti de Franca and Guilherme Seidyo Imai Aldeia. Interaction-transformation evolutionary algorithm for symbolic regression. Evolutionary computation , 29(3):367-390, 2021.                                                                                |
| Brenden K Petersen, Mikel Landajuela Larma, T Nathan Mundhenk, Claudio P Santiago, Soo K Kim, and Joanne T Kim. Deep symbolic regression: Recovering mathematical expressions from data via risk-seeking policy gradients. arXiv preprint arXiv:1912.04871 , 2019.      |
| James Enouen and Yan Liu. Sparse interaction additive networks via feature interaction detection and sparse selection. Advances in Neural Information Processing Systems , 35:13908-13920, 2022.                                                                        |
| Louis Albert Pipes. Car following models and the fundamental diagram of road traffic. Transportation Research/UK/ 1966.                                                                                                                                                 |
| Robert Herman, Elliott WMontroll, Renfrey B Potts, and Richard WRothery. Traffic dynamics: analysis of stability in car following. Operations research , 7(1):86-106, 1959.                                                                                             |
| Leslie C Edie. Car-following and steady-state theory for noncongested traffic. Operations research , 9(1):66-76, 1961.                                                                                                                                                  |
| Reiter Wiedemann and U Reiter. Microscopic traffic simulation: the simulation system mission, background and actual state. Project ICARUS (V1052) Final Report , 2:1-53, 1992.                                                                                          |

| RMMichaels. Perceptual factors in car-following. Proc. of 2nd ISTTF (London) , pages 44-59, 1963.                                                                                                                                                                              |
|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Gerald O Burnham and George A Bekey. A heuristic finite-state model of the human driver in a car-following situation. IEEE Transactions on Systems, Man, and Cybernetics , (8):554-562, 1976.                                                                                  |
| David N Lee. A theory of visual control of braking based on information about time-to-collision. Perception , 5(4): 437-459, 1976.                                                                                                                                             |
| H Kumamoto, K Nishi, K Tenmoku, and H Shimoura. Rule-based cognitive animation simulator for current lane and lane change drivers. Steps Forward. Intelligent Transport Systems World Congress , 1995.                                                                         |
| Stefan Krauß. Microscopic modeling of traffic flow: Investigation of collision free vehicle dynamics, 1998.                                                                                                                                                                    |
| Martin Treiber, Ansgar Hennecke, and Dirk Helbing. Congested traffic states in empirical observations and microscopic simulations. Physical review E , 62(2):1805, 2000.                                                                                                       |
| Venkatesan Kanagaraj, Gowri Asaithambi, CH Naveen Kumar, Karthik K Srinivasan, and R Sivanandan. Evaluation of different vehicle following models under mixed traffic conditions. Procedia-Social and Behavioral Sciences , 104: 390-401, 2013.                                |
| CH Naveen Kumar. Microscopic analysis of vehicle behaviour in mixed traffic. M. Tech Project Report , 2013.                                                                                                                                                                    |
| Martin Treiber and Arne Kesting. Traffic flow dynamics . Springer, 2013.                                                                                                                                                                                                       |
| Michael O'Neill, Leonardo Vanneschi, Steven Gustafson, and Wolfgang Banzhaf. Open issues in genetic programming. Genetic Programming and Evolvable Machines , 11:339-363, 2010.                                                                                                |
| Kathleen Champion, Bethany Lusch, J Nathan Kutz, and Steven L Brunton. Data-driven discovery of coordinates and governing equations. Proceedings of the National Academy of Sciences , 116(45):22445-22451, 2019.                                                              |
| L Gao and J Nathan Kutz. Bayesian autoencoders for data-driven discovery of coordinates, governing equations and fundamental constants. arXiv preprint arXiv:2211.10575 , 2022.                                                                                                |
| Ying Jin, Weilin Fu, Jian Kang, Jiadong Guo, and Jian Guo. Bayesian symbolic regression. arXiv preprint arXiv:1910.08892 , 2019.                                                                                                                                               |
| Silviu-Marian Udrescu, Andrew Tan, Jiahai Feng, Orisvaldo Neto, Tailin Wu, and Max Tegmark. Ai feynman 2.0: Pareto-optimal symbolic regression exploiting graph modularity. Advances in Neural Information Processing Systems 33:4860-4871, 2020.                              |
| Miles D. Cranmer, Alvaro Sanchez-Gonzalez, Peter W. Battaglia, Rui Xu, Kyle Cranmer, David N. Spergel, and Shirley Ho. Discovering symbolic models from deep learning with inductive biases. CoRR , abs/2006.11287, 2020. URL https://arxiv.org/abs/2006.11287 .               |
| Aaron Meurer, Christopher P Smith, Mateusz Paprocki, Ondˇej r ˇ Certík, Sergey B Kirpichev, Matthew Rocklin, AMiT Kumar, Sergiu Ivanov, Jason K Moore, Sartaj Singh, et al. Sympy: symbolic computing in python. PeerJ Computer Science , 3:e103, 2017.                        |
| Mikel Landajuela Larma, Brenden K Petersen, Soo K Kim, Claudio P Santiago, Ruben Glatt, T Nathan Mundhenk, Jacob F Pettit, and Daniel MFaissol. Improving exploration in policy gradient search: Application to symbolic optimization. arXiv preprint arXiv:2107.09158 , 2021. |
| Michael Tsang, Sirisha Rambhatla, and Yan Liu. How does this interaction affect me? interpretable attribution for feature interactions. Advances in neural information processing systems , 33:6147-6159, 2020.                                                                |
| Daniel Krajzewicz, Georg Hertkorn, Christian Rössel, and Peter Wagner. Sumo (simulation of urban mobility)-an open-source traffic simulation. In Proceedings of the 4th middle East Symposium on Simulation and Modelling (MESM20002) , pages 183-187, 2002.                   |
| Martin Treiber, Arne Kesting, Martin Treiber, and Arne Kesting. Car-following models based on driving strategies. Traffic Flow Dynamics: Data, Models and Simulation , pages 181-204, 2013.                                                                                    |
| Hamid Al-Jameel et al. Examining and improving the limitations of the gazis-herman-rothery car-following model. In Salford Postgraduate Annual Research Conference , volume 26, pages 3-4. Citeseer, 2009.                                                                     |
| Adolf D May and Hartmut EMKeller. Non-integer car-following models. Highway Research Record ,                                                                                                                                                                                  |
| 199(1):19-32, 1967.                                                                                                                                                                                                                                                            |
| WJames Murdoch, Chandan Singh, Karl Kumbier, Reza Abbasi-Asl, and Bin Yu. Definitions, methods, and applications                                                                                                                                                               |
| in interpretable machine learning. Proceedings of the National Academy of Sciences , 116(44):22071-22080, 2019.                                                                                                                                                                |

Alejandro Barredo Arrieta, Natalia Díaz-Rodríguez, Javier Del Ser, Adrien Bennetot, Siham Tabik, Alberto Barbado, Salvador García, Sergio Gil-López, Daniel Molina, Richard Benjamins, et al. Explainable artificial intelligence (xai): Concepts, taxonomies, opportunities and challenges toward responsible ai. Information fusion , 58:82-115, 2020.

Tambiama Madiega. Artificial intelligence act. European Parliament: European Parliamentary Research Service , 2021.

Hugh G. Gauch Jr. PARSIMONY AND EFFICIENCY , page 269-326. Cambridge University Press, 2002.

William G Lycan. Occam's razor. Metaphilosophy , 6(3/4):223-237, 1975.

Jon L Bentley and M Douglas McIlroy. Engineering a sort function. Software: Practice and Experience , 23(11): 1249-1265, 1993.
Ohay Angah James Enouen Xuegang Ban Yan Liu
Discovering Car-following Dynamics from Trajectory Data through Deep Learning
This study aims to discover the governing mathematical expressions of
car-following dynamics from trajectory data directly using deep learning
techniques. We propose an expression exploration framework based on deep
symbolic regression (DSR) integrated with a variable intersection selection
(VIS) method to find variable combinations that encourage interpretable and
parsimonious mathematical expressions. In the exploration learning process, two
penalty terms are added to improve the reward function: (i) a complexity
penalty to regulate the complexity of the explored expressions to be
parsimonious, and (ii) a variable interaction penalty to encourage the
expression exploration to focus on variable combinations that can best describe
the data. We show the performance of the proposed method to learn several
car-following dynamics models and discuss its limitations and future research
directions.