• <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ù)的重要不等式
    男女无遮挡免费网站观看| 一区二区三区四区激情视频| 亚洲精品aⅴ在线观看| 成年人免费黄色播放视频 | av有码第一页| 伦理电影大哥的女人| 精品国产一区二区三区久久久樱花| 看免费成人av毛片| 中文字幕久久专区| 欧美丝袜亚洲另类| 美女视频免费永久观看网站| 久久99热6这里只有精品| 国产精品熟女久久久久浪| 国产黄片视频在线免费观看| 久久6这里有精品| 国语对白做爰xxxⅹ性视频网站| 51国产日韩欧美| 国产一级毛片在线| 久久久久久久亚洲中文字幕| 成人午夜精彩视频在线观看| 午夜福利,免费看| 久久女婷五月综合色啪小说| av在线观看视频网站免费| 国产精品秋霞免费鲁丝片| 亚洲美女黄色视频免费看| 亚洲怡红院男人天堂| av又黄又爽大尺度在线免费看| 天堂8中文在线网| 国产成人免费无遮挡视频| 婷婷色综合www| 国产亚洲精品久久久com| 最近2019中文字幕mv第一页| 性色avwww在线观看| 校园人妻丝袜中文字幕| 在线观看美女被高潮喷水网站| 午夜免费鲁丝| 免费观看性生交大片5| 国产av国产精品国产| 亚洲精品日本国产第一区| 极品教师在线视频| 国产欧美日韩精品一区二区| 亚洲成色77777| 久久人人爽人人片av| 亚洲一区二区三区欧美精品| 色吧在线观看| 亚洲内射少妇av| 久久久久久久精品精品| 亚洲,一卡二卡三卡| 女性被躁到高潮视频| 丝袜在线中文字幕| 日韩成人伦理影院| 亚洲成人av在线免费| 两个人的视频大全免费| 日本vs欧美在线观看视频 | av又黄又爽大尺度在线免费看| 日日爽夜夜爽网站| 中文字幕精品免费在线观看视频 | 久久久久久久久久久久大奶| 亚洲精品亚洲一区二区| 免费人妻精品一区二区三区视频| 欧美日韩视频精品一区| 亚洲精品乱久久久久久| 97超视频在线观看视频| 18禁在线播放成人免费| 久久久久久久精品精品| 高清毛片免费看| av有码第一页| 亚洲av福利一区| 国产精品嫩草影院av在线观看| 在线观看一区二区三区激情| 三级国产精品片| 少妇被粗大的猛进出69影院 | 一级毛片aaaaaa免费看小| 在线天堂最新版资源| 日韩中文字幕视频在线看片| 少妇 在线观看| 亚洲精品一二三| 自拍偷自拍亚洲精品老妇| 久久久久人妻精品一区果冻| 国产欧美日韩一区二区三区在线 | 国内揄拍国产精品人妻在线| 2021少妇久久久久久久久久久| 老女人水多毛片| 日本免费在线观看一区| 亚洲av福利一区| 亚洲国产精品成人久久小说| av女优亚洲男人天堂| 久久午夜综合久久蜜桃| 不卡视频在线观看欧美| 久久久国产精品麻豆| 久久国产精品大桥未久av | 国产成人免费无遮挡视频| 成年人午夜在线观看视频| 亚洲国产成人一精品久久久| 另类亚洲欧美激情| 国产黄频视频在线观看| 久久国产亚洲av麻豆专区| 如何舔出高潮| 天堂中文最新版在线下载| 国产国拍精品亚洲av在线观看| 人妻少妇偷人精品九色| 伊人久久精品亚洲午夜| 国产精品国产三级国产av玫瑰| 高清欧美精品videossex| 国产淫语在线视频| 哪个播放器可以免费观看大片| 免费观看无遮挡的男女| 极品少妇高潮喷水抽搐| 亚洲精品自拍成人| 在线观看免费高清a一片| 亚洲国产精品国产精品| 日韩一本色道免费dvd| av天堂久久9| 特大巨黑吊av在线直播| 色94色欧美一区二区| 最近最新中文字幕免费大全7| 99久国产av精品国产电影| 青春草视频在线免费观看| 久久免费观看电影| 免费观看无遮挡的男女| 在线观看免费视频网站a站| 亚洲真实伦在线观看| 久久青草综合色| 97在线人人人人妻| freevideosex欧美| 好男人视频免费观看在线| 国产高清三级在线| 黄片无遮挡物在线观看| 亚洲va在线va天堂va国产| 日韩电影二区| 精华霜和精华液先用哪个| 美女内射精品一级片tv| 日韩视频在线欧美| 久久国产乱子免费精品| 18禁在线无遮挡免费观看视频| 久久精品国产鲁丝片午夜精品| 国产日韩欧美在线精品| 婷婷色综合www| 色视频www国产| 成人综合一区亚洲| 日韩视频在线欧美| 黄色日韩在线| 国产免费一级a男人的天堂| 国产男人的电影天堂91| 精品卡一卡二卡四卡免费| 特大巨黑吊av在线直播| 欧美97在线视频| 亚洲欧美一区二区三区黑人 | 欧美bdsm另类| 成人18禁高潮啪啪吃奶动态图 | 国产精品国产三级国产av玫瑰| 精品久久久久久电影网| 少妇猛男粗大的猛烈进出视频| 99热这里只有是精品在线观看| 国产午夜精品一二区理论片| 男的添女的下面高潮视频| 爱豆传媒免费全集在线观看| 国产免费视频播放在线视频| 五月天丁香电影| 在线亚洲精品国产二区图片欧美 | 国产精品不卡视频一区二区| 免费观看在线日韩| 亚洲怡红院男人天堂| 久久影院123| av女优亚洲男人天堂| 嘟嘟电影网在线观看| av国产久精品久网站免费入址| 亚洲国产精品国产精品| 肉色欧美久久久久久久蜜桃| 在线观看三级黄色| 日韩亚洲欧美综合| 亚洲av在线观看美女高潮| 女人久久www免费人成看片| av福利片在线| 国精品久久久久久国模美| 国产乱人偷精品视频| 高清午夜精品一区二区三区| 高清av免费在线| 日韩一区二区三区影片| 国产日韩欧美亚洲二区| 欧美区成人在线视频| 免费不卡的大黄色大毛片视频在线观看| 国产在视频线精品| 在线观看一区二区三区激情| 边亲边吃奶的免费视频| 一本—道久久a久久精品蜜桃钙片| 卡戴珊不雅视频在线播放| 妹子高潮喷水视频| 久久久久久久久久成人| 久久热精品热| 久久久久久伊人网av| 亚洲情色 制服丝袜| 简卡轻食公司| 尾随美女入室| 中国三级夫妇交换| 香蕉精品网在线| 成人亚洲精品一区在线观看| 免费观看在线日韩| 熟女av电影| 亚洲成人av在线免费| 草草在线视频免费看| 男女国产视频网站| 中文乱码字字幕精品一区二区三区| 精品亚洲乱码少妇综合久久| 精品少妇黑人巨大在线播放| 欧美亚洲 丝袜 人妻 在线| 噜噜噜噜噜久久久久久91| 中文乱码字字幕精品一区二区三区| 国产色爽女视频免费观看| 久久久久久久久久久免费av| 国产精品人妻久久久久久| 亚洲国产精品国产精品| 欧美老熟妇乱子伦牲交| 午夜福利在线观看免费完整高清在| 国产一区亚洲一区在线观看| 亚洲国产精品成人久久小说| 亚洲三级黄色毛片| 国产免费福利视频在线观看| 欧美精品国产亚洲| 亚洲国产色片| 国产精品人妻久久久久久| 亚洲综合色惰| 国产日韩欧美视频二区| 午夜福利网站1000一区二区三区| 午夜老司机福利剧场| kizo精华| 精品久久久久久久久av| 精品亚洲成a人片在线观看| av播播在线观看一区| 久久国内精品自在自线图片| 免费av不卡在线播放| 亚洲精品国产成人久久av| 97在线人人人人妻| 久久精品久久久久久噜噜老黄| 熟女人妻精品中文字幕| 中文字幕人妻丝袜制服| 六月丁香七月| 国产探花极品一区二区| 大片电影免费在线观看免费| 高清在线视频一区二区三区| 婷婷色综合大香蕉| 欧美丝袜亚洲另类| 2022亚洲国产成人精品| 精品人妻熟女av久视频| 精品熟女少妇av免费看| 18禁在线无遮挡免费观看视频| a级片在线免费高清观看视频| 亚洲一级一片aⅴ在线观看| 热re99久久精品国产66热6| 中文字幕人妻丝袜制服| av国产久精品久网站免费入址| 人妻一区二区av| 欧美精品国产亚洲| 伊人久久国产一区二区| 街头女战士在线观看网站| 亚洲精品国产av成人精品| 成年人免费黄色播放视频 | 亚洲国产精品成人久久小说| 亚洲av二区三区四区| 久久影院123| 夜夜爽夜夜爽视频| 亚洲综合色惰| 国产精品久久久久久久久免| 国产精品福利在线免费观看| 精品一区二区三区视频在线| 久久人妻熟女aⅴ| 日日摸夜夜添夜夜爱| 精品亚洲成a人片在线观看| 国产av精品麻豆| 国产色婷婷99| 国产成人freesex在线| av免费观看日本| 青青草视频在线视频观看| 国产亚洲一区二区精品| 色5月婷婷丁香| 日韩精品有码人妻一区| 热99国产精品久久久久久7| 精品一区二区三卡| 男的添女的下面高潮视频| 久久久久人妻精品一区果冻| 欧美成人午夜免费资源| 美女大奶头黄色视频| 99九九在线精品视频 | 一级a做视频免费观看| 欧美3d第一页| 自线自在国产av| 曰老女人黄片| 精品久久久噜噜| 免费观看av网站的网址| 亚洲av.av天堂| 丰满少妇做爰视频| 又爽又黄a免费视频| 日本-黄色视频高清免费观看| 亚州av有码| 丝袜脚勾引网站| 日韩成人伦理影院| 日本wwww免费看| 中文字幕av电影在线播放| 一本一本综合久久| 欧美一级a爱片免费观看看| 成年美女黄网站色视频大全免费 | 亚洲综合色惰| 亚洲欧美中文字幕日韩二区| 国产精品人妻久久久久久| 久久久久视频综合| 丝袜喷水一区| 十分钟在线观看高清视频www | 蜜臀久久99精品久久宅男| 久久久久网色| 久久99一区二区三区| 人妻人人澡人人爽人人| 国产av码专区亚洲av| 国产免费视频播放在线视频| 亚洲中文av在线| 天堂8中文在线网| 国产精品人妻久久久久久| 九色成人免费人妻av| 精品亚洲成国产av| 各种免费的搞黄视频| 久久久久久久大尺度免费视频| 久久久久久久久久久免费av| 一级毛片aaaaaa免费看小| 免费观看av网站的网址| 国产成人免费观看mmmm| 国产精品久久久久久精品电影小说| 午夜免费观看性视频| 看非洲黑人一级黄片| 黄色毛片三级朝国网站 | 在线观看av片永久免费下载| 一级毛片电影观看| 18禁动态无遮挡网站| 插逼视频在线观看| 国产视频首页在线观看| av不卡在线播放| 97超碰精品成人国产| 极品人妻少妇av视频| 麻豆成人午夜福利视频| 国产高清有码在线观看视频| 99热6这里只有精品| 亚洲成人av在线免费| 日韩av免费高清视频| 中文字幕制服av| 国产免费一级a男人的天堂| 又大又黄又爽视频免费| 六月丁香七月| 男人添女人高潮全过程视频| av在线观看视频网站免费| 王馨瑶露胸无遮挡在线观看| 成人午夜精彩视频在线观看| 一区二区av电影网| 97精品久久久久久久久久精品| 老司机影院毛片| 国产亚洲最大av| tube8黄色片| 一区在线观看完整版| 亚洲国产色片| 免费看光身美女| 国产视频首页在线观看| 国产精品麻豆人妻色哟哟久久| 交换朋友夫妻互换小说| 国产成人a∨麻豆精品| 亚洲精品一二三| 九色成人免费人妻av| 精品亚洲乱码少妇综合久久| 日韩一本色道免费dvd| 六月丁香七月| 久久久欧美国产精品| 一区二区av电影网| 亚洲精品中文字幕在线视频 | 欧美精品一区二区免费开放| 男人爽女人下面视频在线观看| 成人亚洲欧美一区二区av| 9色porny在线观看| 中文字幕久久专区| 日韩制服骚丝袜av| 狂野欧美激情性bbbbbb| 国产av码专区亚洲av| 偷拍熟女少妇极品色| 高清毛片免费看| 日韩一本色道免费dvd| 成人无遮挡网站| 中国美白少妇内射xxxbb| 日韩大片免费观看网站| 亚洲av综合色区一区| 国产精品一区二区性色av| 一区二区三区四区激情视频| 午夜视频国产福利| 免费av中文字幕在线| 80岁老熟妇乱子伦牲交| 我要看日韩黄色一级片| 高清视频免费观看一区二区| 91久久精品国产一区二区成人| www.av在线官网国产| 日日摸夜夜添夜夜爱| 国产黄片美女视频| 中文字幕精品免费在线观看视频 | 日韩熟女老妇一区二区性免费视频| 日日爽夜夜爽网站| 九九在线视频观看精品| 久久久久网色| 汤姆久久久久久久影院中文字幕| 91精品国产九色| 免费播放大片免费观看视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品乱久久久久久| 久热这里只有精品99| 中国国产av一级| 一个人看视频在线观看www免费| 黄色怎么调成土黄色| 精品酒店卫生间| 美女福利国产在线| 国产精品久久久久久久电影| 在线 av 中文字幕| 免费高清在线观看视频在线观看| 女的被弄到高潮叫床怎么办| 九九久久精品国产亚洲av麻豆| 精品亚洲成a人片在线观看| 国产精品人妻久久久影院| 亚洲国产毛片av蜜桃av| 亚洲欧美成人综合另类久久久| 国产在视频线精品| 最近手机中文字幕大全| 亚洲欧美成人精品一区二区| av又黄又爽大尺度在线免费看| 亚洲精品456在线播放app| 国产免费福利视频在线观看| 好男人视频免费观看在线| 国产精品99久久久久久久久| 国产亚洲午夜精品一区二区久久| 观看免费一级毛片| 日本wwww免费看| 欧美日韩亚洲高清精品| 国产精品国产三级专区第一集| 97精品久久久久久久久久精品| 黄片无遮挡物在线观看| 亚洲国产精品国产精品| 国产精品免费大片| 99re6热这里在线精品视频| 国产成人a∨麻豆精品| 日本免费在线观看一区| 免费黄色在线免费观看| 免费观看a级毛片全部| 亚洲精品一区蜜桃| 一级a做视频免费观看| 亚洲三级黄色毛片| 国产亚洲欧美精品永久| 国产亚洲av片在线观看秒播厂| 美女国产视频在线观看| 久久精品国产鲁丝片午夜精品| 久久久亚洲精品成人影院| 日本vs欧美在线观看视频 | 丁香六月天网| 桃花免费在线播放| 日本91视频免费播放| 亚洲综合色惰| 99热这里只有是精品50| 日韩三级伦理在线观看| 日本91视频免费播放| 日韩三级伦理在线观看| 国产色婷婷99| 国产成人精品婷婷| 国产在线一区二区三区精| 搡女人真爽免费视频火全软件| 在线亚洲精品国产二区图片欧美 | 日韩三级伦理在线观看| 肉色欧美久久久久久久蜜桃| 少妇裸体淫交视频免费看高清| 国产精品国产av在线观看| 精品一区在线观看国产| 少妇的逼水好多| 久久ye,这里只有精品| 久久久欧美国产精品| 狂野欧美激情性xxxx在线观看| 在线播放无遮挡| 亚洲情色 制服丝袜| 麻豆成人av视频| 免费大片18禁| 国产探花极品一区二区| 观看美女的网站| 欧美激情国产日韩精品一区| 国产日韩一区二区三区精品不卡 | 如何舔出高潮| 91aial.com中文字幕在线观看| 建设人人有责人人尽责人人享有的| 日韩av在线免费看完整版不卡| 精品人妻熟女av久视频| 午夜福利网站1000一区二区三区| 黑丝袜美女国产一区| 婷婷色综合www| 午夜日本视频在线| 高清黄色对白视频在线免费看 | 久久精品国产自在天天线| a级毛片在线看网站| 亚洲av成人精品一二三区| 国产成人精品无人区| 精品国产露脸久久av麻豆| 色吧在线观看| 色婷婷久久久亚洲欧美| 久久久欧美国产精品| 国产亚洲5aaaaa淫片| 免费观看性生交大片5| 黑人巨大精品欧美一区二区蜜桃 | 菩萨蛮人人尽说江南好唐韦庄| 久久 成人 亚洲| 日韩精品免费视频一区二区三区 | 精品久久久久久久久av| 亚洲成人一二三区av| 亚洲精品亚洲一区二区| 2021少妇久久久久久久久久久| 国产又色又爽无遮挡免| 亚洲国产精品专区欧美| 国产黄色视频一区二区在线观看| 国产黄片美女视频| 欧美性感艳星| 中文字幕免费在线视频6| 午夜日本视频在线| 一区二区三区乱码不卡18| 五月开心婷婷网| 最近中文字幕高清免费大全6| 男人添女人高潮全过程视频| 又大又黄又爽视频免费| 青春草国产在线视频| 人妻夜夜爽99麻豆av| 美女福利国产在线| 欧美老熟妇乱子伦牲交| 中文精品一卡2卡3卡4更新| 亚洲精品乱码久久久v下载方式| 日本与韩国留学比较| 亚洲欧美清纯卡通| 大香蕉久久网| 日韩强制内射视频| 亚洲丝袜综合中文字幕| 亚洲精品乱久久久久久| 九九爱精品视频在线观看| 一区在线观看完整版| 狂野欧美白嫩少妇大欣赏| 国产成人精品无人区| 久久久久网色| 免费观看av网站的网址| 国产成人a∨麻豆精品| 在线精品无人区一区二区三| 日韩中文字幕视频在线看片| 欧美成人午夜免费资源| 日韩中字成人| 亚洲精品第二区| 欧美精品一区二区大全| 亚洲av免费高清在线观看| 国产黄片美女视频| 91在线精品国自产拍蜜月| 欧美精品高潮呻吟av久久| 在线观看av片永久免费下载| 插阴视频在线观看视频| 日本黄色日本黄色录像| 日韩成人av中文字幕在线观看| 午夜免费观看性视频| 久久久久视频综合| 国产乱人偷精品视频| 婷婷色综合大香蕉| 国产视频内射| 91午夜精品亚洲一区二区三区| av不卡在线播放| 高清午夜精品一区二区三区| 在线观看免费高清a一片| 伊人久久精品亚洲午夜| 综合色丁香网| 精品少妇黑人巨大在线播放| 高清不卡的av网站| 伊人久久精品亚洲午夜| 免费观看性生交大片5| 看免费成人av毛片| 9色porny在线观看| 亚洲性久久影院| 天堂8中文在线网| 久久久久久久久久成人| 亚洲欧洲日产国产| 少妇的逼水好多| 午夜激情久久久久久久| 精品亚洲成a人片在线观看| av免费观看日本| 在线观看www视频免费| av播播在线观看一区| 人妻系列 视频| 边亲边吃奶的免费视频| 曰老女人黄片| 男人舔奶头视频| 建设人人有责人人尽责人人享有的| 亚洲精品亚洲一区二区| 亚洲人与动物交配视频| 免费高清在线观看视频在线观看| 97精品久久久久久久久久精品| 日韩人妻高清精品专区| 国产精品嫩草影院av在线观看| 成人毛片60女人毛片免费| 国产精品嫩草影院av在线观看| 男的添女的下面高潮视频| 亚洲国产精品成人久久小说| 色网站视频免费| 久久午夜福利片| .国产精品久久| av卡一久久| 老司机影院毛片| 亚洲四区av| 22中文网久久字幕| 国产日韩欧美视频二区| 亚洲激情五月婷婷啪啪| 日本黄色片子视频| 日本与韩国留学比较| 亚洲成人一二三区av| 免费观看a级毛片全部| 国产欧美日韩一区二区三区在线 | 亚洲av中文av极速乱| 国产黄色视频一区二区在线观看| 这个男人来自地球电影免费观看 | 3wmmmm亚洲av在线观看| 晚上一个人看的免费电影|