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

    ARIMA與GM(1,1)模型對我國肺結(jié)核年發(fā)病人數(shù)預(yù)測情況的比較*

    2019-03-18 03:31:20濰坊醫(yī)學(xué)院261053
    中國衛(wèi)生統(tǒng)計(jì) 2019年1期
    關(guān)鍵詞:肺結(jié)核病乘積差分

    濰坊醫(yī)學(xué)院(261053)

    孫 娜 許小珊 馮佳寧 馬 潔 黃 璐 田 野 杜澤玉 孟維靜 王素珍△ 石福艷△

    s

    國家掌握肺結(jié)核的發(fā)展情況,有利于及時制定方針政策,更好的做到“早發(fā)現(xiàn)、早隔離、早治療”,從而及時做好肺結(jié)核的防治工作[1]。本次研究利用2011年至2016年我國肺結(jié)核的發(fā)病人數(shù),分別構(gòu)建ARIMA乘積季節(jié)模型和GM(1,1)模型,采用平均絕對百分誤差(MAPE)來評價模型的擬合效果,并給出肺結(jié)核發(fā)病人數(shù)的預(yù)測值,為我國傳染病應(yīng)急預(yù)案的制定提供依據(jù)。

    資料與方法

    1.資料來源

    研究中使用的2011年至2016年肺結(jié)核的發(fā)病人數(shù)資料來源于中國疾病預(yù)防控制中心發(fā)布的中華人民共和國國家衛(wèi)生和計(jì)劃生育委員會公報(bào)。利用2011年至2016年肺結(jié)核的月發(fā)病人數(shù)擬合ARIMA乘積季節(jié)模型,利用2011年至2016年肺結(jié)核的年發(fā)病人數(shù)擬合GM(1,1)模型。

    2.分析方法

    分別利用ARIMA乘積季節(jié)模型和GM(1,1)模型擬合我國肺結(jié)核的發(fā)病人數(shù),用于我國肺結(jié)核發(fā)病人數(shù)的預(yù)測,并通過比較平均絕對百分誤差(MAPE)得到擬合效果較好的模型。

    (1)ARIMA乘積季節(jié)模型

    當(dāng)序列既具有短期相關(guān)性還存在季節(jié)效應(yīng),同時短期相關(guān)性和季節(jié)效應(yīng)通過加法模型無法充分、有效提取時,就可以使用乘積季節(jié)模型來進(jìn)行擬合。

    (2)ARIMA乘積季節(jié)模型的構(gòu)造原理

    當(dāng)原始序列具有季節(jié)效應(yīng),季節(jié)相關(guān)性可以通過以周期步長為單位的ARMA(P,Q)模型提取;當(dāng)原始序列具有短期相關(guān)性時,一般可通過低階ARMA(p,q)模型提取。當(dāng)短期相關(guān)性和季節(jié)相關(guān)性之間具有乘積關(guān)系時,我們要擬合模型的模型則為ARMA(p,q)模型和ARMA(P,Q)模型的乘積。

    (3)ARIMA乘積季節(jié)模型的建模流程

    ①平穩(wěn)性檢驗(yàn):根據(jù)時序圖顯示的特征與ADF單位根檢驗(yàn)結(jié)果判斷序列的平穩(wěn)性。②白噪聲檢驗(yàn):對平穩(wěn)序列進(jìn)行白噪聲檢驗(yàn),若序列為白噪聲序列則分析結(jié)束,若序列不是白噪聲序列則進(jìn)一步擬合ARMA模型。③擬合ARMA模型。

    其中擬合ARMA模型的過程主要包括:①計(jì)算得到平穩(wěn)后序列的樣本自相關(guān)系數(shù)(ACF)與樣本偏自相關(guān)系數(shù)(PACF)的值。②根據(jù)樣本ACF和PACF,選擇合適的ARMA(p,q)模型進(jìn)行擬合。③運(yùn)用極大似然估計(jì)法估計(jì)模型參數(shù)。④模型的檢驗(yàn)。當(dāng)擬合模型不能通過檢驗(yàn)時,則重新進(jìn)行步驟(2)。⑤模型的優(yōu)化。當(dāng)擬合模型檢驗(yàn)通過了,仍然重新進(jìn)行步驟(2),擬合新的模型,并且通過AIC和SBC最小原則選擇最優(yōu)模型。⑥預(yù)測。利用所選的最優(yōu)模型進(jìn)行預(yù)測。

    (4)GM(1,1)模型的建立

    利用灰色預(yù)測法來擬合GM(1,1)模型[2]。設(shè)時間序列x(1)有n個觀察值,x(1)={X(0)(1),X(0)(2),…,X(0)(n)},通過一次累加生成以弱化其隨機(jī)性,強(qiáng)化其規(guī)律性的序列X(1)={X(1)(1),X(1)(2),…,X(1)(n)},則GM(1,1)模型相應(yīng)的微分方程為:

    式中,a稱為發(fā)展灰數(shù),μ稱為內(nèi)生控制灰數(shù)。利用最小二乘法求解a和μ,可以得到預(yù)測模型:

    (5)GM(1,1)模型的檢驗(yàn)

    (6)外推預(yù)測

    若所建模型的殘差檢驗(yàn)、后驗(yàn)差檢驗(yàn)和關(guān)聯(lián)度檢驗(yàn)均符合條件,我們就可以用該模型進(jìn)行預(yù)測,預(yù)測公式為:

    以2011-2016年我國肺結(jié)核發(fā)病人數(shù)為原始序列,利用SAS 9.4調(diào)用IML宏程序建立肺結(jié)核發(fā)病人數(shù)預(yù)測模型[4],然后對模型進(jìn)行檢驗(yàn),檢驗(yàn)通過后,運(yùn)用該模型預(yù)測我國2017-2018年肺結(jié)核的發(fā)病人數(shù)。

    (7)ARIMA與GM(1,1)模型效果評價

    本次研究采用平均絕對百分誤差(MAPE),比較ARIMA乘積季節(jié)模型和GM(1,1)模型兩種模型在預(yù)測我國肺結(jié)核年發(fā)病人數(shù)中的效果。

    MAPE=(∑|(xi-x)/x|*100%)/n

    式中xi與x分別代表為預(yù)測值和實(shí)際值。

    結(jié) 果

    1.ARIMA乘積季節(jié)模型建立

    (1)判斷序列的平穩(wěn)性

    圖1 2011-2016年我國肺結(jié)核發(fā)病人數(shù)時序圖

    時序圖顯示,該序列具有明顯的單調(diào)遞減的趨勢,同時也具有周期變化規(guī)律。我們嘗試進(jìn)行1階差分,實(shí)現(xiàn)其趨勢平穩(wěn);進(jìn)行12步差分,提取其中的周期信息。

    時序圖顯示,差分后序列在0附近比較穩(wěn)定地波動,可以初步認(rèn)為差分運(yùn)算后的序列平穩(wěn)。ADF單位根檢驗(yàn)結(jié)果如表1所示。

    由于ADF單位根檢驗(yàn)結(jié)果顯示檢驗(yàn)統(tǒng)計(jì)量所對應(yīng)的P值均小于0.05。因此,可以認(rèn)為經(jīng)過差分后的序列已經(jīng)平穩(wěn)。

    (2)白噪聲檢驗(yàn)

    圖2 差分后序列時序圖

    表1 增廣 Dickey-Fuller 單位根檢驗(yàn)

    表2 差分運(yùn)算后序列白噪聲檢驗(yàn)

    在檢驗(yàn)的顯著性水平取0.05的條件下,由于延遲6階的卡方統(tǒng)計(jì)量的值為23.16,P值為0.0007,小于0.05,所以該差分后的序列不能視為白噪聲序列,即差分后序列還蘊(yùn)含著相關(guān)信息可供提取。

    (3)ARIMA模型識別

    圖3 差分后序列自相關(guān)圖

    樣本自相關(guān)圖顯示延遲12階自相關(guān)系數(shù)顯著大于2倍標(biāo)準(zhǔn)差范圍,這說明差分后序列中仍蘊(yùn)含著非常顯著的季節(jié)效應(yīng)。而且延遲1階的自相關(guān)系數(shù)也大于2倍標(biāo)準(zhǔn)差,這說明差分后序列還具有短期相關(guān)性,考慮擬合乘積模型ARIMA(p,d,q)(P,D,Q)s。根據(jù)ACF圖和PACF圖所顯示的特征,認(rèn)為12階以內(nèi)的自相關(guān)圖1階截尾,偏自相關(guān)圖拖尾,嘗試ARMA(0,1)模型提取成分后序列的短期相關(guān)信息。一般P、Q階值比較難判斷,根據(jù)有關(guān)文獻(xiàn)[5]提示,P、Q階值通常不會超過2階,所以我們對P、Q分別取0、1、2,由低階到高階逐個實(shí)驗(yàn),根據(jù)AIC和SBC最小原則選擇最優(yōu)模型[6]。最終得到的可以充分提取序列信息的模型有以下兩個:

    圖4 差分后序列偏自相關(guān)圖

    表3 ARIMA(p,d,q)(P,D,Q)s模型的組合與最優(yōu)模型的選取準(zhǔn)則檢驗(yàn)

    根據(jù)AIC和SBC最小原則確定最優(yōu)模型為ARIMA(0,1,1)(0,1,1)12。

    (4)ARIMA模型參數(shù)估計(jì)與模型檢驗(yàn)

    運(yùn)用極大似然估計(jì)法估計(jì)模型參數(shù)。

    表4 ARIMA(0,1,1)(0,1,1)12模型參數(shù)顯著性檢驗(yàn)

    模型參數(shù)檢驗(yàn)的檢驗(yàn)統(tǒng)計(jì)量的P值都小于0.05,因此,在顯著性水平為0.05條件下,可以認(rèn)為模型參數(shù)值與0之間的差異具有統(tǒng)計(jì)學(xué)意義。

    表5 ARIMA(0,1,1)(0,1,1)12模型殘差白噪聲檢驗(yàn)

    延遲各階的LB統(tǒng)計(jì)量的P值均大于0.05,所以擬合的ARIMA(0,1,1)(0,1,1)12模型成立。擬合模型的具體形式為▽▽12xt=(1-0.83148B)(1-0.73177B12)εt

    根據(jù)擬合出的ARIMA(0,1,1)(0,1,1)12模型得到2017-2018年我國肺結(jié)核的月發(fā)病人數(shù)依次為82561、78303、109046、103622、99448、92873、93436、90499、86277、79478、82595、83082、78809、74552、105295、99871、95697、89122、89685、86748、82525、75726、78844、79330人。因此,2017-2018年我國肺結(jié)核的年發(fā)病人數(shù)依次為1081220、1036204人。

    2.GM(1,1)模型

    以我國2011-2016年肺結(jié)核發(fā)病人數(shù)為原始序列,運(yùn)用SAS 9.4計(jì)算得出模型參數(shù)a=0.053161,μ=1469122.88,因此,應(yīng)用GM(1,1)模型建立我國肺結(jié)核發(fā)病人數(shù)預(yù)測模型:

    (1)殘差檢驗(yàn)

    表6 預(yù)測值與真實(shí)值之間的比較

    結(jié)果顯示,所有的相對誤差均小于0.05,所以該模型可以通過殘差檢驗(yàn)。

    (2)后驗(yàn)差檢驗(yàn)

    經(jīng)檢驗(yàn),后驗(yàn)差C值為0.17385,小概率P值為1.00,模型預(yù)測精度等級為好[3],可以將該模型用于外推預(yù)測。

    (3)關(guān)聯(lián)度檢驗(yàn)

    運(yùn)用SAS 9.4進(jìn)行計(jì)算,結(jié)果顯示,關(guān)聯(lián)度r=0.63042>0.6,滿足檢驗(yàn)準(zhǔn)則。因此,模型通過了關(guān)聯(lián)度檢驗(yàn)。

    (4)肺結(jié)核發(fā)病人數(shù)的預(yù)測

    根據(jù)擬合的GM(1,1)模型進(jìn)行肺結(jié)核發(fā)病人數(shù)的預(yù)測。預(yù)測我國2017-2018年肺結(jié)核的年發(fā)病人數(shù)依次為1042909、988915人。

    3.模型效果評價

    根據(jù)模型預(yù)測結(jié)果,計(jì)算2011-2016年兩模型的MAPE值。肺結(jié)核年發(fā)病人數(shù)ARIMA乘積季節(jié)模型和GM(1,1)模型MAPE分別為1.70%和1.45%。通過比較MAPE可以發(fā)現(xiàn)GM(1,1)模型擬合效果要好于ARIMA乘積季節(jié)模型。

    討 論

    目前在時間序列預(yù)測中應(yīng)用比較多的模型就是ARIMA乘積季節(jié)模型和GM(1,1)模型。本研究中運(yùn)用GM(1,1)模型和ARIMA乘積季節(jié)模型對我們肺結(jié)核發(fā)病人數(shù)進(jìn)行擬合,用于比較ARIMA乘積季節(jié)模型和GM(1,1)模型在我們年肺結(jié)核發(fā)病人數(shù)預(yù)測中的準(zhǔn)確性。兩種模型擬合均通過了模型檢驗(yàn),均可用于我國肺結(jié)核發(fā)病人數(shù)的外推預(yù)測,通過比較兩種模型的MAPE可以發(fā)現(xiàn),GM(1,1)模型的MAPE小于ARIMA乘積季節(jié)模型的MAPE。因此GM(1,1)模型對我國肺結(jié)核年發(fā)病人數(shù)的擬合效果高于ARIMA乘積季節(jié)模型。可能是因?yàn)槲覀冊谶M(jìn)行年發(fā)病人數(shù)預(yù)測時,體現(xiàn)不出ARIMA乘積季節(jié)模型的預(yù)測優(yōu)勢。但是,通過ARIMA乘積季節(jié)模型可以得到我國肺結(jié)核的月發(fā)病人數(shù),也具有其獨(dú)特的優(yōu)勢。

    運(yùn)用所得到的GM(1,1)灰色預(yù)測模型預(yù)測我國未來兩年肺結(jié)核的發(fā)病人數(shù)依次為1042909、988915人。預(yù)測結(jié)果顯示,我國2017-2018年肺結(jié)核的發(fā)病人數(shù)將呈現(xiàn)出下降的趨勢,但是其發(fā)病數(shù)仍然處于較高的水平。

    肺結(jié)核是慢性傳染性病,其流行趨勢是有規(guī)律可循得,因此疾病預(yù)防控制中心應(yīng)積極采取綜合措施,聚焦重點(diǎn)人群、強(qiáng)化業(yè)務(wù)知識培訓(xùn)、全面深入推進(jìn)全民肺結(jié)核預(yù)防策略與措施的宣傳,提高公眾對肺結(jié)核預(yù)防知識的知曉率[7],更好的做到“早發(fā)現(xiàn)、早隔離、早治療”,從而提高肺結(jié)核病人的發(fā)現(xiàn)率。疾病預(yù)防控制中心還應(yīng)該采取措施強(qiáng)化聯(lián)防聯(lián)控,在宣傳肺結(jié)核預(yù)防策略與措施的同時宣傳現(xiàn)代肺結(jié)核治療策略與措施,讓公眾更好地了解到國家對肺結(jié)核病人所采用的政策與措施,使肺結(jié)核的病人可以得到規(guī)范的治療和管理,從而使肺結(jié)核病人的治愈率得到提高。對于確診為肺結(jié)核的患者,我們應(yīng)該根據(jù)他們的實(shí)際病情,及時的采取有效的治療措施,對患者的病情進(jìn)行控制,才能有效控制傳染源,從而降低肺結(jié)核病人的傳染率。對于肺結(jié)核的病人,我們還應(yīng)該及時的掌握病人的病情資料以及藥物的使用情況,督促病人及時用藥,加強(qiáng)肺結(jié)核病人規(guī)律治療的意識,避免不規(guī)律服藥造成的治療失敗等后果。對肺結(jié)核病人進(jìn)行長期的追蹤治療,還可以及時獲得他們的病情詳細(xì)資料,為將來對肺結(jié)核的觀察性研究提供資料。

    猜你喜歡
    肺結(jié)核病乘積差分
    健康教育對于肺結(jié)核病防控管理的效果研究
    數(shù)列與差分
    乘積最大
    Dirichlet級數(shù)及其Dirichlet-Hadamard乘積的增長性
    莫西沙星和左氧氟沙星在耐多藥肺結(jié)核病治療方案中的效果比較
    110例耐藥肺結(jié)核病臨床特征分析
    復(fù)變?nèi)呛瘮?shù)無窮乘積的若干應(yīng)用
    基于差分隱私的大數(shù)據(jù)隱私保護(hù)
    牛肺結(jié)核病的病理形態(tài)學(xué)探討
    Dirichlet級數(shù)的Dirichlet-Hadamard乘積
    中文欧美无线码| 日本欧美国产在线视频| 黄色毛片三级朝国网站| 精品一区二区三区av网在线观看 | 一级毛片我不卡| 大片免费播放器 马上看| 大片免费播放器 马上看| 精品一区二区三区av网在线观看 | 国产一级毛片在线| 午夜两性在线视频| 老司机在亚洲福利影院| 亚洲午夜精品一区,二区,三区| 最近手机中文字幕大全| 水蜜桃什么品种好| 婷婷色综合www| 色视频在线一区二区三区| 日本av手机在线免费观看| 最近手机中文字幕大全| 热99国产精品久久久久久7| 国产高清videossex| 大片电影免费在线观看免费| 1024视频免费在线观看| 在线天堂中文资源库| 亚洲精品成人av观看孕妇| 欧美人与性动交α欧美软件| 亚洲国产最新在线播放| 国产片特级美女逼逼视频| 日本欧美国产在线视频| 日本av手机在线免费观看| 观看av在线不卡| 老司机影院毛片| 国产精品香港三级国产av潘金莲 | 亚洲七黄色美女视频| 麻豆av在线久日| 亚洲av成人精品一二三区| 色综合欧美亚洲国产小说| 国产高清国产精品国产三级| 亚洲av男天堂| 欧美人与性动交α欧美精品济南到| 亚洲av欧美aⅴ国产| 一区二区av电影网| 黄色a级毛片大全视频| 日韩 亚洲 欧美在线| 午夜影院在线不卡| 天天操日日干夜夜撸| 亚洲伊人久久精品综合| 一本综合久久免费| 久久精品国产a三级三级三级| 日本黄色日本黄色录像| 日本一区二区免费在线视频| 啦啦啦在线免费观看视频4| 啦啦啦啦在线视频资源| 在线天堂中文资源库| 久久九九热精品免费| 一边亲一边摸免费视频| 亚洲欧美精品综合一区二区三区| 国产欧美日韩精品亚洲av| 国产精品av久久久久免费| 两性夫妻黄色片| 国产精品亚洲av一区麻豆| 热re99久久国产66热| 每晚都被弄得嗷嗷叫到高潮| av视频免费观看在线观看| 久久99一区二区三区| 18禁裸乳无遮挡动漫免费视频| 精品卡一卡二卡四卡免费| 久久久久久久国产电影| 亚洲人成77777在线视频| 麻豆av在线久日| 波多野结衣一区麻豆| 人成视频在线观看免费观看| 亚洲精品久久成人aⅴ小说| 中文字幕制服av| 国产日韩欧美视频二区| 青春草亚洲视频在线观看| 精品国产乱码久久久久久小说| 国产免费现黄频在线看| 亚洲中文日韩欧美视频| 成人黄色视频免费在线看| 日本午夜av视频| 国产熟女欧美一区二区| 黄片小视频在线播放| 好男人电影高清在线观看| 精品人妻1区二区| 欧美精品高潮呻吟av久久| 巨乳人妻的诱惑在线观看| 国产真人三级小视频在线观看| 高清欧美精品videossex| 国产精品一区二区免费欧美 | 另类亚洲欧美激情| 精品少妇一区二区三区视频日本电影| 国产99久久九九免费精品| 国产1区2区3区精品| 久久国产精品影院| 大陆偷拍与自拍| 亚洲熟女毛片儿| 日韩av在线免费看完整版不卡| 大陆偷拍与自拍| 亚洲成国产人片在线观看| 成人影院久久| 丝袜在线中文字幕| 亚洲精品乱久久久久久| 黄色一级大片看看| 国产一区二区 视频在线| 国产成人一区二区在线| a 毛片基地| av片东京热男人的天堂| 日韩一区二区三区影片| 电影成人av| 视频在线观看一区二区三区| 午夜日韩欧美国产| 极品人妻少妇av视频| 国产精品麻豆人妻色哟哟久久| av一本久久久久| 欧美精品高潮呻吟av久久| √禁漫天堂资源中文www| 精品人妻在线不人妻| 国产人伦9x9x在线观看| 亚洲专区中文字幕在线| 免费看不卡的av| 一级毛片我不卡| 免费观看av网站的网址| 久久热在线av| 热re99久久国产66热| 国产成人av教育| 最新的欧美精品一区二区| 视频区欧美日本亚洲| 亚洲人成网站在线观看播放| 国产日韩欧美视频二区| 日韩视频在线欧美| 亚洲av国产av综合av卡| 我的亚洲天堂| av线在线观看网站| 你懂的网址亚洲精品在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧洲国产日韩| 中文欧美无线码| 美女脱内裤让男人舔精品视频| 天天躁狠狠躁夜夜躁狠狠躁| 欧美中文综合在线视频| 久久人妻福利社区极品人妻图片 | 国产成人精品无人区| 久久久久精品人妻al黑| 国产午夜精品一二区理论片| 9热在线视频观看99| 91精品伊人久久大香线蕉| 2018国产大陆天天弄谢| 大话2 男鬼变身卡| 欧美人与善性xxx| 热re99久久国产66热| 99久久99久久久精品蜜桃| 日本a在线网址| 国产av国产精品国产| 肉色欧美久久久久久久蜜桃| 精品久久久久久久毛片微露脸 | 国产成人影院久久av| 亚洲,欧美,日韩| 中文字幕亚洲精品专区| 国产极品粉嫩免费观看在线| 在线精品无人区一区二区三| 999精品在线视频| 国产黄频视频在线观看| 丝袜脚勾引网站| 老司机影院毛片| 欧美精品亚洲一区二区| 老司机亚洲免费影院| 国产一区二区在线观看av| 日韩熟女老妇一区二区性免费视频| 狠狠精品人妻久久久久久综合| 久久久久网色| 日日摸夜夜添夜夜爱| 日韩一卡2卡3卡4卡2021年| 视频区图区小说| 超色免费av| 国产黄频视频在线观看| 欧美精品一区二区免费开放| 久久精品亚洲av国产电影网| 久久99精品国语久久久| 色网站视频免费| 日本午夜av视频| 真人做人爱边吃奶动态| 日韩视频在线欧美| 欧美人与性动交α欧美精品济南到| av国产精品久久久久影院| 老汉色∧v一级毛片| 五月开心婷婷网| 亚洲欧洲国产日韩| 亚洲精品一卡2卡三卡4卡5卡 | 天天操日日干夜夜撸| 在线观看人妻少妇| 黄色 视频免费看| 青青草视频在线视频观看| 自拍欧美九色日韩亚洲蝌蚪91| 大香蕉久久成人网| 成人黄色视频免费在线看| 免费高清在线观看视频在线观看| a级毛片在线看网站| av网站免费在线观看视频| 天堂中文最新版在线下载| www.av在线官网国产| www.av在线官网国产| 高清黄色对白视频在线免费看| 男女免费视频国产| 亚洲欧洲国产日韩| 久久精品久久精品一区二区三区| 国产精品一区二区精品视频观看| 搡老岳熟女国产| 99国产综合亚洲精品| 久久这里只有精品19| 久久99热这里只频精品6学生| 黑人欧美特级aaaaaa片| 一二三四在线观看免费中文在| 国产又色又爽无遮挡免| 国产精品偷伦视频观看了| 亚洲精品日韩在线中文字幕| 天天躁夜夜躁狠狠久久av| 亚洲av日韩精品久久久久久密 | 又大又爽又粗| 首页视频小说图片口味搜索 | 美女高潮到喷水免费观看| cao死你这个sao货| 成人国产av品久久久| www日本在线高清视频| 91老司机精品| 热99国产精品久久久久久7| xxxhd国产人妻xxx| 丝袜脚勾引网站| 中文字幕制服av| 亚洲少妇的诱惑av| 成人亚洲精品一区在线观看| 国产精品一国产av| 青青草视频在线视频观看| 亚洲欧美激情在线| 亚洲成人国产一区在线观看 | 欧美 亚洲 国产 日韩一| 欧美成人精品欧美一级黄| 深夜精品福利| 视频在线观看一区二区三区| 久久久精品94久久精品| 久久国产精品影院| 国产午夜精品一二区理论片| 母亲3免费完整高清在线观看| 极品人妻少妇av视频| 日韩中文字幕欧美一区二区 | 电影成人av| 久久青草综合色| 伊人久久大香线蕉亚洲五| 久久中文字幕一级| 日日摸夜夜添夜夜爱| 少妇精品久久久久久久| 亚洲av成人精品一二三区| 美女国产高潮福利片在线看| 大陆偷拍与自拍| 国产激情久久老熟女| 国产高清视频在线播放一区 | 成人国语在线视频| 在线亚洲精品国产二区图片欧美| 久久99精品国语久久久| 精品福利观看| 久久久国产精品麻豆| 黑人猛操日本美女一级片| 99久久人妻综合| 久久精品成人免费网站| 国产精品熟女久久久久浪| 久久久久精品人妻al黑| 男女高潮啪啪啪动态图| 精品免费久久久久久久清纯 | 亚洲精品国产av蜜桃| 好男人电影高清在线观看| 精品一品国产午夜福利视频| 亚洲av电影在线进入| 少妇人妻久久综合中文| 国产精品一区二区免费欧美 | 久久久久久久大尺度免费视频| 在线观看免费视频网站a站| 亚洲,欧美,日韩| 精品久久久久久久毛片微露脸 | 色94色欧美一区二区| 人人妻人人澡人人看| 美女福利国产在线| 日本五十路高清| 欧美日韩综合久久久久久| 性少妇av在线| 热99国产精品久久久久久7| 在线观看免费视频网站a站| 18在线观看网站| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品中文字幕在线视频| 国产在线免费精品| 菩萨蛮人人尽说江南好唐韦庄| 五月天丁香电影| 亚洲第一青青草原| 少妇人妻 视频| 国产精品免费视频内射| 久久精品国产亚洲av涩爱| 免费人妻精品一区二区三区视频| 久久久久久久久久久久大奶| 亚洲国产毛片av蜜桃av| 黄片小视频在线播放| av视频免费观看在线观看| 国产熟女欧美一区二区| 日本wwww免费看| 99久久99久久久精品蜜桃| 欧美精品高潮呻吟av久久| 欧美精品人与动牲交sv欧美| 久久国产精品人妻蜜桃| 蜜桃国产av成人99| www.av在线官网国产| 免费不卡黄色视频| 女性被躁到高潮视频| 蜜桃在线观看..| 后天国语完整版免费观看| 美女福利国产在线| tube8黄色片| 母亲3免费完整高清在线观看| 精品第一国产精品| 久久久久久久国产电影| 日韩av不卡免费在线播放| 老司机深夜福利视频在线观看 | 成年美女黄网站色视频大全免费| av国产精品久久久久影院| kizo精华| 欧美国产精品一级二级三级| 美女中出高潮动态图| 国产一区二区在线观看av| 一边摸一边抽搐一进一出视频| 国产欧美日韩一区二区三 | 人人妻人人澡人人爽人人夜夜| 免费观看人在逋| 国产精品秋霞免费鲁丝片| 国产在视频线精品| 伊人久久大香线蕉亚洲五| 欧美激情极品国产一区二区三区| 色精品久久人妻99蜜桃| 免费av中文字幕在线| 久久久久久免费高清国产稀缺| 岛国毛片在线播放| 又大又黄又爽视频免费| 男女边吃奶边做爰视频| 亚洲专区国产一区二区| 黄色视频不卡| 夫妻午夜视频| 欧美+亚洲+日韩+国产| 日本a在线网址| 亚洲国产av新网站| 国产99久久九九免费精品| 丝袜喷水一区| 亚洲美女黄色视频免费看| 脱女人内裤的视频| 亚洲精品一区蜜桃| av有码第一页| 日韩大片免费观看网站| 欧美精品亚洲一区二区| 国产真人三级小视频在线观看| 欧美日韩视频高清一区二区三区二| 亚洲欧美精品综合一区二区三区| 一个人免费看片子| av有码第一页| 亚洲国产中文字幕在线视频| 久久久久久免费高清国产稀缺| 老司机影院毛片| 丰满人妻熟妇乱又伦精品不卡| 999精品在线视频| 亚洲图色成人| 国产精品一区二区在线观看99| 狂野欧美激情性xxxx| 99热全是精品| 成人亚洲欧美一区二区av| 日韩大片免费观看网站| 狂野欧美激情性bbbbbb| 亚洲成色77777| 久久精品久久精品一区二区三区| www.av在线官网国产| 国产99久久九九免费精品| 99国产精品免费福利视频| 在线观看国产h片| 久热这里只有精品99| 成人国语在线视频| 久久国产精品男人的天堂亚洲| 日韩大片免费观看网站| 国产又爽黄色视频| 一边亲一边摸免费视频| 人妻 亚洲 视频| 久久久国产一区二区| 亚洲欧美色中文字幕在线| 爱豆传媒免费全集在线观看| 国产日韩欧美在线精品| 亚洲国产中文字幕在线视频| 午夜免费男女啪啪视频观看| 久久ye,这里只有精品| 少妇被粗大的猛进出69影院| 久久久久精品国产欧美久久久 | 久久久精品区二区三区| 亚洲熟女精品中文字幕| 中文字幕亚洲精品专区| 久久精品国产亚洲av高清一级| 看免费av毛片| 免费人妻精品一区二区三区视频| 国产精品免费大片| 亚洲精品在线美女| 日本a在线网址| 久久青草综合色| 无遮挡黄片免费观看| 国产精品久久久久久人妻精品电影 | 久久免费观看电影| 久久 成人 亚洲| 国产精品久久久人人做人人爽| 国产国语露脸激情在线看| 欧美精品啪啪一区二区三区 | 国产伦理片在线播放av一区| 肉色欧美久久久久久久蜜桃| 国产成人影院久久av| 夫妻午夜视频| 新久久久久国产一级毛片| 青草久久国产| 99国产精品99久久久久| 亚洲国产成人一精品久久久| 母亲3免费完整高清在线观看| 国产精品免费大片| 亚洲国产欧美日韩在线播放| 免费人妻精品一区二区三区视频| 日本欧美国产在线视频| 国产一区二区激情短视频 | 国产欧美日韩综合在线一区二区| 脱女人内裤的视频| av线在线观看网站| 国产高清不卡午夜福利| 欧美大码av| 久久久国产欧美日韩av| 亚洲色图 男人天堂 中文字幕| 亚洲精品日本国产第一区| 19禁男女啪啪无遮挡网站| 中文字幕最新亚洲高清| av在线app专区| 婷婷色综合大香蕉| 国产黄色视频一区二区在线观看| 亚洲中文av在线| 亚洲av电影在线观看一区二区三区| 无遮挡黄片免费观看| 七月丁香在线播放| 日韩精品免费视频一区二区三区| 亚洲精品国产色婷婷电影| 亚洲av美国av| 久久国产亚洲av麻豆专区| 色94色欧美一区二区| 亚洲欧美一区二区三区久久| 久久精品国产亚洲av高清一级| 欧美日韩精品网址| 久久九九热精品免费| 美女视频免费永久观看网站| 极品人妻少妇av视频| 午夜影院在线不卡| 久久ye,这里只有精品| 91精品三级在线观看| 啦啦啦视频在线资源免费观看| 午夜福利免费观看在线| 免费在线观看视频国产中文字幕亚洲 | 欧美日韩综合久久久久久| 人体艺术视频欧美日本| 精品一区二区三卡| 婷婷色综合大香蕉| 国产视频一区二区在线看| 精品免费久久久久久久清纯 | 咕卡用的链子| 日本黄色日本黄色录像| 成人18禁高潮啪啪吃奶动态图| 80岁老熟妇乱子伦牲交| 精品少妇久久久久久888优播| 性高湖久久久久久久久免费观看| av不卡在线播放| 男女边摸边吃奶| 国产亚洲一区二区精品| 精品国产一区二区久久| 国产在线观看jvid| 国产片内射在线| 欧美人与善性xxx| 中文字幕av电影在线播放| 性高湖久久久久久久久免费观看| 不卡av一区二区三区| 久久精品成人免费网站| 三上悠亚av全集在线观看| 七月丁香在线播放| 国产精品久久久人人做人人爽| 国产成人av激情在线播放| 美女高潮到喷水免费观看| 又黄又粗又硬又大视频| 日韩av不卡免费在线播放| 欧美亚洲日本最大视频资源| av福利片在线| 巨乳人妻的诱惑在线观看| 国产高清国产精品国产三级| 国产在线免费精品| 免费日韩欧美在线观看| 日韩大码丰满熟妇| 美女主播在线视频| 国产精品麻豆人妻色哟哟久久| 男女高潮啪啪啪动态图| av有码第一页| 色婷婷久久久亚洲欧美| 成年人黄色毛片网站| 亚洲精品一区蜜桃| 免费av中文字幕在线| 十八禁网站网址无遮挡| 夜夜骑夜夜射夜夜干| 色精品久久人妻99蜜桃| 国产爽快片一区二区三区| 各种免费的搞黄视频| 亚洲人成网站在线观看播放| 国产1区2区3区精品| 你懂的网址亚洲精品在线观看| 久热爱精品视频在线9| 69精品国产乱码久久久| 波多野结衣av一区二区av| 亚洲国产中文字幕在线视频| 十八禁网站网址无遮挡| 丰满少妇做爰视频| 亚洲国产精品一区二区三区在线| www日本在线高清视频| 国产成人精品无人区| 久久久久久免费高清国产稀缺| av福利片在线| 欧美+亚洲+日韩+国产| 国产成人av教育| 国产精品香港三级国产av潘金莲 | 国产男女超爽视频在线观看| 黄色视频在线播放观看不卡| 久9热在线精品视频| 精品亚洲成国产av| 亚洲精品国产色婷婷电影| 女人久久www免费人成看片| 热99国产精品久久久久久7| 91国产中文字幕| av天堂在线播放| 夫妻午夜视频| 视频区欧美日本亚洲| 久久久久国产精品人妻一区二区| 人人妻,人人澡人人爽秒播 | 亚洲伊人久久精品综合| 亚洲少妇的诱惑av| 国产97色在线日韩免费| 免费黄频网站在线观看国产| 亚洲精品乱久久久久久| 美女午夜性视频免费| 精品欧美一区二区三区在线| 欧美av亚洲av综合av国产av| xxx大片免费视频| 久久狼人影院| 国产亚洲午夜精品一区二区久久| 精品视频人人做人人爽| 欧美国产精品一级二级三级| 国产又色又爽无遮挡免| 青春草亚洲视频在线观看| 一区二区日韩欧美中文字幕| 大片电影免费在线观看免费| 一区二区日韩欧美中文字幕| 激情五月婷婷亚洲| 考比视频在线观看| 久久午夜综合久久蜜桃| 男女床上黄色一级片免费看| 亚洲精品久久午夜乱码| 亚洲av日韩精品久久久久久密 | 国产日韩欧美亚洲二区| av片东京热男人的天堂| 国产亚洲av片在线观看秒播厂| av网站在线播放免费| a级片在线免费高清观看视频| 国产成人av激情在线播放| 看免费成人av毛片| 啦啦啦视频在线资源免费观看| 看十八女毛片水多多多| 精品高清国产在线一区| av有码第一页| 免费在线观看日本一区| 新久久久久国产一级毛片| 在线观看一区二区三区激情| 色综合欧美亚洲国产小说| 宅男免费午夜| 高潮久久久久久久久久久不卡| 久久久精品免费免费高清| 丰满人妻熟妇乱又伦精品不卡| 午夜免费观看性视频| 久久久精品94久久精品| av视频免费观看在线观看| 亚洲国产精品一区三区| 老汉色av国产亚洲站长工具| 午夜福利视频精品| 丝袜美腿诱惑在线| 人妻人人澡人人爽人人| 国产成人精品久久久久久| 国产黄色视频一区二区在线观看| 汤姆久久久久久久影院中文字幕| 啦啦啦啦在线视频资源| 欧美日韩福利视频一区二区| 亚洲少妇的诱惑av| 一级片'在线观看视频| 在线观看免费视频网站a站| 成人18禁高潮啪啪吃奶动态图| 日韩,欧美,国产一区二区三区| 色播在线永久视频| 51午夜福利影视在线观看| 欧美精品亚洲一区二区| 香蕉国产在线看| av线在线观看网站| 国产男人的电影天堂91| 午夜精品国产一区二区电影| 欧美黄色淫秽网站| 宅男免费午夜| 色婷婷久久久亚洲欧美| 亚洲av成人精品一二三区| 午夜免费男女啪啪视频观看| 一区二区三区乱码不卡18| 亚洲美女黄色视频免费看| 午夜影院在线不卡| 青春草亚洲视频在线观看| 久久久国产精品麻豆|