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

    干預(yù)分析在腎綜合出血熱發(fā)病過程模擬中的應(yīng)用

    2016-10-13 02:06:07程幸福
    關(guān)鍵詞:新息分析模型均值

    程幸福

    (武漢大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,湖北武漢430072)

    干預(yù)分析在腎綜合出血熱發(fā)病過程模擬中的應(yīng)用

    程幸福

    (武漢大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,湖北武漢430072)

    為探討應(yīng)用干預(yù)分析模型模擬腎綜合出血熱發(fā)病數(shù)的可行性,對(duì)武漢市蔡甸區(qū)1980~2009年的月度發(fā)病數(shù)據(jù)建立干預(yù)分析模型,并將模擬效果與單純使用ARIMA模型的模擬效果進(jìn)行對(duì)比,發(fā)現(xiàn)干預(yù)分析模型能夠很好地?cái)M合HFRS的月發(fā)病數(shù),給疾病的預(yù)防和控制提供科學(xué)依據(jù)。

    干預(yù)分析;ARIMA模型;腎綜合征出血熱

    干預(yù)分析的研究始于威斯康辛大學(xué)統(tǒng)計(jì)系教授Box和Tiao于1975年聯(lián)合發(fā)表《Intervention Analysis with Applications to Economic and Environ?mental Problems》,并迅速引起人們的關(guān)注,用于描述經(jīng)濟(jì)政策的變化及其給經(jīng)濟(jì)帶來的影響。研究干預(yù)分析模型就是從定量分析的角度來評(píng)估政策干預(yù)或突發(fā)事件對(duì)經(jīng)濟(jì)過程的具體影響。

    1 干預(yù)分析模型

    1.1干預(yù)分析簡(jiǎn)介

    時(shí)間序列常受到諸如節(jié)假日、罷工、促銷和其他政策變化之類的外部事件的影響,這類外部事件稱為干預(yù)。干預(yù)分析就是在時(shí)間序列建模中納入干預(yù)信息的建模方法。

    1.2常見的干預(yù)變量

    為簡(jiǎn)單起見,考慮單一干預(yù)的簡(jiǎn)單情況,假設(shè)…,Yt-1,Yt,Yt+1,…是得到的時(shí)間序列,經(jīng)過適當(dāng)變換后,可表示為

    式中:yt是由Yt經(jīng)某個(gè)變換得到的;mt是均值函數(shù)的變化;Nt是某個(gè)自回歸滑動(dòng)平均求和過程,代表未受干預(yù)的時(shí)間序列所適合的模型,可能是平穩(wěn)的也可能是非平穩(wěn)的、季節(jié)的或非季節(jié)的。假設(shè)Nt可以表示為ARIMA(p,d,q)模型:

    式中:B是向后推移算子,假設(shè)這個(gè)時(shí)間序列在T時(shí)刻受到某個(gè)干擾,即在時(shí)刻T之前,mt=0。

    干預(yù)分析模型的建立基于干預(yù)變量的引入,通過若干參數(shù)來說明干預(yù)對(duì)均值函數(shù)的影響,例如階梯函數(shù)函數(shù):階梯函數(shù)通常代表一種干預(yù)一旦發(fā)生將永久持續(xù)下去的情形。類似用脈沖函數(shù)作為干預(yù)發(fā)生時(shí)間的指示器,即Pt(T)是干預(yù)發(fā)生的標(biāo)記:

    1.3干預(yù)模型的表示

    干預(yù)事件的影響多種多樣,但幾乎所有的影響都可以用以上2種干預(yù)變量的線性組合來表示:

    1)如果干預(yù)導(dǎo)致了均值函數(shù)發(fā)生了即時(shí)且永久性的偏移,那么,偏移可以建模為

    式中:ω表示干預(yù)導(dǎo)致均值未知且永久的變化。如果干預(yù)經(jīng)過d個(gè)時(shí)間單位的延遲后作用才顯現(xiàn),并且d已知,那么可以設(shè)定:

    2)如果干預(yù)只是逐漸地影響均值函數(shù),其全部影響作用只有經(jīng)歷很長(zhǎng)時(shí)間才能充分體現(xiàn)出來,這時(shí)可以用AR(1)模型來表示mt,其誤差項(xiàng)由滯后1的St

    T的倍數(shù)充任:

    3)如果干預(yù)僅在t=T時(shí)刻產(chǎn)生影響,則令

    4)如果干預(yù)影響逐漸消失,通過AR(1)類型的模型表示均值函數(shù)的變化:

    如果均值的變化發(fā)生于d個(gè)時(shí)間單位的延遲之后,且該干預(yù)效應(yīng)逐漸消失,則可規(guī)定:也可以用延遲算子B來表示上述模型,Bmt=mt-1且

    1.4異常值

    異常值是一些不規(guī)則的觀測(cè)值,其出現(xiàn)可能源自測(cè)量誤差與復(fù)制誤差其中之一,或者兩者皆有,也可能是源于基礎(chǔ)過程發(fā)生了突發(fā)的短期性變化。對(duì)于時(shí)間序列而言,可識(shí)別的異常值,即可加異常值與新息異常值,通常分別記為AO與IO。

    1)可加異常值A(chǔ)O如果基礎(chǔ)過程yt在時(shí)刻T時(shí)受到了可疊加性的擾動(dòng)而得到y(tǒng)t,此時(shí)就會(huì)出現(xiàn)可加異常值,它們之間有下面的關(guān)系:

    僅在T時(shí)刻有yT′=yT+ωA,其余情況下yt′=yt,即如果時(shí)間序列在T時(shí)刻存在一個(gè)可加異常值,意味著它只在時(shí)刻T受到影響。可加異常值也可被看作一種在時(shí)刻T上具有脈沖響應(yīng)的干預(yù),因此均值函數(shù)建模為

    2)新息異常值IO如果處于時(shí)間T時(shí)的誤差(也稱為新息)受到了干擾(即誤差變?yōu)閑t′=et+ωI,其中et是零均值白噪聲)。因此僅在時(shí)刻T有eT′=eT+ωI,其余情況et′=et。假設(shè)無擾過程是平穩(wěn)的,且可表示為一個(gè)MA(∞)過程:

    則受到擾動(dòng)的過程可以表示為

    式中:ψ0=1且當(dāng) j<0時(shí)有ψj=0。因此,一旦發(fā)生新息異常值,則新息異常值將從該異常值發(fā)生的時(shí)刻T起對(duì)之后所有的觀測(cè)產(chǎn)生擾動(dòng)。若Nt滿足φ(B) Nt=θ(B) et,則在此情況下,模型可表示為

    2 腎綜合出血熱發(fā)病過程的模擬

    2.1數(shù)據(jù)來源

    選擇武漢市蔡甸區(qū)HFRS月發(fā)病數(shù)序列(已作開方處理)作為分析對(duì)象,記錄了1980年1月~2009年12月的HFRS發(fā)病數(shù)共360個(gè)觀測(cè),如圖1所示。

    2.2對(duì)基礎(chǔ)無擾過程建模

    圖1 HFRS序列圖

    從圖1發(fā)現(xiàn)1982年12月(即第36個(gè)數(shù)據(jù))迅速增長(zhǎng),可以認(rèn)為在1982年12月有一個(gè)對(duì)HFRS發(fā)病數(shù)序列產(chǎn)生干擾的事件發(fā)生。在進(jìn)行干擾分析之前,檢查時(shí)間序列并且建立基礎(chǔ)模型,為此選取干預(yù)發(fā)生之前的數(shù)據(jù)(即HFRS序列的前35個(gè)數(shù)據(jù))進(jìn)行建模。通過分析發(fā)現(xiàn)基礎(chǔ)無擾過程可以用季模型ARIMA(0,1,1)×(0,1,1)12表示。即如果用Nt表示基礎(chǔ)無擾過程,則有

    2.3對(duì)干擾過程建模

    在確定了基礎(chǔ)無擾過程的模型后,模型的殘差即干擾過程mt,通過對(duì)mt進(jìn)行簡(jiǎn)單自相關(guān)分析發(fā)現(xiàn)起始于1982年12月份的干擾過程可以用一個(gè)AR(2)過程來描述,因此對(duì)干擾過程可以建模為

    式中:T表示1982年12月。另外,模型診斷表明:在1981年2~3月有可加異常值,1981年12月、1987 年1月和1990年12月存在新息異常值。因此,最終確定的模型為

    2.4參數(shù)估計(jì)

    用極大似然估計(jì)法估計(jì)模型的參數(shù),基礎(chǔ)無擾過程N(yùn)t的參數(shù)估計(jì)結(jié)果如表1所示。由表1所示參數(shù)得到干預(yù)模型:它反映了干預(yù)事件對(duì)HFRS序列的影響。

    2.5與ARIMA模型比較

    為了考察干預(yù)分析模型的模擬效果,考慮與單純的使用ARIMA模型對(duì)比,兩模型的參數(shù)估計(jì)結(jié)果及檢驗(yàn)結(jié)果見表2。由表2可見干預(yù)分析殘差平方和與赤池信息量(AIC)都較小,且極大似然估計(jì)量也有所提高。因此選擇干預(yù)分析模型比較合適。

    表1 參數(shù)估計(jì)結(jié)果

    表2 干預(yù)分析模型與ARIMA模型比較

    3 總結(jié)

    利用帶干預(yù)的ARIMA模型預(yù)測(cè)武漢蔡甸區(qū)腎綜合征出血熱月發(fā)病數(shù)的變化趨勢(shì)時(shí),由于考慮了干預(yù)的影響,預(yù)測(cè)結(jié)果一般要比傳統(tǒng)ARIMA模型預(yù)測(cè)結(jié)果準(zhǔn)確,當(dāng)然也相對(duì)正確客觀。本文中采用干預(yù)分析模型對(duì)腎綜合征出血熱月發(fā)病數(shù)發(fā)病趨勢(shì)的變化進(jìn)行建模,得到了較精確的模型,為武漢市蔡甸區(qū)腎綜合征出血熱的預(yù)防和防治提供了一定的理論依據(jù)。

    [1]馮文權(quán),蔡基棟.干預(yù)分析及其應(yīng)用[J].武漢大學(xué)學(xué)報(bào),1994(6):1-7.

    [2]G E P Box,G C Tiao.Intervention Analysis with Applica?tions to Economic and Environmental Problems[J].Jour?nal of the American Statistical Association,1975,70 (349):70-79.

    [3]George E P Box,Gwilym M Jenkins,Gregory C Reinsel. Time Series Analysis Forecasting and Control[M].3版.北京:人民郵電出版社,2005:462-471.

    [4]Jonathan D.Cryer,Kung-Sik Chan.時(shí)間序列分析及應(yīng)用[M].潘紅宇,譯.2版.北京:機(jī)械工業(yè)出版社,2011.

    [5]馬亮亮,田富鵬.ARIMA模型在膽結(jié)石病發(fā)病率預(yù)測(cè)中的應(yīng)用[J].軍事醫(yī)學(xué),2010(5):469-472.

    [6]孫玉環(huán).ARIMA模型在測(cè)算重大突發(fā)事件影響中的應(yīng)用[J].統(tǒng)計(jì)與決策,2006(14):24-26.

    Application of Intervention Analysis in Simulation of Hemorrhagic Fever with Renal Syndrome

    Cheng Xingfu
    (School of Mathematics and Statistics,Wuhan University,Wuhan 430072,China)

    The intervention analysis model was established based on the monthly incidence data of Caid?ian District of Wuhan from 2009 to 1980 in order to simulate the feasibility of comprehensive kidney hemorrhagic fever incidence.And the simulation results were compared with the simulation results of the ARIMA model.It is found that the intervention analysis model can fit well the monthly incidence of HFRS.And this result provides scientific basis for the prevention and control of diseases.

    intervention analysis;ARIMA model;HFRS

    O212.1

    A

    1008-5483(2016)02-0074-03

    10.3969/j.issn.1008-5483.2016.02.017

    2016-04-27

    程幸福(1990-),男,安徽渦陽人,碩士生,從事隨機(jī)分析方面的研究。E-mail:1318247907@qq.com

    猜你喜歡
    新息分析模型均值
    基于BERT-VGG16的多模態(tài)情感分析模型
    傳遞函數(shù)辨識(shí)(21):線性回歸系統(tǒng)的遞階遞推參數(shù)估計(jì)
    M估計(jì)的強(qiáng)跟蹤SVD-UKF算法在組合導(dǎo)航中的應(yīng)用
    電子科技(2018年7期)2018-07-23 05:30:32
    層次分析模型在結(jié)核疾病預(yù)防控制系統(tǒng)中的應(yīng)用
    自適應(yīng)卡爾曼濾波在航空重力異常解算的應(yīng)用研究
    均值不等式失效時(shí)的解決方法
    均值與方差在生活中的應(yīng)用
    全啟發(fā)式語言分析模型
    基于新息正交性自適應(yīng)濾波的慣性/地磁組合導(dǎo)航方法
    關(guān)于均值有界變差函數(shù)的重要不等式
    18禁美女被吸乳视频| 精品久久久久久,| 久久精品91无色码中文字幕| 亚洲国产高清在线一区二区三 | 777久久人妻少妇嫩草av网站| 成人国产一区最新在线观看| 久久久久久久午夜电影| 欧美乱妇无乱码| 国内久久婷婷六月综合欲色啪| 无限看片的www在线观看| 午夜视频精品福利| 1024香蕉在线观看| 久久久久久久久久久久大奶| 成人三级做爰电影| 国产精品乱码一区二三区的特点 | 非洲黑人性xxxx精品又粗又长| 中文字幕久久专区| 欧美+亚洲+日韩+国产| 久久精品影院6| 黄片大片在线免费观看| 久久青草综合色| 亚洲熟女毛片儿| 美女午夜性视频免费| av在线天堂中文字幕| 热99re8久久精品国产| 亚洲欧美日韩高清在线视频| 欧美亚洲日本最大视频资源| 欧美乱妇无乱码| 精品国产一区二区久久| 波多野结衣巨乳人妻| 亚洲成人久久性| 99精品欧美一区二区三区四区| 日韩高清综合在线| 久久青草综合色| 麻豆久久精品国产亚洲av| 午夜免费观看网址| 真人一进一出gif抽搐免费| 国产日韩一区二区三区精品不卡| 亚洲精品中文字幕在线视频| 琪琪午夜伦伦电影理论片6080| 午夜精品在线福利| 看片在线看免费视频| 一进一出抽搐gif免费好疼| 午夜福利高清视频| 看黄色毛片网站| 精品电影一区二区在线| 久久久精品欧美日韩精品| 中文字幕人妻熟女乱码| 久久精品国产综合久久久| ponron亚洲| 日韩三级视频一区二区三区| 久久久国产成人免费| 亚洲一区高清亚洲精品| 日韩一卡2卡3卡4卡2021年| 婷婷丁香在线五月| 国产成人欧美在线观看| 日本三级黄在线观看| 男男h啪啪无遮挡| 欧美乱色亚洲激情| 亚洲av五月六月丁香网| 12—13女人毛片做爰片一| 不卡一级毛片| 午夜久久久久精精品| 亚洲在线自拍视频| 亚洲一区二区三区不卡视频| 久热爱精品视频在线9| 激情视频va一区二区三区| 麻豆成人av在线观看| 欧美乱码精品一区二区三区| 国产精品免费视频内射| 在线观看午夜福利视频| 亚洲精品粉嫩美女一区| 久久久久九九精品影院| 免费不卡黄色视频| 黑人欧美特级aaaaaa片| 亚洲人成网站在线播放欧美日韩| 午夜福利欧美成人| 亚洲熟妇熟女久久| 中亚洲国语对白在线视频| 久久中文字幕一级| 亚洲黑人精品在线| 免费不卡黄色视频| 桃色一区二区三区在线观看| 日韩有码中文字幕| 午夜福利成人在线免费观看| 国产99久久九九免费精品| 成人国产综合亚洲| 男女午夜视频在线观看| 亚洲熟妇中文字幕五十中出| 亚洲人成电影观看| √禁漫天堂资源中文www| 人人澡人人妻人| www.999成人在线观看| 国产三级黄色录像| 色播亚洲综合网| 亚洲成a人片在线一区二区| 久久国产亚洲av麻豆专区| 女人被躁到高潮嗷嗷叫费观| 巨乳人妻的诱惑在线观看| 亚洲无线在线观看| 在线观看免费午夜福利视频| 国产成人av激情在线播放| 亚洲一区中文字幕在线| 亚洲成人精品中文字幕电影| 国产伦人伦偷精品视频| 日韩有码中文字幕| 国产精品 欧美亚洲| 国产亚洲精品av在线| 一本大道久久a久久精品| netflix在线观看网站| 国产视频一区二区在线看| 国产高清激情床上av| 欧美日韩中文字幕国产精品一区二区三区 | 欧美中文日本在线观看视频| 国产激情欧美一区二区| 又黄又爽又免费观看的视频| 亚洲第一av免费看| 俄罗斯特黄特色一大片| 老司机午夜十八禁免费视频| 丁香六月欧美| ponron亚洲| 他把我摸到了高潮在线观看| 日韩欧美免费精品| 很黄的视频免费| 麻豆一二三区av精品| 少妇粗大呻吟视频| 欧美在线黄色| 美女高潮喷水抽搐中文字幕| 免费人成视频x8x8入口观看| 日本撒尿小便嘘嘘汇集6| 国产成人免费无遮挡视频| 久久国产亚洲av麻豆专区| 九色亚洲精品在线播放| 国产精品久久久久久亚洲av鲁大| 日韩成人在线观看一区二区三区| 亚洲av美国av| 国产成人免费无遮挡视频| 青草久久国产| 美女 人体艺术 gogo| 熟女少妇亚洲综合色aaa.| 国产亚洲av高清不卡| 亚洲人成网站在线播放欧美日韩| 999久久久国产精品视频| 一区二区日韩欧美中文字幕| 视频区欧美日本亚洲| 国产精品久久视频播放| 免费在线观看黄色视频的| 成人亚洲精品av一区二区| 欧美一级a爱片免费观看看 | 亚洲av成人av| 男女之事视频高清在线观看| 成人手机av| 亚洲第一欧美日韩一区二区三区| 久久人人爽av亚洲精品天堂| 午夜福利影视在线免费观看| 黑人巨大精品欧美一区二区mp4| 亚洲熟女毛片儿| 国产精品1区2区在线观看.| 老汉色∧v一级毛片| av视频免费观看在线观看| 日韩欧美免费精品| 午夜老司机福利片| 夜夜躁狠狠躁天天躁| 亚洲午夜精品一区,二区,三区| 乱人伦中国视频| 久久久水蜜桃国产精品网| 制服诱惑二区| 丰满人妻熟妇乱又伦精品不卡| 成人永久免费在线观看视频| 亚洲 国产 在线| 变态另类丝袜制服| 18禁国产床啪视频网站| 99在线视频只有这里精品首页| 黄色丝袜av网址大全| 国产免费男女视频| 亚洲久久久国产精品| 岛国在线观看网站| 好男人在线观看高清免费视频 | 身体一侧抽搐| 亚洲av电影在线进入| 精品熟女少妇八av免费久了| 美女 人体艺术 gogo| 九色国产91popny在线| 中文字幕av电影在线播放| 在线观看免费视频日本深夜| 国产在线观看jvid| 午夜福利免费观看在线| 久久九九热精品免费| 国内久久婷婷六月综合欲色啪| 精品久久久久久久人妻蜜臀av | 99热只有精品国产| 一本久久中文字幕| 亚洲欧洲精品一区二区精品久久久| 国产成人免费无遮挡视频| 免费在线观看影片大全网站| 一级a爱视频在线免费观看| 搡老岳熟女国产| 少妇被粗大的猛进出69影院| 日日干狠狠操夜夜爽| 国产97色在线日韩免费| 香蕉久久夜色| 国产男靠女视频免费网站| 午夜影院日韩av| 国产又爽黄色视频| 久久久久久国产a免费观看| 中文字幕人妻熟女乱码| 亚洲欧美日韩高清在线视频| 亚洲 欧美 日韩 在线 免费| 亚洲男人天堂网一区| 男人操女人黄网站| 欧美日韩精品网址| 大陆偷拍与自拍| 欧美大码av| 最好的美女福利视频网| 大香蕉久久成人网| 免费无遮挡裸体视频| 免费在线观看完整版高清| 欧美黄色淫秽网站| 欧美一区二区精品小视频在线| 黄频高清免费视频| 看黄色毛片网站| av视频免费观看在线观看| 亚洲人成伊人成综合网2020| 999久久久精品免费观看国产| 午夜亚洲福利在线播放| 亚洲 欧美 日韩 在线 免费| av在线天堂中文字幕| 黑人欧美特级aaaaaa片| 免费少妇av软件| 国产熟女午夜一区二区三区| 亚洲人成网站在线播放欧美日韩| 日韩av在线大香蕉| 欧美乱码精品一区二区三区| 国产欧美日韩精品亚洲av| 免费搜索国产男女视频| 亚洲精品中文字幕一二三四区| 久久久久久久久免费视频了| 亚洲欧美日韩无卡精品| 久久久精品国产亚洲av高清涩受| 亚洲精品中文字幕在线视频| 亚洲精品久久国产高清桃花| 国产精品国产高清国产av| www日本在线高清视频| 日韩三级视频一区二区三区| 精品第一国产精品| 免费少妇av软件| 操出白浆在线播放| 变态另类成人亚洲欧美熟女 | 日韩精品青青久久久久久| 日本精品一区二区三区蜜桃| 日韩欧美在线二视频| 好男人电影高清在线观看| 色播亚洲综合网| 热99re8久久精品国产| 美女免费视频网站| 精品久久久久久,| 少妇裸体淫交视频免费看高清 | 免费看a级黄色片| 国产成人精品在线电影| av福利片在线| 亚洲 欧美一区二区三区| 老司机在亚洲福利影院| 久久亚洲精品不卡| 最近最新免费中文字幕在线| 中文字幕最新亚洲高清| 中文亚洲av片在线观看爽| 国产精品一区二区三区四区久久 | 淫秽高清视频在线观看| 午夜久久久久精精品| 国产在线观看jvid| 在线av久久热| 亚洲av日韩精品久久久久久密| 久热这里只有精品99| 国产精品秋霞免费鲁丝片| 精品国产超薄肉色丝袜足j| а√天堂www在线а√下载| 日日爽夜夜爽网站| 精品乱码久久久久久99久播| 亚洲欧美精品综合久久99| 成人国产一区最新在线观看| 91国产中文字幕| 精品福利观看| 一个人观看的视频www高清免费观看 | 两个人看的免费小视频| 少妇被粗大的猛进出69影院| 亚洲第一欧美日韩一区二区三区| 亚洲精品久久成人aⅴ小说| 制服诱惑二区| 又紧又爽又黄一区二区| 日韩欧美免费精品| 久久 成人 亚洲| 亚洲,欧美精品.| 国产精品久久电影中文字幕| 美女大奶头视频| 国产成人av教育| 亚洲欧美激情在线| www.www免费av| 曰老女人黄片| 亚洲av第一区精品v没综合| 波多野结衣巨乳人妻| 99国产精品免费福利视频| 高清在线国产一区| 一区二区日韩欧美中文字幕| 欧美另类亚洲清纯唯美| 中文字幕人成人乱码亚洲影| 99久久综合精品五月天人人| 亚洲精品一区av在线观看| 后天国语完整版免费观看| 精品不卡国产一区二区三区| 视频在线观看一区二区三区| 麻豆成人av在线观看| 亚洲一码二码三码区别大吗| 天天一区二区日本电影三级 | 免费无遮挡裸体视频| 一区在线观看完整版| 中文字幕精品免费在线观看视频| 欧美日韩亚洲综合一区二区三区_| 国产1区2区3区精品| 亚洲av成人av| 窝窝影院91人妻| 中文字幕人妻丝袜一区二区| 欧美一级a爱片免费观看看 | 操美女的视频在线观看| 亚洲五月色婷婷综合| 大型av网站在线播放| 成人永久免费在线观看视频| 日韩欧美免费精品| 国产伦一二天堂av在线观看| 国产精品av久久久久免费| 亚洲精品av麻豆狂野| www.999成人在线观看| 少妇粗大呻吟视频| 男人的好看免费观看在线视频 | 国产亚洲精品综合一区在线观看 | 国产一区在线观看成人免费| 国产精品亚洲av一区麻豆| 神马国产精品三级电影在线观看 | 大香蕉久久成人网| 亚洲 欧美 日韩 在线 免费| 日韩大码丰满熟妇| 亚洲成人免费电影在线观看| 日韩高清综合在线| 夜夜看夜夜爽夜夜摸| 最近最新免费中文字幕在线| 一边摸一边做爽爽视频免费| a在线观看视频网站| 亚洲欧美精品综合一区二区三区| 久久久久久久精品吃奶| 午夜激情av网站| 亚洲av熟女| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲专区中文字幕在线| 国产欧美日韩综合在线一区二区| 波多野结衣高清无吗| 热99re8久久精品国产| 999久久久国产精品视频| 亚洲人成电影免费在线| 成人国产一区最新在线观看| 国产国语露脸激情在线看| ponron亚洲| 精品熟女少妇八av免费久了| 久久人人97超碰香蕉20202| 日韩有码中文字幕| 少妇粗大呻吟视频| 18禁黄网站禁片午夜丰满| 天天躁夜夜躁狠狠躁躁| 精品国产一区二区久久| 午夜影院日韩av| 久久久久久亚洲精品国产蜜桃av| 国产精品av久久久久免费| 黄色a级毛片大全视频| 国产精品一区二区三区四区久久 | 窝窝影院91人妻| 99在线人妻在线中文字幕| 亚洲欧美日韩高清在线视频| 国产精品九九99| 免费一级毛片在线播放高清视频 | 岛国在线观看网站| 一级毛片女人18水好多| 99国产精品免费福利视频| 午夜免费鲁丝| 欧美激情极品国产一区二区三区| 99久久99久久久精品蜜桃| 亚洲av电影在线进入| 99久久99久久久精品蜜桃| 黄色视频,在线免费观看| 在线观看舔阴道视频| 国内毛片毛片毛片毛片毛片| 无限看片的www在线观看| 亚洲av美国av| 变态另类丝袜制服| 高潮久久久久久久久久久不卡| 欧美黑人精品巨大| 久久久久精品国产欧美久久久| 国产精品久久久久久精品电影 | 国产精品99久久99久久久不卡| 久久久久久久久久久久大奶| 一级a爱视频在线免费观看| 一边摸一边抽搐一进一小说| 国产午夜福利久久久久久| 十八禁人妻一区二区| 精品国产亚洲在线| 操出白浆在线播放| 亚洲精品一区av在线观看| 国产精品亚洲av一区麻豆| 久久国产精品男人的天堂亚洲| 成人三级做爰电影| 国产欧美日韩一区二区三| 久久人妻av系列| 深夜精品福利| 啦啦啦观看免费观看视频高清 | 欧美黑人精品巨大| 9色porny在线观看| 色综合站精品国产| 又黄又粗又硬又大视频| 亚洲国产精品成人综合色| 国产精品免费视频内射| 男人的好看免费观看在线视频 | 手机成人av网站| 免费搜索国产男女视频| 成人特级黄色片久久久久久久| 国产av又大| 日韩免费av在线播放| 国产麻豆69| 人人妻人人澡欧美一区二区 | 亚洲,欧美精品.| 午夜免费鲁丝| 精品久久久久久久毛片微露脸| 中文字幕精品免费在线观看视频| 最近最新中文字幕大全免费视频| 一边摸一边做爽爽视频免费| 在线观看日韩欧美| 国产av在哪里看| 精品欧美一区二区三区在线| 一进一出好大好爽视频| 欧美一级a爱片免费观看看 | 中文字幕另类日韩欧美亚洲嫩草| 亚洲va日本ⅴa欧美va伊人久久| 一个人观看的视频www高清免费观看 | 在线观看一区二区三区| 午夜久久久久精精品| 久久精品影院6| 色综合婷婷激情| 久久中文字幕一级| 97人妻天天添夜夜摸| 久久久久九九精品影院| 伊人久久大香线蕉亚洲五| 欧美黄色淫秽网站| 精品一品国产午夜福利视频| 亚洲精品中文字幕在线视频| 国产麻豆成人av免费视频| 无遮挡黄片免费观看| 精品日产1卡2卡| 777久久人妻少妇嫩草av网站| 日本 欧美在线| 夜夜看夜夜爽夜夜摸| 在线观看免费午夜福利视频| 久久国产乱子伦精品免费另类| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av成人av| 亚洲国产日韩欧美精品在线观看 | 亚洲在线自拍视频| 国产亚洲欧美在线一区二区| 不卡av一区二区三区| 亚洲国产精品成人综合色| 亚洲欧美激情在线| bbb黄色大片| 亚洲自拍偷在线| 色尼玛亚洲综合影院| 国产精品一区二区三区四区久久 | 欧美激情久久久久久爽电影 | 日本五十路高清| 免费久久久久久久精品成人欧美视频| 午夜福利,免费看| 看片在线看免费视频| 久久中文字幕一级| 国产一卡二卡三卡精品| 人人澡人人妻人| 久久久国产精品麻豆| 不卡av一区二区三区| 99riav亚洲国产免费| 午夜福利高清视频| 亚洲av电影不卡..在线观看| 亚洲成av片中文字幕在线观看| 9191精品国产免费久久| 亚洲五月天丁香| 国产高清视频在线播放一区| 中出人妻视频一区二区| 男女床上黄色一级片免费看| 搡老熟女国产l中国老女人| 母亲3免费完整高清在线观看| 动漫黄色视频在线观看| 桃色一区二区三区在线观看| 国产欧美日韩一区二区三| 免费一级毛片在线播放高清视频 | 国产精品亚洲一级av第二区| 国产国语露脸激情在线看| 美女免费视频网站| 黄色a级毛片大全视频| 亚洲人成77777在线视频| 午夜激情av网站| 午夜精品久久久久久毛片777| 色在线成人网| 午夜福利视频1000在线观看 | 满18在线观看网站| 欧美成人午夜精品| 精品一区二区三区四区五区乱码| 男女下面进入的视频免费午夜 | 久热爱精品视频在线9| 99在线视频只有这里精品首页| 女人高潮潮喷娇喘18禁视频| 正在播放国产对白刺激| 黑丝袜美女国产一区| 97人妻天天添夜夜摸| 亚洲人成电影免费在线| 一a级毛片在线观看| 中文字幕久久专区| 精品免费久久久久久久清纯| 亚洲精品国产区一区二| 中文字幕色久视频| 十八禁网站免费在线| 国产精品亚洲美女久久久| 男人的好看免费观看在线视频 | 久久国产精品影院| 久久久精品国产亚洲av高清涩受| 一区二区日韩欧美中文字幕| 国产1区2区3区精品| 国产精品久久久久久人妻精品电影| 久久婷婷成人综合色麻豆| 大码成人一级视频| www.999成人在线观看| 老司机深夜福利视频在线观看| 国产亚洲av嫩草精品影院| av中文乱码字幕在线| 一边摸一边抽搐一进一出视频| 后天国语完整版免费观看| 日本欧美视频一区| 99精品在免费线老司机午夜| 亚洲av熟女| 女人被躁到高潮嗷嗷叫费观| 成人国语在线视频| 99精品欧美一区二区三区四区| 久久精品国产综合久久久| www.www免费av| 国产精品影院久久| 最新美女视频免费是黄的| 午夜精品国产一区二区电影| 十分钟在线观看高清视频www| 岛国视频午夜一区免费看| 好看av亚洲va欧美ⅴa在| 在线观看免费日韩欧美大片| av中文乱码字幕在线| 九色国产91popny在线| 久久久精品欧美日韩精品| 少妇 在线观看| 国产麻豆成人av免费视频| 男女做爰动态图高潮gif福利片 | 欧美黑人欧美精品刺激| 超碰成人久久| 久9热在线精品视频| 国产不卡一卡二| 真人做人爱边吃奶动态| 中亚洲国语对白在线视频| 久久欧美精品欧美久久欧美| 天堂影院成人在线观看| 成在线人永久免费视频| 一区二区三区激情视频| 欧美成狂野欧美在线观看| 欧美在线黄色| 好男人在线观看高清免费视频 | 男女下面插进去视频免费观看| 黄色成人免费大全| 午夜老司机福利片| 久久精品国产清高在天天线| 熟妇人妻久久中文字幕3abv| 国产熟女午夜一区二区三区| 欧美黑人精品巨大| 亚洲午夜理论影院| 欧美成狂野欧美在线观看| 啦啦啦免费观看视频1| 国产成人免费无遮挡视频| av视频在线观看入口| 91字幕亚洲| 午夜久久久久精精品| 午夜亚洲福利在线播放| 99久久久亚洲精品蜜臀av| svipshipincom国产片| 天天添夜夜摸| 国产99白浆流出| 亚洲精品国产一区二区精华液| 90打野战视频偷拍视频| 亚洲中文字幕日韩| 成熟少妇高潮喷水视频| 午夜福利免费观看在线| 亚洲 欧美 日韩 在线 免费| 亚洲欧美精品综合久久99| 久久精品人人爽人人爽视色| 88av欧美| 国产精品乱码一区二三区的特点 | 精品久久久久久,| 亚洲色图综合在线观看| 亚洲精品久久国产高清桃花| 天天躁狠狠躁夜夜躁狠狠躁| 国产99久久九九免费精品| 一级,二级,三级黄色视频| 在线视频色国产色| 国产成人免费无遮挡视频| 中国美女看黄片| 美女大奶头视频| 热re99久久国产66热| 18禁裸乳无遮挡免费网站照片 | 在线免费观看的www视频| 久久久久国产精品人妻aⅴ院| 久久久精品国产亚洲av高清涩受| 国产成年人精品一区二区|