# A Logistic Regression Model with a Hierarchical Random Error Term for Analyzing the Utilization of Public Transport

Authors: Chong Wei; Tingting Lu; Xuedong Yan
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Logistic regression models have been widely used in previous studies to analyze public transport utilization. These studies have shown travel time to be an indispensable variable for such analysis and usually consider it to be a deterministic variable. This formulation does not allow us to capture travelers’ perception error regarding travel time, and recent studies have indicated that this error can have a significant effect on modal choice behavior. In this study, we propose a logistic regression model with a hierarchical random error term. The proposed model adds a new random error term for the travel time variable. This term structure enables us to investigate travelers’ perception error regarding travel time from a given choice behavior dataset. We also propose an extended model that allows constraining the sign of this error in the model. We develop two Gibbs samplers to estimate the basic hierarchical model and the extended model. The performance of the proposed models is examined using a well-known dataset.


---

## Body

## 1. Introduction

Understanding the utilization of public transport is important for policy design and urban traffic planning. From a behavior analysis perspective, such utilization can be analyzed in terms of a binary choice problem in which the traveler must choose between public transit and a private mode of transport. Previous studies usually employ logistic regression models to discuss this binary choice problem. These models can be used to predict choice probability and to evaluate the effect of various attitudes on the utilization of public transport. The effects of demography, travel cost, travel time, and accessibility on such utilization can be analyzed through estimating the parameters of these models.McGillivray [1] proposed the binary choice problem in the two-modal (public and private) case and applied a logistic model to investigate the dependence of modal choice on the individual values of time and cost by modal. Gebeyehu and Takano [2] showed that the logistic model is a useful tool for evaluating transit systems; their study developed a logistic model to investigate citizens’ perceptions of bus conditions in Addis Ababa. Hess [3] attempted to apply a logistic regression model to assist planners and policymakers in expanding the mobility and accessibility of public transport for older adults by analyzing the influence of the accessibility of public transport on ridership for older people in California and New York. Johansson et al. [4] further formulated the binary logit model with latent variables; this structure can examine the effects of attitudes and personality traits on mode choice. Muley et al. [5] and Badland et al. [6] focused on microscopic behavior in the binary choice problem. Muley et al. [5] employed a binary logistic regression model to explore the impact of personal and transit characteristics on the utilization of public transport. Badland et al. [6] adopted a logistic regression model to discuss how parking availability and public transport accessibility influenced the split between uptakes of the two modes. However, Buehler and Pucher [7] studied the binary choice problem from a macroscopic perspective and employed a logistic model to compare utilization of public transport in the US and Germany.One benefit of using a logistic regression model to analyze the utilization of public transport is that the model can consider the combined effects of attributes through a linear combination and can easily evaluate the contribution of various attributes to such utilization. Previous studies agree that travel time is an important variable for formulating a logistic regression model to analyze this binary choice problem. Although previous studies usually prefer to treat travel time as a deterministic variable in the model, one cannot neglect that travelers’ perception error can prevent travelers from accurately evaluating their actual travel time. On the other hand, travelers also cannot say exactly how long a given travel time, for instance, 10 minutes, actually is. Carrion [8] provided a case analysis for this error using self-reported as well as GPS measured travel times and suggested that perception error regarding travel time can influence travel behavior. Cheng and Tsai [9] studied travelers’ perceptions of travel time through a questionnaire survey. They found that travel time perceptions can be influenced by personal characteristics such as gender and age.As aforementioned, researchers have recognized the effect of perception errors regarding travel time on travelers’ choice behavior. However, the classic formulation of a random error term in logistic regression models cannot reflect this perception error appropriately (Chen et al. [10]). Analyzing travel behavior using the classic logistic regression model might therefore prevent us from looking for further insights into the utilization of public transport. Hence, this study proposes a hierarchical logistic regression model to fill the gap. Previous studies have successfully applied hierarchical error terms to model heterogeneous unobservable utility in logistic regression models (e.g., Tilahun et al. [11] and Czado and Prokopenko [12]). Inspired by these studies, we adopt hierarchical error terms to solve the problem. However, we add these terms after important attributes such as travel time rather than after the entire utility function in order to capture travelers’ perception error regarding the attributes. This study is not merely an effort to improve the performance of the regression model; more importantly, the proposed model offers an alternative way to estimate the statistical characteristics of perception error regarding travel time and allows us to explore the property of this error. Although this study prefers to focus on modeling perception error regarding travel time, we believe that individuals might have perception errors on other attributes such as travel cost. The proposed model can be easily used to estimate other such errors as well.We first propose a basic hierarchical model and then develop an extended model in this study. The extended model allows us to constrain the sign of perception error regarding travel time. Correspondingly, this study also develops two Gibbs samplers to estimate the parameters of the proposed models. We evaluate the performance of the proposed models using a well-known dataset provided by Horowitz [13, 14].


## 2. The Hierarchical Logistic Regression Model

For simplicity, we describe the proposed model based on a binary choice problem that was provided by Horowitz [13]. In this choice problem, an agent faces making a choice decision between using a private car or public transport. The attributes considered in this choice problem are CARS, DCOST, DOVTT, and DIVTT whereDCOST is “public transport fare minus private car travel cost,”CARS is “private cars owned by the traveler’s household,”DOVTT is “public transport out-of-vehicle time minus private car out-of-vehicle time,”DIVTT is “public transport in-vehicle time minus private car in-vehicle time.”We formulate the choice problem through a logistic regression model. If we present this logistic regression model as a latent-variable model, then the model can be obtained as follows:(1)Yik=1ifUik>00otherwiseUik=β0+β1·CARSi+β2·DCOSTik+β3·DOVTTik+β4·DIVTTik+τikτik~Logistic0,1,where i and k are the indexes of individual and survey questions, respectively. Yik=1 denotes that the choice result is private car and is otherwise public transport. The logistic regression model uses a random error term τik to present stochastic behavior; however, this formulation does not allow us to investigate the perception error of attributes. In this study, we introduce a random variable ΔTi into the logistic regression model. Here ΔTi denotes the random perception error on DIVTT for individual i. We consider that the value of ΔTi can be different among individuals and that ΔTi follows a probability distribution. The logistic regression model with a hierarchical random error term can be obtained as follows:(2)Yik=1ifUik>00otherwiseUik=β0+β1·CARSi+β2·DCOSTik+β3·DOVTTik+β4·DIVTTik+ΔTi+τikτik~logistic0,1ΔTi~Nμ,1μ~πμβ~ζβ.In the above equation, we use β to denote the parameter vector [β0,β1,β2,β3,β4]. In the proposed model, ΔTi is a latent random variable, which serves as the random error term for DIVTT rather than the random term of the utility function. This structure enables us to investigate the perception error of DIVTT. We assume that ΔTi follows the normal distribution with unknown mean μ and specify the variance of the normal distribution as 1 for identification purposes. π(μ) and ζ(β) denote the prior distributions of μ and β, respectively. We wish to estimate the value of μ along with β. In this hierarchical model, DIVTT is divided between DIVTTik and ΔTi; this leads to DIVTT no longer being a deterministic variable in the model.The logistic regression model shown by (2) can be further extended. Here, we consider that ΔTi<0 indicates that individual i subjectively weakens the difference of in-vehicle travel time when the agent makes the choice decision; on the other hand, ΔTi>0 suggests that the individual subjectively enlarges the difference of in-vehicle travel time. In addition, we impose the value of DIVTTik+ΔTi to be a nonnegative (nonpositive) value if DIVTTik>0 (≤0). Following these considerations, we rewrite Uik of (2) as(3)Uik=Vik+τik=β0+β1·CARSi+β2·DCOSTik+β3·DOVTTik+β4·DIVTTik∗+τik,where(4)DIVTTik∗=ϕDIVTTikmin⁡0,DIVTTik-ΔTi+1-ϕDIVTTikmax⁡0,DIVTTik+ΔTi.In (4), if DIVTTik<0 then the value of ϕDIVTTik is 1 or else 0.As shown by (2) and (3), the combination of ΔTi and τik leads to a hierarchical random error term structure; this structure can be presented as a Directed Arc Graph (see Figure 1).Figure 1
The DAG of the parameters of the hierarchical logistic regression model.


## 3. Estimation Methods

First, we discuss how to estimate parameters of the model shown by (2). Stefanski [15] indicated that the logistic distribution can be represented as a normal scale mixture. Accordingly, Holmes and Held [16] suggested an auxiliary variable method to present the logistic regression model. Along the same lines, the regression model presented by (2) can also be derived in the following form:(5)Yik=1ifUik>00otherwiseUik=β0+β1·CARSi+β2·DCOSTik+β3·DOVTTik+β4·DIVTTik+ΔTi+εikεik~N0,λikλik=4ψik2ψ~KSΔTi~Nμ,1μ~πμβ~ζβ.μ and β are the parameters of the proposed model.U=[Uik∀i,k], ΔT=[ΔTi∀i], and λ=[λik∀i,k] are the latent parameters of the proposed model. P(β,U,ΔT,μ,λ∣X,Y) is the joint distribution of these parameters, which can be uniquely identified by (2). π(μ) and ζ(β) are the prior distributions for μ and β, respectively. KS denotes the Kolmogorov-Smirnov distribution.We develop a Gibbs sampler to draw the samples ofβ, U, λ, ΔT, and μ from the joint distribution and calculate the estimates of the parameters through aggregating the samples. Let Xik, X, and Y denote the vector [CARSi,DCOSTik,DOVTTik,DIVTTik], X=[Xik∀i,k], and Y=[Yik∀i,k], respectively. We use t to denote the indexes of the sample. The outline of the Gibbs sampler can then be given as follows:Algorithm 1.
Step  0 (initialization). Set initial values: β0=0, U0=0, λ0=1, ΔT0=0, and μ0=0.
Sett←1.
Step  1. Draw tth sample of β(6)βt~Pβ∣Ut-1,λt-1,ΔTt-1,μt-1,X,Y.
Step  2. Draw tth sampler of U(7)Ut~PU∣βt,λt-1,ΔTt-1,μt-1,X,Y.
Step  3. Draw tth sampler of λ(8)λt~Pλ∣βt,Ut,ΔTt-1,μt-1,X,Y.
Step  4. Draw tth sample of ΔT(9)ΔTt~PΔT∣βt,Ut,λt,μt-1,X,Y.
Step  5. Draw tth sample of μ(10)μt~Pμ∣βt,Ut,λt,ΔTt,X,Y.If t<T then go toStep  1; otherwise, stop the sampling procedure.The set of conditional distributionsβ, U, λ, ΔT, and μ is referred to as the full set of conditional distributions of P(β,U,ΔT,μ,λ∣X,Y). Sampling β, U, λ, ΔT, and μ in turn from the conditional distributions is equivalent to sampling these parameters from P(β,U,ΔT,μ,λ∣X,Y) simultaneously. We then illustrate how the Gibbs sampler draws tth samplers for β, U, λ, ΔT, and μ in detail.Algorithm 2.
Step  1. Update Xik′ for drawing tth samplers fori=1 to I dofork=1 to Ki doDIVTTik,t-1′←DIVTTik+ΔTit-1Xik,t-1′←[CARSi,DCOSTik,DOVTTik,DIVTTik′]endend
Step  2. Draw tth sample for β.
The conditional distributionP(β∣Ut-1,λt-1,ΔTt-1,μt-1,X,Y) is equivalent to P(β∣Ut-1,λt-1,X′,Y). On replacing X in Albert and Chib [17] by X′, the sample of β can be generated through the sampling scheme proposed by Albert and Chib [17]. As shown by them, ζ(β) can be specified as a diffuse distribution, and P(β∣Ut-1,λt-1,X′,Y) can be derived as a normal distribution.
Step  3. Draw tth sample for U.
The conditional distributionPU∣βt,λt-1,ΔTt-1,μt-1,X,Y is equivalent to P(U∣βt,λt,X′,Y). Albert and Chib [17] indicated that P(U∣βt,λt,X′,Y) is a truncated normal distribution.
Step  4. Draw tth sample for λ.
The conditional distributionλt~Pλ∣βt,Ut,ΔTt-1,μt-1,X,Y is equivalent to λt~Pλ∣βt,Ut,X′. Holmes and Held [16] proposed a rejection sampling scheme that can be used to draw samples from Pλ∣βt,Ut,X′.
Step  5. Draw tth sample for ΔT.
Δ
T
t can be generated through sampling ΔTiti=1,…,I in turn. The conditional distribution of ΔTit is equivalent to PΔTi∣βt,X,Y,Ut,λt. The probability density function of ΔTi is a normal distribution and is defined by (11)–(13).
Step  6. Draw tth sample for μ.
The conditional distributionPμ∣βt,Ut,λt,ΔTt,X,Y is equivalent to Pμ∣ΔTt. The probability density function of Pμ∣ΔTt is a normal distribution and is given by (16).Applying the Bayesian theorem, we obtainPΔTi∣βt,X,Y,Zt,λt as a normal distribution. Without loss of generality, we only present how to obtain PΔTi∣βt,X,Y,Ut,λt for the case Ki=1 for i=1,…,I: (11)PΔTi∣βt,X,Y,Ut,λt=1a2πe-ΔTi-Uik∗/β4t2/2λikt/β4t212πe-ΔTi-μt-12/2∝e-ΔTi-A2/2B.The above equation shows that the conditional probability of ΔTi is a normal distribution. The mean of the distribution is(12)A=Uik∗/β4t+μtλikt/β4t21+λikt/β4t2and the variance of the distribution is(13)B=λikt/β4t21+λikt/β4t2.In (12), Uik∗ represents(14)Uik∗=Uikt-β0t-β1t·CARSi-β2t·DCOSTik-β3t·DOVTTik-β4t·DIVTTik.We consider the prior distribution of μ as a normal distribution with mean 0 and variance σ2. If we further consider the prior distribution as a diffuse prior then this probability function of the prior can be given as(15)Pμ=lim⁡σ→+∞⁡1σ2πeμ-02/2σ2.Accordingly, the conditional distribution Pμ∣ΔTiti=1,…,I can be obtained as(16)Pμ∣ΔTiti=1,…,I=11/I2πeμ-∑i=1IΔTit/I2/2/I.This conditional distribution is a normal distribution too; the mean of the distribution is ∑i=1IΔTit/I and the variance of the distribution is 1/I.Now, we discuss how to estimate the parameters of the extended model described by (3). In this case, we can also employ Steps 1 to 3 and 5 of Algorithm 2 to draw the samples of β, U, λ, and μ. We need only to use DIVTTik,t-1∗ and Xik,t-1∗ to replace DIVTTik,t-1′ and Xik,t-1′ in Step  1, respectively, where DIVTTik,t-1∗ is defined by (4) and Xik,t-1∗ is defined as [CARSi,DCOSTik,DOVTTik,DIVTTik∗]. On the other hand, instead of Step  4, we develop an M-H sampling scheme to drawΔTii=1,…,I.We derive the formulation ofPΔTi∣βt,Xi,Uit,λit,μt-1 as(17)PΔTi∣βt,Xi,Uit,λit,μt-1=PΔTi∣μt-1∏k=1KiPYik∣Xik,βt,ΔTi,where PΔTi∣μt-1 is the probability density function of the normal distribution with mean μt-1 and variance 1. Following (3), we get PYik∣Xik,βt,ΔTi as(18)PYik∣Xik,βt,ΔTi=exp⁡Vik1+exp⁡Vik,where Vik is defined by (3). The H-M sampling scheme is as follows.Algorithm 3.
Step  0 (initialization). Set i←1.
Step  1. Generate a candidate value for ΔTit.
DrawΔTi∗ from a normal distribution: ΔTi∗~N(ΔTit-1,1), where ΔTi∗ denotes the candidate value forΔTit.
Step  2. Calculate the value of η(19)η=min⁡PΔTi∗∣μt-1∏k=1KiPYik∣Xik,βt,ΔTi∗PΔTit-1∣μt-1∏k=1KiPYik∣Xik,βt,ΔTit-1,1.
Step  3. Get ΔTitγ~U(0,1), if η>γ then ΔTit=ΔTi∗ else ΔTit=ΔTit-1.Ifi<I then i←i+1 and go toStep  1.The proposed model structure also allows us to further investigate traveler’s perception error regarding both DIVTT and DOVTT. To do this, we just need to modify (3) as follows:(20)Uik=β0+β1·CARSi+β2·DCOSTik+β3·DOVTTik∗+β4·DIVTTik∗+τik,where(21)DOVTTik∗=ϕDOVTTikmin⁡0,DOVTTik-ΔSi+1-ϕDOVTTikmax⁡0,DOVTTik+ΔSi.We use ΔSi to denote the traveler’s perception error on DOVTT and use ν to denote the mean of ΔSi. Let ΔS be [ΔSi∀i]; the outline of the sampling scheme for tth sample of β, U, λ, ΔT, ΔS, μ, and ν can be described as follows.Algorithm 4.
Step  1. Consider βt~Pβ∣Ut-1,λt-1,ΔTt-1,μt-1,ΔSt-1,νt-1,X,Y.
Step  2. Consider Ut~PU∣βt,λt-1,ΔTt-1,μt-1,ΔSt-1,νt-1,X,Y.
Step  3. Consider λt~Pλ∣βt,Ut,ΔTt-1,μt-1,ΔSt-1,νt-1,X,Y.
Step  4. Consider ΔTt~PΔT∣βt,Ut,λt,μt-1,ΔSt-1,νt-1,X,Y.
Step  5. Consider ΔSt~PΔS∣βt,Ut,λt,ΔTt,μt-1,νt-1,X,Y.
Step  6. Consider μt~Pμ∣βt,Ut,λt,ΔTt,ΔSt,νt-1,X,Y.
Step  7. Consider νt~Pμ∣βt,Ut,λt,ΔTt,μt,ΔSt,X,Y.β, U, and λ can be generated along the same line of Algorithm 2. Since DOVTTik∗ and DIVTTik∗ hold the same formulation, therefore, ΔSt and νt can be sampled the same as ΔTt and μt (see Algorithm 2 and (16)).To estimate the parameters of the proposed model, we can apply the sampling algorithm to draw random samples forβ, U, λ, ΔS, ν, ΔT, and μ. These samples can be denoted as βt, Ut, λt, ΔSt, νt, ΔTt, and μt for t=1 to M, where M is the number of drawings. The estimates of the parameters can be obtained through averaging the samples. For example, the estimates of ν and μ can be obtained as follows:(22)ν^=∑t=1MνtM,μ^=∑t=1MμtM.ν^ and μ^ reflect the perception error estimate regarding DOVTT and DIVTT, respectively.


## 4. Numerical Example

In this section, we use a modal choice dataset provided by Horowitz [13] to examine the performance of the proposed models. This dataset was collected in Washington D.C. and contains 842 persons’ modal choice results (private car or public transport) for the daily trip from home to work. Table 1 provides the data structure of the dataset.Table 1
Data structure of the dataset from Horowitz (1993) [13].
Attribute
Scale
CARS
0~7
DCOST
−111~89
DOVTT
−6~48
D
V
I
T
T
−59~102
CHOICE
CAR = 1, TRANSIT = 0
Number of samples
842We first estimate the hierarchical model defined by (2) using the dataset. We conduct the estimation using the Gibbs samplers. The Gibbs samplers generate 20,000 samples for each parameter. The first 5,000 samples are dropped as the burn-in procedure. The remaining samples are used to aggregate the estimates of the parameters. We use Table 2 to show the estimation results. The sign of the estimates of β2, β3, and β4 suggests that DCOST, DOVTT, and DIVTT have a positive effect on the choice probability of car (and have a negative effect on the utilization of public transport). The value of DIC of the model is 465.764.Table 2
Estimates of the parameters of the basic model defined by (2).
Parameter
Value
Lower 95% CI
Upper 95% CI
SD
p value
β
0
−1.1731
−1.179
−1.1671
0.3023
<0.0001
β
1 (CARS)
2.3149
2.3105
2.3193
0.2233
<0.0001
β
2 (DCOST)
0.0172
0.0171
0.0173
0.0038
<0.0001
β
3 (DOVTT)
0.0620
0.0617
0.0624
0.0182
<0.0001
β
4 (DIVTT)
0.0096
0.0095
0.0098
0.0096
<0.0001
μ
−3.5245
−3.5521
−3.4970
1.9869
<0.0001
DIC
465.764
Number of observations
894We estimate the parameter of the extended model through the sampling scheme with H-M step. We draw 20,000 samples of the parameter and also treat the first 5,000 samples as the burn-in procedure. The DIC of the extended model is 461.687. This result indicates that the performance of the extended model is better than that of the basic model defined by (2). We also estimate the parameters of the logistic regression model defined by (1). The value of DIC of this model is 465.848. This value is slightly higher than the DIC of the basic model defined by (2) and clearly higher than the DIC of the extended model defined by (3).The samples of the parameters are used to investigate the shapes of the distributions of the parameters. Table3 reports the estimates of the parameters of the extended model. Figures 2 and 3 provide the histogram for β4 and μ of the extended model. The distributions of the parameters have symmetrical sides and one single peak.Table 3
Estimates of the parameters of the extended hierarchical model defined by (3).
Parameter
Value
Lower 95% CI
Upper 95% CI
SD
p value
β
0
−1.2490
−1.2531
−1.2448
0.3012
<0.0001
β
1 (CARS)
2.3322
2.3292
2.3352
0.2161
<0.0001
β
2 (DCOST)
0.0174
0.0174
0.0175
0.0038
<0.0001
β
3 (DOVTT)
0.0590
0.0587
0.0592
0.0190
<0.0001
β
4 (DIVTT)
0.0156
0.0155
0.0158
0.0101
<0.0001
μ
−1.022
−1.0234
−1.0206
0.1026
<0.0001
DIC
461.687
Number of observations
894Figure 2
The histogram forβ4.Figure 3
The histogram forμ.In addition to the improvement of the performance, the main contribution of the proposed model is that it can allow us to capture the perception error on travel time and analyze the characteristics of the perception error. Let us look at the extended model. As shown by Table3, the estimate of μ is −1.022; this result implies that people tend to weaken the difference of in-vehicle travel time when making the choice decision. The sampling schemes also draw samples of ΔTi for i=1 to I. The latent parameter ΔTi directly reflects the perception error on DIVTT of individual i.To further investigate the property ofΔTi of the extended model, we group these samples according to the value of DCARS, which is related to the demography of an individual. We calculate the mean and variance of the samples of ΔTi for each group. Figure 4 plots the mean of the samples, and Figure 5 shows the variances of the samples versus the value of DCARS. Since the values of DCARS in the dataset do not contain 6 and only one-record DCARS is equal to 7, we study the case where the values of DCARS range from 1 to 5 to avoid bias.Figure 4
Mean of the samples ofΔTi grouped by the value of CARS.Figure 5
Variance of the samples ofΔTi grouped by the value of CARS.One can find that the mean of the samples ofΔTi tends to be close to 0 as DCARS increases. On the other hand, the variance decreases with the value of DCARS. This result implies that travelers who hold more cars will be more sensitive to the difference of travel time between public transport and private car. For example, the mean of ΔTi for travelers who own 5 cars is −0.357; this value can change to −1.5623 for travelers who do not own any cars. As mentioned above, ΔTi<0 indicates that individual i subjectively weakens the difference of travel time when he/she makes the choice decision; therefore, if ΔTi<0, a smaller ΔTi can reduce the contribution of the difference of travel time on the travel mode choice decision even more.To investigate traveler’s perception error regarding both DIVTT and DOVTT, we also use Algorithm4 to estimate the model defined by (20). Table 4 shows that the DIC of the model is significantly smaller than the standard logistic model defined by (1). In the model defined by (3), traveler’s perception error regarding travel time can be characterized by the estimate ofμ. On the other hand, in the model defined by (20), traveler’s perception error regarding travel time should be characterized by the summation of the estimates of μ and ν. Looking at Table 4, one can find that the summation of the estimates of μ and ν is −1.011. This value is closed to the estimate of μ (−1.022) in the model defined by (3).Table 4
Estimates of the parameters of the extended hierarchical model defined by (20).
Parameter
Value
Lower 95% CI
Upper 95% CI
SD
p value
β
0
−1.3607
−1.8495
−0.8590
0.3005
<0.0001
β
1 (CARS)
2.2569
1.9034
2.6305
0.2131
<0.0001
β
2 (DCOST)
0.0149
0.0086
0.0213
0.0039
<0.0001
β
3 (DOVTT)
0.0907
0.0569
0.1250
0.0209
<0.0001
β
4 (DIVTT)
0.0082
−0.0072
0.0247
0.0097
0.1973
μ
−0.4310
−0.5229
−0.3363
0.0559
<0.0001
ν
−0.5801
−0.6804
−0.4785
0.0558
<0.0001
DIC
436.240
Number of observations
894


## 5. Conclusions

This study proposes a logistic regression model with a hierarchical random error term to analyze the binary choice problem. The proposed model can account for travelers’ perception errors regarding attributes. Since a number of studies have shown perception error regarding travel time to have a significant impact on modal choice, this study focuses in particular on how to capture this error from behavior data.We construct a hierarchical random error term structure in the logistic regression model though introducing a random error term for the travel time variable. In the proposed model, travel time is no longer a deterministic variable. To make the model more sensible, we also propose an extended model, where the sign of the DIVTT and DOVTT variables can be constrained. We develop a Gibbs sample to estimate the basic hierarchical model, while developing a Gibbs sampler with M-H step to estimate the parameters of the extended model. The binary choice dataset provided by Horowitz [13] is used to examine model performance and demonstrate the proposed model’s benefit. The numerical example shows that the proposed model can estimate the perception error regarding travel time exactly and that the samples of μ and ν can be used to further analyze the characteristics of this error. Although this study prefers to look at the perception error regarding travel time, the proposed methodology can also be used to investigate travelers’ perception error regarding travel cost, among others.



---

*Source: 430926-2015-08-12.xml*