基于D-optimal試驗設(shè)計的公共交通出行l(wèi)ogit模型

李 信，楊橋東，曾大堃

(西南交通大學(xué) 交通運輸與物流學(xué)院，四川成都 610031)

基于D-optimal試驗設(shè)計的公共交通出行l(wèi)ogit模型

李 信，楊橋東，曾大堃

(西南交通大學(xué) 交通運輸與物流學(xué)院，四川成都 610031)

以成都到龍泉片區(qū)的交通通道為研究對象，假設(shè)通道開通地鐵的情況下，采用D-optimal試驗設(shè)計意愿調(diào)查問卷，選擇費用、時間等變量以及個人屬性作為效用變量，建立公共交通出行的多項logit模型，結(jié)合意愿調(diào)查數(shù)據(jù)對模型進(jìn)行求解，分析通道中公共交通出行選擇行為，并預(yù)測成都到龍泉片區(qū)主要交通方式的客流分擔(dān)率。分析表明，地鐵將承擔(dān)通道中大部分的客流。

公共交通出行；D-optimal試驗設(shè)計；意愿調(diào)查；多項logit模型

離散選擇模型是描述人們選擇行為的模型，廣泛應(yīng)用于具有選擇行為的社會各個方面，在交通工程界的應(yīng)用主要有需求量預(yù)測、政策評價、交通服務(wù)效益評價等。離散選擇模型是基于效用函數(shù)最大化的假說，以Logit模型應(yīng)用最為廣泛，假設(shè)效用函數(shù)的概率項服從二重指數(shù)分布。Logit模型中需要的數(shù)據(jù)主要采用意愿調(diào)查(Stated Preference，SP)得到，意愿調(diào)查獲取假設(shè)情景下出行者的偏好或者選擇結(jié)果[1-2]。意愿調(diào)查的問卷設(shè)計方法有多種，應(yīng)用比較多的是正交設(shè)計法，但本文采用一種新的試驗設(shè)計方法——D-optimal試驗設(shè)計法進(jìn)行問卷設(shè)計，該方法的本質(zhì)是使模型的漸進(jìn)協(xié)方差矩陣的行列式最小，本質(zhì)是使模型參數(shù)的標(biāo)準(zhǔn)誤差最小，從而得到更可靠的參數(shù)估計結(jié)果，進(jìn)而能夠更精確的分析出行者的出行偏好。

1 D-optimal試驗設(shè)計

1.1 D-optimal試驗設(shè)計原理

D-optimal試驗設(shè)計的目的是使模型漸進(jìn)協(xié)方差矩陣的行列式Derror最小，則模型參數(shù)估計的標(biāo)準(zhǔn)差也就越小，參數(shù)估計的精度就越高，確保所選屬性水平的組合方案能夠最大程度體現(xiàn)受調(diào)查者出行選擇時的偏好信息最大化[3]。

多項Logit(MNL)模型的對數(shù)似然函數(shù)L為

(1)

式中yni為列矩陣，選擇枝i在選擇情景中被回答者n選擇時，yni值為1，反之為0；Pni為選擇模型中選擇枝i被出行者n選擇的概率；I為選擇枝總數(shù)；N為出行者總數(shù)。

定義的選擇模型考慮一個特殊的選擇數(shù)據(jù)集，使式(1)最大化，得到參數(shù)的最大似然估計量。

漸進(jìn)協(xié)方差矩陣Ω可以計算為Logit模型的對數(shù)似然函數(shù)L的海賽矩陣逆矩陣的負(fù)數(shù)[4-6]。海賽矩陣H中的元素計算式為

式中Xikn為選擇枝i中屬性k的值，該值不隨出行者n改變；βik為選擇枝i中屬性k的系數(shù)；K為試驗設(shè)計中需要標(biāo)定的一般參數(shù)的總數(shù)量。

1.2 D-optimal試驗設(shè)計

1)定義選擇枝集合、選擇枝屬性以及屬性水平

假設(shè)成都地鐵2#線東沿線開通，成都到龍泉片區(qū)交通出行的選擇枝集合為：大巴、公交車、出租車、私家車、地鐵。利用SP調(diào)查問卷采集出行者從成都到龍泉片區(qū)的選擇偏好的數(shù)據(jù)。選擇枝屬性定義時，應(yīng)選能夠有效反應(yīng)選擇枝特性且不具有相關(guān)性的屬性。在定義選擇枝屬性水平時，每個屬性的水平個數(shù)不需要相等，屬性水平定義越多，在效用空間中能夠得到的信息量越大，但屬性水平的個數(shù)會影響試驗組合個數(shù)，因此，應(yīng)適當(dāng)選擇[7]。選擇枝屬性水平的選取應(yīng)結(jié)合實際情況，適當(dāng)增加和減少相應(yīng)的比例，結(jié)合從成都到龍泉片區(qū)各種交通方式的實際出行時間，根據(jù)兩地距離和地鐵的運行速度計算出車內(nèi)時間，最后確定選擇枝的屬性及屬性水平定義如表1所示。

表1 選擇枝屬性及屬性水平定義

注：到站時間，指從出發(fā)點到乘車站點的時間。

表2 選擇枝屬性預(yù)先值

2)計算參數(shù)的預(yù)先值

D-optimal試驗設(shè)計法需要先計算參數(shù)的預(yù)先值。調(diào)查問卷采用正交設(shè)計，采集192組觀測結(jié)果，利用采集的數(shù)據(jù)對MNL模型進(jìn)行參數(shù)估計，參數(shù)估計結(jié)果為屬性參數(shù)的預(yù)先值，從而得到5個選擇枝中屬性參數(shù)的預(yù)先值，如表2所示。

3)設(shè)計試驗組合

根據(jù)表1定義的選擇枝屬性以及屬性水平，利用全因子設(shè)計法可以有648個試驗組合，把這648個試驗組合作為候選試驗組合。本文使 D-optimal試驗設(shè)計的試驗組合數(shù)為32個。根據(jù)D-optimal試驗設(shè)計的步驟，利用Matlab軟件編程，從具有648個組合的候選矩陣中任意挑選32個試驗組合，計算Derror，重復(fù)迭代，得出使Derror最小的32個試驗組合，該Derror為0.065。

同樣，為了減少被調(diào)查者的調(diào)查情景數(shù)，得到更可靠的調(diào)查數(shù)據(jù)，把32個組合隨機(jī)分成8塊，由8個被調(diào)查者來完成，每個被調(diào)查者完成4個組合情景，再加上個人屬性等問題的選擇，構(gòu)成完整的問卷。試驗組合示例如表3所示。

表3 D-optimal試驗設(shè)計組合示例

2 城市居民公共交通出行l(wèi)ogit模型建立

2.1模型建立

Logit模型是基于隨機(jī)效用最大化理論的假設(shè)。根據(jù)效用最大化理論，出行者總是選擇具有最大出行效用的出行方式，出行者n選擇第i種出行方式的概率Pn(i)為

Pn(i)=P(Uni≥Unj),?j∈Cn,j≠i，

式中Uni為出行方式i對于出行者n的效用；Unj為出行方式j(luò)對于出行者n的效用；Cn為選擇枝集合。

效用Uni、Unj與出行方式的屬性和出行者本身屬性相關(guān)。隨機(jī)效用理論認(rèn)為效用是一個隨機(jī)變量，并將效用函數(shù)U分為可以觀測的影響因素所構(gòu)成的固定項(V)和不可觀測的影響因素所構(gòu)成的隨機(jī)項(ε)兩大部分，并假設(shè)它們呈線性關(guān)系。因此交通方式i對出行者n的效用函數(shù)可以表示為

Uni=Vni+εni，

Vni=α+β1Xi1n+β2Xi2n+…+βkXikn+λ1Si1n+λ2Si2n+…+λmSimn，

(2)

式中α為選擇枝i的常數(shù)項；βk為每個選擇枝的效用中屬性k的權(quán)重，即選擇枝屬性的參數(shù)值；λm為每個選擇枝的效用中出行者屬性m的權(quán)重，即出行者屬性的參數(shù)值；Sinm為出行者n屬性m的值。

當(dāng)εni服從相互獨立的Gumbel分布時，可以得出Multinomial logit(MNL)模型的一般式，logit模型中出行者n選擇第i種出行方式的概率[8]為

(3)

2.2數(shù)據(jù)分析

基于D-optimal試驗設(shè)計法的一份問卷包括32個選擇情景，每個被調(diào)查者完成4個選擇情景。在實際調(diào)查中不僅采集被調(diào)查者的年齡、收入、性別等個人信息，而且被調(diào)查者必須回答在不同的費用與出行時間組合情況下首選的交通方式。每份問卷包含多個完全獨立的情景，被調(diào)查者需要對其中調(diào)查內(nèi)容做出回答。此次，調(diào)查對象為1 344人，男性占50.2%。樣本中擁有私家車的出行者占41%，樣本平均年齡為32.8歲，年齡的分布如圖1如示，其中男性和女性的平均年齡差距比較小，總體平均年齡分別是33.4歲和32.2歲。樣本總體月平均收入為3 788元，與成都市居民的月平均收入相差不大，樣本的收入分布比例如圖2所示。

圖1 總體樣本年齡分布圖

圖2 總體樣本收入分布圖

3 模型標(biāo)定結(jié)果分析

MNL模型效用參數(shù)的標(biāo)定通過軟件NLOGIT 3.0實現(xiàn)。NLOGIT軟件可以實現(xiàn)多項式選擇數(shù)據(jù)的估測、模擬和分析，已經(jīng)成為估計和模擬多項離散選擇模型的首要軟件。把調(diào)查的數(shù)據(jù)樣本數(shù)量1 344輸入標(biāo)定軟件，NLOGIT軟件利用最大似然估計法得到模型參數(shù)的標(biāo)定結(jié)果，如表4所示。由表4可得，似然函數(shù)值為1 004.118，擬含優(yōu)度為0.37。

表4 MNL模型標(biāo)定結(jié)果

注：性別、年齡、月收入分別為對應(yīng)車型駕駛員的性別、年齡、月收入。

一般離散選擇模型的擬合優(yōu)度為0.2～0.4時，代表模型的擬合度較好，所建模型的擬合優(yōu)度為0.37，模型擬合度較好，說明所建模型正確。由表4可知，廣義變量“費用”、“到站時間”、“候車時間”、“車內(nèi)時間”以及“下車到目的地時間”的參數(shù)值均為負(fù)值，表示費用越高，選擇該種交通方式的概率越小，出行時間越長，選擇該種交通方式的概率越小，與實際相符。

通過t檢驗值的絕對值判斷該參數(shù)值所對應(yīng)的屬性是否能夠顯著影響出行者的選擇行為。通常，在5%置信度區(qū)間下，如果t檢驗的絕對值>1.96，則認(rèn)為該屬性能夠顯著影響出行者的選擇行為，反之則認(rèn)為該屬性對出行者的選擇行為不會造成顯著影響。根據(jù)模型結(jié)果，變量“費用”、“到站時間”、“車內(nèi)時間”t檢驗的絕對值均>1.96，表明這3個屬性對出行者的選擇行為具有顯著影響。

精神衛(wèi)生既是全球性的重大公共衛(wèi)生問題,也是較為嚴(yán)重的社會問題。由于精神障礙患者病情波動的不可預(yù)見性,在預(yù)防控制精神障礙患者突發(fā)肇事肇禍?zhǔn)录洗嬖谥喈?dāng)困難,因此嚴(yán)重精神障礙應(yīng)急處置就顯得異常重要。應(yīng)急處置主要是通過有效的應(yīng)急救援行動,盡可能地降低事故的后果,包括人員傷亡、財產(chǎn)損失和環(huán)境破壞等。在應(yīng)急處置過程中,需要包括綜治、公安、衛(wèi)生、街道、居委會、家屬等多部門的參與,各部門要分工明確,職責(zé)分明,流程順暢,應(yīng)對及時;應(yīng)急處置后還有可能需要民政、司法、殘聯(lián)等相關(guān)部門的后續(xù)支持,因此掌握嚴(yán)重精神障礙患者應(yīng)急處置流程,明確各部門職責(zé)就顯得異常重要。

出行者的個人屬性也作為解釋變量進(jìn)入模型，由表4可知，選擇枝大巴、公交車和地鐵對應(yīng)變量“收入”的參數(shù)值均為負(fù)數(shù)，說明收入越高，選擇這幾種交通方式的概率越小。出租車對應(yīng)變量“收入”的參數(shù)值為正數(shù)，表明收入越高，選擇出租車的概率就會越大?！笆欠駬碛兴郊臆嚒睂λ郊臆嚨氖褂糜兄@著地影響，擁有私家車的出行者在出行上非常依賴于私家車。

將表4的模型標(biāo)定參數(shù)值代入式(2)可以得出5種出行方式對出行者的效用。

如：一年齡28歲的男性出行者，月收入5 600元，無私家車，針對表3中的第4種組合的出行方式，分別計算選擇5種出行方式的概率。

當(dāng)出行者性別為男時，屬性值為1；性別為女時，屬性值為0；擁有私家車時，屬性值為1，否則為0；年齡為25～30歲時，屬性值為3；收入為4 000～6 000元時，屬性值為3。將表4中相應(yīng)的數(shù)據(jù)代入式(2)，計算大巴對該出行者的效用為

V1=3.328-0.028 9×10-0.108 9×20-0.020 9×15-0.034 2×40-0.031 7×15+0.351 6×1+0.036 9×3-0.292 6×3=-1.713 0.

同理計算出公交車、出租車、私家車、地鐵對該出行者的效用分別為：-1.249 8、-2.396 5、-1.946 1、1.628 9。

同理可計算出該出行者選擇公交車、出租車、私家車、地鐵的概率分別為：0.49，0.16，0.25，0.879。這樣，計算出每個選擇情景中出行者對每種交通方式的選擇概率，再計算每種交通方式選擇概率的平均值，得到通道中大巴、公交車、出租車、私家車、地鐵5種交通方式的客流分擔(dān)率分別為2.5%、22.6%、1.2%、25.6%、48.1%?？梢姷罔F對旅客具有強(qiáng)大的吸引力，將承擔(dān)通道內(nèi)大部分的客流。

4 結(jié)論

基于D-optimal 試驗設(shè)計的公共交通logit模型的擬合優(yōu)度為0.37，擬合優(yōu)度較好，說明所建模型的正確性。模型求解結(jié)果表明：

1)費用、到站時間以及車內(nèi)時間顯著影響出行者的出行行為，通過提高公交車的站點覆蓋率、縮短公交線路的非直線系數(shù)等手段來縮短出行者的到站時間和車內(nèi)時間，可以提高公交車的選擇概率。

2)擁有私家車的出行者對私家車依賴度非常高，因此可以采用一些交通需求管理政策控制小汽車的出行量，以達(dá)到緩解交通擁擠的目的。

3)成都到龍泉片區(qū)地鐵建成后，成都到龍泉片區(qū)大巴、公交車、出租車、私家車、地鐵5種交通方式的客流分擔(dān)率分別為2.5%、22.6%、1.2%、25.6%、48.1%，地鐵將承擔(dān)成都到龍泉片區(qū)的大部分客流。

[1]趙鵬,藤原章正.SP調(diào)查方法在交通預(yù)測中的應(yīng)用[J].北京交通大學(xué)學(xué)報,2000,24(6):29-32.

[2]焦朋朋,陸化普.基于意向調(diào)查數(shù)據(jù)的非集計模型研究[J].公路交通科技,2005,22(6):114-116.

[3]Michiel C J Bliemer, John M Rose. Experimental Design Influences on Stated Choice Outputs: An Empirical Study in Air Travel Choice[J].Transportation Research: Part A , 2011,45:63-79.

[4]John M Rose, Michiel C J Bliemer, David A Hensher, et al.Designing E Cient Stated Choice Experiments in the Presence of Reference Alternatives[J].Transportation Research: Part B,2008,42(4):395-406.

[5]Rose J M, Bliemer M C J, Hensher D A ,er al. Handbook of Transport Modelling[K].Oxford:Oxford University Press,2008:151-180.

[6]Kenneth Train, Wesley W Wilson. Estimation on Stated-Preference Experiments Constructed from Revealed-Preference Choices[J].Transportation Research： Part B, 2008,42(3): 191-203.

[7]David Alan Hensher, William Greene. Applied Choice Analysis [M].England：Cambridge University Press,2005.

[8]關(guān)宏志.非集計模型——交通行為分析的工具[M].北京:人民交通出版社,2004.

(責(zé)任編輯：楊秀紅)

ApplicationofD-OptimalDesignMethodforLogitModelConstructionofPublicTransportTravelBehavior

LIXin,YANGQiao-dong，ZENGDa-kun

(CollegeofTransportationandLogistics,SouthwestJiaotongUniversity,Chengdu610031,China)

A traffic corridor between Chengdu and Longquan district is selected to analyze the public transport travel behavior. Assuming that the subway is opened in the traffic corridor, the D-optimal design is adopted in the stated preference questionnaires. Then, such utility variables as time, cost, sex, age and income are selected to construct the multinomial Logit (MNL) model and the stated preference data are used to calibrate the Logit model to analyze peoples′ public transport behavior. Finally, the sharing rate of main traffic modes is forecast. As a result, the analysis shows that the subway will undertake the main part of passenger flow in the traffic corridor.

public transport travel behavior; D-optimal design;stated preference; multinomial logit model

2013-03-03

李 信(1990—)，女，云南彌勒人，西南交通大學(xué)碩士研究生，主要研究方向為交通運輸規(guī)劃與管理、交通規(guī)劃與交通行為理論與實證.

10.3969/j.issn.1672-0032.2014.01.008

U491

A

1672-0032(2014)01-0033-06