Simulation Research on Safe Flow Rate of Bidirectional Crowds Using Bayesian-Nash Equilibrium

Current research on pedestrian flows has mainly focused on evacuation optimization during or after emergencies; however, crowd management before emergencies has received little attention.This paper examines the management of a Safe Pedestrian Flow Rate, in which the Bayesian-Nash Equilibriummimics pedestrians’ decision-making, and a multiagent system is employed to reproduce pedestrians’ interactions. In the model, the pedestrian tunnel is divided into cells, with each pedestrian in a cell receiving a utility depending on the distance to the exit and the number of pedestrians in the cell. Then, each pedestrian uses the Bayesian-Nash Equilibrium to search for the target cell with maximum expected utility, moves in, andmakes next decision until exiting the tunnel. The simulation model is calibrated and validated from a real scenario. Finally, from the experimental data collected from different simulation scenarios, this research reaches the conclusion that the Safe Pedestrian Flow Rate increases by about 2.96ped/s as the tunnel width expanded by 1m.This paper offers a novel method for reducing potential losses caused by crowd emergencies and can be a valuable reference for managing pedestrian flows and designing public places.


Introduction
As urbanization increases, urban populations worldwide have grown exponentially. As a result, stampede accidents such as the Shanghai New Year's Eve stampede in 2015, the Phnom Penh stampede in 2010, the Beijing Lantern Festival stampede in 2004, and the Mecca Hajj stampedes have become more frequent. According to incomplete statistics, over 100 stampede accidents have occurred since 2001, causing more than 10,000 deaths [1]. Therefore, research on pedestrian flows, especially in dense situations, has become a more important research focus.
However, as current crowd research mainly focused on evacuation optimization after an emergency, the management of pedestrian flows before an emergency received little attention. From the two photographs in Figure 1, it can be seen that as roads or venues are always crowded with pilgrims or pedestrians moving in opposite directions during ceremonies and holidays, accidents or near accidents can easily occur. Therefore, studies on managing pedestrian flows to prevent or control stampedes or other emergencies are required.
In this paper, a model is developed to mimic pedestrian crowd decision-making and determine a Safe Pedestrian Flow Rate (SPFR) for different situations. The remainder of this paper is organized as follows. Section 2 gives a brief review of previous research on pedestrian flows. Section 3 introduces the proposed simulation model in detail, including the decision-making processes and the collision avoidance strategies. Section 4 calibrates and validates the model with data collected from an actual scene. Section 5 conducts simulation experiments that consider pedestrians' walking preferences under different circumstances and reveals the relationships between the SPFR and the tunnel width. Section 6 concludes the work and discusses prospects for future research.  and after emergency. This paper classifies current emergency research into two pre-emergency research and postemergency research, as shown in Figure 2. Generally, preemergency research focuses on the work before emergencies happen, such as crowd management, risk analysis, and public space design, while post-emergency research keeps an eye on the work during or after emergencies happen, such as evacuation planning and optimizing and loss prevention. Present pedestrian flow models can be divided into macroscopic and microscopic models. Macroscopic models examine overall pedestrian flows and research pedestrian characteristics using qualitative theories or methods from other disciplines, such as the Fluid Dynamics Model proposed by Huges [2,3] and Wave Theory introduced by Lu [4]. Microscopic models, on the other hand, such as the Lattice Gas Model [5], the Social Force Model [6,7], the Cellular Automata Model [8,9], and the Agent-based Model [10,11], focus on the influence of individual behavior on pedestrian flows. As macroscopic models study overall pedestrian flows and ignore the heterogeneity and complexity of individuals, they are unable to accurately analyze the mechanisms between the microscopic pedestrian behavior and the macroscopic phenomena of the crowd. Therefore, microscopic models have been the main methods used for pedestrian flow research [12][13][14][15]. Table 1 summarizes recent research on pedestrian flows.
From Table 1, it can be seen that nearly 75% of recent research has been focused on post-emergency, such as reproducing evacuation scenes and optimizing evacuations, with very few examining situations before the crises occurred. What is more, as most pre-emergency research has focused on the reproduction of pedestrian flow phenomena and the optimization of emergency plan, crowd management to avoid accidents has received very little research attention. Further, while many studies have examined the factors that influence egress efficiency, neither the quantitative nor concrete relationships between these factors have been explored. The common characteristic in current research has been that almost all authors have employed a methodology that combined at least two modeling methods, with one being a decision or behavioral model that mimics pedestrians' decision-making or behaviors and the other being an action model that motivates pedestrians' movement.
The main purpose of this research, therefore, is to determine a scientific, practical method for managing crowds at a safe level to prevent possible emergencies. Game theory has been found to be able to properly capture a pedestrian's decision-making processes and interactions, and agent-based  modeling methods have been proven to effectively mimic pedestrian movements [13]. This paper, therefore, adopts a multiagent system and the Bayesian-Nash Equilibrium to formulate the simulation model. In the model, one or more pedestrians are permitted to stay in a cell, with each pedestrian assigned a utility that depends on their distance to the exit and the number of other pedestrians sharing the cell. The Bayesian-Nash Equilibrium is an equilibrium achieved in incomplete information games in which each player chooses a strategy to gain a maximum expected utility based on the probability distributions of the others' decisions. The simulation experiment in this model is selforganized once it begins to run, as the pedestrian does not know the exact choices their neighbors are going to make.
They choose a target cell to gain as much utility as possible based on the probability distributions of the neighboring pedestrians' choices. Hughes pointed out that emergencies can easily occur when pedestrian density is greater than 4ped/m 2 [2]; therefore, this paper regards the density of 4ped/m 2 to be the Critical Emergency Density and defines the 4 Complexity  maximal pedestrian flow rate before the critical density as the Safe Pedestrian Flow Rate (SPFR). In other words, if timely measures are not taken when the pedestrian flow exceeds the SPFR, it is more likely that an emergency could occur. To be clear, the main pedestrian flow measures in this paper are as shown in Table 2.

Simulation Model Framework
. . Initial Model Environment. This paper studies pedestrian flows using Agent-Based Modeling and Simulation (ABMS). Pedestrians enter the pedestrian tunnel from both ends and walk in two opposite directions: left to right and right to left; as shown in Figure 3, and as with most models, the tunnel is divided into cells. In most existing microscopic models, the evacuation venues are partitioned into square (diamond, or hexagon) cells with 0.4m sides, which are based on pedestrian projections. In other words, these models cannot simulate pedestrians density more than 1ped/(0.4m * 0.4m), which is about 6ped/m 2 . What is more, when the pedestrian moves from one cell to another, the distance is not in accordance with the actual average pedestrian walking step length. When the pedestrian moves more than one cell in each time step, their decision-making processes become rather complex as many factors need to be considered, such as decision-making efficiency and the target cell's availability.
Based on experimental data and statistical analyses, Zeng determined that in a normal condition a pedestrian walks with a step length of 0.7m and at a speed of 1.4m/s [40]. This paper, therefore, sets the side length of the cell at 0.7m based on the average pedestrian step length, with pedestrians making decisions every 0.5s. In other words, at each second the pedestrian moves twice and moves no more than one cell each time. The probability that a pedestrian chooses each target cell is shown in Figure 3.
In Figure 3, Pedestrian A and Pedestrian B are pedestrians walking from left to right and from right to left, 1 , 2 , 3 , 4 , and 5 are the successive probabilities of the adjacent cells being chosen by the pedestrian from their left side to right side, and 0 is the probability that a pedestrian remains stationary.
. . Utility Function. The utility that a pedestrian gains while walking depends on their desire to achieve their destination [35]. In this research, the utility is divided into movement utility V and comfort utility , as follows.
The movement utility depends on the distance the pedestrian moves toward to the exit and is defined as (1) The Numbering Rule: Cell (r, c) means the cell in Row r and Column c (2) The decision-maker and ambient pedestrians in the game where is the distance the pedestrian moves each time and is the angle between the movement direction and the exit direction. For example, when a pedestrian walks straight to the exit, equals 0 and when they walk in a vertical direction to the exit, equals 90.
A relationship has been found between the pedestrian's speed and crowd density; however, few studies have developed a formula for this relationship. Zeng summarized former studies and provided a function for pedestrian velocity: pedestrians are able to move freely if the density is below 4ped/m 2 , have severely restricted movements when the density exceeds 7ped/m 2 , and when the density is between 4ped/m 2 and 7ped/m 2 , the pedestrian has restricted velocity, which is related to the density [40]. Based on the function, this paper defines 1.00 as the pedestrian comfort utility when they can move freely, with the ratio of restricted speed to free speed being their comfort utility value in other cases, as follows.
where denotes the number of pedestrians in the same cell.
. . Expected Comfort Utility. Pedestrians in the model do not know the other pedestrians' decisions, which could affect their own comfort utility, and therefore make independent decisions when walking. That is, before moving, each pedestrian has to calculate the expected utility of every possible decision based on the probability distribution of the other pedestrians' choices, after which they choose the cell with the maximum utility as the target cell. In this paper, a pedestrian is assumed to choose the target cell from six possible cells, as shown in Figure 3, and each of the six alternative cells has eight neighboring cells, which means that many pedestrians may choose the same target cell, as shown in Figure 4.
Subgraph (1) in Figure 4 shows the cell numbering rule. Specifically, cell (r, c) denotes the cell in Row r Column c. Subgraph (2) illustrates the game interactions between the decision-maker and their neighboring pedestrians. The circles represent the pedestrians who may potentially choose Cell (r, c). Further, blue circles indicate pedestrians walking from left to right, and green ones denote pedestrians walking from right to left. It can be seen in subgraph (2) that the ambient pedestrians impact the expected comfort utility of the decision-maker in Cell (r, c-1) who chooses Cell (r, c) as their target cell. This means the decision-maker must consider the choices made by its surrounding pedestrians. Expressions in each circle represent the possibilities of pedestrians moving into the Cell (r, c) separately, according to Figure 3.
Take the decision-maker in Figure 4 as an example. As in the simulation model, a pedestrian should take all the situations into consideration before making a decision. Thus, a pedestrian's expected utility of the target cell can be divided into four parts, which are utility of one pedestrian, of two pedestrians, of three pedestrians, and of more than three pedestrians. The expected comfort utility is where (0), (1), (2) represent the possibility of no pedestrians, one pedestrian, or two pedestrians in Cell (r, c), respectively.
The probability that no pedestrians choose Cell (r, c) is and the probability that only one pedestrian chooses Cell (r, c) is When two pedestrians choose the target cell, they may come from the same group or two different groups, in which case, it could be calculated separately.
(2) = (2) + (2) If the two pedestrians come from the same group, and if they come from two different groups, When a pedestrian chooses other cells or is at other positions, the method used to calculate their expected comfort utility is the same. Further, to keep pedestrians in the tunnel, the expected comfort utility is defined as negative infinity when the target cell is beyond the upper or the lower boundaries.
. . Collision Avoidance Action. In most microscopic models, conflicts occur when pedestrians make decisions or move at the same time, as each cell can contain only one pedestrian. If conflicts are not properly dealt with in the model, pedestrians will go directly through other pedestrians or some other uncontrollable phenomena will occur, which inevitably leads to a false conclusion. This paper allows two or more pedestrians to share a cell and as a result avoids conflicts caused by pedestrians making simultaneous decisions; however, when  two or more pedestrians move at the same time, once their routes intersected, these collisions are still unavoidable in the model, as shown in Figure 5.
As we all know, collisions seldom happen even in extremely dense crowds because pedestrians can generally perceive possible collisions before happening and immediately change their routes. This paper, therefore, makes a rule that if a collision is likely to happen, each of the involved pedestrians changes their decisions to an adjacent alternative cell with a chance of V . Based on the above collision definition, it can be concluded that collisions can only occur between pedestrians moving in opposite directions, as shown in Figure 5. For instance, a pedestrian in Cell (r, c-1) walking from left to right predicts a collision only when they choose Cell (r-1, c) or Cell (r+1, c) and if there are pedestrians in Cell (r, c) walking from right to left. It is important to note that even if the pedestrian in Cell (r, c-1) does not change their decision, there is still a chance that a collision does not occur as pedestrians in Cell (r, c) may not choose the same cell. From this point of view, the chance that a pedestrian changes their minds to avoid a collision should be approximated to the possibility that a collision may occur. In this research, V is defined as follows: where , , 2 , and 4 are in accordance with previous sections.
. . Decision-Making Process. The detailed process for the pedestrians' decision-making is shown in Figure 6. In each time step, a pedestrian calculates its utilities of all available choices. After that, it chooses the cell with maximal expected utility as its target cell. Then it takes an action to avoid possible collisions. Finally, it moves into the chosen cell and begins the next decision-making round until passing through the tunnel. Before simulation, the length and the width of the tunnel, the pedestrian flow rate, and the probabilities from 0 to 5 are put into the model to initialize the experimental environment. During simulation, each pedestrian agent makes decisions according to the flow chart. After simulation, experimental data, such as the real-time velocity and realtime pedestrian distribution, are collected and analyzed.

Running the Simulation
. . Parameters Initialization. The simulation is based on data collected from a real tunnel about 300m long and 10m wide with a 12ped/s peak pedestrian flow. The peak of pedestrian flow lasts about 15 minutes. Over 10 gigabytes of videos on pedestrian movements is collected from the monitoring system. Through statistical analyzing, the possibilities of each alternative cell being chosen by the pedestrians are found to be 0 = 4 = 0.1, 1 = 3 = 0.2, 2 = 0.35, 5 = 0.05 on average. After the simulation begins, the pedestrians walk into the tunnel from either of the two ends, and experimental data are collected and analyzed after 1200 seconds (20 minutes) to identify characteristics of the pedestrian flow in the tunnel. The detailed input parameter is listed in Table 3.

. . Model Calibration and
Validation. An important pedestrian flow indicator, the Velocity of Pedestrian Flow (VPF), is defined as the average speed of all pedestrians in the tunnel at the same time. In Figure 7, the red dashed line stands for the actual average speed of pedestrians passing through tunnel, which is about 1.125m/s, while the blue fitting curve represents the VPF transformed over time in simulation. Though the model assigned an initial speed of 1.4m/s to each pedestrian (0.7m each step and two steps each second), it can be found from the experimental data that overall the pedestrians pass through the tunnel at a speed of around 1.11m/s. Considering the measurement error, it can be concluded that the VPF in the simulation experiment agrees with the real situation. Figure 8 demonstrates the number of pedestrians getting out of tunnel over time, with circles representing data of real situation and the blue curve representing data of simulation. Each circle stands for the average number of pedestrians in a minute. It can be found that the peak number of pedestrians getting out of tunnel highly accords with the real situation. From Figures 7 and 8, it can be concluded that the proposed model can effectively mimic the bidirectional pedestrian flow. To study the factors that influence the AVPF, experiments in scenarios with different tunnel lengths, different tunnel widths, and different pedestrian flow rates are conducted, and the results are shown in Figure 9. There are nine subgraphs in Figure 9, each with a layout of three rows and three columns. The left to right subgraphs in the three columns are the results for tunnel widths of 3m, 5m, and 10m, while the top to bottom subgraphs in the three rows show the experimental results for tunnel lengths of 100m, 200m, and 300m (this research has conducted experiments in various scenarios, such as different tunnel lengths (under 100m, exceeding 300m), different tunnel widths (under 3m, exceeding 10m), different rates of pedestrian flow (under 3ped/s, exceeding 12ped/s), and has gotten similar results. This paper chooses data collected in some typical scenarios to display). The line in each subgraph indicates the AVPF changing with the pedestrian flow increasing from 3ped/s to 12ped/s in each scenario.
It can be seen from Figure 9 that subgraphs in the same column have few differences, indicating that the tunnel length has limited effect on the AVPF; however, as lines varied across the same row, the tunnel width is found to have a significant influence on the AVPF; that is, the AVPF increases as the width broadened. Curves in each subgraph decrease as the pedestrian flow rate increases; therefore, it can be surmised that the wider the tunnel is, the weaker the influence of pedestrian flow rate has on the AVPF. In summary, the simulations indicate that the tunnel width and the pedestrian flow rate significantly affect the AVPF, while tunnel length has limited influence.
It is a common sense that the greater the pedestrian density is, the more dangerous the space will be. As the model allows more than one pedestrian to stand in a cell, potential risk will increase with the number of pedestrians in a cell growing. This research collects real-time data of pedestrian distribution in the tunnel throughout experiments and the average data for the pedestrian distributions in different scenarios during experiments are displayed in Figure 10.
In accordance with Figure 9, subgraphs in Figure 10 show average pedestrian distributions with different tunnel lengths, different tunnel widths, and different pedestrian rates. In each subgraph, there are five differently colored bars, which represent the percentages of cells containing zero pedestrians, one pedestrian, two pedestrians, three pedestrians, and over three pedestrians. From a comparative analysis, it is concluded that tunnel length has a limited impact on pedestrian distribution, but tunnel width has a significant influence. Concretely, the proportion of multipedestrian cells increases while the tunnel width narrows down, and in certain tunnels, the ratio of cells containing multiple pedestrians increases with the pedestrian rate. Besides, it can be found that, except for the extreme conditions (3m wide tunnel with 12ped/s), the percentage of cells containing more than three pedestrians is below 5%, which indicates that pedestrians in the simulation are under tolerable situations (pedestrian density under 9ped/m 2 ) [2].
As tunnel length is found to have limited impact on pedestrian flow movements, experiments are conducted using a fixed tunnel length of 100m. Experimental data are collected in scenarios with different tunnel widths and different pedestrian flow rates. Figure 11 shows the average pedestrian density from the start to the completion of each experiment. The three lines in the figure, respectively, illustrate the average density transferring with the pedestrian rate when the tunnel widths are, respectively, 3m, 5m, and 10m. From the experimental data, it is found that tunnel width and pedestrian flow rate have significant influences on average density. Therefore, based on the critical density definition (4ped/m 2 ), the SPFR is a function of tunnel width; that is, under certain circumstances, the SPFR is only related to tunnel width.
Through the many simulation experiments, SPFRs are collected under different scenarios, the fitting curve for which is shown in Figure 12, and the function is as follows: where y is the SPFR and x is the tunnel width. We thus conclude that the SPFR has a linear relationship with tunnel width. Specifically, the SPFR increases by 2.96ped/s as the tunnel width is expanded by 1m, which can be useful information for managing pedestrian flows and designing public places.

Simulation Experiments with Walking Preferences
. . Improved Simulation Model. In the former experiments, pedestrians use the Bayesian-Nash Equilibrium to estimate their expected utility for each choice. After that they take collision avoidance actions before moving into their chosen cells. However, pedestrians usually have preferences when making their own decisions; for example, pedestrians in China usually walk in the right while in Japan, they walk in the left. Further, adults tend to move forward, while children are more likely to move around. Therefore, even in the same situations, different pedestrians will make different decisions based on their preferences. As walking preferences affect pedestrian decisions, which may affect the pedestrian flow, the model is improved to study in depth the influence of microscopic walking preferences on macroscopic pedestrian flows.
A weight coefficient set is employed to illustrate pedestrian walking preferences; when a pedestrian makes their preferred choice, their comfort utility is times the previous. In the improved model, the total utility that a pedestrian expects to gain is defined as where is the serial number of the six choices, in which the numbering rule is in accordance with the possibility in Figure 3, denotes the total expected utility of choice , and , respectively, are the expected movement utility and the comfort utility, and denotes the pedestrian preference Table 4: Weight coefficient set in each scenario.

Scenario
The degree for choice . The two subgraphs in Figure 13 show the expected comfort utilities with and without walking preferences for pedestrians walking from left to right. By adjusting the weight coefficients, the walking preferences that affect pedestrian flow are examined in this research.
. . Simulation with Walking Preference. Based on the improved model, experiments are conducted with forward walking preferences and right walking preferences; the tunnel in each experiment is 100m, the width is 10m, and the rate of pedestrian flow entering the tunnel from both ends is 12ped/s in total.
( ) Forward Walking Preference. When a pedestrian is eager to move out of a tunnel, they tend to prefer cells in front. Therefore, to examine effects of forward walking preferences on pedestrian flow, comparative experiments in different scenarios are conducted. The weight coefficient set combinations for each scenario are shown in Table 4.
As different pedestrians may have different forward walking preferences, in the experimental scenarios the pedestrian preference degrees for cells in front are expected to obey a normal distribution with a mean value and a standard deviation 0.1. It can be found that the pedestrian preference degree for the cells in front increases with , and the number of pedestrians choosing the front cells increases as well. Figure 14 demonstrates the relationship between the AVPF and .
From Figure 14, it can be observed that the AVPF is about 1.1m/s when pedestrians have no forward walking preferences; however, when ≤ 1.2, the AVPF increases with and decreases when 1.2 ≤ ≤ 1.6. This is mainly because the pedestrian's desire to exit the tunnel increases, resulting in the pedestrians' distribution becoming uneven, and more pedestrians have to take collision avoidance actions, which in turn slow down the AVPF. When ≥ 1.6, the AVPF is stable at around 1.0m/s, which indicates that excessive forward walking preferences decrease the AVPF. The experiment indicates that a slight forward walking tendency promotes the pedestrian flow's movement, which can be used to improve pedestrian flow's efficiency in tunnels.
( ) Right Walking Preference. Note that the right in the walking direction of "from left to right" is relative to the tunnel, while in the "right-moving preference" it is relative to the pedestrians themselves. As mentioned, pedestrians in some areas prefer walking in the right, while in other areas they prefer walking in the left. The right walking pedestrian is taken as an example to examine the effect of side preferences on pedestrian flow. Note that naturally a pedestrian's right walking preference is no more than a forward walking preference; that is, pedestrians choose vacant cells in front  first. Therefore, the weight coefficient set combinations in each right walking preference experiment are shown in Table 5.
Similarly, the pedestrian's right moving preference degree is defined to obey a normal distribution with a mean value and a standard deviation 0.1. It can be found that the pedestrian preference degree to right hand cells increases with . Because of the right walking preference, which keeps bidirectional pedestrians in good order, pedestrians in two directions are divided into two lanes immediately. As a result, the right walking preference increases AVPF overall, which is different from forward walking preferences. Further, when ≥ 1.4, the AVPF fluctuated around 1.26m/s, which indicates that a right walking preference has a limited effect on the AVPF. This is mainly because, with increasing, the right hand side becomes crowded, and the speed is restrained. Figure 15 shows the AVPF changes over .
We thus conclude that it is beneficial for pedestrians to walk on the right (or left) as this separates the bidirectional Complexity 13 Table 5: Weight coefficient set in each scenario.

Scenario
The pedestrians and improves walking speeds. Therefore, by comparing Figures 14 and 15, it can be seen that side walking preferences better improve pedestrian flow efficiency than forward walking preferences.
. . Safe Rate of Pedestrian Flow considering Walking Preferences. Therefore, it has been proven that both forward walking and side walking preferences influence average pedestrian flow velocity, which can affect the SPFR. To examine this further, simulation experiments over seven scenarios are conducted to explore this effect. The seven scenarios are a comparison scenario (no walking preferences), three forward walking preference scenarios ( = 1.2, 1.4, 1.6), and three side walking scenarios ( = 1.2, 1.4, 1.6). In each scenario, many experiments are conducted to assess the SPFR effect, as shown in Figure 16.
As can be seen in Figure 16, both forward walking and side walking preferences increase the SPFR to some degree. The side walking preference (right) always increased the SPFR; however, the increase is not very obvious when ≥ 1.4. In the forward walking preference, the SPFR first increases and then decreases, with the conclusion being that, in consideration of public safety, the SPFR should not exceed that in the comparison scenario; that is, = 2.9643 − 1.0333 (13) where y is the SPFR and x is the tunnel width. In other words, the total flow rate of pedestrians entering a tunnel from both ends should satisfy this relationship for security reasons; that is, the SPFR can increase by 2.96ped/s as the tunnel width expands by 1m.

Conclusion
This paper studies crowd management by combining a multiagent system with the Bayesian-Nash Equilibrium. In the model, pedestrians make decisions based on their expected utility, which is calculated through the Bayesian-Nash Equilibrium and their own preferences, and then take collision avoidance actions before moving. The proposed model is calibrated and validated using data collected from a real tunnel. Simulation experiments over multiple scenarios are then conducted. From the experimental data, we find that the Safe Pedestrian Flow Rate increased by 2.96ped/s as the tunnel width expanded by 1m. We also find that both forward walking preferences and side (left or right) walking preferences are able to improve average pedestrian flow velocity, which can be beneficial for the crowd management.
Different from other researches, this paper focuses on mitigating the harmful effects caused by crowd accidents through managing safe pedestrian flows in advance. That is, according to the conclusion, crowd managers can decide how many pedestrians are allowed to go through a tunnel at the same time. Since the emergency management for crowds is a social science, which is mainly based on experience, our work can give references for the mangers and be beneficial for the public safety.
However, even at a low pedestrian flow rate, pedestrians may panic in an emergency; therefore, in future research, pedestrian psychology, social relationships, information transmission, and obstacles should be considered for results that are more accurate.

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.