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

    BBN模型在江蘇居民氰戊菊酯長期膳食暴露評估中的應(yīng)用*

    2012-07-27 09:30:06東南大學(xué)公共衛(wèi)生學(xué)院流行病與衛(wèi)生統(tǒng)計教研室210009王燦楠
    中國衛(wèi)生統(tǒng)計 2012年5期
    關(guān)鍵詞:氰戊菊酯貝塔正態(tài)分布

    東南大學(xué)公共衛(wèi)生學(xué)院流行病與衛(wèi)生統(tǒng)計教研室(210009) 田 野 閔 捷 劉 沛 王燦楠

    在進行污染物的長期膳食暴露評估時,合理地估計每人每天長期攝入某種物質(zhì)平均量(usual intake,簡稱日常攝入量)往往是風(fēng)險評估過程中的重點。日常攝入量的計算涉及該物質(zhì)在各種食物中的濃度和人群消費各種食物的量,前者可用各種食物污染物檢測的平均濃度估計,不考慮其變異;后者需要采集人群長期消費各種食物的信息,實際工作中是難以實現(xiàn)的,如何利用短期2~7天的24小時膳食回顧調(diào)查數(shù)據(jù)來估計人群長期暴露量就顯得尤為重要。許多專家及機構(gòu)提出了估計及改進的方法〔1-3〕,目前大家普遍接受的是貝塔二項-正態(tài)分布模型(Beta binomial-normal,BBN),該方法是2005年由De Boer和Van der Voe等人提出的〔4〕。其基本思想是考慮調(diào)查人群中,個體間變異可近似反映人群的長期飲食情況,個體內(nèi)變異則被認為僅反映調(diào)查個體在調(diào)查的n天內(nèi)短期膳食結(jié)構(gòu)的波動,采用合適的統(tǒng)計模型去除個體內(nèi)變異,保留個體間變異的暴露量,近似地作為人群的長期暴露量。本文利用BBN模型對江蘇省各年齡、性別組人群氰戊菊酯的長期膳食安全風(fēng)險水平進行評估,以期為農(nóng)藥膳食長期暴露風(fēng)險評估提供科學(xué)依據(jù)。

    BBN模型原理

    BBN模型包括攝入頻率和陽性暴露量兩部分。其中,貝塔二項分布被用來構(gòu)建攝入頻率模型,對數(shù)正態(tài)分布被用來構(gòu)建陽性暴露量模型。

    1.攝入頻率模型

    假設(shè)n和npos分別代表每個調(diào)查對象的總調(diào)查天數(shù)和攝入食物中含有該污染物的天數(shù)(即陽性攝入天數(shù))。由于每人每天攝入的食物中可能包含污染物,也可能不包含污染物,因此陽性攝入天數(shù)npos服從陽性攝入概率為p,重復(fù)數(shù)為n的二項分布;其中,陽性攝入概率p的取值在0到1之間,服從貝塔分布;那么,陽性攝入天數(shù)就服從貝塔二項分布:

    貝塔分布的均數(shù)和方差分別為

    令π=α/(α+β),φ=1/(α+β+1),則該分布的均數(shù)和方差可表示為:

    通過最大似然估計法擬合貝塔二項分布模型??梢缘玫溅泻挺盏墓烙嬛郸泻挺眨瑥亩赐频玫溅梁挺碌墓啦閷ο笤谀骋徽{(diào)查日食用某種事物的陽性攝入頻率分布即可用貝塔分布

    2.陽性暴露量模型

    長期膳食暴露評估主要反應(yīng)人群中某污染物日常攝入情況,假設(shè)陽性暴露量即調(diào)查對象每人每天膳食污染物暴露量不為零,經(jīng)過對數(shù)變換后近似服從正態(tài)分布,變換后的陽性暴露量,通過最大似然估計法擬合線性隨機效應(yīng)模型(linear random effects model):

    其中,transf(yposij)表示經(jīng)對數(shù)變換后的陽性暴露量,ci和uij表示人間效應(yīng)和人內(nèi)效應(yīng),其二者分布分別滿足正態(tài)分布 N(0,σ2B)和N(0,σ2w)。通過最大似然估計法可以得到三個參數(shù)的估計值m,分別代表個體平均值、人間效應(yīng)估計值及人內(nèi)效應(yīng)估計值。

    3.計算長期暴露量分布

    長期暴露量是通過對陽性暴露量分布和攝入頻率分布進行Monte Carlo抽樣,然后進行整合獲得。其具體過程為:

    首先通過模擬獲得對數(shù)轉(zhuǎn)換后的陽性暴露量m+zi,其中zi從正態(tài)分布N(0)中隨機抽取10萬次,然后反對數(shù)變換獲得陽性暴露量:

    實例分析

    1.資料情況主要涉及的數(shù)據(jù)有膳食調(diào)查數(shù)據(jù)、污染物常規(guī)監(jiān)測數(shù)據(jù)、研究人群人口學(xué)基本情況數(shù)據(jù)。膳食調(diào)查及人口學(xué)數(shù)據(jù):來源于2002年全國居民營養(yǎng)與健康狀況調(diào)查中的江蘇省數(shù)據(jù);污染物監(jiān)測數(shù)據(jù):來源于2000~2006年全國14個省/地區(qū)食品污染物監(jiān)測網(wǎng)數(shù)據(jù)及2005~2006年海關(guān)出口農(nóng)產(chǎn)品監(jiān)測數(shù)據(jù)中的江蘇省數(shù)據(jù)。本研究資料中氰戊菊酯檢出率為3.10%,按照WHO關(guān)于污染物未檢出值處理的原則,本文中未檢出值均采用0替換〔5〕。目前,F(xiàn)AO/WHO食品添加劑聯(lián)合專家委員會(JECFA)制定的食品中每日允許攝入氰戊菊酯量(ADI)為2μg/kg bw/day。將模型結(jié)果與之進行比較,評價人群的氰戊菊酯攝入水平。

    2.結(jié)果與分析

    應(yīng)用SAS9.1統(tǒng)計軟件進行數(shù)據(jù)整理和膳食暴露量的計算,結(jié)果如下。

    (1)BBN模型參數(shù)

    利用SAS的Proc nlmixed過程的最大似然估計法迭代獲得貝塔分布參數(shù)和的估計值。在對陽性暴露量進行方差分解時,首先將陽性暴露量進行對數(shù)變換,使其近似正態(tài)分布,然后利用SAS的Proc mixed過程將殘差分解為個體間效應(yīng)和個體內(nèi)效應(yīng),通過最大似然估計法獲得個體平均值m、人間效應(yīng)估計值及人內(nèi)效應(yīng)估計值。計算的氰戊菊酯膳食暴露BBN模型參數(shù)見下表。

    表1 氰戊菊酯膳食暴露BBN評估模型參數(shù)

    (2)氰戊菊酯長期膳食暴露評估結(jié)果

    如表2可見,氰戊菊酯的長期暴露量在平均值及各百分位數(shù)均低于其每日允許攝入量ADI(2 μg/kg bw/day),說明氰戊菊酯殘留在長期膳食暴露評估中是相對安全的。

    表2 不同人群氰戊菊酯長期膳食暴露量(單位μg/kg bw/day)

    (3)氰戊菊酯膳食暴露量的不確定性分析

    為了解人群氰戊菊酯暴露變異的各百分位數(shù)的可信區(qū)間,采用Bootstrap方法有放回的從原樣本中進行抽樣,并對抽取的每一個Bootstrap樣本進行Monte Carlo模擬以計算暴露量百分位數(shù),經(jīng)過n次Bootstrap抽樣計算百分位數(shù)分布及可信區(qū)間,代表百分位數(shù)的不確定性。本文取200次Bootstrap抽樣進行不確定性分析。

    由表3可見,各暴露水平的暴露量均位于95%可信區(qū)間內(nèi),說明暴露量百分位數(shù)能夠達到穩(wěn)定。由于本次研究采用的是全比例抽樣,其不確定性主要存在于攝入頻率p~Beta(α,β)和人間效應(yīng)zi~N(0)的抽樣中,同時Monte Carlo模擬較好的量化了變異性,故在可信區(qū)間內(nèi)變異較小。

    討 論

    1.長期膳食暴露評估模型

    評價污染物長期膳食暴露有諸多不同的數(shù)學(xué)模型,其共同思想是考慮在大規(guī)模人群短期橫斷面膳食調(diào)查情況下,數(shù)據(jù)資料中變異主要包括個體間變異和個體內(nèi)變異。不同個體間的變異可以近似反映人群的長期膳食狀況,而個體內(nèi)變異則被看做僅反映調(diào)查個體在橫斷面調(diào)查中短期的膳食結(jié)構(gòu)波動,因此,通過采用合適的數(shù)學(xué)模型來去除個體內(nèi)變異,而保留個體間變異,從而達到用污染物短期的膳食暴露量近似代表長期的膳食暴露量的目的。本研究在計算長期膳食暴露評估時采用歐盟MCRA(Monte Carlo Risk Assessment)推薦的BBN模型。將長期的污染物膳食陽性攝入頻率用短期陽性攝入頻率的貝塔二項分布模型來表示;同時,用正態(tài)分布構(gòu)建長期陽性暴露量模型。從分布中抽樣,進行Monte Carlo模擬獲得污染物的長期暴露量。利用BBN模型進行長期膳食暴露評估的優(yōu)點是將個體作為研究對象,通過統(tǒng)計方法進行模擬,使結(jié)果更符合實際并獲得人群的暴露量分布,其對膳食調(diào)查數(shù)據(jù)及污染物殘留數(shù)據(jù)的信息利用程度遠大于點評估模型。國外一些研究進行了模型的比較也證明了這一結(jié)論〔6〕。當然,此種暴露評估方法也有其局限性,如模擬次數(shù)足夠多才能使最終結(jié)果穩(wěn)定及對數(shù)據(jù)質(zhì)量要求較高等。另外該模型適用于污染物檢出率不高的暴露評估,如果檢出率高,例如鉛、鎘等的污染物,則BBN中可不考慮Beta分布。

    表3 不同人群氰戊菊酯長期膳食暴露量不確定性分析(單位μg/kg bw/day)

    2.關(guān)于陽性暴露量的正態(tài)性轉(zhuǎn)換

    BBN模型中,需要利用正態(tài)分布構(gòu)建陽性攝入量模型。人群的陽性攝入量的原始數(shù)據(jù)往往不滿足正態(tài)性,常常需要用變量變換的方法將偏態(tài)資料正態(tài)化。在BBN模型中常用的變量變換有兩種:對數(shù)變換和冪變換〔7〕。經(jīng)過變換后的變量,一般需對其分布進行擬合優(yōu)度檢驗,常用的檢驗有Anderson-Darling檢驗、Kolmogorov-Smirnov檢驗、Chisquare檢驗等。本研究中,陽性暴露量的原始數(shù)據(jù)表現(xiàn)為正偏態(tài)的基本特征,本研究采用對數(shù)變換對陽性攝入量的原始數(shù)據(jù)進行轉(zhuǎn)換。

    3.Bootstrap抽樣次數(shù)的確定

    關(guān)于Bootstrap抽樣次數(shù)的確定,一般根據(jù)研究目的的不同取不同的值,如果用來計算均值、方差、標準誤等統(tǒng)計量,200次抽樣即可,但如果計算極端百分位數(shù)如99.9th的可信區(qū)間,Bootstrap的模擬次數(shù)應(yīng)當增加到至少500次〔8〕。本研究Monte Carlo模擬的10萬個人天數(shù)作為一個Bootstrap樣本含量,根據(jù)需要估計的百分位數(shù)情況,將Bootstrap樣本數(shù)定為200次、500次、1000次分別進行模擬,比較抽樣結(jié)果發(fā)現(xiàn)可信區(qū)間寬度沒有明顯差異,最終確定本研究中Bootstrap模擬次數(shù)為200次比較合適。需要說明的是,對不確定度更多的應(yīng)從專業(yè)上加以評價,如果風(fēng)險管理者認為數(shù)據(jù)的不確定度過大,可以拒絕概率估計的結(jié)論,或者采取一定措施(如提高數(shù)據(jù)質(zhì)量)來減少不確定度〔9〕。

    1.Nusser SM,Carriquiry AL,Dodd KW,et al.A semi-parametric transformation approach to estimating usual intake daily intake distributions.J Am Stat Assoc,1996,91:1440-1449.

    2.Slob W.Modeling Long-term exposure of the whole population to chemicals in food.RISK ANAL,1993,13(5):525-530.

    3.Slob W.Probabilistic dietary exposure assessment taking into account variability in both amount and frequency of consumption.Food and Chemical Toxicology,2006,44(7):933-951.

    4.De Boer WJ,van der Voet H.MCRA,Release 5,a web-based program for Monte Carlo Risk Assessment.Biometris and RIKILT,Wageningen University and Research Centre,2006.

    5.王緒卿,吳永寧.食品污染監(jiān)測低水平數(shù)據(jù)處理問題.中華預(yù)防醫(yī)學(xué)雜志,2002:837-852.

    6.Boon PE,Bonthuis M.Comparison of different exposure assessment methods to estimate the long-term dietary exposure to dioxins and ochratoxin A.Food and Chemical Toxicology,2011,49(9):1979-1988.

    7.De Boer WJ,van der Voet H.Comparison of two models for the estimation of usual intake addressing zero consumption and non-normality.Food Additives and Contaminants,2009,26(11):1433-1449.

    8.Cullen AC,F(xiàn)rey HC.Probabilistic Techniques in Exposure Assessment:A Handbook for Dealing with Variability and Uncertainty in Models and Inputs.Society for Risk Analysis,1999.

    9.Thompson KM.Variability and uncertainty meet risk management and risk communication.Risk Anal,2002,22:647-654.

    猜你喜歡
    氰戊菊酯貝塔正態(tài)分布
    用于MEMS電源的貝塔輻伏電池研究進展
    云南化工(2021年8期)2021-12-21 06:37:12
    貝塔魚
    下雨啦!
    氰戊菊酯超標茶葉
    基于對數(shù)正態(tài)分布的出行時長可靠性計算
    正態(tài)分布及其應(yīng)用
    氰戊菊酯降解菌的篩選與鑒定及其降解條件優(yōu)化
    正態(tài)分布題型剖析
    χ2分布、t 分布、F 分布與正態(tài)分布間的關(guān)系
    這件事我做對了
    丁香六月欧美| av欧美777| 国产午夜精品久久久久久| 午夜福利在线观看吧| 自拍欧美九色日韩亚洲蝌蚪91| 真人做人爱边吃奶动态| 好男人电影高清在线观看| 99国产精品免费福利视频| 国产精品 国内视频| 午夜a级毛片| 国产三级黄色录像| 在线观看免费午夜福利视频| 中文字幕久久专区| 黄色 视频免费看| 亚洲成人国产一区在线观看| 一区二区三区精品91| √禁漫天堂资源中文www| 搞女人的毛片| 久久精品aⅴ一区二区三区四区| 久久人妻熟女aⅴ| 久久这里只有精品19| 亚洲色图综合在线观看| 韩国精品一区二区三区| 精品久久久久久久毛片微露脸| 国产熟女午夜一区二区三区| 黄色毛片三级朝国网站| 国产欧美日韩一区二区三区在线| 久久香蕉国产精品| 老司机午夜十八禁免费视频| 中文字幕av电影在线播放| 亚洲精品在线美女| 19禁男女啪啪无遮挡网站| 一区二区三区精品91| 91大片在线观看| 久热爱精品视频在线9| 精品国产一区二区三区四区第35| 亚洲五月天丁香| 精品午夜福利视频在线观看一区| av电影中文网址| 女人高潮潮喷娇喘18禁视频| 制服丝袜大香蕉在线| 桃色一区二区三区在线观看| 欧美精品啪啪一区二区三区| 久久狼人影院| a级毛片在线看网站| 校园春色视频在线观看| 激情在线观看视频在线高清| 免费av毛片视频| 亚洲精华国产精华精| 美女扒开内裤让男人捅视频| 国产精品亚洲一级av第二区| 变态另类丝袜制服| 免费在线观看亚洲国产| 天天一区二区日本电影三级 | 国产精品 国内视频| 国产亚洲欧美精品永久| 精品少妇一区二区三区视频日本电影| 男女床上黄色一级片免费看| 久99久视频精品免费| 成人欧美大片| 99国产极品粉嫩在线观看| 国产成人精品在线电影| 亚洲免费av在线视频| 亚洲国产日韩欧美精品在线观看 | 国产av精品麻豆| 日本撒尿小便嘘嘘汇集6| 午夜日韩欧美国产| 精品久久久久久,| 免费女性裸体啪啪无遮挡网站| 十分钟在线观看高清视频www| 久久精品国产综合久久久| 国产野战对白在线观看| 久久久久久久久久久久大奶| 9色porny在线观看| 精品国产超薄肉色丝袜足j| 天堂动漫精品| 久99久视频精品免费| 精品高清国产在线一区| 亚洲美女黄片视频| 日本 欧美在线| 亚洲精品国产色婷婷电影| 国产亚洲精品久久久久久毛片| 亚洲视频免费观看视频| 亚洲人成电影免费在线| 老司机深夜福利视频在线观看| 国产成人系列免费观看| 麻豆av在线久日| 91大片在线观看| 久久久久久免费高清国产稀缺| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区四区五区乱码| 曰老女人黄片| 一进一出好大好爽视频| 亚洲人成电影免费在线| 亚洲国产精品合色在线| 日韩欧美三级三区| 日韩有码中文字幕| 少妇熟女aⅴ在线视频| 午夜福利欧美成人| 欧美成狂野欧美在线观看| 国产99久久九九免费精品| 91字幕亚洲| 亚洲狠狠婷婷综合久久图片| 欧美黄色淫秽网站| 国产精品98久久久久久宅男小说| 黄片小视频在线播放| 好男人在线观看高清免费视频 | 天天躁夜夜躁狠狠躁躁| 久久久久久久精品吃奶| 夜夜躁狠狠躁天天躁| 丝袜在线中文字幕| 黑丝袜美女国产一区| 中出人妻视频一区二区| 一级黄色大片毛片| 国产黄a三级三级三级人| 一进一出抽搐动态| 国产乱人伦免费视频| 亚洲,欧美精品.| 一级作爱视频免费观看| 757午夜福利合集在线观看| 黄色女人牲交| 老汉色∧v一级毛片| 美女高潮到喷水免费观看| 久久久精品国产亚洲av高清涩受| 手机成人av网站| 色综合亚洲欧美另类图片| 免费搜索国产男女视频| 久久午夜综合久久蜜桃| 亚洲精品在线观看二区| 1024视频免费在线观看| 一级毛片高清免费大全| 久久精品aⅴ一区二区三区四区| 一区二区三区精品91| 中文字幕av电影在线播放| 亚洲美女黄片视频| 成人手机av| 日日摸夜夜添夜夜添小说| 一进一出抽搐gif免费好疼| 久久久久久大精品| www.精华液| 亚洲 欧美 日韩 在线 免费| 国内毛片毛片毛片毛片毛片| 国产日韩一区二区三区精品不卡| 999久久久国产精品视频| 亚洲人成77777在线视频| 久久草成人影院| 亚洲精品一卡2卡三卡4卡5卡| 亚洲男人天堂网一区| 国产国语露脸激情在线看| 午夜福利影视在线免费观看| 午夜福利,免费看| 人人妻人人澡人人看| 一本大道久久a久久精品| 国产精品香港三级国产av潘金莲| 可以在线观看的亚洲视频| 国产成+人综合+亚洲专区| 日韩精品免费视频一区二区三区| 国产欧美日韩综合在线一区二区| 欧美不卡视频在线免费观看 | 久久亚洲精品不卡| bbb黄色大片| 无限看片的www在线观看| 在线观看免费视频网站a站| 大香蕉久久成人网| 欧美 亚洲 国产 日韩一| 在线永久观看黄色视频| 国产精品美女特级片免费视频播放器 | 国产一区在线观看成人免费| 久久香蕉激情| 欧美日韩精品网址| 一个人免费在线观看的高清视频| 欧美成人免费av一区二区三区| 国产精品亚洲一级av第二区| 国产伦人伦偷精品视频| 日韩大尺度精品在线看网址 | 伦理电影免费视频| 国产男靠女视频免费网站| 免费人成视频x8x8入口观看| 欧美一级a爱片免费观看看 | 纯流量卡能插随身wifi吗| 91av网站免费观看| 国产aⅴ精品一区二区三区波| 日韩精品青青久久久久久| 满18在线观看网站| 丝袜美足系列| 少妇被粗大的猛进出69影院| 亚洲人成电影免费在线| 一边摸一边做爽爽视频免费| 国产亚洲欧美在线一区二区| 法律面前人人平等表现在哪些方面| 亚洲精品中文字幕在线视频| 夜夜爽天天搞| 欧美乱码精品一区二区三区| 一二三四社区在线视频社区8| 女人精品久久久久毛片| 国产亚洲精品久久久久久毛片| or卡值多少钱| 成人18禁在线播放| 国产成人系列免费观看| 一区在线观看完整版| 搞女人的毛片| 一二三四社区在线视频社区8| 成年女人毛片免费观看观看9| 亚洲 欧美 日韩 在线 免费| 亚洲黑人精品在线| 午夜福利影视在线免费观看| 国产精品九九99| 人妻久久中文字幕网| 欧美亚洲日本最大视频资源| 一级,二级,三级黄色视频| 一级毛片高清免费大全| 日本精品一区二区三区蜜桃| 熟女少妇亚洲综合色aaa.| 悠悠久久av| 亚洲精品中文字幕一二三四区| 精品一区二区三区av网在线观看| 男人舔女人下体高潮全视频| 无人区码免费观看不卡| 亚洲精品国产区一区二| 免费人成视频x8x8入口观看| 国产精品亚洲美女久久久| 在线免费观看的www视频| 97超级碰碰碰精品色视频在线观看| 少妇 在线观看| 久久久久久大精品| 九色国产91popny在线| 久久久国产欧美日韩av| 操出白浆在线播放| 精品一区二区三区视频在线观看免费| 99国产精品一区二区蜜桃av| 国产伦一二天堂av在线观看| 女生性感内裤真人,穿戴方法视频| 校园春色视频在线观看| 久久人人精品亚洲av| 午夜福利视频1000在线观看 | www.www免费av| 女人爽到高潮嗷嗷叫在线视频| 亚洲,欧美精品.| 久久草成人影院| 12—13女人毛片做爰片一| 亚洲av熟女| 怎么达到女性高潮| 国产91精品成人一区二区三区| 亚洲一区高清亚洲精品| 97碰自拍视频| 一本久久中文字幕| 国产精品99久久99久久久不卡| 亚洲狠狠婷婷综合久久图片| 中文字幕av电影在线播放| 少妇被粗大的猛进出69影院| 午夜福利成人在线免费观看| 天天躁夜夜躁狠狠躁躁| 黄片大片在线免费观看| 成年版毛片免费区| 日本免费a在线| 亚洲 欧美一区二区三区| 欧美黑人精品巨大| 狂野欧美激情性xxxx| 悠悠久久av| 日本 欧美在线| 国产亚洲精品一区二区www| 国产成年人精品一区二区| 免费久久久久久久精品成人欧美视频| 成人三级黄色视频| 亚洲免费av在线视频| 黄色 视频免费看| 久久天躁狠狠躁夜夜2o2o| 女人爽到高潮嗷嗷叫在线视频| 亚洲成av片中文字幕在线观看| 亚洲avbb在线观看| 成人三级黄色视频| 操美女的视频在线观看| 999久久久国产精品视频| 十八禁人妻一区二区| 99国产精品一区二区三区| 亚洲精品国产精品久久久不卡| 亚洲精品av麻豆狂野| 久久精品aⅴ一区二区三区四区| e午夜精品久久久久久久| 亚洲精品国产区一区二| 禁无遮挡网站| 国产视频一区二区在线看| а√天堂www在线а√下载| 波多野结衣巨乳人妻| 在线免费观看的www视频| 免费久久久久久久精品成人欧美视频| 国产一区在线观看成人免费| 精品久久蜜臀av无| 亚洲国产精品sss在线观看| 中文字幕人妻熟女乱码| 999精品在线视频| 视频在线观看一区二区三区| 手机成人av网站| 中文字幕另类日韩欧美亚洲嫩草| 咕卡用的链子| 久久久精品欧美日韩精品| 高潮久久久久久久久久久不卡| 婷婷六月久久综合丁香| 校园春色视频在线观看| 精品国产一区二区三区四区第35| 91麻豆av在线| 国产真人三级小视频在线观看| 亚洲人成电影免费在线| 悠悠久久av| 桃红色精品国产亚洲av| 黄网站色视频无遮挡免费观看| av网站免费在线观看视频| 欧美乱色亚洲激情| 免费女性裸体啪啪无遮挡网站| 丝袜美足系列| 一个人免费在线观看的高清视频| 午夜日韩欧美国产| 国产精品自产拍在线观看55亚洲| 亚洲中文日韩欧美视频| 久久香蕉国产精品| 成人精品一区二区免费| 国产成人精品在线电影| 99久久99久久久精品蜜桃| 啦啦啦 在线观看视频| av天堂久久9| 久久狼人影院| www.精华液| 啦啦啦免费观看视频1| 老司机午夜十八禁免费视频| 日本五十路高清| 亚洲人成电影免费在线| 亚洲 国产 在线| 后天国语完整版免费观看| 少妇 在线观看| 亚洲国产精品成人综合色| 欧美乱色亚洲激情| 亚洲欧美一区二区三区黑人| 欧美国产精品va在线观看不卡| 国产精品自产拍在线观看55亚洲| 女人爽到高潮嗷嗷叫在线视频| 露出奶头的视频| 久久 成人 亚洲| 国产精品一区二区三区四区久久 | 一级毛片高清免费大全| 露出奶头的视频| 一区二区三区国产精品乱码| 90打野战视频偷拍视频| 国产亚洲精品综合一区在线观看 | 99热只有精品国产| 亚洲精品国产色婷婷电影| 久热这里只有精品99| 成人三级做爰电影| 在线观看一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 一进一出好大好爽视频| 国产在线精品亚洲第一网站| 国产欧美日韩精品亚洲av| 精品久久久久久久久久免费视频| 一进一出好大好爽视频| 黄片播放在线免费| 十八禁网站免费在线| 午夜免费鲁丝| 国产精品九九99| 黄色毛片三级朝国网站| 亚洲狠狠婷婷综合久久图片| 无限看片的www在线观看| 久久人妻福利社区极品人妻图片| 欧美乱色亚洲激情| 满18在线观看网站| 国产精品免费视频内射| 一区二区三区精品91| 亚洲熟女毛片儿| 一本久久中文字幕| 老司机靠b影院| 欧美成狂野欧美在线观看| 国产一级毛片七仙女欲春2 | 成人精品一区二区免费| 国产精品精品国产色婷婷| 咕卡用的链子| 狠狠狠狠99中文字幕| www日本在线高清视频| 搡老妇女老女人老熟妇| 国产乱人伦免费视频| av超薄肉色丝袜交足视频| 9热在线视频观看99| 欧美在线一区亚洲| 黄色成人免费大全| 一本久久中文字幕| 国产高清视频在线播放一区| 亚洲人成伊人成综合网2020| 亚洲国产毛片av蜜桃av| 成人亚洲精品av一区二区| avwww免费| www.精华液| 亚洲片人在线观看| 熟女少妇亚洲综合色aaa.| 欧美黑人精品巨大| 久久人妻av系列| 国产精品,欧美在线| 久久影院123| 欧美激情 高清一区二区三区| 中文字幕最新亚洲高清| 91麻豆精品激情在线观看国产| 国产精品一区二区在线不卡| АⅤ资源中文在线天堂| 精品国产国语对白av| 欧美日韩亚洲综合一区二区三区_| 咕卡用的链子| 日韩av在线大香蕉| 在线观看免费视频日本深夜| 99在线视频只有这里精品首页| 女人被躁到高潮嗷嗷叫费观| x7x7x7水蜜桃| 巨乳人妻的诱惑在线观看| www.精华液| 黄色女人牲交| 免费观看精品视频网站| 一个人免费在线观看的高清视频| 欧美性长视频在线观看| 亚洲精华国产精华精| 夜夜看夜夜爽夜夜摸| 久9热在线精品视频| 午夜精品国产一区二区电影| 日本撒尿小便嘘嘘汇集6| 成人三级黄色视频| 国产片内射在线| 亚洲av美国av| 一区二区三区国产精品乱码| 搡老妇女老女人老熟妇| 精品少妇一区二区三区视频日本电影| 一区在线观看完整版| 亚洲av成人不卡在线观看播放网| 久久久国产成人精品二区| 男女下面进入的视频免费午夜 | 91精品国产国语对白视频| 成人亚洲精品一区在线观看| 日韩国内少妇激情av| 99国产精品一区二区蜜桃av| 一区二区三区激情视频| 一进一出抽搐动态| 国产又爽黄色视频| 国产aⅴ精品一区二区三区波| 中文字幕人妻丝袜一区二区| 丝袜美足系列| 欧美在线一区亚洲| 日韩精品中文字幕看吧| 午夜亚洲福利在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 制服人妻中文乱码| 亚洲成人国产一区在线观看| 真人做人爱边吃奶动态| 亚洲精品中文字幕一二三四区| av免费在线观看网站| 青草久久国产| av中文乱码字幕在线| 老司机福利观看| 宅男免费午夜| 窝窝影院91人妻| 日韩视频一区二区在线观看| 女性被躁到高潮视频| 国产成人影院久久av| 黄色 视频免费看| 性少妇av在线| 一a级毛片在线观看| 成人18禁高潮啪啪吃奶动态图| 国产三级黄色录像| www日本在线高清视频| 日韩高清综合在线| 男人舔女人的私密视频| 欧美日韩亚洲综合一区二区三区_| 91精品国产国语对白视频| 天堂√8在线中文| 91大片在线观看| 在线永久观看黄色视频| 成人亚洲精品av一区二区| 亚洲国产看品久久| 国产人伦9x9x在线观看| 国产午夜精品久久久久久| 久久精品亚洲精品国产色婷小说| 一边摸一边抽搐一进一小说| 精品电影一区二区在线| 久久久久国产一级毛片高清牌| 黄色女人牲交| 91精品国产国语对白视频| 51午夜福利影视在线观看| 亚洲五月色婷婷综合| 国产精品美女特级片免费视频播放器 | 国产熟女午夜一区二区三区| 亚洲成人精品中文字幕电影| 最近最新中文字幕大全电影3 | 夜夜爽天天搞| 国产精品综合久久久久久久免费 | 嫩草影视91久久| 老司机福利观看| 天堂动漫精品| 老司机深夜福利视频在线观看| www.www免费av| 韩国精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区激情短视频| 欧美乱色亚洲激情| 大香蕉久久成人网| 国产野战对白在线观看| 最新美女视频免费是黄的| 亚洲av电影在线进入| 97人妻天天添夜夜摸| 久9热在线精品视频| 国产精品久久久av美女十八| 国产aⅴ精品一区二区三区波| 国产精品爽爽va在线观看网站 | 久久婷婷成人综合色麻豆| 欧美黄色淫秽网站| 免费在线观看完整版高清| 午夜日韩欧美国产| 欧美成人一区二区免费高清观看 | 欧洲精品卡2卡3卡4卡5卡区| 久久香蕉国产精品| 亚洲av五月六月丁香网| 精品熟女少妇八av免费久了| 在线观看一区二区三区| 国产精品野战在线观看| 国产成人精品久久二区二区91| 欧美亚洲日本最大视频资源| 亚洲欧美精品综合久久99| 免费在线观看完整版高清| 欧美乱码精品一区二区三区| 十八禁人妻一区二区| 女人高潮潮喷娇喘18禁视频| 制服人妻中文乱码| 19禁男女啪啪无遮挡网站| 免费一级毛片在线播放高清视频 | 不卡一级毛片| 91麻豆av在线| 亚洲精品美女久久av网站| 亚洲va日本ⅴa欧美va伊人久久| 给我免费播放毛片高清在线观看| 在线视频色国产色| 欧美久久黑人一区二区| 国产精品av久久久久免费| 亚洲欧美日韩另类电影网站| 精品一区二区三区四区五区乱码| 他把我摸到了高潮在线观看| 91在线观看av| 91国产中文字幕| 亚洲国产精品sss在线观看| 我的亚洲天堂| 亚洲精品中文字幕一二三四区| 欧美亚洲日本最大视频资源| 老司机福利观看| 欧美激情极品国产一区二区三区| 激情在线观看视频在线高清| 一进一出抽搐动态| 少妇被粗大的猛进出69影院| 欧美日本视频| 久久九九热精品免费| 国产成人系列免费观看| 欧美最黄视频在线播放免费| 国产精品久久久久久人妻精品电影| 女人被躁到高潮嗷嗷叫费观| 日韩欧美在线二视频| 国产成人精品在线电影| 中文亚洲av片在线观看爽| 首页视频小说图片口味搜索| 91精品三级在线观看| 不卡av一区二区三区| 欧美色视频一区免费| 国产精品香港三级国产av潘金莲| 成人精品一区二区免费| 99riav亚洲国产免费| 午夜福利18| 啦啦啦观看免费观看视频高清 | 亚洲精品国产区一区二| 欧美成人免费av一区二区三区| 亚洲精品粉嫩美女一区| 欧美日韩精品网址| 性少妇av在线| 亚洲久久久国产精品| 欧美国产精品va在线观看不卡| 在线观看免费午夜福利视频| 国产精品久久久久久精品电影 | 精品久久久久久,| 免费观看精品视频网站| 久久亚洲真实| 最好的美女福利视频网| 久久久久久久久免费视频了| 大型av网站在线播放| 18禁观看日本| 制服人妻中文乱码| 亚洲欧美一区二区三区黑人| 日韩免费av在线播放| 色尼玛亚洲综合影院| 99国产综合亚洲精品| 亚洲中文字幕一区二区三区有码在线看 | 国产精品综合久久久久久久免费 | 久久国产精品男人的天堂亚洲| 亚洲 欧美一区二区三区| 色精品久久人妻99蜜桃| 亚洲自拍偷在线| 亚洲精品av麻豆狂野| 这个男人来自地球电影免费观看| 久久久国产欧美日韩av| 国产激情欧美一区二区| 最近最新免费中文字幕在线| av中文乱码字幕在线| 日韩大尺度精品在线看网址 | 亚洲全国av大片| 久久久久精品国产欧美久久久| 午夜亚洲福利在线播放| 大型黄色视频在线免费观看| 国产日韩一区二区三区精品不卡| 最近最新中文字幕大全免费视频| 亚洲av成人不卡在线观看播放网| 久久久久久久久免费视频了| 一个人观看的视频www高清免费观看 | 国产一区二区三区在线臀色熟女| 在线观看www视频免费| 色av中文字幕| 亚洲人成77777在线视频| 老汉色∧v一级毛片| 在线播放国产精品三级|