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

    具有相關(guān)關(guān)系的靈敏度和特異度的Monte Carlo模擬方法研究*

    2015-01-27 10:31:04李劍龍陳方堯李丹玲周基元陳平雁
    中國衛(wèi)生統(tǒng)計(jì) 2015年3期
    關(guān)鍵詞:四格錯(cuò)誤率樣本量

    李劍龍 陳方堯 李丹玲 周基元 陳平雁

    南方醫(yī)科大學(xué)公共衛(wèi)生與熱帶醫(yī)學(xué)學(xué)院“器官衰竭防治”國家重點(diǎn)實(shí)驗(yàn)室(510515)

    具有相關(guān)關(guān)系的靈敏度和特異度的Monte Carlo模擬方法研究*

    李劍龍#陳方堯#李丹玲 周基元△陳平雁

    南方醫(yī)科大學(xué)公共衛(wèi)生與熱帶醫(yī)學(xué)學(xué)院“器官衰竭防治”國家重點(diǎn)實(shí)驗(yàn)室(510515)

    目的 旨在提出一種Monte Carlo模擬方法,對(duì)具有相關(guān)關(guān)系的靈敏度和特異度進(jìn)行計(jì)算機(jī)模擬,并運(yùn)用其相關(guān)系數(shù)估計(jì)值的均方根誤及偏差對(duì)該模擬方法的有效性進(jìn)行評(píng)價(jià);進(jìn)而探討分析具有相關(guān)關(guān)系的靈敏度和特異度數(shù)據(jù)時(shí)單樣本Youden指數(shù)u檢驗(yàn)方法的可靠性。方法 將診斷試驗(yàn)的四格表資料看成來自服從四元正態(tài)分布總體的樣本,借助R語言的mvrnorm函數(shù)產(chǎn)生四格表資料,其滿足靈敏度和特異度之間的總體相關(guān)系數(shù)。計(jì)算相關(guān)系數(shù)估計(jì)值的均方根誤及偏差,并計(jì)算Youden指數(shù)u檢驗(yàn)的第I類錯(cuò)誤率與檢驗(yàn)效能。結(jié)果 相關(guān)系數(shù)估計(jì)值的均方根誤及偏差均較小。當(dāng)相關(guān)系數(shù)ρ為0時(shí),u檢驗(yàn)?zāi)芎芎玫乜刂品傅贗類錯(cuò)誤的概率;然而,當(dāng)ρ不是0時(shí),u檢驗(yàn)的結(jié)果偏保守,并且隨著|ρ|的增加,其結(jié)果偏離0.05更甚。結(jié)論 本文提出的具有相關(guān)關(guān)系靈敏度和特異度的Monte Carlo模擬方法能真實(shí)地反應(yīng)靈敏度和特異度之間的相關(guān)性;當(dāng)相關(guān)性強(qiáng)時(shí),Youden指數(shù)u檢驗(yàn)的前提假設(shè)不正確。

    相關(guān)系數(shù) Monte Carlo模擬 靈敏度 特異度

    科學(xué)的評(píng)價(jià)診斷試驗(yàn)是臨床醫(yī)師選擇診斷試驗(yàn)的基礎(chǔ)[1]。為了合理選用診斷試驗(yàn)以避免盲目性,臨床醫(yī)師就需要對(duì)診斷試驗(yàn)進(jìn)行科學(xué)研究并做出科學(xué)的評(píng)價(jià)。在進(jìn)行兩個(gè)診斷系統(tǒng)比較時(shí),單獨(dú)使用靈敏度與特異度作為評(píng)價(jià)指標(biāo)不能完美判斷哪個(gè)診斷系統(tǒng)更好。由此,Youden建議將靈敏度與特異度綜合起來提出了Youden指數(shù)來評(píng)價(jià)診斷試驗(yàn)的準(zhǔn)確度[2]。Youden指數(shù)(Youden index,YI)是評(píng)價(jià)診斷方法診斷能力的綜合性指標(biāo)[3]。YI的值介于-1與1之間,其值越大,說明診斷試驗(yàn)的真實(shí)性越好;當(dāng)YI小于等于0時(shí),該診斷試驗(yàn)無任何臨床應(yīng)用價(jià)值[1]。

    然而,對(duì)于現(xiàn)有診斷試驗(yàn)的評(píng)價(jià)方法能否正確地對(duì)診斷試驗(yàn)進(jìn)行評(píng)價(jià)都是需要進(jìn)行檢驗(yàn)的。單樣本Youden指數(shù)u檢驗(yàn)可以用于檢驗(yàn)單樣本的Youden指數(shù)是否為某一常數(shù)。注意到其需要一個(gè)前提假設(shè),即要求靈敏度和特異度之間是相互獨(dú)立的。然而,現(xiàn)實(shí)中靈敏度和特異度通常是相關(guān)的[4-5]。有研究曾提出采用多項(xiàng)分布來模擬產(chǎn)生兩個(gè)相關(guān)的二項(xiàng)分布數(shù)據(jù)(四格表數(shù)據(jù)),然而這種方法早前曾被研究人員證實(shí)所產(chǎn)生的數(shù)據(jù)并不能良好地反應(yīng)出所設(shè)定的參數(shù)[6-7]。因此,本研究旨在提出一種Monte Carlo模擬方法[8],對(duì)具有相關(guān)關(guān)系的靈敏度和特異度進(jìn)行計(jì)算機(jī)模擬,并運(yùn)用其相關(guān)系數(shù)估計(jì)值的均方根誤(root mean squared error,RMSE)及偏差對(duì)該模擬方法的有效性進(jìn)行評(píng)價(jià);進(jìn)而計(jì)算出其樣本Youden指數(shù),運(yùn)用現(xiàn)有的單樣本Youden指數(shù)u檢驗(yàn)方法對(duì)模擬數(shù)據(jù)進(jìn)行分析,重復(fù)多次模擬時(shí)計(jì)算其第I類錯(cuò)誤率與檢驗(yàn)效能,以探討運(yùn)用u檢驗(yàn)方法分析具有相關(guān)關(guān)系的靈敏度和特異度數(shù)據(jù)時(shí)的可靠性。

    對(duì)象與方法

    1.理論背景

    (1)靈敏度與特異度之間的相關(guān)系數(shù)

    在實(shí)際應(yīng)用過程中,靈敏度與特異度之間存在著一定的相關(guān)關(guān)系[4-5],體現(xiàn)二者相關(guān)關(guān)系的方法很多[9]。本文只考慮將二者之間的相關(guān)關(guān)系與二者之間的線性回歸聯(lián)系起來的方法,其滿足r=bSpe.Sen(SSen/SSpe)=bSen.Spe(SSpe/SSen),其中r為靈敏度與特異度之間的樣本相關(guān)系數(shù),bSpe.Sen與bSen.Spe分別表示由靈敏度到特異度的回歸系數(shù)和由特異度到靈敏度的回歸系數(shù),SSen與SSpe分別為靈敏度與特異度的標(biāo)準(zhǔn)誤。

    (2)靈敏度Sen、特異度Spe及Youden指數(shù)的u檢驗(yàn)[1]

    (1)

    給定一個(gè)四格表診斷試驗(yàn)數(shù)據(jù)與三個(gè)常數(shù)Sen0、Spe0與YI0,對(duì)靈敏度、特異度與Youden指數(shù)可以做單樣本的u檢驗(yàn),其檢驗(yàn)統(tǒng)計(jì)量分別為

    (2)

    2.模擬方法

    (1)模擬第I類錯(cuò)誤率時(shí)的參數(shù)設(shè)置

    設(shè)p1、p2、p3與p4分別為四格表中四個(gè)實(shí)際頻數(shù)a、b、c與d的理論頻率,t=p1+p3=0.5,1-t=p2+p4=0.5,給定靈敏度總體均數(shù)μ1=0.65,特異度總體均數(shù)μ2=0.75,那么p1=0.325,p2=0.175,p3=0.125,p4=0.375。總樣本量n為100、200、300、500與1000,靈敏度和特異度之間的總體相關(guān)系數(shù)ρ為0、-0.3與-0.4。從而,三個(gè)檢驗(yàn)的原假設(shè)分別為Sen0=0.65、Spe0=0.75、YI0=0.4。

    (2)模擬靈敏度和特異度之間的相關(guān)關(guān)系與四格表數(shù)據(jù)

    (3)靈敏度和特異度之間相關(guān)系數(shù)的估計(jì)

    (4)模擬檢驗(yàn)效能時(shí)的參數(shù)設(shè)置

    三個(gè)檢驗(yàn)的原假設(shè)仍然分別取為Sen0=0.65、Spe0=0.75、YI0=0.4。而在備擇假設(shè)條件下Sen、Spe和YI的設(shè)置如表1,其中D=YI-YI0。另外,總樣本量n、靈敏度和特異度之間總體相關(guān)系數(shù)ρ的取法與第I類錯(cuò)誤率模擬中的參數(shù)設(shè)置相同。

    (5)模擬第I類錯(cuò)誤率與檢驗(yàn)效能

    結(jié) 果

    1.靈敏度和特異度之間相關(guān)系數(shù)的均方根誤與偏差

    2.單樣本靈敏度、特異度與Youden指數(shù)u檢驗(yàn)的第I類錯(cuò)誤率

    在不同ρ值與不同樣本量條件下,單樣本靈敏度、特異度與Youden指數(shù)u檢驗(yàn)的第I類錯(cuò)誤率結(jié)果見表3。無論ρ為0還是不為0,靈敏度與特異度的u檢驗(yàn)第I類錯(cuò)誤率基本控制在0.05的附近。因此,本文新提出的模擬方法能準(zhǔn)確地模擬出所設(shè)定的參數(shù)。另一方面,由于單樣本Youden指數(shù)u檢驗(yàn)有一個(gè)前提假設(shè),即要求靈敏度和特異度之間是相互獨(dú)立的。因此,當(dāng)靈敏度和特異度的相關(guān)系數(shù)為0時(shí),其u檢驗(yàn)?zāi)芎芎玫乜刂品傅贗類錯(cuò)誤的概率。然而,當(dāng)相關(guān)系數(shù)ρ不是0時(shí),單樣本Youden指數(shù)u檢驗(yàn)的結(jié)果偏保守,并且隨著相關(guān)性的增強(qiáng),其結(jié)果偏離0.05更甚。

    3.單樣本Youden指數(shù)u檢驗(yàn)的檢驗(yàn)效能

    在不同ρ值與不同樣本量條件下,單樣本Youden指數(shù)u檢驗(yàn)的檢驗(yàn)效能結(jié)果見圖1。注意,當(dāng)YI與YI0之間的差值D為0時(shí),其為u檢驗(yàn)的第I類錯(cuò)誤率。從圖形可以看出,當(dāng)相關(guān)系數(shù)與樣本量固定時(shí),YI與YI0之間差值D的絕對(duì)值越大,單樣本Youden指數(shù)u檢驗(yàn)的檢驗(yàn)效能越大。另外,隨著樣本量的增加,其檢驗(yàn)效能逐漸增大。

    討 論

    本文提出了一種基于四元正態(tài)分布總體、對(duì)具有相關(guān)關(guān)系的靈敏度和特異度進(jìn)行Monte Carlo模擬的方法。從模擬結(jié)果可以看出,模擬方法能真實(shí)地反應(yīng)靈敏度和特異度之間的相關(guān)系數(shù)。運(yùn)用現(xiàn)有的單樣本Youden指數(shù)u檢驗(yàn)方法對(duì)模擬數(shù)據(jù)進(jìn)行分析,重復(fù)多次模擬時(shí)計(jì)算其第I類錯(cuò)誤率。從模擬結(jié)果可以看出,當(dāng)靈敏度和特異度的相關(guān)系數(shù)為0時(shí),u檢驗(yàn)?zāi)芎芎玫乜刂品傅贗類錯(cuò)誤的概率。然而,當(dāng)相關(guān)系數(shù)ρ不是0時(shí),單樣本Youden指數(shù)u檢驗(yàn)的結(jié)果偏保守,并且隨著相關(guān)性的增強(qiáng),其結(jié)果偏離0.05更甚。因此,當(dāng)靈敏度和特異度之間的相關(guān)性強(qiáng)時(shí),反應(yīng)診斷試驗(yàn)真實(shí)性的綜合評(píng)價(jià)指標(biāo)Youden指數(shù)單樣本u檢驗(yàn)的前提假設(shè)不正確。

    本文只模擬了總樣本量大于等于100的情形,將來要針對(duì)小樣本做適當(dāng)?shù)哪M。另外,本文中靈敏度與特異度之間的相關(guān)關(guān)系采用了二者標(biāo)準(zhǔn)化后將二者之間的回歸與相關(guān)聯(lián)系起來的方法,但是靈敏度與特異度之間的相關(guān)關(guān)系并不僅限于此種形式,在將來的研究中要選取更多體現(xiàn)相關(guān)關(guān)系的情況進(jìn)行模擬分析。

    [1]Galen RS,Brennan L.Laboratory Diagnosis and Patient Monitoring:Clinical Chemistry.Oradell,New Jersey:Medical Economics Company,1981.

    [2]Yerushalmy J.Statistical problems in assessing methods of medical diagnosis,with special reference to X-ray techniques.Public Health Reports (1896-1970),1947,62(40):1432-1449.

    [3]Youden WJ.Index for rating diagnostic tests.Cancer,1950,3(1):32-35.

    [4]Zhou XH,Obuchowski NA,McClish DK.Statistical Method in Diagnostic Medicine.New York:John Wiley & Sons,2002.

    [5]Fleiss JL,Levin B,Paik MC.Statistical Methods for Rates and Proportions.New York:John Wiley & Sons,2013.

    [6]Johnson NL,Kotz S,Balakrishnan N.Discrete Multivariate Distributions.New York:Wiley,1997.

    [7]Mardia KV.Families of Bivariate Distributions.London:Griffin,1970.

    [8]文德智,卓仁鴻,丁大杰,等.蒙特卡羅模擬中相關(guān)變量隨機(jī)數(shù)序列的產(chǎn)生方法.物理學(xué)報(bào),2012,61(22):220204-220204.

    [9]Rodgers JL,Nicewander WA.Thirteen ways to look at the correlation coefficient.The American Statistician,1988,42(1):59-66.

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

    A Monte Carlo Simulation Method for Generating the Correlated Sensitivity and Specificity

    Li Jianlong,Chen Fangyao,Li Danling,et al

    (State Key Laboratory of Organ Failure Research,School of Public Health and Tropical Medicine,Southern Medical University (510515),Guangzhou)

    Objective This study aims to propose a Monte Carlo simulation method to generate the correlated sensitivity and specificity in diagnostic trials,and then calculate the root mean squared error (RMSE) and bias of the estimate of the correlation coefficient between them so as to assess the validity of the proposed simulation method.Further,we study the reliability of the existingutest for single sampleYoudenindex based on the data where the sensitivity and specificity are correlated.Methods In the simulation study,the diagnostic test data are considered to be generated from a four-dimensional normal distribution and the four-fold table data in the diagnostic test are generated from the mvrnorm function in R,where the sample correlation coefficient between the sensitivity and specificity is kept close to the population correlation coefficientρ.As such,the corresponding RMSE and bias,the type I error rate and power can be calculated.Results Both the RMSE and bias are very small.Whenρis 0,theutest can control the type I error rate well.However,whenρis not 0,theutest is conservative.Further,with the increase of |ρ|,the type I error rate is much more faraway from 0.05.Conclusion The Monte Carlo simulation method proposed in this article can effectively simulate the correlation relationship between the sensitivity and specificity.When there is a strong correlation between the sensitivity and specificity,the independence assumption between them is not so reasonable for theutest of single sampleYoudenindex in diagnostic test.

    Correlation coefficient;Monte Carlo simulation;Sensitivity;Specificity

    國家自然科學(xué)基金項(xiàng)目(81373098,81072386,81402759);廣東省醫(yī)學(xué)科研基金(B2013222)和全國統(tǒng)計(jì)科學(xué)研究計(jì)劃項(xiàng)目(2013LY059)共同資助。

    # 李劍龍與陳方堯?yàn)楣餐谝蛔髡?/p>

    △通信作者:周基元,E-mail:zhoujiyuan5460@hotmail.com

    猜你喜歡
    四格錯(cuò)誤率樣本量
    油爆四格
    限制性隨機(jī)試驗(yàn)中選擇偏倚導(dǎo)致的一類錯(cuò)誤率膨脹*
    醫(yī)學(xué)研究中樣本量的選擇
    爆笑四格
    爆笑四格
    爆笑四格
    航空裝備測試性試驗(yàn)樣本量確定方法
    Sample Size Calculations for Comparing Groups with Binary Outcomes
    正視錯(cuò)誤,尋求策略
    教師·中(2017年3期)2017-04-20 21:49:49
    解析小學(xué)高段學(xué)生英語單詞抄寫作業(yè)錯(cuò)誤原因
    精品国产乱码久久久久久小说| 激情 狠狠 欧美| 国产精品久久久久久av不卡| 国产v大片淫在线免费观看| 久久久欧美国产精品| 制服丝袜香蕉在线| 国产精品免费大片| 国产人妻一区二区三区在| 国产欧美日韩一区二区三区在线 | 成人无遮挡网站| 建设人人有责人人尽责人人享有的 | 青春草视频在线免费观看| 成人亚洲欧美一区二区av| 亚洲一区二区三区欧美精品| 亚洲av电影在线观看一区二区三区| 国产精品久久久久久精品电影小说 | 黄片无遮挡物在线观看| 久久ye,这里只有精品| 在线 av 中文字幕| 蜜桃亚洲精品一区二区三区| 精品一品国产午夜福利视频| 国产高清不卡午夜福利| 久久久亚洲精品成人影院| 国产男人的电影天堂91| 春色校园在线视频观看| 精品一区二区三卡| 两个人的视频大全免费| 久久ye,这里只有精品| 少妇 在线观看| 秋霞伦理黄片| 国内少妇人妻偷人精品xxx网站| 亚洲自偷自拍三级| 最近2019中文字幕mv第一页| 欧美日本视频| 精华霜和精华液先用哪个| 久久综合国产亚洲精品| 在线亚洲精品国产二区图片欧美 | 精品人妻熟女av久视频| 国产精品99久久99久久久不卡 | 啦啦啦中文免费视频观看日本| 久久99热6这里只有精品| 欧美变态另类bdsm刘玥| 这个男人来自地球电影免费观看 | 日韩三级伦理在线观看| 简卡轻食公司| 91狼人影院| 久久99精品国语久久久| 国产高清有码在线观看视频| av线在线观看网站| 制服丝袜香蕉在线| 国产精品一区二区在线不卡| 亚洲精品中文字幕在线视频 | 成人美女网站在线观看视频| 成人毛片a级毛片在线播放| 久久久国产一区二区| 亚洲精品日本国产第一区| 久久久久久久亚洲中文字幕| 免费观看在线日韩| 黑人高潮一二区| 国产亚洲精品久久久com| 国产大屁股一区二区在线视频| 波野结衣二区三区在线| 亚洲av福利一区| 热99国产精品久久久久久7| 一级爰片在线观看| 亚洲精品自拍成人| 中国国产av一级| 性色av一级| 在线精品无人区一区二区三 | av在线老鸭窝| 日韩av在线免费看完整版不卡| 夫妻午夜视频| 蜜桃亚洲精品一区二区三区| 插逼视频在线观看| 男女边吃奶边做爰视频| 久热久热在线精品观看| 国产精品一区二区三区四区免费观看| 99热网站在线观看| 热99国产精品久久久久久7| 成人综合一区亚洲| 黄色配什么色好看| 内射极品少妇av片p| 国产精品久久久久久精品电影小说 | 天堂8中文在线网| 美女内射精品一级片tv| av女优亚洲男人天堂| 男女无遮挡免费网站观看| av网站免费在线观看视频| 黑人猛操日本美女一级片| 国产中年淑女户外野战色| 91狼人影院| 亚洲成人中文字幕在线播放| 成人无遮挡网站| 亚洲av成人精品一区久久| 大话2 男鬼变身卡| 美女xxoo啪啪120秒动态图| 干丝袜人妻中文字幕| 国产探花极品一区二区| 日韩中字成人| 色婷婷久久久亚洲欧美| 色5月婷婷丁香| 成人黄色视频免费在线看| 久久97久久精品| 99热这里只有精品一区| av在线蜜桃| 欧美国产精品一级二级三级 | 日韩成人伦理影院| 哪个播放器可以免费观看大片| 久久97久久精品| 亚洲图色成人| 天天躁夜夜躁狠狠久久av| 噜噜噜噜噜久久久久久91| 国产无遮挡羞羞视频在线观看| 免费大片黄手机在线观看| 亚洲国产av新网站| 欧美精品亚洲一区二区| 看免费成人av毛片| 黑人高潮一二区| 亚洲欧美一区二区三区国产| 黄色配什么色好看| 亚洲内射少妇av| 国产精品国产三级国产专区5o| 高清视频免费观看一区二区| 亚洲av.av天堂| 色婷婷久久久亚洲欧美| 伊人久久精品亚洲午夜| 免费不卡的大黄色大毛片视频在线观看| 国产老妇伦熟女老妇高清| 视频中文字幕在线观看| 在线观看免费高清a一片| 亚洲婷婷狠狠爱综合网| 亚洲精品成人av观看孕妇| 男女无遮挡免费网站观看| 中国美白少妇内射xxxbb| 亚洲真实伦在线观看| 国国产精品蜜臀av免费| 人妻少妇偷人精品九色| 久久久久久久国产电影| 男女无遮挡免费网站观看| 午夜福利高清视频| 男女国产视频网站| 国产精品一二三区在线看| 亚洲性久久影院| 深夜a级毛片| 亚洲精品久久午夜乱码| 亚洲av不卡在线观看| 国产在线免费精品| 99热国产这里只有精品6| 国产成人91sexporn| 精品久久久久久久末码| 在线免费十八禁| 国产av一区二区精品久久 | 97超视频在线观看视频| 中文乱码字字幕精品一区二区三区| 欧美97在线视频| 男女免费视频国产| 国产成人精品婷婷| 在现免费观看毛片| 久久这里有精品视频免费| 这个男人来自地球电影免费观看 | 亚洲成人一二三区av| 国产大屁股一区二区在线视频| 国产白丝娇喘喷水9色精品| www.色视频.com| 99久久综合免费| 国产 一区 欧美 日韩| 黄色日韩在线| 国产黄色免费在线视频| 天美传媒精品一区二区| 内地一区二区视频在线| 国产在线男女| 街头女战士在线观看网站| 中国国产av一级| 亚洲成人一二三区av| 免费看光身美女| 少妇被粗大猛烈的视频| 天天躁日日操中文字幕| 一个人免费看片子| 我的老师免费观看完整版| www.色视频.com| 老司机影院成人| 午夜免费男女啪啪视频观看| 亚洲精品日韩av片在线观看| 超碰av人人做人人爽久久| 免费黄色在线免费观看| 亚洲精华国产精华液的使用体验| 日本与韩国留学比较| 欧美一区二区亚洲| 特大巨黑吊av在线直播| 五月开心婷婷网| 深夜a级毛片| 夫妻午夜视频| 欧美xxxx性猛交bbbb| 18+在线观看网站| 日韩欧美 国产精品| 久久人人爽人人爽人人片va| 久久久久网色| a级一级毛片免费在线观看| 久久久久久久久久人人人人人人| 国产精品偷伦视频观看了| 日本猛色少妇xxxxx猛交久久| 欧美xxxx黑人xx丫x性爽| 欧美极品一区二区三区四区| 国产又色又爽无遮挡免| 午夜福利网站1000一区二区三区| 久久久久久久久大av| 狂野欧美激情性xxxx在线观看| 中文乱码字字幕精品一区二区三区| 日韩视频在线欧美| 欧美xxⅹ黑人| 又粗又硬又长又爽又黄的视频| 久久久久久久精品精品| 日本黄色片子视频| 在线播放无遮挡| 国产极品天堂在线| 激情 狠狠 欧美| a 毛片基地| 另类亚洲欧美激情| 男女边摸边吃奶| 成人18禁高潮啪啪吃奶动态图 | 欧美三级亚洲精品| 日本vs欧美在线观看视频 | 老司机影院成人| 亚洲精品一区蜜桃| 美女脱内裤让男人舔精品视频| 看非洲黑人一级黄片| 卡戴珊不雅视频在线播放| 中文字幕亚洲精品专区| 国产高潮美女av| av国产精品久久久久影院| 一级毛片aaaaaa免费看小| 少妇丰满av| 麻豆成人午夜福利视频| 你懂的网址亚洲精品在线观看| 国产综合精华液| 亚洲国产毛片av蜜桃av| 看非洲黑人一级黄片| 狂野欧美激情性bbbbbb| av福利片在线观看| 亚洲,一卡二卡三卡| 亚洲av中文字字幕乱码综合| 中文精品一卡2卡3卡4更新| 一级av片app| 高清黄色对白视频在线免费看 | 亚洲av在线观看美女高潮| 这个男人来自地球电影免费观看 | 啦啦啦视频在线资源免费观看| 日韩一区二区视频免费看| 搡女人真爽免费视频火全软件| 少妇人妻精品综合一区二区| 三级国产精品欧美在线观看| 丰满少妇做爰视频| 久久久久久久久久久免费av| 校园人妻丝袜中文字幕| 欧美老熟妇乱子伦牲交| 中文精品一卡2卡3卡4更新| 蜜桃久久精品国产亚洲av| 亚洲精品久久久久久婷婷小说| 亚洲综合精品二区| 看十八女毛片水多多多| 深爱激情五月婷婷| 亚洲人成网站在线观看播放| 国产亚洲91精品色在线| 夜夜爽夜夜爽视频| 亚洲精品aⅴ在线观看| 只有这里有精品99| 三级经典国产精品| 只有这里有精品99| 在线观看美女被高潮喷水网站| 久久精品夜色国产| 国产成人精品一,二区| 夫妻性生交免费视频一级片| 亚洲怡红院男人天堂| 99久久综合免费| 高清在线视频一区二区三区| 人体艺术视频欧美日本| 国产精品国产三级国产av玫瑰| 男人和女人高潮做爰伦理| 熟妇人妻不卡中文字幕| 大陆偷拍与自拍| 久久久久精品久久久久真实原创| 久久精品久久精品一区二区三区| 亚洲成人一二三区av| 国产成人精品福利久久| 亚洲国产av新网站| 国内精品宾馆在线| 久久婷婷青草| 国产极品天堂在线| 超碰av人人做人人爽久久| 超碰97精品在线观看| 少妇高潮的动态图| 亚洲怡红院男人天堂| 丰满少妇做爰视频| 一级毛片aaaaaa免费看小| 只有这里有精品99| 亚洲精品久久午夜乱码| 国产成人a区在线观看| 欧美亚洲 丝袜 人妻 在线| 18禁在线播放成人免费| 国产成人精品一,二区| 日韩制服骚丝袜av| 午夜精品国产一区二区电影| 久久国产乱子免费精品| 嘟嘟电影网在线观看| 蜜臀久久99精品久久宅男| 在线 av 中文字幕| 亚洲精品aⅴ在线观看| 婷婷色综合www| 女人久久www免费人成看片| 视频中文字幕在线观看| 不卡视频在线观看欧美| 久久青草综合色| 亚洲国产精品999| 欧美另类一区| 欧美变态另类bdsm刘玥| 伦理电影免费视频| 亚洲欧美成人精品一区二区| 搡老乐熟女国产| 成人一区二区视频在线观看| 久久人人爽人人爽人人片va| 久久婷婷青草| 免费av中文字幕在线| 大又大粗又爽又黄少妇毛片口| 中国三级夫妇交换| 国产一区有黄有色的免费视频| 精品国产三级普通话版| 内地一区二区视频在线| 日韩制服骚丝袜av| 极品少妇高潮喷水抽搐| 丰满人妻一区二区三区视频av| 亚洲图色成人| 亚洲综合色惰| 国产日韩欧美亚洲二区| 亚洲美女视频黄频| av在线app专区| 最新中文字幕久久久久| 亚洲不卡免费看| 免费播放大片免费观看视频在线观看| 久久久久久久国产电影| 亚洲怡红院男人天堂| 蜜桃久久精品国产亚洲av| 美女主播在线视频| 免费黄网站久久成人精品| 哪个播放器可以免费观看大片| 啦啦啦啦在线视频资源| 哪个播放器可以免费观看大片| 热re99久久精品国产66热6| 国产极品天堂在线| 三级国产精品片| 国产极品天堂在线| 干丝袜人妻中文字幕| 国产精品福利在线免费观看| 国产精品久久久久久精品古装| 中文欧美无线码| 三级国产精品片| 汤姆久久久久久久影院中文字幕| 精品人妻熟女av久视频| 久久久午夜欧美精品| 少妇 在线观看| 亚洲精品自拍成人| 色视频www国产| 久久99精品国语久久久| 亚洲成人一二三区av| 街头女战士在线观看网站| 性色av一级| 亚洲美女黄色视频免费看| 欧美丝袜亚洲另类| 久久热精品热| 建设人人有责人人尽责人人享有的 | av在线观看视频网站免费| 狂野欧美激情性bbbbbb| 亚洲精品,欧美精品| 99热全是精品| 中文欧美无线码| 亚洲av综合色区一区| 欧美亚洲 丝袜 人妻 在线| 青春草国产在线视频| 亚洲精品久久久久久婷婷小说| 麻豆精品久久久久久蜜桃| 久久久午夜欧美精品| 青春草国产在线视频| 中文欧美无线码| 成人二区视频| 亚洲av成人精品一二三区| 欧美极品一区二区三区四区| 亚洲精品,欧美精品| 国产成人aa在线观看| 久久女婷五月综合色啪小说| 国产国拍精品亚洲av在线观看| 少妇熟女欧美另类| 性高湖久久久久久久久免费观看| 国产成人免费观看mmmm| 国产精品国产三级专区第一集| 自拍偷自拍亚洲精品老妇| 99热6这里只有精品| 只有这里有精品99| 亚洲av不卡在线观看| 精品国产三级普通话版| 又爽又黄a免费视频| 亚洲一级一片aⅴ在线观看| 亚洲人成网站高清观看| 高清av免费在线| 国产精品麻豆人妻色哟哟久久| 久久久久视频综合| 80岁老熟妇乱子伦牲交| 青青草视频在线视频观看| 亚洲欧美日韩东京热| 美女福利国产在线 | 国产精品三级大全| 色5月婷婷丁香| 免费观看性生交大片5| 国产在线免费精品| 久久国产精品男人的天堂亚洲 | 尾随美女入室| 亚洲欧美日韩卡通动漫| 插逼视频在线观看| 精品人妻熟女av久视频| 91aial.com中文字幕在线观看| 嘟嘟电影网在线观看| 视频区图区小说| 中文字幕人妻熟人妻熟丝袜美| 亚州av有码| 国产淫语在线视频| 深爱激情五月婷婷| 成人亚洲精品一区在线观看 | 国产精品.久久久| 久久午夜福利片| 日本-黄色视频高清免费观看| 国产精品久久久久久精品电影小说 | 好男人视频免费观看在线| 人体艺术视频欧美日本| 国产 精品1| 日日啪夜夜撸| 中文字幕制服av| 亚洲av男天堂| 久久精品熟女亚洲av麻豆精品| 免费av中文字幕在线| 高清视频免费观看一区二区| 在线观看国产h片| 99国产精品免费福利视频| 91久久精品国产一区二区成人| 精品国产一区二区三区久久久樱花 | 三级国产精品片| freevideosex欧美| av国产精品久久久久影院| 天堂俺去俺来也www色官网| 纵有疾风起免费观看全集完整版| 久久青草综合色| 亚洲人成网站在线观看播放| 中文精品一卡2卡3卡4更新| 久久久久久久大尺度免费视频| 两个人的视频大全免费| 亚洲四区av| 欧美xxⅹ黑人| 丰满少妇做爰视频| 99国产精品免费福利视频| 男女边吃奶边做爰视频| 最近的中文字幕免费完整| 看免费成人av毛片| 在线观看三级黄色| 超碰97精品在线观看| 精品酒店卫生间| 国产精品人妻久久久影院| 极品教师在线视频| 国产乱人视频| 国内精品宾馆在线| 国产av精品麻豆| 国产成人a∨麻豆精品| 99国产精品免费福利视频| 久久久久人妻精品一区果冻| 小蜜桃在线观看免费完整版高清| 777米奇影视久久| 国产精品麻豆人妻色哟哟久久| 少妇人妻久久综合中文| av专区在线播放| 亚洲综合色惰| 尤物成人国产欧美一区二区三区| 视频区图区小说| 亚洲精品一区蜜桃| 夫妻午夜视频| 精品酒店卫生间| 男人和女人高潮做爰伦理| 久久99精品国语久久久| 又爽又黄a免费视频| 亚洲欧洲日产国产| 一个人看视频在线观看www免费| 亚洲高清免费不卡视频| 亚洲精品aⅴ在线观看| 免费看av在线观看网站| 最近最新中文字幕免费大全7| 日韩欧美一区视频在线观看 | 亚洲精品第二区| 国产真实伦视频高清在线观看| 国产深夜福利视频在线观看| 少妇的逼好多水| 人人妻人人爽人人添夜夜欢视频 | 最近2019中文字幕mv第一页| 久久精品国产自在天天线| 免费少妇av软件| 久久人人爽av亚洲精品天堂 | 欧美精品亚洲一区二区| 精品国产乱码久久久久久小说| 成人美女网站在线观看视频| 综合色丁香网| 天堂中文最新版在线下载| 国产亚洲一区二区精品| 日本vs欧美在线观看视频 | 婷婷色麻豆天堂久久| 精品99又大又爽又粗少妇毛片| 国产男人的电影天堂91| 欧美+日韩+精品| 婷婷色综合www| 免费人妻精品一区二区三区视频| 新久久久久国产一级毛片| 亚洲国产欧美人成| 少妇高潮的动态图| 人妻少妇偷人精品九色| 免费播放大片免费观看视频在线观看| h视频一区二区三区| 亚洲精品久久午夜乱码| 三级国产精品欧美在线观看| 国产v大片淫在线免费观看| av又黄又爽大尺度在线免费看| 中国美白少妇内射xxxbb| 能在线免费看毛片的网站| 亚洲精品日韩在线中文字幕| 直男gayav资源| 男女国产视频网站| 国产男女超爽视频在线观看| 少妇熟女欧美另类| 国产成人免费观看mmmm| 久久av网站| 一本久久精品| 91aial.com中文字幕在线观看| 国产高清三级在线| 久久久久久久精品精品| 1000部很黄的大片| 欧美日韩亚洲高清精品| 一级片'在线观看视频| 亚洲欧美日韩东京热| 成年人午夜在线观看视频| 久久国产乱子免费精品| 精品一区二区三卡| 啦啦啦视频在线资源免费观看| 欧美精品人与动牲交sv欧美| 亚洲精品乱码久久久v下载方式| 免费av中文字幕在线| 波野结衣二区三区在线| 在线免费观看不下载黄p国产| 视频区图区小说| 赤兔流量卡办理| 亚洲婷婷狠狠爱综合网| 日产精品乱码卡一卡2卡三| 伊人久久国产一区二区| 18禁在线播放成人免费| 国产成人精品久久久久久| 欧美日韩亚洲高清精品| 97在线视频观看| 国产精品一区二区在线观看99| 国产成人aa在线观看| 久久久成人免费电影| 日韩成人伦理影院| 一本一本综合久久| 又大又黄又爽视频免费| 国产精品人妻久久久久久| 午夜福利影视在线免费观看| 欧美日韩视频精品一区| 成年女人在线观看亚洲视频| 91精品国产国语对白视频| 在线观看av片永久免费下载| 久久精品国产亚洲av天美| 舔av片在线| 午夜福利在线观看免费完整高清在| 极品教师在线视频| 黑人猛操日本美女一级片| 美女脱内裤让男人舔精品视频| 妹子高潮喷水视频| 国产国拍精品亚洲av在线观看| 国产白丝娇喘喷水9色精品| 秋霞在线观看毛片| 97超视频在线观看视频| 国产欧美亚洲国产| 国产成人aa在线观看| 亚洲综合色惰| 国产探花极品一区二区| 国产美女午夜福利| 我的女老师完整版在线观看| 伊人久久精品亚洲午夜| 国产黄色视频一区二区在线观看| 日韩制服骚丝袜av| 国产精品99久久久久久久久| 丝袜脚勾引网站| 亚洲伊人久久精品综合| 99re6热这里在线精品视频| 国产 精品1| 日日撸夜夜添| 国产成人freesex在线| 狂野欧美激情性bbbbbb| 18禁裸乳无遮挡动漫免费视频| 国产欧美亚洲国产| 99热全是精品| 又粗又硬又长又爽又黄的视频| 亚洲成色77777| 热re99久久精品国产66热6| 国产精品免费大片| 精品国产露脸久久av麻豆| 女的被弄到高潮叫床怎么办| 国产精品欧美亚洲77777| 日本免费在线观看一区| 国产精品蜜桃在线观看| 一级毛片电影观看| 成年美女黄网站色视频大全免费 | 妹子高潮喷水视频| 精品久久国产蜜桃| 国产av精品麻豆|