Why Do Students Choose Buses over Private Motorcycles and Motorcycle-Based Ride-Sourcing? A Hybrid Choice Approach

: Improving public transport (PT) services is critical to maintaining resilient and sustainable transportation. While the investment in PT (i.e., buses) has not been able to shift private vehicle users, especially motorcyclists, to buses, the emergence of motorcycle-based ride-sourcing (MBRS) has become a popular travel mode that threatens the existing bus demand. Due to this, this study aimed to investigate the factors affecting the choice of buses over private motorcycles and MBRS among high school and university students in Yogyakarta, Indonesia, through the hybrid choice model. The explored variables include sociodemographic characteristics, travel time, travel cost, and latent variables of travel mode perceptions. By involving 527 respondents, the analysis showed that travel time, travel cost, and sociodemographic factors, including education level, income, and motorcycle ownership, signiﬁcantly impacted travel mode choices. The ﬁndings also revealed that convenience causes buses to be the chosen transport means among students, while lifestyle generates students to choose MBRS services. Based on ﬁndings, this study proposed some policies to increase bus demand, especially related to the identiﬁed convenience factors and travel time.


Introduction
In many cities worldwide, public transportation is heavily relied upon to take people to work, study, or other activities, especially for low-income households [1], students [2], the elderly [3], and disabled people [4]. Public transportation is beneficial because it uses road space much more efficiently than private transport modes [5]. Thus, public transport is a pivotal tool to overcome the dependence on cars, urban congestion problems, and poor air quality [6]. As traffic congestion is getting worse in urban areas, more cities are becoming aware that public transportation investments, such as subways, bus rapid transit, or buses, should be one of the top priorities [7].
Amid the efforts of many cities to promote the sharing of urban public transport and reduce the use of private transport, the use of online transportation is increasing. The emergence of online transportation brings major changes in the usage patterns of transportation modes. In the United States, the emergence of two online transportation services, namely Uber and Lyft, has replaced the role of public transportation in urban areas [8]. Online transportation services have-unavoidably-changed the preference of some public transportation passengers [9]. Approximately half of the trips using ridesourcing (online taxis) have replaced the preference for using conventional taxis or public transportation, going on foot, cycling, and driving a private car [10]. Moreover, Uber pool can replace urban public transportation between neighborhoods that require transfer [11]. countries due to the benefits it offers, such as requiring small investments, its increased flexibility, its ability to be established rapidly, and its offer of a larger capacity compared to traditional bus services [31]. A study in Multan, Pakistan, showed that approximately 54% of passengers were very satisfied and selected the BRT due to its convenience [32]. Furthermore, BRT-like (without separate bus lanes) transport modes in Yogyakarta consist of three LVs that most significantly affect the selection of bus services over MBRS services, namely being forced to use buses, the quality of bus services, and the benefits of using buses. These three variables have a more significant effect than the potential cost of and time saving offered when using MBRS. Additionally, the decision to use the BRT-like in Yogyakarta is also significantly influenced by users' age and income [16].

Methods Used in Determining Mode Choice Decision
Previous studies concluded that there are three main factors that affect transport mode choices [33]. The first is the characteristics of travelers, such as their background, household structures and income, vehicle ownership, and the availability of transport mode options. The second is the characteristics of the trip, such as the trip's purpose, travel time, and travel distance. The third is the characteristics of transportation facilities, such as travel time, cost, service quality, and availability of a parking lot.
However, previous studies applied different methods in exploring the factors affecting mode choice decisions. The binary and multinomial logit model is the most frequently used model in exploring the factors influencing mode choice decisions [34][35][36][37][38]. For example, a study involving 450 university students and 450 employees implemented a binary logit model in analyzing the choice between motorcycles and buses in Ho Chi Minh City, Vietnam [34]. A study in the Greater Toronto and Hamilton areas in Canada expanded the model by combining independent availability logit and constrained multinomial logit model formulation to explore the demand for ride-sourcing (e.g., Uber) services [39]. Meanwhile, Si et al. integrated the binary logit model with structural equation modeling to investigate the factors influencing the choice decision between regular taxis and online taxis in China [40]. Another study also attempted to integrate the nested logit model with structural equation modeling in the case of travel mode choice between bicycle, bus, and online taxi in Jinan City, China [41]. Additionally, Aguilera-Garca et al. used a generalized structural equation model (GSEM) to examine ride-sharing usage and personal choices between conventional taxi and ride-sharing services in Spain [42].
Out of all of these methods, the latest development related to the mode choice model is the inclusion of LVs. The integrated choice and latent variable (ICLV) models, also referred to as the HCM, was initially proposed by McFadden in 1986, and further developed by Ben-Akiva et al. in 2002 [43]. The ICLV model involves, simultaneously, survey theory, econometric theory, logical psychology theory, and economic theory. Furthermore, it involves and integrates various dimensions. Integrating LVs improves the understanding of transport mode options. To illustrate, a study discovered that the LVs enrich the model of discrete choices, outperforming the model of traditional discrete options [44]. Another study shows that the LVs impact the behavior of transport mode choice, and the integrated model of structural equation modeling-the discrete choice model-has higher accuracy and a better explanation of the variables [45]. Other studies have also involved LVs, for example, a study investigating 8284 respondents and discovering that a lack of of comfort and convenience in the public transport system was the main factor causing people to avoid using public transport [46]. The use of the ICLV increased the value of McFadden's R 2 from the MNL condition by 7% to 11.5% [47]. The goodness of the ICLV fit model was 0.201 higher than the traditional logit, demonstrating that the LVs impacted transport mode choice [45]. A study of transport mode choice for going to work found that the prediction level of transport mode choice models involving LVs increased by 7.48%, as compared to the models not involving any LVs [48].
To date, investigations of the LVs of mode choices have only focused on one mode [16,47]. By involving the LVs for all travel mode options (bus, MBRS, and private motorcycle), it is expected that this study could address a gap in previous studies. By this, the study could explore the contribution of each LV for every travel mode, which is hypothesized to have different impacts in each mode of travel. However, we note that recently there have been new developments in modeling choice decisions, such as the use of the deep neural network (DNN), a machine learning method. For example, a study in Seoul, Korea, analyzed travel mode choice by applying an interpretable machine learning approach [49]. Zheng et al. stated that DNN, a machine learning method, can produce lower prediction bias than the conventional discrete choice model [50]. However, for our study, it was challenging to integrate latent psychological construct into a DNN based model. Therefore, the present study utilized the HCM, an integrated latent variable model, to determine which factors affect the choice of transport modes among high school and university students in Yogyakarta, Indonesia. It should be noted that HCM has been applied in a significant amount of transportation literature, for example, in the field of electric vehicle choice [51][52][53], autonomous vehicle [54,55], bikesharing [56], mobility as a service [57], toll and residential choice [58,59], MBRS use [60], and BRT and motorcycle taxi choice [61].

Study Variables
This study discussed the travel choices between a private motorcycle, MBRS, and buses. The study used revealed preference data, namely measurement based on actual market behavior [62]. The involved variables included sociodemographic characteristics [63], trip characteristics, and LVs. The sociodemographic variables included gender, age, level of education, income, vehicle ownership, and driving license ownership. The trip characteristics comprised travel time and cost of any transport mode. This study used an approach based on the respondents' perceptions of the services of a transportation mode [64]. This study focused on the students' trip to their school/college. The target group was only high school and university students, and not all Yogyakarta city residents, because bus passengers and MBRS consumers in Yogyakarta are dominated by high school and university students [14,65]. Table 1 presents the questions related to the respondents' sociodemographic variables.

Survey Instrument
The respondents also answered questions related to their perceptions of three travel modes: private motorcycle, MBRS, and city bus. Each transport mode had 5 LVs, and 16 indicators. These 16 indicators were answered on a 5-point Likert scale ranging from 1 = strongly disagree to 5 = strongly agree. The 16 indicators can be categorized into five LV groups: flexibility [66][67][68][69], security and safety [70][71][72][73][74] convenience/comfort [70,75,76], pro-environment [67,75,77], and images/lifestyles [48,78]. These chosen variables were based on prior literature. To illustrate, transportation research has widely examined the variables of flexibility, security, and safety. The pro-environment variable was included because transportation contributes to deteriorating environmental conditions, such as global warming [79]. Meanwhile, the images or lifestyles variable was included because young people or students usually prefer to follow ongoing trends. Table 2 shows the questions regarding perceptions of transport modes.
The survey was conducted online using the probability sampling method distributed via a WhatsApp group. This type of survey was conducted for two reasons. First, since the COVID-19 pandemic was still ongoing, it was not possible to conduct the survey through face-to-face interviews. Second, the online survey method was more cost-and time-efficient. Although the survey was conducted during the COVID-19 pandemic, the respondents were asked about their travel behavior before the pandemic. Eleventh-grade high school students 3 Twelfth-grade high school students 4 First-year university students 5 Second-year university students or more

Sample and Data Collection
The sample included senior high school and university students in Yogyakarta. The inclusion criteria were students who had lived in Yogyakarta for at least one year. This aimed to ensure that they understood the condition of the three transportation modes. This study also required that the respondents lived in a place where the distance to the bus stop was no more than 400 m. It considers that people's ability to walk to the bus stop is about 400 m [80], and therefore, respondents have an opportunity to use bus mode instead of MBRS and private motorcycle. Furthermore, it should be noted that even though we acknowledge that there is another method to collect data using big data provided by a navigation service provider, such as Baidu Maps in the China case [81], the present study incorporates latent psychological construct which can not be collected using navigation service provider app. Therefore, the questionnaire survey keeps used for this study. The respondents were randomly obtained from WhatsApp groups to obtain a random sample. Since the research subjects were school and college students, WhatsApp groups in which the members were predominantly students were selected to maximize the response rates. However, public WhatsApp groups could still receive the online survey forms if the number of WhatsApp group members who were students could be identified. In addition, each of the group members were not allowed to send the Google Form to other groups, making it possible to identify the response rate. To ensure the respondents' seriousness in completing the online form, the survey used experimental control sentences; these sentences contained questions that could easily be answered. When the respondents gave an incorrect answer in an experimental control sentence, the form could not be included in the following analysis phase.
Data were collected in November and December 2020. Overall, 1250 respondents completed the survey. However, only 603 respondents passed the initial stage of the survey; 23 respondents answered the experimental control questions incorrectly, and 622 respondents did not meet the criteria because the distance of their houses to the Trans Jogja bus stop was more than 400 m.
Then, the survey results of the 603 respondents who passed the initial stage were analyzed, particularly concerning to questions related to travel time and cost when using motorcycles, MBRS, and Trans Jogja buses. Unfortunately, 76 respondents answered with illogical data. For example, they mentioned that the travel time using the Trans Jogja bus was much faster than riding a private motorcycle. In addition, some respondents answered that it took less than 10 min to reach a destination when using the Trans Jogja bus, including the time they needed to walk to the bus stop. All the forms containing illogical data were excluded from the following stage of analysis. Thus, there were only 527 respondents for which the data were used for the next stage of analysis. These 527 respondents met the minimum requirement of 400 [84] or 500 people [43,85]. Table 3 shows the survey results for respondents' socio-economic characteristics, travel time, and cost of each transport mode. It shows that most of the respondents (62.3%) were female, and were aged 15 to 24 years old. The higher number of females is similar to previous studies, where bus passengers and MBRS customers in Yogyakarta are dominated by females, accounting for 61.3% and 57.56%, respectively [14,65]. A total of 50% of them were above 23 years old. More than 50% of the participants had an income of IDR 500,000 or more. A total of 62% of them had a driving license, and only 16% did not have a private motorcycle. Furthermore, the relationship between sociodemographic conditions of the respondents and the use of transportation modes is presented in Figure 1. It shows that the proportion of male and female users of motorcycles, MBRS, and buses was relatively balanced. Most motorcycle riders were aged 25 years old and above, while 100% of MBRS users were 15 years old. The proportion of the use of motorcycles increased, along with a decrease in the use of MBRS and buses when the education level was higher. Regarding income, the proportion of the use of motorcycles, MBRS, and buses fluctuated, although the proportion of motorcycle use in various income levels was always above 50%. The proportion of the use of motorcycles was lower among college students who did not have a private motorcycle. Moreover, the proportion of the use of transport modes (motorcycle, MBRS, and buses) was relatively evenly divided among the college students who did not have a private motorcycle.

Perception Variables
The reliability of the perception variables constructing the LVs was tested to investigate the consistency of the respondents' answers. The results of the reliability test are presented in Table 4. It shows that all of the variables had a Cronbach's Alpha greater than 0.6. This figure indicates that the consistency of the respondents' answers met the criteria. Then, the validity test was conducted to determine the level of accuracy of the respondents' answers. The corrected item-total correlation showed that the r-count was comparable with the significance level. All of the answers were valid because the r-count was greater than the r-table (0.085). Meanwhile, the details of the results of respondents' perceptions of the five LVs for each transport mode is presented in Table 5 below. Based on Table 5, the highest LV assessment of each transport mode is on the flexibility of motorbikes/MBRS and the pro-environment factor for the Trans-Jogja bus. The lowest rating is on the pro-environment for motorbikes/MBRS and the flexibility of the Trans-Jogja bus. Furthermore, the relationship between sociodemographic conditions of the respondents and the use of transportation modes is presented in Figure 1. It shows that the proportion of male and female users of motorcycles, MBRS, and buses was relatively balanced. Most motorcycle riders were aged 25 years old and above, while 100% of MBRS users were 15 years old. The proportion of the use of motorcycles increased, along with a decrease in the use of MBRS and buses when the education level was higher. Regarding income, the proportion of the use of motorcycles, MBRS, and buses fluctuated, although the proportion of motorcycle use in various income levels was always above 50%. The proportion of the use of motorcycles was lower among college students who did not have a private motorcycle. Moreover, the proportion of the use of transport modes (motorcycle, MBRS, and buses) was relatively evenly divided among the college students who did not have a private motorcycle.

Perception Variables
The reliability of the perception variables constructing the LVs was tested to investigate the consistency of the respondents' answers. The results of the reliability test are presented in Table 4. It shows that all of the variables had a Cronbach's Alpha greater than 0.6. This figure indicates that the consistency of the respondents' answers met the criteria. Then, the validity test was conducted to determine the level of accuracy of the respondents' answers. The corrected item-total correlation showed that the r-count was comparable with the significance level. All of the answers were valid because the r-count

Hybrid Choice Model
The HCM or integrated choice and latent variable model is an extension of a discrete choice model that involves LVs of psychological constructs. To estimate the parameter in the HCM, this study used Biogeme, an open source Python software, to maximize the likelihood estimation [86]. However, given that it takes a long time to run the Biogeme program, LVs that significantly affect the mode choice of the five LVs for each transportation mode were initially selected. By using principal component analysis, the results showed that LV 2 (security and safety) and LV 5 (image/lifestyles) did not significantly influence the choice of motorcycle. Meanwhile, LV1 (flexibility), LV 3 (convenience/ comfort), and LV 4 (Pro-environment) did not significantly influence the choice of MBRS and LV 1 (flexibility), LV 2 (security and safety), and LV 4 (pro-environment) did not significantly influence the choice of bus mode. Thus, these insignificant LVs were excluded from the HCM. The remaining LVs, together with the sociodemographic variables, travel time, and travel costs, were then examined to find the significant parameter influencing transport mode choice. Figure 2 shows the proposed modeling framework in this study.   There are two equations in the HCM. The first is a structural equation exploring the effect of respondents' sociodemographic characteristics on the LVs. For example, the structural equation exploring the effect of respondents' sociodemographic characteristics on the LV of flexibility (LV 1 ) is formulated in Equation (1) below. LV 1.n = ζ 1.n + λ 1.1.n .gender + λ 2.1.n .age + λ 3.1.n .education + λ 4.1.n .income+ λ 5.1.n .vehicle ownership + λ 6.1.n .driving licence + ε 1.n with n = 1, 2, 3 (1) where n = 1 for motorcycles, n = 2 for MBRS, and n = 3 for bus. The notation of λ 1.1.n to λ 6.1.n is the estimated parameter for gender, age, education, income, vehicle ownership, and driving license, respectively. The notation of ζ 1.n is the intercept, and ε 1.n is the error terms, which is normally distributed with a zero mean and a standard deviation of σ. A similar equation was also applied to the structural equation for other LVs (LV 2 to LV 5 ) The second equation is a measurement equation, which expresses the relationship between the LVs and their indicators (observed variables). The measurement equation for the LV of flexibility for motorcycles (LV 1.1 ) is shown in Equations (2)-(4). (2) In the equation, I 1.1 shows an observed indicator (i.e., being able to go anytime) for the LV of flexibility for motorcycles (LV 1.1 ). The notation of δ 1.1 and η 1.1 is the intercept and error term for the observed variable of being able to go anytime, respectively. The error term is normally distributed with a zero mean and a standard deviation of σ. A similar explanation also occurs for I 2.1 (an observed indicator of being able to go anywhere) and I 3.1 (an observed indicator for not requiring a transfer).
Furthermore, the formulation of the utility (U) of the HCM in each transport mode is written according to Equation (5) for motorcycle, Equation (6) for the MBRS, and Equation (7) In the equation, α 1 and α 2 are the parameter of preference for travel time and cost, while β m.n is the parameter of the LVs. Meanwhile, ω is the intercept of each mode, and ν is the error factor.

Results and Discussion
Tables 6-12 show the model results. There were three models in this study. The first was a measurement model explaining the relationship between latent and observed variables constructing the LVs. The second was a structural model revealing the relationship between sociodemographic characteristics and LVs. The last model was the HCM, showing the factors influencing transport mode choices, consisting of LVs, travel time, and travel cost. Each of the models consisted of three transport mode choices: MBRS, motorcycle, and bus. Tables 6-8 shows the measurement model for motorcycle, MBRS, and bus, respectively.      Based on the Tables 6-8 of the measurement model, the parameter estimation related to the LV was positive, as expected. When an LV is positive, an increase in the values of the LV increases the possibility for respondents to "strongly agree" with the questions containing the indicators of the LV. Furthermore, Tables 9-11 shows the structural model for motorcycles, MBRS, and buses, respectively. Table 9 denotes that the sociodemographic conditions affecting the LV of motorcycles were income and vehicle ownership. The higher the respondents' income, the higher the perceived flexibility of motorcycles. Comparatively, the higher the income, the lower the perceived pro-environment factor and convenience of motorcycles. When the respondents did not have full ownership of a motorcycle, the flexibility, pro-environment factor, and convenience of a motorcycle decreased. Meanwhile, Table 10 shows that only age had a significant correlation with the lifestyle and security and safety variables of MBRS. That is, the older the student, the lower the lifestyle and perceived safety of MBRS.
Additionally, Table 11 shows that the LVs of perceived travel mode performances were influenced by education, vehicle ownership, and driving license ownership. A previous study in Bandung, Indonesia, also revealed that sociodemograhic factors, including age, gender, education, and income, significantly influenced the perceived performance of travel mode [60]. Meanwhile, this study finding shows that the higher the education, the lower the lifestyle and convenience. Similarly, when the respondents had full ownership of a motorcycle, they perceived lower flexibility, pro-environment, and convenience of buses. Meanwhile, when the respondents did not have full ownership of a driving license, they perceived a lower lifestyle and convenience of buses.
Finally, Table 12 shows the model result of the HCM. When the value of the r-square is 0.502, it means that the considered variables can explain 50.2% of the factors influencing transport mode choices. As expected, travel time and travel cost were negative, meaning that the longer the travel time or the more expensive the travel cost, the lower the utility of the transport mode. The travel time variable with a coefficient of −0.523 was more important than the travel cost variable (coefficient of −0.0619), because the coefficient was greater. That is, the respondents preferred a shorter travel time to more affordable travel costs regarding the trip from home to school or college.
The three LVs included as independent variables in the HCM were also significant and positive, meaning that the higher the assessment of the LVs, the higher the probability of the transport mode to be selected. These three LVs had a greater coefficient/parameter than the coefficient of travel time or cost. This means that the LVs had a more significant effect on transport mode choices than travel time and cost. This finding is consistent with another study focused only on the competition between public transport and MBRS in Yogyakarta [16].
Of the three LVs, the highest coefficient was obtained by the convenience of motorcycles, at 14.4, meaning that motorcycles had a much greater convenience/comfort than buses (2.74). The LV of lifestyle became the determinant for the transportation mode choice of MBRS with a coefficient of 5.73. This suggests that, for students, lifestyles are an important variable to consider, because they must maintain their prestige by following the trend according to the era. A study on MBRS in Yogyakarta also shows a similar result, where feeling in touch with modernity becomes a reason why millennials in that city use MBRS services [14].

Conclusions and Policy Recommendations
The purpose of this study was to investigate the travel mode choice for high school and university students in Yogyakarta by including sociodemographic characteristics, travel mode characteristics, and the LV of travel modes' performances. By applying an HCM and involving 527 respondents, this study produced three important findings. Those findings could be used as a reference in determining the appropriate policies, especially to increase public transport demand in Yogyakarta. First, out of the considered variables, the LV of travel modes' performances had the highest effect on students' decisions in using a travel mode. In this case, the LV of perceived convenience influenced the use of buses, and the LV of lifestyle affected the use of MBRS services. Meanwhile, LV of perceived convenience also significantly influenced the choice of private motorcycles. Second, both travel time and cost significantly affected the choice of buses, MBRS services, and private motorcycles, in which the travel times more affected the travel mode decision rather than travel costs. Third, sociodemographic characteristics influenced the LV of transport means' performances, except for gender. The model results revealed that gender was not correlated with the LVs of all transport means. The variable of age affected the LV of lifestyle in using MBRS. However, the age factor did not affect the perceived convenience of city bus use. Meanwhile, education, income, and vehicle ownership significantly affected the LV of lifestyle in using MBRS, as well as the LV of perceived convenience in using city buses and private motorcycles. Additionally, the factor of driving license ownership also affected the LV of convenience for bus use and the LV of lifestyle for MBRS use, but did not affect the convenience LV of motorcycle use.
Based on the findings, some policy recommendations are proposed to the government to increase passengers' preference of using city buses in Yogyakarta. First, decreasing the travel time by decreasing the waiting time, decreasing the distance from home to a bus stop, and decreasing in-vehicle travel time should be implemented. The waiting time can be reduced in two ways: developing a real-time bus arrival information system with intelligent transportation systems [87], and adding to the number of buses. Decreasing the distance from home to a bus stop can be achieved by expanding the coverage of bus services in urban areas. In addition, creating a convenient pedestrian pavement to reach the bus stop is also necessary. Meanwhile, decreasing in-vehicle travel time can be achieved by building a bus priority section in intersections or building a special lane for buses. However, choosing which strategies should be prioritized to reduce the travel time of city buses should be determined by conducting a more in-depth study.
Second, competitive prices for city buses should be created to not burden school or college students. At the time when this study was conducted, the ticket price of the Trans Jogja bus was IDR 3500 for a single trip, while the subscription ticket price was IDR 2700 for a single trip for non-students and IDR 1800 for students. These prices are still quite competitive for passengers. The ticket prices of new buses can be raised if passengers gain other new advantages which exceed the extra money they spend on the ticket. The advantages can be in the form of significant reduction of travel time. However, further studies regarding reduced travel time that is worth a price increase should be conducted.
Third, it is important to maintain the convenience of city buses, so passengers do not choose other transport modes. The convenience variable in this study was defined using three indicators, namely big seats, comfortable suspension, and a chance to relax during the trip. Big seats are maintained by focusing on the load factor [88], so the Trans Jogja buses will not be overloaded. Comfortable suspension is ensured by periodical bus maintenance.
In addition, bus passengers should be able to relax and perform other activities during the trip, such as using a smartphone. Fourth, since MBRS has become a part of students' lifestyles, the government needs to ensure the safety of MBRS, such as by creating a safety campaign. Furthermore, for MBRS to contribute to the use of public transport, they should act as a feeder service for city buses if the travel distance is far enough. If the travel distance is higher, the cost of MBRS increase exponentially; therefore, it is better to act as a feeder for long-distance travel. Last, based on the survey results, most motorcycle users did not have a driving license. Therefore, it is crucial to create a campaign for driving license ownership among students. In addition, controlling the number of motorcycle users can lower the rate of accidents and air pollution. These objectives can be achieved by strictly regulating driver's license tests and creating a cycling campaign. Strictly regulating driving license tests will ensure students' ability to drive and ride safely on the highways, thus reducing the rate of traffic accidents.
This study has several limitations that need to be considered and followed up for future studies. First, this study used the revealed preference approach; thus, it cannot be used as a simulation of various scenarios. Therefore, future studies should be conducted which use the stated preference approach by providing various scenarios for travel time and cost. Second, we realized that our two-mode alternatives in this study, motorcycle and MBRS, might not be safe, and lead to high vulnerability on the road. Therefore, future studies might also consider the safety aspect, especially considering health care costs, as discussed in Corazza et al. [89]. Third, this study needs to be carried out in other Indonesian cities, since each city has specific characteristics resulting in different travel behavior and travel mode decisions, as discussed and compared by Nugroho et al. [13]. Although this study had several limitations, overall, this study's results offer important insights as an initial step to improving public transportation in the future, and can guide future research on the topic.

Funding:
This study was self-funded by Universitas Janabadra.

Institutional Review Board Statement:
The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of Universitas Islam Indonesia.