• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    城市購物出發(fā)時間的非參數(shù)生存分析

    2014-06-05 15:27:39李明楊小寶環(huán)梅董苓
    山東科學(xué) 2014年4期
    關(guān)鍵詞:生存率購物變量

    李明,楊小寶,環(huán)梅,董苓

    (北京交通大學(xué)交通運輸學(xué)院,北京 100044)

    城市購物出發(fā)時間的非參數(shù)生存分析

    李明,楊小寶,環(huán)梅,董苓

    (北京交通大學(xué)交通運輸學(xué)院,北京 100044)

    根據(jù)濟(jì)南市2011年的出行數(shù)據(jù),研究了該市城市購物出發(fā)時間的分布特性。采用生存分析方法,建立城市購物出發(fā)時間的持續(xù)時間模型。運用非參數(shù)方法,對購物出發(fā)的持續(xù)時間進(jìn)行估計,討論了城市購物出行的時間分布和家庭社會經(jīng)濟(jì)屬性變量對購物出發(fā)時間的影響。結(jié)果表明,總體樣本的5.4%在7:00之前出發(fā)購物,69.7%在7:00~10:00出發(fā)購物,在10:00之后出發(fā)購物的少于24.9%。性別、年齡和家庭中兒童的數(shù)量對購物出發(fā)時間的選擇有著顯著的影響。女性比男性購物出發(fā)早,老年人更易在非高峰期進(jìn)行購物出行,沒有兒童的家庭購物出行在時間上更加靈活。本研究為購物出行的定量研究提供了準(zhǔn)確、有效的分析工具,也對通過動態(tài)價格機(jī)制控制交通擁堵有重要意義。

    城市購物出發(fā)時間;連續(xù)時間模型;非參數(shù)生存分析

    在城市發(fā)展過程中,交通系統(tǒng)和城市發(fā)展相互影響,相互制約。而個人的出行選擇對城市交通有非常重要的影響,對個體出行行為的研究,是進(jìn)行城市交通規(guī)劃、建設(shè)和管理的一項不可缺少的基礎(chǔ)工作[1]。研究個人出行的目的是為了改善個人出行計劃和提出有效的交通管理措施。根據(jù)基于活動的出行理論,出行可以看作為到達(dá)活動空間的必要方式,這意味著出行時間的重要性[2]。因此,在構(gòu)造與出行有關(guān)的模型時,活動的出發(fā)時間是一個重要的因素。

    現(xiàn)有對出行的研究主要集中在對通勤者的研究方面,通勤者的出發(fā)時間一般是在一天中的高峰擁堵時段。通勤出行是城市居民最基本和最重要的出行,影響到城市生活的方方面面,特別是早、晚高峰通勤時段的交通擁堵,已成為城市交通最為突出的問題之一,并制約著城市社會和經(jīng)濟(jì)的進(jìn)一步發(fā)展。然而,目前對非通勤者出行研究較少,由于非通勤出行在城市出行中所占比重越來越大,且非通勤出行比通勤出行時間上更加靈活,這意味著社會人口學(xué)變量的變化和交通控制措施對非通勤出行時間的影響要多于不太靈活的通勤出行,所以有必要對非通勤者出行進(jìn)行深入研究。以前的研究多采用logit離散模型,將時間分成若干時段[3-6],這樣的方法雖具有一定的科學(xué)性,但是將出發(fā)時間作為離散的變量,而不是連續(xù)的變量。另外,當(dāng)劃分時間的方法不同時,就會出現(xiàn)不同的結(jié)果,很難對交通狀況做出正確的評價。

    在眾多非通勤出行活動中,購物出行是一個重要的出行目的。本文主要對非通勤出行人群中的購物人群的出發(fā)時間進(jìn)行研究?;跐?jì)南市2011年工作日一天中居民出行的調(diào)查數(shù)據(jù),主要探討了購物出發(fā)時間的整體分布特征,利用生存分析方法,構(gòu)造了購物出發(fā)時間的選擇模型,并從變量的顯著性檢驗方面,討論并分析了對購物出發(fā)時間有顯著影響的社會經(jīng)濟(jì)屬性變量。本研究將為制定擁擠收費、錯時上班等交通需求管理措施,提供重要的參考依據(jù)。

    1 模型和方法

    1.1 生存分析風(fēng)險模型

    生存分析(survival analysis)是一種根據(jù)實驗或調(diào)查數(shù)據(jù),將事件的結(jié)果(終點事件)和出現(xiàn)這一結(jié)果所經(jīng)歷的時間結(jié)合起來分析的一種統(tǒng)計分析方法,也稱之為風(fēng)險模型或持續(xù)模型(hazard model/duration model)[7]。目前生存分析已廣泛應(yīng)用于生物工程、社會科學(xué)等多個領(lǐng)域,近年來也逐步被運用到城市交通研究中,部分學(xué)者運用生存分析方法研究了交通事故的清理時間[8]、基于活動的出行行為[9]以及機(jī)非混行的通行時間[10]等問題。

    廣義的生存時間指生物體存活的時間,或所關(guān)心的某種現(xiàn)象的持續(xù)時間。城市購物出發(fā)時間指從凌晨開始,一直到出發(fā)去購物為止的持續(xù)時間[11],屬于廣義生存時間范疇,可運用生存分析方法來對其進(jìn)行研究。設(shè)T是非負(fù)的隨機(jī)變量,表示的是從凌晨到出發(fā)去購物的持續(xù)時間。令f(t)表示T的概率密度函數(shù),則T的分布函數(shù)為

    生存函數(shù)S(t),也叫累積生存率,表示生存時間大于t的概率。表達(dá)式為

    本研究中,由于沒有刪失觀測值,生存函數(shù)可用在時刻t之前仍沒有出去購物的人數(shù)占總樣本的比例來估計,即

    生存分析中,T的風(fēng)險函數(shù)也叫危險率函數(shù),指事件在持續(xù)到t時刻時沒有結(jié)束,但在接下來的一段很小的時間Δt內(nèi)結(jié)束的概率。用公式表示如下

    在本文中,風(fēng)險函數(shù)指在時刻t之前沒有出去購物,但在接下來的一段很小的時間Δt內(nèi)出去購物的概率,可用下式來估計

    其中N[t,t+Δt]表示在區(qū)間[t,t+Δt]內(nèi)出發(fā)購物的人數(shù);Nbefor[t]表示在t時刻前未出發(fā)購物的人數(shù)。

    風(fēng)險函數(shù)與密度函數(shù)和生存函數(shù)之間有如下關(guān)系

    持續(xù)時間模型用風(fēng)險函數(shù)和生存函數(shù)來描述事件的持續(xù)發(fā)展過程。從式(6)可知,風(fēng)險函數(shù)是生存函數(shù)變化率的相反數(shù),反映了事件持續(xù)過程的時間變化規(guī)律。

    1.2 估計生存函數(shù)的非參數(shù)方法

    生存分析通常采用生存函數(shù)、風(fēng)險函數(shù)來描述,估計生存函數(shù)的方法有參數(shù)方法和非參數(shù)方法。當(dāng)生存時間所服從的理論分布已知時用參數(shù)方法。當(dāng)生存時間所服從的理論分布的類型未知時,非參數(shù)方法比參數(shù)方法有較高的效率,對數(shù)據(jù)的分布沒有特殊的要求。參數(shù)模型的一個缺點是當(dāng)參數(shù)形式不正確的時候,估計的風(fēng)險函數(shù)也不正確。非參數(shù)模型的一個優(yōu)點是不論資料是什么樣的分布形式,只需要根據(jù)樣本提供的順序統(tǒng)計量對生存率進(jìn)行估計。本研究中,我們采用非參數(shù)法,即Kaplan-Meier乘積限方法,來估計購物出發(fā)時間的生存函數(shù)。以城市購物出行行為為研究對象,假定共有n個出發(fā)時間樣本,這些樣本的出發(fā)時間共有k(k≤n)個不同的取值,由于沒有刪失數(shù)據(jù),直接將它們從小到大排序t1<t2<…<tk。令dj為tj-1到tj時段內(nèi)出發(fā)去購物的樣本數(shù),nj為各時刻初未去購物的樣本數(shù)。生存函數(shù)S(t)的乘積限估計可用下式表示

    2 數(shù)據(jù)與實證分析

    2.1 數(shù)據(jù)的提取和變量的定義

    本研究以濟(jì)南市2011年居民出行調(diào)查數(shù)據(jù)庫為數(shù)據(jù)來源,調(diào)查的是一個工作日,一天中(0:00~20:00)的活動,調(diào)查包括家庭中所有成員的出行日志。出行日志包括出行的目的、出行的方式、出發(fā)的時刻、出發(fā)地、到達(dá)地、出行方式、所居住小區(qū)的性質(zhì)等。另外,調(diào)查還包括個人和家庭的社會經(jīng)濟(jì)屬性變量,如年齡、性別、月交通支出、月收入和職業(yè)類型等。

    文章重點分析城市購物出行的出發(fā)時間,樣本來自1 246個外出購物者。表1提供的是購物出發(fā)相關(guān)變量及其定義。

    表1 變量的定義及取值Table 1 Definition and value of the variables

    2.2 購物出發(fā)時間的總體分布

    運用1中介紹的生存函數(shù)非參數(shù)估計方法,可求出購物出發(fā)時間的生存函數(shù)及其風(fēng)險函數(shù)。樣本的最早出發(fā)時間為2:00。表2給出了樣本的出發(fā)時間生存函數(shù)和風(fēng)險函數(shù)的估計值,圖1則給出了生存函數(shù)和風(fēng)險函數(shù)隨著購物出發(fā)時間的變化規(guī)律。

    表2 購物出發(fā)時間的分析Table 2 Analysis of shopping departure time

    由圖1a生存函數(shù)曲線可知,在7:00之前,生存率變化較平穩(wěn),占到總體樣本的5.4%,表明在此時段出去購物的人很少。在7:00~10:00之間生存率急劇下降,樣本占到總體的69.7%,說明購物人群中有近70%的比例選擇在這一時段內(nèi)出發(fā)。在10:00之后變化較平穩(wěn),占到總體樣本的24.9%。由圖1b風(fēng)險率函數(shù)可知,該函數(shù)不是單調(diào)的,在8:00~10:00之間出現(xiàn)購物高峰,加入購物的概率比較大。在10:00~16:00風(fēng)險率較平穩(wěn),說明10:00沒有出去購物的人在這一時段內(nèi)出去購物的可能性不大。在16:00之后,隨著時間的延長,風(fēng)險率函數(shù)趨勢是上升的,說明在16:00前仍未去購物的人群,在這之后出去購物的可能性越來越大。

    圖1 購物出發(fā)時間的生存函數(shù)和風(fēng)險率函數(shù)曲線Fig.1 Survival and dangerous functions curves of shopping departure time

    表3給出了樣本出行時間分布情況,T(min)表示0:00到出發(fā)去購物的持續(xù)時間。由表3可知,75%分位點為600 min(10:00),表明很大比例的購物人群是在10:00之前出去的。

    表3 購物出發(fā)時間統(tǒng)計Table 3 Statistics of shopping departure time

    2.3 變量的影響

    除了購物出發(fā)時間的整體分布規(guī)律之外,各屬性變量對購物出發(fā)時間分布也有一定影響。本文運用生存分析的KM方法來分析各變量對購物出發(fā)時間的影響,通過生存分析中的非參數(shù)方法對各變量在不同取值下的生存率和風(fēng)險率函數(shù)進(jìn)行比較,可以得出各變量對出發(fā)時間選擇的影響。由于7:00~10:00的樣本占到總體的將近70%,本文將重點針對這一時間段內(nèi)購物人群的出發(fā)時間及其影響因素進(jìn)行分析,見圖2~4。

    圖2 性別對購物出發(fā)時間的生存函數(shù)和風(fēng)險函數(shù)的影響Fig.2 Impact of gender on survival and dangerous functions of shopping departure time

    圖3 年齡對購物出發(fā)時間的生存函數(shù)和風(fēng)險函數(shù)的影響Fig.3 Impact of age on survival and dangerous functions of shopping departure time

    圖4 家庭中兒童數(shù)量對購物出發(fā)時間的生存函數(shù)和風(fēng)險函數(shù)的影響Fig.4 Impact of children amount on survival and dangerous functions of shopping departure time

    由圖2~4可以看出,性別、年齡和家中有無兒童對購物出發(fā)時間有顯著影響。從不同性別的生存函數(shù)曲線(圖2)可以看出,男性外出購物的生存函數(shù)一般比女性的低,這說明女性的購物出行活動總體比男性更多。另外,從風(fēng)險函數(shù)曲線可以看出,男性外出購物的風(fēng)險率變化較為平緩,說明男性在7:00~10:00之間購物出行的時間分布上比較均勻;而女性外出購物的風(fēng)險函數(shù)的變化較為明顯,在8:50前后有個明顯的高峰值,對應(yīng)的生存函數(shù)曲線上這一時段內(nèi)生存率有急劇的下降趨勢(由50%迅速下降到20%),表明樣本中有近30%的女性選擇在8:50前后出發(fā)去購物。

    從不同年齡段的生存曲線(圖3)可以看出,8:40之前,老年人外出購物的生存率曲線一般低于中青年,表明在這一時段內(nèi)老年人出外購物的比例比中青年高,這可能與更多老年人有逛早市的習(xí)慣有關(guān)。此外,在7:40~8:50之間老年人外出購物的生存率曲線下降得很慢,表明在這一時段內(nèi)老年人出外購物的比例非常低。從風(fēng)險率曲線可以看出,中青年外出購物的曲線變化較為平緩,說明他們購物的出發(fā)時間分布更為均勻,而老年人外出購物的風(fēng)險率函數(shù)中有幾個明顯的峰值,對應(yīng)在生存率曲線上,老年人外出購物在7:40和8:50前后都有急劇的下降趨勢,表明老年人的購物活動大多會選擇在上班高峰期(7:40~8:50)之前或之后出行,以此避開上班高峰期。

    家庭中兒童的數(shù)量也是影響購物出發(fā)時間的另一個重要因素。由圖4可以看出,有兒童的家庭要比沒有兒童的外出購物生存率要平穩(wěn),沒有兒童的家庭外出購物在7:40和8:50出現(xiàn)生存率急劇下降,表明沒有兒童的家庭在這兩個時點前后出外購物的比例較高,以此避開上班高峰期,他們在購物出行時間選擇上比有兒童的家庭更加靈活。家庭中有兒童的會在8:30左右出現(xiàn)第一個購物高峰,可能是由于這段時間內(nèi)老人或家庭婦女帶領(lǐng)兒童出來購物的比較多。

    3 結(jié)論

    本文通過非參數(shù)生存分析的方法分析了城市購物出發(fā)的時間分布,討論了購物出發(fā)時間隨出行者社會經(jīng)濟(jì)屬性的連續(xù)分布規(guī)律,建立的出發(fā)時間的持續(xù)時間模型可以更好地分析購物出行的有效數(shù)據(jù),為出行時間的定量研究提供了更準(zhǔn)確有效的分析工具,對于準(zhǔn)確把握居民購物出行行為有重要意義。利用模型得到的影響購物出發(fā)時間的顯著因素主要是個人和家庭屬性,而未考慮實時的道路狀況[12],今后應(yīng)該對這方面做出更廣泛的估計,包括更廣泛的社會資源和爭議性的變量。

    [1]雋志才,鮮于建川.基于持續(xù)時間模型的日出行時間分析[J].交通運輸系統(tǒng)工程與信息,2010,10(4):62-65.

    [2]陳團(tuán)生.通勤者出行行為特征與分析方法研究[D].北京:北京交通大學(xué),2007.

    [3]ARENTZE,TA,TIMMERMANS,H JP.A Learning based transportation oriented simulation system.[J].Transportation research part B:Methodological,2004,38(7):613-633.

    [4]ETTEMA D,BASTIN F,POLAK J,et al.Modelling the joint choice of activity timing and duration[J].Transportation Research Part A,2007,41(9):827-841.

    [5]ETTEMA D,ASHIRU O,POLAK JW.Modeling tim ing and duration of activities and trips in response to road-pricing policies[J].Transportation Research Record:Journal of the Transportation Research Board,2004,1894:1-10.

    [6]LEE E T,WANG J.Statisticalmethods for survival data analysis[M].New York:John Wiley&Sons,Inc,2003.

    [7]CHUNG Y.Development of an accident duration prediction model on the Korean Freeway Systems[J].Accident Analysis and Prevention,2010,42(1):282-289.

    [8]van den BERG P,ARENTZE T,TIMMERMANSH.A latent class accelerated hazard model of social activity duration[J]. Transportation Research Part A:Policy and Practice,2012,46(1):12-21

    [9]GUO H W,GAO Z Y,YANG X B,et al.Modeling travel time underthe influence of on-street parking[J].Journal of Transportation Engineering,2012,138(2):229-235.

    [10]YANG X B,HUAN M,GUO H W,et al.Car travel time estimation near a bus stop with non-motorized vehicles[J]. International Journal of Computational Intelligence Systems,2011,4(6):1350-1357.

    [11]BHATC R,STEED J L.A continuous-timemodel of departure time choice for urban shopping trips[J].Transportation research part B:Methodological,2002,36(3):207-224.

    [12]GADDA S,KOCKELMAN KM,DAMIEN P,etal.Contimuous departure timemodels[J].Transportation Research Record:Journal of the Transportation Research Board,2009,2132:13-24.

    Nonparametric survival analysis of departure time of urban shopping trips

    LI Ming,YANG Xiao-bao,HUAN Mei,DONG Ling
    (School of Transportation and Traffic,Beijing Jiaotong University,Beijing 100044,China)

    We addressed distribution characteristic of departure time of urban shopping trips based on the travel data of Jinan in 2011.We established its duration model with survival analysis method.We also estimated the duration time of shopping trips with a nonparametric method.We further discussed the impacts of departure time distribution of shopping trips and individual social-demographic properties on departure time of urban shopping trips.Results show that 5.4% of the total samples depart for shopping before 7:00,69.7% between 7:00 and 10:00,and less than 24.9% after 10:00.Gender,age and children amount of a family have significant impacts on departure time:females are early than males in departure time;elders are easier to go shopping at non-peak hours;departure time is more flexible for a no-child family.Our results can provide an accurate and effec tive analysis tool for quantitative investigation of departure time of shopping trips and have guiding significance for the controlof traffic congestion through dynamic pricing mechanism.

    departure time of urban shopping trips;continuous time model;nonparametric survival analysis

    U491

    A

    1002-4026(2014)04-0085-07

    10.3976/j.issn.1002-4026.2014.04.015

    2013-11-01

    國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)(2012CB725400);國家自然科學(xué)基金(70901005;71071016;71131001)

    李明(1990-),女,碩士,研究方向為智能交通工程。Email:12120965@bjtu.edu.cn

    猜你喜歡
    生存率購物變量
    抓住不變量解題
    也談分離變量
    “五年生存率”不等于只能活五年
    人工智能助力卵巢癌生存率預(yù)測
    “五年生存率”≠只能活五年
    圣誕購物季
    快樂六一,開心購物!
    不可錯過的“購物”APP
    Coco薇(2015年5期)2016-03-29 23:34:28
    HER2 表達(dá)強(qiáng)度對三陰性乳腺癌無病生存率的影響
    SL(3,3n)和SU(3,3n)的第一Cartan不變量
    九色成人免费人妻av| 免费观看性生交大片5| 18禁在线无遮挡免费观看视频| 丰满少妇做爰视频| 国产片内射在线| 免费播放大片免费观看视频在线观看| 亚洲人成网站在线观看播放| 亚洲精品第二区| 亚洲国产最新在线播放| 国产免费视频播放在线视频| 女的被弄到高潮叫床怎么办| 免费大片黄手机在线观看| 18在线观看网站| 亚洲综合色惰| 国产一区亚洲一区在线观看| 纵有疾风起免费观看全集完整版| 中国美白少妇内射xxxbb| 超碰97精品在线观看| 在线免费观看不下载黄p国产| 在线 av 中文字幕| 国产高清三级在线| 我的老师免费观看完整版| 色婷婷av一区二区三区视频| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看三级黄色| 日本与韩国留学比较| 精品国产露脸久久av麻豆| 麻豆成人av视频| 蜜桃在线观看..| 国产精品人妻久久久久久| 午夜福利网站1000一区二区三区| 一级片'在线观看视频| 99九九在线精品视频| 最近2019中文字幕mv第一页| 久久精品久久久久久久性| 大香蕉97超碰在线| 777米奇影视久久| 国产 一区精品| 午夜日本视频在线| 亚洲激情五月婷婷啪啪| 男女免费视频国产| 国产一区二区在线观看av| 亚洲,一卡二卡三卡| 亚洲欧美一区二区三区国产| 性色av一级| 欧美日韩国产mv在线观看视频| √禁漫天堂资源中文www| 九九久久精品国产亚洲av麻豆| 日韩不卡一区二区三区视频在线| 亚洲天堂av无毛| 亚洲怡红院男人天堂| 一二三四中文在线观看免费高清| 精品亚洲成a人片在线观看| 国产男女超爽视频在线观看| 久久久久久久久大av| 黑人猛操日本美女一级片| 97在线视频观看| 我的女老师完整版在线观看| 26uuu在线亚洲综合色| 91国产中文字幕| 亚洲精品,欧美精品| 五月伊人婷婷丁香| 亚洲婷婷狠狠爱综合网| 天天影视国产精品| 久久久午夜欧美精品| 国产成人免费观看mmmm| 久久久久视频综合| 国产精品一区二区在线不卡| xxxhd国产人妻xxx| 好男人视频免费观看在线| 国产精品嫩草影院av在线观看| 国产日韩欧美视频二区| 天天影视国产精品| 国产免费一级a男人的天堂| 国产欧美另类精品又又久久亚洲欧美| 99热6这里只有精品| 18+在线观看网站| 只有这里有精品99| 日韩成人av中文字幕在线观看| 妹子高潮喷水视频| 久久久久久久久久人人人人人人| freevideosex欧美| 一本—道久久a久久精品蜜桃钙片| 久久ye,这里只有精品| 老司机影院成人| 青春草视频在线免费观看| 中文天堂在线官网| 色视频在线一区二区三区| 久热久热在线精品观看| 国产熟女欧美一区二区| 男女国产视频网站| 在线免费观看不下载黄p国产| 在线播放无遮挡| 国产一区二区三区av在线| 少妇被粗大猛烈的视频| 国产在视频线精品| 久热这里只有精品99| 亚洲精品视频女| 中国国产av一级| 国产69精品久久久久777片| 亚洲在久久综合| av福利片在线| 五月天丁香电影| 国产淫语在线视频| 制服诱惑二区| 国产精品久久久久久av不卡| 国产色婷婷99| 一级,二级,三级黄色视频| 色婷婷av一区二区三区视频| 中文字幕av电影在线播放| 国内精品宾馆在线| 伊人久久国产一区二区| 久热这里只有精品99| 亚洲成人av在线免费| 男女啪啪激烈高潮av片| 女人精品久久久久毛片| av女优亚洲男人天堂| 久久毛片免费看一区二区三区| 蜜臀久久99精品久久宅男| 人人澡人人妻人| 超碰97精品在线观看| 看免费成人av毛片| 国产男人的电影天堂91| 香蕉精品网在线| 99久久精品一区二区三区| .国产精品久久| 91aial.com中文字幕在线观看| 在线观看三级黄色| 一本久久精品| 国产爽快片一区二区三区| 精品一区二区三卡| 蜜桃在线观看..| 免费av不卡在线播放| 成人手机av| 日日撸夜夜添| 热99国产精品久久久久久7| 久久这里有精品视频免费| 超色免费av| 久久97久久精品| 午夜精品国产一区二区电影| 国产女主播在线喷水免费视频网站| 亚洲精品国产色婷婷电影| 国产精品久久久久久精品古装| 色94色欧美一区二区| 久久久久视频综合| 在线免费观看不下载黄p国产| 91精品三级在线观看| 一区二区av电影网| 国产 精品1| 人人妻人人爽人人添夜夜欢视频| av又黄又爽大尺度在线免费看| 日韩 亚洲 欧美在线| 国产探花极品一区二区| 在线观看三级黄色| 在线亚洲精品国产二区图片欧美 | 多毛熟女@视频| 免费看av在线观看网站| 亚洲人成77777在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 麻豆乱淫一区二区| 国产熟女欧美一区二区| 国产综合精华液| 久久久久久久久久久丰满| 三级国产精品欧美在线观看| 亚洲精品久久午夜乱码| 亚洲丝袜综合中文字幕| 亚洲av成人精品一二三区| 在线 av 中文字幕| 女人久久www免费人成看片| 黑人高潮一二区| 少妇熟女欧美另类| 国产精品久久久久久久电影| 免费av中文字幕在线| 国产精品久久久久久精品古装| 亚洲成人av在线免费| 夫妻午夜视频| 成年av动漫网址| 永久免费av网站大全| 天堂8中文在线网| 伊人久久国产一区二区| 久久99精品国语久久久| 欧美 亚洲 国产 日韩一| 国产又色又爽无遮挡免| 日本午夜av视频| 王馨瑶露胸无遮挡在线观看| 日日摸夜夜添夜夜添av毛片| 精品久久蜜臀av无| 九九爱精品视频在线观看| 亚洲av日韩在线播放| 久久综合国产亚洲精品| 成人综合一区亚洲| 日本wwww免费看| 国产精品久久久久久精品电影小说| 欧美激情国产日韩精品一区| 少妇的逼好多水| 99热全是精品| 少妇被粗大的猛进出69影院 | 在线观看免费视频网站a站| 亚洲国产日韩一区二区| 天天操日日干夜夜撸| 少妇被粗大猛烈的视频| 9色porny在线观看| 热99国产精品久久久久久7| 在线亚洲精品国产二区图片欧美 | 免费大片18禁| 日韩一区二区三区影片| 国产av一区二区精品久久| 国产永久视频网站| 午夜久久久在线观看| 美女cb高潮喷水在线观看| 丝瓜视频免费看黄片| 九草在线视频观看| av播播在线观看一区| 亚洲精品国产色婷婷电影| 超碰97精品在线观看| 久久这里有精品视频免费| 欧美精品一区二区大全| 国国产精品蜜臀av免费| 狂野欧美激情性bbbbbb| 考比视频在线观看| 亚洲精品久久久久久婷婷小说| 日本猛色少妇xxxxx猛交久久| 亚洲国产成人一精品久久久| 这个男人来自地球电影免费观看 | 纵有疾风起免费观看全集完整版| 丰满少妇做爰视频| av有码第一页| 人人妻人人添人人爽欧美一区卜| 亚洲丝袜综合中文字幕| 草草在线视频免费看| 水蜜桃什么品种好| 18在线观看网站| 丰满少妇做爰视频| 伦理电影大哥的女人| 欧美少妇被猛烈插入视频| 精品亚洲乱码少妇综合久久| 亚洲一区二区三区欧美精品| 国产片特级美女逼逼视频| 大香蕉久久网| 18禁在线播放成人免费| 国产免费又黄又爽又色| 99热国产这里只有精品6| 99久国产av精品国产电影| 久久久久久久久久成人| 国产精品国产三级国产av玫瑰| 超碰97精品在线观看| 三级国产精品欧美在线观看| 一区二区三区乱码不卡18| 欧美老熟妇乱子伦牲交| 日韩av免费高清视频| 中文字幕免费在线视频6| 22中文网久久字幕| 日韩亚洲欧美综合| 国产精品久久久久久久电影| 大片免费播放器 马上看| 日日摸夜夜添夜夜添av毛片| 久久人人爽人人爽人人片va| av在线app专区| 免费人妻精品一区二区三区视频| 热99国产精品久久久久久7| 伦理电影大哥的女人| av天堂久久9| 欧美精品高潮呻吟av久久| av免费观看日本| 最近2019中文字幕mv第一页| 亚洲精品aⅴ在线观看| 久久韩国三级中文字幕| 亚洲精品美女久久av网站| 日韩av免费高清视频| 十八禁高潮呻吟视频| 成人手机av| 国产精品一区www在线观看| 久久久久久久久久人人人人人人| 久久韩国三级中文字幕| 欧美老熟妇乱子伦牲交| 国产片内射在线| 国产成人精品一,二区| 亚洲欧美成人综合另类久久久| 国产乱来视频区| 国产免费一级a男人的天堂| 国产高清有码在线观看视频| 日日啪夜夜爽| 亚洲欧美一区二区三区黑人 | 女人久久www免费人成看片| 亚洲综合色惰| 高清欧美精品videossex| 日本欧美视频一区| 成年美女黄网站色视频大全免费 | 国产精品麻豆人妻色哟哟久久| 成年美女黄网站色视频大全免费 | 纵有疾风起免费观看全集完整版| 99久国产av精品国产电影| 成人亚洲精品一区在线观看| 啦啦啦视频在线资源免费观看| 国产精品不卡视频一区二区| 一边亲一边摸免费视频| 亚洲第一av免费看| 亚洲精品乱久久久久久| 丝瓜视频免费看黄片| 日本色播在线视频| 亚洲精品456在线播放app| 日韩av在线免费看完整版不卡| 欧美 亚洲 国产 日韩一| 国产成人精品在线电影| 91久久精品国产一区二区成人| 高清视频免费观看一区二区| 久久人人爽av亚洲精品天堂| 91精品一卡2卡3卡4卡| 国产亚洲欧美精品永久| 97在线人人人人妻| 日日撸夜夜添| 国产亚洲av片在线观看秒播厂| 国产国语露脸激情在线看| 国产69精品久久久久777片| 亚洲精品自拍成人| 狂野欧美激情性bbbbbb| 亚洲精品视频女| a 毛片基地| 日本黄大片高清| 能在线免费看毛片的网站| 精品酒店卫生间| 日本av手机在线免费观看| 国产成人午夜福利电影在线观看| 日韩强制内射视频| 日韩成人伦理影院| 人人澡人人妻人| 国产成人午夜福利电影在线观看| 国产高清国产精品国产三级| 免费看不卡的av| 一边亲一边摸免费视频| 午夜免费观看性视频| 午夜激情久久久久久久| 美女中出高潮动态图| 97超视频在线观看视频| 亚洲av成人精品一区久久| 十分钟在线观看高清视频www| 乱人伦中国视频| 又大又黄又爽视频免费| 黑人巨大精品欧美一区二区蜜桃 | 亚洲成人手机| 久久国产精品大桥未久av| 国产成人a∨麻豆精品| 亚洲精品美女久久av网站| 视频中文字幕在线观看| 国产精品99久久久久久久久| 校园人妻丝袜中文字幕| 岛国毛片在线播放| 国产精品久久久久久av不卡| 夫妻性生交免费视频一级片| 午夜老司机福利剧场| 老司机影院毛片| 国产精品久久久久久av不卡| 大陆偷拍与自拍| 精品少妇内射三级| 免费高清在线观看日韩| 亚洲精品国产av成人精品| 国产熟女欧美一区二区| 人妻人人澡人人爽人人| 狠狠婷婷综合久久久久久88av| 最后的刺客免费高清国语| 老司机影院毛片| videossex国产| 亚洲成人av在线免费| 成人无遮挡网站| 免费高清在线观看日韩| 一本大道久久a久久精品| 国产有黄有色有爽视频| av不卡在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av免费高清在线观看| 91精品一卡2卡3卡4卡| 五月伊人婷婷丁香| 免费观看的影片在线观看| 亚洲精华国产精华液的使用体验| 视频区图区小说| 国产黄色视频一区二区在线观看| 国产精品一区二区三区四区免费观看| 校园人妻丝袜中文字幕| 少妇高潮的动态图| 精品久久久精品久久久| 男女免费视频国产| av电影中文网址| 人妻一区二区av| 亚洲一级一片aⅴ在线观看| 欧美丝袜亚洲另类| 国产成人精品久久久久久| 亚洲情色 制服丝袜| 特大巨黑吊av在线直播| 午夜福利在线观看免费完整高清在| 久久人人爽av亚洲精品天堂| 午夜影院在线不卡| 亚洲熟女精品中文字幕| 亚洲欧美一区二区三区国产| 熟女电影av网| 欧美日韩视频高清一区二区三区二| 日本91视频免费播放| 三级国产精品欧美在线观看| 免费观看无遮挡的男女| 热99久久久久精品小说推荐| 国产综合精华液| 成年人免费黄色播放视频| 在线观看免费高清a一片| 两个人免费观看高清视频| 男女边吃奶边做爰视频| 国产成人一区二区在线| 日日摸夜夜添夜夜爱| 伦理电影大哥的女人| 丝袜美足系列| 寂寞人妻少妇视频99o| 久久久久网色| 这个男人来自地球电影免费观看 | 亚洲欧洲国产日韩| 久久久a久久爽久久v久久| 欧美精品国产亚洲| 日韩在线高清观看一区二区三区| 久久精品久久精品一区二区三区| 亚洲成色77777| 国语对白做爰xxxⅹ性视频网站| 中文字幕人妻丝袜制服| 午夜激情av网站| 亚洲av中文av极速乱| 在线观看www视频免费| 婷婷色综合大香蕉| av在线app专区| xxx大片免费视频| 两个人免费观看高清视频| 久久影院123| 亚洲,欧美,日韩| 日日摸夜夜添夜夜爱| 黑丝袜美女国产一区| 国内精品宾馆在线| 婷婷色麻豆天堂久久| 999精品在线视频| 久久久久精品久久久久真实原创| 久久久国产欧美日韩av| 有码 亚洲区| 91精品一卡2卡3卡4卡| 国产高清有码在线观看视频| 熟女电影av网| 亚洲人与动物交配视频| 国产熟女午夜一区二区三区 | 色哟哟·www| a级片在线免费高清观看视频| 黄色怎么调成土黄色| 国产成人精品在线电影| 久久久久网色| av播播在线观看一区| 狠狠精品人妻久久久久久综合| 日韩制服骚丝袜av| 男人添女人高潮全过程视频| 中国国产av一级| 精品酒店卫生间| 最近手机中文字幕大全| 亚洲精品日韩在线中文字幕| 国产欧美日韩一区二区三区在线 | 国模一区二区三区四区视频| 精品国产一区二区三区久久久樱花| 大香蕉97超碰在线| 永久网站在线| 亚洲高清免费不卡视频| 精品人妻熟女毛片av久久网站| 午夜久久久在线观看| 久久国产精品男人的天堂亚洲 | 男女边吃奶边做爰视频| 欧美精品高潮呻吟av久久| 一区在线观看完整版| 亚洲精品成人av观看孕妇| 精品熟女少妇av免费看| 久久国产亚洲av麻豆专区| 97超碰精品成人国产| 大片电影免费在线观看免费| 老女人水多毛片| 女人精品久久久久毛片| 另类亚洲欧美激情| 91在线精品国自产拍蜜月| 精品国产一区二区三区久久久樱花| av女优亚洲男人天堂| 街头女战士在线观看网站| 国产在线一区二区三区精| 2018国产大陆天天弄谢| 精品少妇黑人巨大在线播放| 毛片一级片免费看久久久久| 国产一区二区三区综合在线观看 | 日韩成人av中文字幕在线观看| 大又大粗又爽又黄少妇毛片口| 3wmmmm亚洲av在线观看| 亚洲国产色片| 精品国产露脸久久av麻豆| 国语对白做爰xxxⅹ性视频网站| 在线看a的网站| a级毛色黄片| 国产国语露脸激情在线看| 久久国产精品男人的天堂亚洲 | 亚洲精品中文字幕在线视频| 亚洲av.av天堂| 久久久精品免费免费高清| 国产高清有码在线观看视频| 亚洲国产精品一区二区三区在线| 亚洲成人一二三区av| 亚洲国产av影院在线观看| 精品视频人人做人人爽| 国产黄片视频在线免费观看| 国产成人精品久久久久久| 日本黄色片子视频| 国产精品国产三级国产av玫瑰| kizo精华| 精品久久久精品久久久| 亚洲av国产av综合av卡| 国产深夜福利视频在线观看| 国国产精品蜜臀av免费| 成人免费观看视频高清| 熟女人妻精品中文字幕| 少妇丰满av| 免费人成在线观看视频色| 校园人妻丝袜中文字幕| 亚洲国产精品一区二区三区在线| 999精品在线视频| 黑人高潮一二区| 国产国语露脸激情在线看| 亚洲第一av免费看| 精品少妇久久久久久888优播| 狠狠婷婷综合久久久久久88av| 99热这里只有是精品在线观看| 色婷婷av一区二区三区视频| 成人午夜精彩视频在线观看| 午夜福利视频精品| 一级,二级,三级黄色视频| 内地一区二区视频在线| 多毛熟女@视频| 五月开心婷婷网| 久久久久精品性色| 亚洲精品乱码久久久久久按摩| 妹子高潮喷水视频| 爱豆传媒免费全集在线观看| 99re6热这里在线精品视频| 欧美 亚洲 国产 日韩一| 久久99蜜桃精品久久| 国产在线一区二区三区精| 丰满乱子伦码专区| 成人亚洲欧美一区二区av| 最近2019中文字幕mv第一页| 亚洲欧美成人综合另类久久久| 人妻 亚洲 视频| 嫩草影院入口| 汤姆久久久久久久影院中文字幕| 精品人妻一区二区三区麻豆| 18禁在线播放成人免费| 欧美三级亚洲精品| 如何舔出高潮| 男女边吃奶边做爰视频| 免费高清在线观看日韩| 99热这里只有是精品在线观看| 成人国语在线视频| 精品人妻熟女毛片av久久网站| 亚洲精品日本国产第一区| 欧美bdsm另类| 少妇被粗大的猛进出69影院 | 水蜜桃什么品种好| 日韩制服骚丝袜av| 中国美白少妇内射xxxbb| 春色校园在线视频观看| 日本色播在线视频| 久久婷婷青草| 嫩草影院入口| 亚洲精品,欧美精品| 亚洲av成人精品一区久久| 91在线精品国自产拍蜜月| av网站免费在线观看视频| 自线自在国产av| 中文字幕最新亚洲高清| av.在线天堂| 亚洲国产欧美日韩在线播放| 九色成人免费人妻av| 亚洲成色77777| 国产视频内射| 在线观看三级黄色| 性高湖久久久久久久久免费观看| 少妇 在线观看| 成人国产av品久久久| 国精品久久久久久国模美| 久久久亚洲精品成人影院| av在线app专区| 天美传媒精品一区二区| 交换朋友夫妻互换小说| 久久久久久久久久久丰满| 美女xxoo啪啪120秒动态图| 精品熟女少妇av免费看| 亚洲av不卡在线观看| 夫妻午夜视频| 国产熟女欧美一区二区| 晚上一个人看的免费电影| 激情五月婷婷亚洲| 免费看光身美女| 一区二区三区乱码不卡18| 午夜福利网站1000一区二区三区| 久久精品夜色国产| 国产免费现黄频在线看| 日韩欧美精品免费久久| 亚洲av日韩在线播放| 亚洲五月色婷婷综合| 26uuu在线亚洲综合色| 国产亚洲精品第一综合不卡 | 人人妻人人添人人爽欧美一区卜| 成年人免费黄色播放视频| 午夜福利视频精品| 国产成人精品久久久久久| 久久精品国产亚洲网站| av福利片在线| 少妇精品久久久久久久| 日韩强制内射视频| av有码第一页| 日本欧美国产在线视频| 汤姆久久久久久久影院中文字幕| 纯流量卡能插随身wifi吗| 日韩三级伦理在线观看|