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

    任意缺失模式缺失數(shù)據(jù)不同填補(bǔ)方法效果比較*

    2013-09-07 09:02:12張秋菊劉美娜
    中國衛(wèi)生統(tǒng)計(jì) 2013年5期
    關(guān)鍵詞:均數(shù)單調(diào)骨密度

    張 橋 李 寧 張秋菊 劉美娜△

    任意缺失模式缺失數(shù)據(jù)不同填補(bǔ)方法效果比較*

    張 橋1李 寧2張秋菊1劉美娜1△

    目的 探討任意缺失模式下缺失數(shù)據(jù)的填補(bǔ)方法,并對(duì)不同方法填補(bǔ)效果進(jìn)行比較和評(píng)價(jià)。方法 結(jié)合我國北方絕經(jīng)期婦女鈣需要和膳食評(píng)估應(yīng)用研究課題的數(shù)據(jù),調(diào)用SAS軟件中IML模塊產(chǎn)生任意缺失模式模擬數(shù)據(jù),通過MI和MIANALYZE過程實(shí)現(xiàn)缺失數(shù)據(jù)的填補(bǔ),同時(shí)應(yīng)用準(zhǔn)確度和穩(wěn)定度兩個(gè)評(píng)價(jià)指標(biāo)來評(píng)價(jià)各方法填補(bǔ)的效果。結(jié)果 PS方法填補(bǔ)3次在本文模擬的任意缺失模式的缺失數(shù)據(jù)中填補(bǔ)效果最佳,MCMC方法填補(bǔ)效果并不理想。結(jié)論 在填補(bǔ)任意缺失模式的缺失數(shù)據(jù)時(shí),MCMC并不是唯一的多重填補(bǔ)方法,通過多重填補(bǔ)的PS方法、PMM方法和REG方法把數(shù)據(jù)填補(bǔ)成單調(diào)缺失后,再用相同方法進(jìn)行一次填補(bǔ)也是一種可選擇的填補(bǔ)方法。

    缺失數(shù)據(jù) 任意缺失模式 多重填補(bǔ) 數(shù)據(jù)模擬

    1.哈爾濱醫(yī)科大學(xué)公共衛(wèi)生學(xué)院衛(wèi)生統(tǒng)計(jì)學(xué)教研室(150081)

    2.寧波市疾病預(yù)防控制中心免疫預(yù)防所

    △通信作者:劉美娜,E-mail:liumeina369@163.com

    數(shù)據(jù)缺失是實(shí)驗(yàn)研究和調(diào)查研究中一個(gè)普遍存在的問題〔1〕,如何正確的處理、分析所缺失的數(shù)據(jù)在數(shù)據(jù)分析中占有重要地位。缺失數(shù)據(jù)的類型按照不同的分類方法可劃分不同類別,按缺失機(jī)制分類和按缺失模式分類兩種劃分方法〔2-3〕。

    按照由Little和Rubin在1976年提出的缺失機(jī)制分類,缺失數(shù)據(jù)可以分為完全隨機(jī)缺失(missing completely at random,MCAR)、隨機(jī)缺失(missing at random,MAR)和非隨機(jī)缺失(not missing at random,NMAR)三類〔4〕。如果所缺失的數(shù)據(jù)發(fā)生的概率既與已觀察到的數(shù)據(jù)無關(guān)也與未觀察到的數(shù)據(jù)無關(guān),則該缺失數(shù)據(jù)類型為MCAR;如果缺失數(shù)據(jù)的發(fā)生概率與所觀察到的變量是有關(guān)的,而與未觀察到的數(shù)據(jù)特征無關(guān),則該缺失數(shù)據(jù)類型為MAR;若數(shù)據(jù)既不屬于完全隨機(jī)缺失也不屬于隨機(jī)缺失,那么該缺失數(shù)據(jù)類型就屬于NMAR〔5〕。按照數(shù)據(jù)缺失模式可以分為單調(diào)缺失模式和任意缺失模式兩類〔6-7〕,為了簡(jiǎn)單明了可以通過圖1來形象的理解,其中是5個(gè)變量,1~5是5個(gè)樣本,“×”表示數(shù)據(jù)能觀察到,“.”表示數(shù)據(jù)缺失。

    圖1 數(shù)據(jù)缺失模式

    單調(diào)缺失模式如圖1(a)所示,對(duì)數(shù)據(jù)集進(jìn)行適當(dāng)?shù)男辛凶儞Q后,可以得到這樣一個(gè)矩陣,它呈現(xiàn)出一種層級(jí)缺失的模式,矩陣中的元素yj缺失時(shí),則對(duì)任意的P≥j,元素yp也是缺失的;任意缺失模式如圖1(b)所示,數(shù)據(jù)缺失具有隨意性,沒有任何規(guī)律可循,即使通過行列變換也無法看出任何規(guī)律。

    對(duì)于任意缺失模式的數(shù)據(jù)處理,查閱相關(guān)文獻(xiàn)發(fā)現(xiàn)常用的就是把缺失值直接刪除即Ad Hoc法或多重填補(bǔ)(multiple imputation,MI)中的馬爾科夫鏈蒙特卡洛(markov chain monte carlo,MCMC)方法〔8〕,對(duì)于縱向數(shù)據(jù)有時(shí)也采用單一填補(bǔ)中的LOCF(last observation carried forward)方法〔9〕。本文將探討 Ad Hoc法、LOCF填補(bǔ)、多重填補(bǔ)中的回歸方法、預(yù)測(cè)均數(shù)匹配(predictive mean matching,PMM)方法、趨勢(shì)得分(propensity score,PS)方法、MCMC方法這六種方法對(duì)任意缺失模式下缺失數(shù)據(jù)的填補(bǔ)效果。

    資料與方法

    1.資料來源

    本文所用數(shù)據(jù)來源于國家科技支撐計(jì)劃項(xiàng)目:我國北方絕經(jīng)期婦女鈣需要和膳食評(píng)估應(yīng)用研究。此課題是一個(gè)為期兩年人群干預(yù)研究,研究對(duì)象282名,通過分層隨機(jī)方法分為四組,3個(gè)鈣干預(yù)組和1個(gè)信息干預(yù)組。分別在干預(yù)前、干預(yù)1年后、干預(yù)2年后三個(gè)時(shí)間點(diǎn)對(duì)干預(yù)對(duì)象進(jìn)行調(diào)查和樣品采集,獲得研究對(duì)象的體格檢查、一般情況、飲食情況、體力活動(dòng)情況和心理與應(yīng)對(duì)等信息,同時(shí)對(duì)研究對(duì)象進(jìn)行骨密度檢測(cè),所采用儀器是美國Norland XR-36雙能X線骨密度儀,包括腰椎、髖骨和全身骨三個(gè)部位,獲得相應(yīng)部位的骨密度T值。本文主要選用志愿者的身高、體重、年齡以及三次骨密度檢查的腰椎骨密度T值作為模擬實(shí)驗(yàn)的參考數(shù)據(jù)。

    2.數(shù)據(jù)基本狀況

    參考數(shù)據(jù)中身高、體重、年齡和第一次腰椎骨密度T值為完整數(shù)據(jù),共282例,第二次和第三次腰椎骨密度T值分別缺失63人和80人,因此剩余人數(shù)分別是219和202例。參考數(shù)據(jù)中各變量的均數(shù)和標(biāo)準(zhǔn)差見表1。

    表1 參考數(shù)據(jù)各變量的均數(shù)和標(biāo)準(zhǔn)

    表2是參考數(shù)據(jù)中各變量間的相關(guān)系數(shù)矩陣。

    表2 參考數(shù)據(jù)各變量的相關(guān)系數(shù)矩陣

    3.分析方法及評(píng)價(jià)標(biāo)準(zhǔn)

    本文的數(shù)據(jù)分析思路為:根據(jù)實(shí)際研究所獲數(shù)據(jù)模擬出100個(gè)完整數(shù)據(jù)集,在此基礎(chǔ)上,分別根據(jù)完整數(shù)據(jù)中第二次和第三次腰椎骨密度T值的數(shù)據(jù)缺失率(分別為22.34%和28.37%)生成100個(gè)有數(shù)據(jù)缺失的數(shù)據(jù)集,然后再用各種缺失數(shù)據(jù)填補(bǔ)方法對(duì)缺失數(shù)據(jù)集進(jìn)行填補(bǔ),最后根據(jù)評(píng)價(jià)指標(biāo)來評(píng)價(jià)各填補(bǔ)方法的優(yōu)劣。

    數(shù)據(jù)分析軟件為SAS 9.1,模擬數(shù)據(jù)集采用IML模塊和SAS宏程序,缺失數(shù)據(jù)的處理和分析主要采用了PROC MI和PROC MIANALYZE過程。由于REG方法、PMM方法和PS方法只能對(duì)單調(diào)缺失模式的數(shù)據(jù)進(jìn)行填補(bǔ),所以在用如上三種方法進(jìn)行缺失數(shù)據(jù)填補(bǔ)時(shí),本文首先對(duì)第二次腰椎骨密度T值填補(bǔ)N(N=3、5、10、15、20)次,使數(shù)據(jù)變成單調(diào)缺失后,再用相應(yīng)的填補(bǔ)方法對(duì)第三次腰椎骨密度T值填補(bǔ)1次。

    針對(duì)缺失數(shù)據(jù)填補(bǔ)效果優(yōu)劣的評(píng)價(jià)指標(biāo)本文采用準(zhǔn)確度和穩(wěn)定度〔10〕。對(duì)于變量Y,100個(gè)完整數(shù)據(jù)集有100個(gè)均數(shù)Y1,Y2,…,Y100,這 100 個(gè)均數(shù)的平均值為Ymean,缺失數(shù)據(jù)經(jīng)過處理后也會(huì)有100個(gè)均數(shù),…,均數(shù)的平均值為,則準(zhǔn)確度指標(biāo)定義為:BIASmeanmean

    BIASmean指標(biāo)的絕對(duì)值越小說明估計(jì)均數(shù)時(shí)偏差越小,準(zhǔn)確度越高。

    MSEmean指標(biāo)越大說明估計(jì)均數(shù)時(shí)穩(wěn)定度越好〔11〕。

    同理可以計(jì)算100個(gè)標(biāo)準(zhǔn)誤的BIASstderr和MSEstderr。

    結(jié) 果

    用不同填補(bǔ)方法對(duì)模擬的缺失數(shù)據(jù)集進(jìn)行填補(bǔ),第二次和第三次腰椎骨密度T值填補(bǔ)效果較好的前五位的評(píng)價(jià)指標(biāo)結(jié)果分別如表3和表4所示:

    表3 不同填補(bǔ)方法對(duì)第二次腰椎骨密度T值填補(bǔ)效果

    從表3中可以看出對(duì)于第二次腰椎骨密度T值均數(shù)準(zhǔn)確性的評(píng)價(jià)指標(biāo)BIASmean絕對(duì)值最小的前五位為:0.0005、0.0006、0.0009、0.0009、0.0012、0.0012 分別為Ad Hoc方法、MCMC方法填補(bǔ)10次、MCMC方法填補(bǔ)5次、PS方法填補(bǔ)3次、MCMC方法填補(bǔ)3次、PMM方法填補(bǔ)15次。對(duì)于第二次腰椎骨密度T值均數(shù)穩(wěn)定性的評(píng)價(jià)指標(biāo)MSEmean最大的前五位為:2.4732、1.9634、1.9466、1.9107、1.9023 分別為 Ad Hoc方法、PS方法填補(bǔ)15次、PS方法填補(bǔ)10次、PS方法填補(bǔ)20次、PS方法填補(bǔ)3次。

    對(duì)于第二次腰椎骨密度T值標(biāo)準(zhǔn)誤準(zhǔn)確性的評(píng)價(jià)指標(biāo)BIASstderr絕對(duì)值最小的前五位為:0.0001、0.0007、0.0008、0.0009、0.0012 分別為 LOCF 方法、PMM方法填補(bǔ)20次、PMM方法填補(bǔ)15次、PMM方法填補(bǔ)10次、PMM方法填補(bǔ)3次。對(duì)于第二次腰椎骨密度T值標(biāo)準(zhǔn)誤穩(wěn)定性的評(píng)價(jià)指標(biāo)MSEstderr最大的前五位為:0.0362、0.0280、0.0136、0.0118、0.0112 分別為PS方法填補(bǔ)3次、PS方法填補(bǔ)5次、PS方法填補(bǔ)10次、Ad Hoc方法、PS方法填補(bǔ)15次。

    表4 不同填補(bǔ)方法對(duì)第三次腰椎骨密度T值填補(bǔ)效果

    從表4中可以看出對(duì)于第三次腰椎骨密度T值均數(shù)準(zhǔn)確性的評(píng)價(jià)指標(biāo)BIASmean絕對(duì)值最小的前五位為:0.0014、0.0020、0.0030、0.0032、0.0039 分別為 PS方法填補(bǔ)3次、REG方法填補(bǔ)3次、PS方法填補(bǔ)15次、PS方法填補(bǔ)5次、REG方法填補(bǔ)5次。對(duì)于第三次腰椎骨密度T值均數(shù)穩(wěn)定性的評(píng)價(jià)指標(biāo)MSEmean最大的 前 五 位 為:2.0351、1.8190、1.8099、1.7976、1.7753分別為Ad Hoc方法、PS方法填補(bǔ)5次、PS方法填補(bǔ)10次、PS方法填補(bǔ)3次、PS方法填補(bǔ)15次。

    從如上的結(jié)果綜合來看,PS方法填補(bǔ)3次在本文模擬的數(shù)據(jù)中填補(bǔ)效果最佳,而MCMC方法除在第二次腰椎骨密度T值的BIASmean指標(biāo)上表現(xiàn)較好外,在其他指標(biāo)中都沒有進(jìn)入填補(bǔ)效果最好的前五位。

    討 論

    在多重填補(bǔ)的四種方法里,PS方法在第二次和第三次腰椎骨密度T值的MSEmean指標(biāo)和MSEstderr指標(biāo)上都有很好的效果,REG方法在第三次腰椎骨密度T值的BIASstderr指標(biāo)上有很好的效果,PMM方法在第二次腰椎骨密度T值的BIASstderr指標(biāo)上有很好的效果,而MCMC方法只在第二次腰椎骨密度T值的BIASmean指標(biāo)上有較好的效果。填補(bǔ)次數(shù)越多填補(bǔ)效果不一定越好。

    因此從本文可以看出,對(duì)于任意缺失模式的缺失數(shù)據(jù)集,多重填補(bǔ)的MCMC并不是唯一的多重填補(bǔ)方法,采用單調(diào)缺失模式下的多重填補(bǔ)方法把任意缺失數(shù)據(jù)填補(bǔ)成單調(diào)缺失,在此基礎(chǔ)上再進(jìn)行一次該方法的填補(bǔ),在某些條件下比MCMC填補(bǔ)的效果好。對(duì)于填補(bǔ)的次數(shù)并不是越多越好,而是要根據(jù)實(shí)際情況,進(jìn)行數(shù)據(jù)模擬,從而找出最佳的填補(bǔ)次數(shù)。

    1.Amold AM,Kronmal RA.Multiple imputation of baseline data in the cardiovascular health study.American Journal of Epidemiology,2003,157(1):74-84.

    2.Abraham,Todd W,Russell,et al.Missing data:a review of current methods and applications in epidemiological research.Current Opinion in Psychiatry,2004,17(4):315-321.

    3.James M,Robins,Wang N.Inference for imputation estimators.Biometrika,2000,87(1):113-124.

    4.Little RJ,Rubin DB.Statistical Analysis with Missing Data.New York:John Wiley&Sons,1987.

    5.Little RJ,Rubin DB.Statistical Analysis with Missing Data.2nd ed.Hoboken,NJ:John Wiley&Sons,2002.

    6.曹陽,謝萬軍,張羅漫.多重填補(bǔ)的方法及其統(tǒng)計(jì)推斷原理.中國醫(yī)院統(tǒng)計(jì),2003,10(2):77-81.

    7.李新華,夏結(jié)來.多重填補(bǔ)處理有缺失數(shù)據(jù)的2×2交叉設(shè)計(jì)資料的應(yīng)用.2004中國衛(wèi)生統(tǒng)計(jì)學(xué)術(shù)會(huì)議論文集,2004:181-187.

    8.張熙,林燧恒.多重填補(bǔ)在隨機(jī)干預(yù)實(shí)驗(yàn)研究中的應(yīng)用.中國衛(wèi)生統(tǒng)計(jì),2011,28(5):537-539.

    9.茅群霞.缺失值處理統(tǒng)計(jì)方法的模擬比較研究及應(yīng)用.四川大學(xué)碩士畢業(yè)論文,2005.

    10.Collins LM,Schafer JL,Kam CM.A comparison of inclusive and restrictive strategies in modern missing data procedures.Psychol Methods,2001,6(4):330-351.

    11.李寧.鈣干預(yù)試驗(yàn)骨密度缺失值的填補(bǔ)研究.哈爾濱醫(yī)科大學(xué)碩士畢業(yè)論文,2010.

    A Simulated Comparison between Different Imputation Meth-ods in Arbitrary Missing Data

    Zhang Qiao,Li Ning,Zhang Qiuju,et al.Department of Health Statistics,Harbin Medical University(150086),Harbin

    ObjectiveTo evaluate the imputation effect of different imputation methods in arbitrary missing data.MethodsFirst of all,we use the IML model in SAS software to simulate arbitrary missing data,which is about the calcium requirements and dietary evaluation of postmenopausal women in the north of china.Imputing the missing data through the MI and MIANALYZE processes.Accuracy and stability were used for the evaluation indices to compare the imputation effect of different methods.ResultsThe effect of PS method when imputing 3 times is the best in this data,while the effect of MCMC method is not ideal.Conclusion The MCMC is not the unique multiple imputation method when imput arbitrary missing data.The PS,PMM,REG methods could turn the arbitrary missingness pattern into monotone missingness pattern,then we use the same method to imput once again.It is also an alternative imputation method.

    Missing data;Arbitrary missingness pattern;Multiple imputation;Data simulation

    國家科技支撐計(jì)劃(2011BAI09B02)

    (責(zé)任編輯:郭海強(qiáng))

    猜你喜歡
    均數(shù)單調(diào)骨密度
    預(yù)防骨質(zhì)疏松,運(yùn)動(dòng)提高骨密度
    中老年保健(2022年3期)2022-08-24 03:00:58
    天天喝牛奶,為什么骨密度還偏低
    數(shù)列的單調(diào)性
    數(shù)列的單調(diào)性
    對(duì)數(shù)函數(shù)單調(diào)性的應(yīng)用知多少
    不要輕易給兒童做骨密度檢查
    關(guān)于均數(shù)與偏差
    OSTEOSPACE型超聲骨密度儀故障案例解析
    關(guān)于均數(shù)與偏差
    關(guān)于均數(shù)與偏差
    日本a在线网址| 久久久精品大字幕| 九九久久精品国产亚洲av麻豆 | 真实男女啪啪啪动态图| 99国产精品一区二区蜜桃av| 精品不卡国产一区二区三区| 99热这里只有是精品50| 99久久无色码亚洲精品果冻| 亚洲熟妇中文字幕五十中出| 中文资源天堂在线| 在线看三级毛片| 最近最新中文字幕大全免费视频| 啦啦啦免费观看视频1| 国内精品一区二区在线观看| 国内久久婷婷六月综合欲色啪| 欧美乱妇无乱码| 久久久精品欧美日韩精品| 久久人人精品亚洲av| 搞女人的毛片| xxxwww97欧美| 中亚洲国语对白在线视频| 床上黄色一级片| 日韩欧美国产一区二区入口| 看黄色毛片网站| 国产v大片淫在线免费观看| 男插女下体视频免费在线播放| 亚洲欧美日韩无卡精品| 黑人巨大精品欧美一区二区mp4| 日日摸夜夜添夜夜添小说| 国产激情久久老熟女| 别揉我奶头~嗯~啊~动态视频| 91在线精品国自产拍蜜月 | 精品无人区乱码1区二区| 啪啪无遮挡十八禁网站| 亚洲国产欧美人成| 男女做爰动态图高潮gif福利片| 午夜福利在线观看免费完整高清在 | 三级毛片av免费| 免费av毛片视频| 国产一区二区三区视频了| 91九色精品人成在线观看| 熟女人妻精品中文字幕| 黄频高清免费视频| 99精品欧美一区二区三区四区| 女人高潮潮喷娇喘18禁视频| 日韩欧美在线乱码| 精品一区二区三区av网在线观看| 亚洲成人久久爱视频| 日韩三级视频一区二区三区| 国产三级中文精品| 老熟妇乱子伦视频在线观看| 丁香六月欧美| 人妻久久中文字幕网| 午夜亚洲福利在线播放| 久久亚洲精品不卡| 每晚都被弄得嗷嗷叫到高潮| 亚洲中文字幕一区二区三区有码在线看 | 久久精品国产99精品国产亚洲性色| 两性夫妻黄色片| 亚洲成av人片免费观看| 日本a在线网址| 亚洲av熟女| 成人午夜高清在线视频| 淫秽高清视频在线观看| 在线免费观看的www视频| 久久这里只有精品中国| 婷婷亚洲欧美| 午夜日韩欧美国产| 久久午夜综合久久蜜桃| 久久久久国产一级毛片高清牌| 免费在线观看成人毛片| 国模一区二区三区四区视频 | 黑人巨大精品欧美一区二区mp4| 午夜精品一区二区三区免费看| 在线观看日韩欧美| АⅤ资源中文在线天堂| 国产精品国产高清国产av| 色视频www国产| 哪里可以看免费的av片| 久久国产精品影院| 俺也久久电影网| 亚洲18禁久久av| 99精品欧美一区二区三区四区| 日本免费一区二区三区高清不卡| 香蕉久久夜色| 久久午夜亚洲精品久久| 国产一区二区在线av高清观看| 亚洲片人在线观看| 黄色 视频免费看| av福利片在线观看| 国产探花在线观看一区二区| 很黄的视频免费| 亚洲人成电影免费在线| 舔av片在线| 亚洲精品粉嫩美女一区| 琪琪午夜伦伦电影理论片6080| 69av精品久久久久久| 国产淫片久久久久久久久 | x7x7x7水蜜桃| 啦啦啦观看免费观看视频高清| 岛国在线观看网站| 精品电影一区二区在线| 国内揄拍国产精品人妻在线| 免费高清视频大片| 天天躁日日操中文字幕| 性欧美人与动物交配| 男女那种视频在线观看| 不卡av一区二区三区| 亚洲熟女毛片儿| 午夜福利在线在线| 男女床上黄色一级片免费看| 母亲3免费完整高清在线观看| 欧美色欧美亚洲另类二区| xxx96com| 国产亚洲精品一区二区www| 日日摸夜夜添夜夜添小说| 老司机在亚洲福利影院| av欧美777| 色综合站精品国产| 毛片女人毛片| 免费看日本二区| 国内精品久久久久精免费| 淫妇啪啪啪对白视频| 两性午夜刺激爽爽歪歪视频在线观看| 伊人久久大香线蕉亚洲五| 国产91精品成人一区二区三区| 国产精品精品国产色婷婷| 麻豆国产97在线/欧美| 一a级毛片在线观看| 99在线人妻在线中文字幕| 性色av乱码一区二区三区2| 国产欧美日韩一区二区精品| 中文字幕av在线有码专区| 2021天堂中文幕一二区在线观| 少妇裸体淫交视频免费看高清| www.精华液| 欧美三级亚洲精品| 久久久久久国产a免费观看| 国产淫片久久久久久久久 | 成人av在线播放网站| 国产单亲对白刺激| 一区二区三区国产精品乱码| 不卡一级毛片| 黄色女人牲交| 日本 欧美在线| 99久久综合精品五月天人人| 无限看片的www在线观看| 特级一级黄色大片| 久久国产精品影院| 1024香蕉在线观看| 久久久久国内视频| 成人性生交大片免费视频hd| 精品国产乱子伦一区二区三区| 欧美激情久久久久久爽电影| 亚洲精品在线观看二区| 综合色av麻豆| 别揉我奶头~嗯~啊~动态视频| 最新美女视频免费是黄的| 一本综合久久免费| 国产三级中文精品| 国产成年人精品一区二区| 露出奶头的视频| 欧美成狂野欧美在线观看| 国产久久久一区二区三区| 中文资源天堂在线| www日本在线高清视频| 麻豆一二三区av精品| 999久久久精品免费观看国产| 草草在线视频免费看| 久久久精品欧美日韩精品| 一个人免费在线观看电影 | 看免费av毛片| 三级国产精品欧美在线观看 | 俺也久久电影网| 久久中文字幕人妻熟女| 国产av在哪里看| 婷婷精品国产亚洲av| 日日摸夜夜添夜夜添小说| 国产欧美日韩精品一区二区| 啦啦啦韩国在线观看视频| 成熟少妇高潮喷水视频| 国产高清视频在线播放一区| 久久久久亚洲av毛片大全| svipshipincom国产片| 国产三级黄色录像| 精品乱码久久久久久99久播| 亚洲中文av在线| 亚洲欧美日韩无卡精品| 精品国产亚洲在线| 亚洲中文字幕日韩| 亚洲专区国产一区二区| 国产亚洲av高清不卡| 亚洲人成网站在线播放欧美日韩| 熟女电影av网| 不卡一级毛片| 日韩欧美国产一区二区入口| 99视频精品全部免费 在线 | 大型黄色视频在线免费观看| 成年女人永久免费观看视频| 99热6这里只有精品| 精品久久久久久久久久免费视频| 老熟妇乱子伦视频在线观看| 久久精品国产清高在天天线| 色播亚洲综合网| 狂野欧美激情性xxxx| 亚洲欧美日韩高清在线视频| 欧美三级亚洲精品| 亚洲精品美女久久久久99蜜臀| 男插女下体视频免费在线播放| 老司机福利观看| 搡老熟女国产l中国老女人| 久久久久亚洲av毛片大全| 亚洲av日韩精品久久久久久密| 在线观看日韩欧美| 亚洲av成人一区二区三| 人人妻人人澡欧美一区二区| 国产一区二区激情短视频| 99久久国产精品久久久| 成年免费大片在线观看| 亚洲一区高清亚洲精品| 亚洲午夜理论影院| 精华霜和精华液先用哪个| 偷拍熟女少妇极品色| 国产成人aa在线观看| 国产精品一区二区免费欧美| 国产午夜福利久久久久久| 美女高潮喷水抽搐中文字幕| 两性夫妻黄色片| 亚洲av成人精品一区久久| 亚洲成人中文字幕在线播放| 香蕉丝袜av| 99久久精品国产亚洲精品| 亚洲熟女毛片儿| 一个人免费在线观看电影 | 亚洲国产欧美人成| 亚洲av电影不卡..在线观看| 日本五十路高清| 十八禁人妻一区二区| 欧美日本视频| 亚洲熟妇熟女久久| 婷婷丁香在线五月| 大型黄色视频在线免费观看| 国模一区二区三区四区视频 | 这个男人来自地球电影免费观看| 免费在线观看视频国产中文字幕亚洲| 国产亚洲欧美在线一区二区| 国产久久久一区二区三区| 欧美+亚洲+日韩+国产| 中亚洲国语对白在线视频| 国产欧美日韩精品一区二区| 欧美黑人欧美精品刺激| 精品久久蜜臀av无| 人妻夜夜爽99麻豆av| 欧美黑人巨大hd| 久久久久免费精品人妻一区二区| 床上黄色一级片| 高清在线国产一区| 久久精品夜夜夜夜夜久久蜜豆| 国产精品美女特级片免费视频播放器 | 国产精品一及| 久久久国产成人精品二区| 99久久久亚洲精品蜜臀av| 69av精品久久久久久| 亚洲av成人一区二区三| 琪琪午夜伦伦电影理论片6080| 亚洲欧洲精品一区二区精品久久久| 国产男靠女视频免费网站| 国产亚洲av嫩草精品影院| 天天一区二区日本电影三级| 日韩欧美免费精品| 久久精品国产综合久久久| 免费高清视频大片| 亚洲在线观看片| 99久久精品热视频| 欧美一级毛片孕妇| 男人舔女人的私密视频| tocl精华| 麻豆成人av在线观看| 国产精品自产拍在线观看55亚洲| 国产激情欧美一区二区| 香蕉av资源在线| 欧美中文日本在线观看视频| 欧美+亚洲+日韩+国产| 成年人黄色毛片网站| 激情在线观看视频在线高清| or卡值多少钱| 国产激情偷乱视频一区二区| 国产精品一区二区三区四区久久| 久久午夜综合久久蜜桃| 亚洲激情在线av| 99久国产av精品| 大型黄色视频在线免费观看| 亚洲成人久久爱视频| 久久精品91蜜桃| 亚洲精华国产精华精| 亚洲最大成人中文| 91在线精品国自产拍蜜月 | 亚洲精品在线观看二区| av女优亚洲男人天堂 | 午夜激情欧美在线| 亚洲中文字幕一区二区三区有码在线看 | 国产亚洲精品综合一区在线观看| 国产午夜精品久久久久久| 级片在线观看| 久久久成人免费电影| 成年女人毛片免费观看观看9| 国产乱人视频| 无人区码免费观看不卡| ponron亚洲| 男人的好看免费观看在线视频| 欧美一级a爱片免费观看看| 国产精品一区二区三区四区免费观看 | 亚洲avbb在线观看| 9191精品国产免费久久| 神马国产精品三级电影在线观看| 国产精品爽爽va在线观看网站| 热99在线观看视频| 中文字幕熟女人妻在线| 日韩欧美在线二视频| 欧美三级亚洲精品| 久久久久久大精品| 成年免费大片在线观看| 不卡av一区二区三区| 18禁国产床啪视频网站| 十八禁人妻一区二区| 校园春色视频在线观看| 精品国产超薄肉色丝袜足j| 国产精品 国内视频| 色老头精品视频在线观看| av在线天堂中文字幕| 国产精品野战在线观看| 99在线人妻在线中文字幕| 中文在线观看免费www的网站| 精华霜和精华液先用哪个| aaaaa片日本免费| 不卡一级毛片| 成人av在线播放网站| 三级国产精品欧美在线观看 | 12—13女人毛片做爰片一| 黄片大片在线免费观看| 欧美日韩综合久久久久久 | 国产欧美日韩精品一区二区| 哪里可以看免费的av片| 国内精品久久久久精免费| 97超视频在线观看视频| 国内久久婷婷六月综合欲色啪| 精品国产三级普通话版| 国产久久久一区二区三区| 99re在线观看精品视频| 久久久国产成人精品二区| 国产野战对白在线观看| 亚洲一区高清亚洲精品| av国产免费在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲国产看品久久| 国产久久久一区二区三区| 日本 av在线| 欧美中文综合在线视频| 精品国产三级普通话版| 天堂影院成人在线观看| 99精品欧美一区二区三区四区| 日本一本二区三区精品| 久久天堂一区二区三区四区| 久久午夜综合久久蜜桃| 欧美黄色片欧美黄色片| 脱女人内裤的视频| 国产一区二区在线av高清观看| 欧美一级a爱片免费观看看| 中文字幕人妻丝袜一区二区| 美女被艹到高潮喷水动态| 国产一区二区在线av高清观看| 丁香六月欧美| 欧美性猛交╳xxx乱大交人| 看黄色毛片网站| a级毛片a级免费在线| 国产1区2区3区精品| 日韩三级视频一区二区三区| 人人妻人人看人人澡| 最近视频中文字幕2019在线8| 免费一级毛片在线播放高清视频| 国产精品日韩av在线免费观看| 久久香蕉国产精品| 看黄色毛片网站| 最近视频中文字幕2019在线8| 在线观看舔阴道视频| 亚洲av电影在线进入| 看黄色毛片网站| av在线天堂中文字幕| 国产免费男女视频| 99国产精品一区二区蜜桃av| 黑人巨大精品欧美一区二区mp4| 18禁裸乳无遮挡免费网站照片| 全区人妻精品视频| tocl精华| 欧美又色又爽又黄视频| 成人鲁丝片一二三区免费| 美女大奶头视频| 成年女人毛片免费观看观看9| 国产真人三级小视频在线观看| 特大巨黑吊av在线直播| 精品国产美女av久久久久小说| 国产av一区在线观看免费| 丰满人妻一区二区三区视频av | 国产在线精品亚洲第一网站| 国产高清视频在线观看网站| 午夜日韩欧美国产| 51午夜福利影视在线观看| 国产精品99久久久久久久久| 国产成人精品久久二区二区91| 婷婷精品国产亚洲av| a在线观看视频网站| 国产一级毛片七仙女欲春2| 午夜福利视频1000在线观看| 国产高潮美女av| 精品一区二区三区视频在线 | 国产午夜精品久久久久久| 中文字幕av在线有码专区| 成年女人永久免费观看视频| 特级一级黄色大片| 午夜福利在线观看吧| 婷婷精品国产亚洲av| 嫩草影院入口| 久久香蕉精品热| 高清毛片免费观看视频网站| 国产爱豆传媒在线观看| 国产精品亚洲av一区麻豆| 少妇丰满av| 国产人伦9x9x在线观看| 亚洲精品粉嫩美女一区| 91九色精品人成在线观看| 日日干狠狠操夜夜爽| 亚洲专区国产一区二区| 亚洲国产色片| 久久国产乱子伦精品免费另类| 婷婷丁香在线五月| 国产精品,欧美在线| 亚洲 欧美一区二区三区| 亚洲狠狠婷婷综合久久图片| 超碰成人久久| 午夜福利在线观看吧| 日本与韩国留学比较| 美女高潮的动态| av欧美777| 久久国产精品人妻蜜桃| 国产精品九九99| 中文资源天堂在线| 欧美黑人巨大hd| 亚洲电影在线观看av| 一进一出抽搐动态| 精品人妻1区二区| 亚洲精品456在线播放app | 成人欧美大片| 免费看美女性在线毛片视频| 欧美日韩精品网址| 欧美成狂野欧美在线观看| 波多野结衣高清无吗| 91麻豆av在线| 亚洲av美国av| 欧美乱码精品一区二区三区| 亚洲国产精品成人综合色| 女同久久另类99精品国产91| 亚洲人与动物交配视频| 成人三级黄色视频| 国产高清videossex| 九色成人免费人妻av| 一本综合久久免费| 丝袜人妻中文字幕| 日韩精品青青久久久久久| 又爽又黄无遮挡网站| 日韩欧美一区二区三区在线观看| 看免费av毛片| 制服丝袜大香蕉在线| h日本视频在线播放| 久久性视频一级片| 人妻夜夜爽99麻豆av| 精品一区二区三区视频在线 | 亚洲专区中文字幕在线| 久久九九热精品免费| 一二三四在线观看免费中文在| 国产99白浆流出| 熟女电影av网| 在线观看午夜福利视频| 不卡av一区二区三区| 亚洲无线在线观看| 久久国产精品人妻蜜桃| 婷婷六月久久综合丁香| 一进一出抽搐gif免费好疼| 岛国在线免费视频观看| 超碰成人久久| 香蕉丝袜av| 国产午夜精品论理片| 亚洲性夜色夜夜综合| 欧美高清成人免费视频www| 成年女人永久免费观看视频| 免费在线观看日本一区| 美女高潮喷水抽搐中文字幕| 欧美大码av| 97超视频在线观看视频| 高潮久久久久久久久久久不卡| av在线天堂中文字幕| 三级男女做爰猛烈吃奶摸视频| 操出白浆在线播放| 最新中文字幕久久久久 | 最新中文字幕久久久久 | 18禁国产床啪视频网站| 久久精品国产综合久久久| 嫩草影院入口| 97碰自拍视频| 国产精品久久电影中文字幕| 女生性感内裤真人,穿戴方法视频| 欧美精品啪啪一区二区三区| 麻豆国产av国片精品| 亚洲美女黄片视频| 级片在线观看| 听说在线观看完整版免费高清| 亚洲av五月六月丁香网| 一进一出抽搐动态| 久久亚洲精品不卡| 国产激情欧美一区二区| 国产精品久久久久久人妻精品电影| 亚洲中文字幕一区二区三区有码在线看 | 精品国产亚洲在线| 欧美3d第一页| 国产91精品成人一区二区三区| 精品久久久久久久人妻蜜臀av| 日日干狠狠操夜夜爽| 老汉色∧v一级毛片| 99久国产av精品| 一二三四社区在线视频社区8| 不卡av一区二区三区| 久久久色成人| 亚洲国产精品999在线| 久久国产精品人妻蜜桃| 欧美日韩乱码在线| 亚洲一区二区三区色噜噜| 欧美黄色淫秽网站| 美女黄网站色视频| 国产成+人综合+亚洲专区| 国产亚洲av嫩草精品影院| 亚洲欧美日韩无卡精品| 久久国产精品人妻蜜桃| 观看美女的网站| 亚洲一区二区三区色噜噜| 欧美一级毛片孕妇| 亚洲精品粉嫩美女一区| 网址你懂的国产日韩在线| 亚洲人成电影免费在线| 国产伦精品一区二区三区视频9 | 老熟妇乱子伦视频在线观看| bbb黄色大片| 在线永久观看黄色视频| 欧美乱码精品一区二区三区| 男插女下体视频免费在线播放| 无遮挡黄片免费观看| xxx96com| 岛国在线免费视频观看| 日韩人妻高清精品专区| 99久久成人亚洲精品观看| a级毛片a级免费在线| 国产精品亚洲一级av第二区| 欧美色欧美亚洲另类二区| 99热只有精品国产| 国产精品99久久久久久久久| 最近在线观看免费完整版| av欧美777| 日韩欧美在线乱码| 日本一二三区视频观看| 巨乳人妻的诱惑在线观看| 国产精品野战在线观看| 在线观看午夜福利视频| 亚洲自拍偷在线| 99热这里只有是精品50| 亚洲欧美日韩卡通动漫| 男人舔奶头视频| 国产 一区 欧美 日韩| 精品免费久久久久久久清纯| ponron亚洲| 99久久成人亚洲精品观看| 久久中文字幕一级| 人妻夜夜爽99麻豆av| www.精华液| 黄色视频,在线免费观看| 午夜福利在线观看免费完整高清在 | 国产成人精品无人区| 欧美+亚洲+日韩+国产| 国产 一区 欧美 日韩| 亚洲欧美一区二区三区黑人| 亚洲专区中文字幕在线| 制服人妻中文乱码| 制服丝袜大香蕉在线| 一进一出抽搐gif免费好疼| www.熟女人妻精品国产| 亚洲中文av在线| 欧美日韩瑟瑟在线播放| 香蕉国产在线看| 悠悠久久av| 亚洲第一欧美日韩一区二区三区| 国产三级黄色录像| 熟女电影av网| 色av中文字幕| 九九热线精品视视频播放| 一区二区三区激情视频| 一卡2卡三卡四卡精品乱码亚洲| 成年女人永久免费观看视频| 全区人妻精品视频| e午夜精品久久久久久久| АⅤ资源中文在线天堂| 三级男女做爰猛烈吃奶摸视频| 久久精品国产99精品国产亚洲性色| 亚洲中文av在线| 成人特级黄色片久久久久久久| 久久草成人影院| 国产三级黄色录像| 亚洲国产高清在线一区二区三| 亚洲avbb在线观看|