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

    無金標(biāo)準(zhǔn)部分核實(shí)數(shù)據(jù)下基于風(fēng)險(xiǎn)差的等價(jià)性檢驗(yàn)

    2019-11-15 06:54:42邱世芳
    關(guān)鍵詞:假設(shè)檢驗(yàn)樣本量分類器

    邱世芳,何 杰

    (重慶理工大學(xué) 理學(xué)院, 重慶 400054)

    在流行病學(xué)研究中,為了估計(jì)某種疾病的流行率,常使用價(jià)格便宜的篩檢方法對(duì)個(gè)體進(jìn)行第一次診斷。由于篩檢方法常存在誤判,從而導(dǎo)致基于有誤判數(shù)據(jù)的估計(jì)常常是有偏的[1]。在實(shí)際中,完全的金標(biāo)準(zhǔn)是不存在的,同時(shí)其花費(fèi)很高且非常耗時(shí)。為此,Tenenbein[2]提出了二重抽樣方法,即從感興趣的總體中隨機(jī)抽取N個(gè)個(gè)體接受易犯錯(cuò)的分類器分類,再從中隨機(jī)抽取n個(gè)個(gè)體接受金標(biāo)準(zhǔn)檢驗(yàn)。這種方法叫作二重抽樣方法,得到的數(shù)據(jù)稱為部分核實(shí)數(shù)據(jù)[3]。

    基于Tenenbein[2]提出的二重抽樣方法得到的部分核實(shí)數(shù)據(jù),Tang等[3]研究了單個(gè)總體下的疾病流行率的假設(shè)問題,考慮了基于Score檢驗(yàn)、似然比檢驗(yàn)以及2種Wald型檢驗(yàn)的大樣本漸近的檢驗(yàn)過程以及小樣本近似非條件檢驗(yàn)過程;Tang等[4]基于上述檢驗(yàn)統(tǒng)計(jì)量提出了關(guān)于疾病流行率的12種區(qū)間估計(jì)方法;并從顯著性檢驗(yàn)的角度和區(qū)間估計(jì)的角度考慮了樣本量的確定方法[5];對(duì)于2組樣本情形,Tang等[6]基于比例差考慮了2組疾病流行率的顯著性差異的假設(shè)檢驗(yàn)過程,以及從檢驗(yàn)功效的角度研究了樣本量的確定問題;同時(shí),Qiu等[7]考慮了2組樣本下關(guān)于疾病流行率差的區(qū)間估計(jì)方法,從置信區(qū)間的角度研究了2組樣本的疾病流行率的顯著性差異問題。

    以上研究中有一個(gè)重要的假定是存在一種完全無誤判的金標(biāo)準(zhǔn)。但在實(shí)際應(yīng)用中,完全準(zhǔn)確的金標(biāo)準(zhǔn)檢驗(yàn)也許根本不存在。為此,安金兵等[8]考慮了在無金標(biāo)準(zhǔn)條件下的診斷試驗(yàn);Qiu等[9]在2種分類器都有誤判的情況下,基于單組樣本下的部分核實(shí)數(shù)據(jù)考慮了基于疾病流行率的假設(shè)檢驗(yàn)過程以及樣本量的確定,以及邱世芳等[10]從區(qū)間估計(jì)的角度研究了在給定置信水平下區(qū)間寬度控制下的樣本量的估計(jì)問題。但Qiu等[9-10]只考慮單組樣本下的部分核實(shí)數(shù)據(jù)對(duì)疾病流行率的顯著性檢驗(yàn)過程。而實(shí)際生活中,疾病流行率可能會(huì)受到年齡、生活習(xí)慣等的影響。為此,本文在此基礎(chǔ)上進(jìn)一步考慮2組獨(dú)立的無金標(biāo)準(zhǔn)部分核實(shí)數(shù)據(jù)下基于風(fēng)險(xiǎn)差的顯著性檢驗(yàn)問題。

    1 數(shù)據(jù)結(jié)構(gòu)

    在2種分類器都有誤判的情況下,假設(shè)從第j組感興趣總體中隨機(jī)抽取了Nj個(gè)個(gè)體進(jìn)行初級(jí)分類器(J)分類,再從Nj個(gè)個(gè)體中隨機(jī)抽取nj個(gè)個(gè)體進(jìn)行高級(jí)分類器(S)分類,得到頻數(shù)為nikj(i,k=0,1,j=1,2)。其中,Nj-nj個(gè)個(gè)體只接受了初級(jí)分類器的分類;xj表示檢驗(yàn)為陽性的個(gè)體數(shù);yj表示檢驗(yàn)為陰性的個(gè)體數(shù),數(shù)據(jù)結(jié)構(gòu)由表1所示。

    表1 第j組的數(shù)據(jù)結(jié)構(gòu)

    令Dj=1表示患病,Dj=0表示不患病,Jj=1表示初級(jí)分類器診斷個(gè)體呈陽性,反之Jj=0;Sj=1表示高級(jí)分類器診斷個(gè)體呈陽性,反之Sj=0。令πj=P(Dj=1)為第j組感興趣總體真實(shí)的患病率,ηj=P(Jj=1|Dj=1)和θj=P(Sj=1|Dj=1)分別表示第j組樣本下初級(jí)分類器和高級(jí)分類器的敏感度。本文假定兩種分類器都不存在假陽性。

    2 假設(shè)檢驗(yàn)

    本文感興趣的問題是如下的假設(shè)檢驗(yàn):

    H0:π1=π2?H1:π1≠π2

    令δ=π2-π1,上述的假設(shè)檢驗(yàn)就可變?yōu)椋?/p>

    H0:δ=0?H1:δ≠0

    2.1 模型1和檢驗(yàn)統(tǒng)計(jì)量

    假定初級(jí)分類器和高級(jí)分類器的分類結(jié)果是條件獨(dú)立的,即滿足P(Jj,Sj|Dj)=P(Jj|Dj)P(Sj|Dj),j=1,2。在此假定下,可得如下的概率模型(Nedelman[11]),其概率分布如表2所示。

    表2 第j組的概率分布

    令m={(n11j,n10j,n01j,n00j,xj,yj):j=1,2}為觀測數(shù)據(jù),則觀測數(shù)據(jù)m的對(duì)數(shù)似然函數(shù)為:

    l1(m;π1,δ,η1,θ1,η2,θ2)=

    n01j+xj)logηj+n10jlog(1-ηj)+

    (n11j+n10j)logθj+n01jlog(1-θj)+

    n00jlog[1-πj(ηj+θj-ηjθj)]+

    yjlog(1-πjηj)}

    (1)

    其中,π2=π1+δ;C1是一個(gè)與參數(shù)無關(guān)的常數(shù)。

    簡單計(jì)算可得π1、δ、ηj、θj(j=1,2)的非限制性極大似然估計(jì),分別為:

    (2)

    (3)

    在H0:δ=0下,π1的限制性極大似然估計(jì)為:

    (4)

    (5)

    此方程組的解沒有顯表達(dá)式,可用迭代方法如牛頓迭代法求解。

    2.1.1Wald型檢驗(yàn)統(tǒng)計(jì)量

    在原假設(shè)H0:δ=0下,當(dāng)N1和N2趨近無窮大時(shí),Tw1和Tw2的漸近分布為標(biāo)準(zhǔn)正態(tài)分布。

    2.1.2似然比檢驗(yàn)統(tǒng)計(jì)量

    根據(jù)式(1)給出的對(duì)數(shù)似然函數(shù),對(duì)于假設(shè)檢驗(yàn)H0:δ=0?H1:δ≠0有如下的似然比檢驗(yàn)統(tǒng)計(jì)量:

    在原假設(shè)H0:δ=0下,當(dāng)N1和N2趨近無窮大時(shí),Tl的漸近分布為自由度為1的卡方分布。

    2.1.3Score檢驗(yàn)統(tǒng)計(jì)量

    根據(jù)Rao[12]提出的score理論,對(duì)于假設(shè)檢驗(yàn)H0:δ=0?H1:δ≠0,經(jīng)過推導(dǎo)得到如下的score檢驗(yàn)統(tǒng)計(jì)量:

    2.1.4基于反雙曲正切變換的檢驗(yàn)統(tǒng)計(jì)量

    在原假設(shè)H0:δ=0下,當(dāng)N1和N2趨近無窮大時(shí),Ttan1和Ttan2的漸近分布為標(biāo)準(zhǔn)正態(tài)分布。

    2.2 模型2和檢驗(yàn)統(tǒng)計(jì)量

    假定初級(jí)分類器和高級(jí)分類器存在假陽性且條件不獨(dú)立時(shí),根據(jù)Lie[13]提出的模型,其概率分布如表3所示。

    表3 第j組的概率分布

    注:j=1,2

    觀測樣本m的對(duì)數(shù)似然函數(shù)為:

    l2(m;π1,δ,η1,θ1,η2,θ2)=

    C2+(nj-n00j+xj)log(πj)+

    n00jlog(1-πj)+xjlogηj+

    n10jlog(1-ηj)+n01jlog(1-θj)+

    yjlog(1-πjηj)+n11jlog(ηj+θj-1)

    (6)

    其中:π2=π1+δ;C2是一個(gè)與參數(shù)無關(guān)的常數(shù)。

    經(jīng)過簡單計(jì)算,可得π1,δ,ηj,θj(j=1,2)的非限制性極大似然估計(jì),分別為:

    (7)

    (8)

    在H0:δ=0下,π1的限制性極大似然估計(jì)為:

    (9)

    由此,得到如下的估計(jì):

    (10)

    2.2.1Wald型檢驗(yàn)統(tǒng)計(jì)量

    2.2.2似然比檢驗(yàn)統(tǒng)計(jì)量

    由式(6)可得到如下的似然比檢驗(yàn)統(tǒng)計(jì)量:

    在原假設(shè)H0:δ=0下,當(dāng)N1和N2趨近無窮大時(shí),Tl的漸近分布為自由度為1的卡方分布。

    2.2.3Score檢驗(yàn)統(tǒng)計(jì)量

    在原假設(shè)H0:δ=0成立下,可以得到Score檢驗(yàn)統(tǒng)計(jì)量:

    2.2.4基于反雙曲正切變換的檢驗(yàn)統(tǒng)計(jì)量

    同模型1,可以得到2種反雙曲正切逆變換的檢驗(yàn)統(tǒng)計(jì)量:

    在原假設(shè)H0:δ=0下,當(dāng)N1和N2趨近無窮大時(shí),Ttan1和Ttan2的漸近分布為標(biāo)準(zhǔn)正態(tài)分布。

    3 模擬研究

    為評(píng)估Tw1、Tw2、Tsc、Tl、Ttan1和Ttan2檢驗(yàn)方法的有效性,考慮了如下樣本量的模擬研究:① 小樣本:(n1,n2,N1,N2)=(20,20,30,30),(30,30,50,50);② 中等樣本(n1,n2,N1,N2)=(50,50,100,100);③ 大樣本:(n1,n2,N1,N2)=(200,200,500,500)。

    在不同的樣本量設(shè)置下,考察犯第一類錯(cuò)誤的概率時(shí)考慮了如下參數(shù)設(shè)置:π1=0.1(0.2)0.5,η1=0.5(0.1)0.7,θ1=0.7(0.1)0.9,η2=η1+0.05,θ2=θ1+0.05,即考慮了3(π1的值)×3(η1(η2)的值)×3(θ1(θ2)的值)=27種參數(shù)組合。對(duì)于檢驗(yàn)功效考慮了:δ=0.1(0.01)0.3,π1=0.1(0.2)0.5,(η1,θ1)=(0.5,0.7),(0.6,0.8),(0.7,0.9),η2,θ2同以上設(shè)置,即考慮了3(δ的值)×3(π1的值)×3((η1,θ1)的值)=27種參數(shù)組合。這里的a(b)c表示取值是從a以步長b變化到c。在每個(gè)樣本量設(shè)置和每種參數(shù)設(shè)置下,隨機(jī)產(chǎn)生5 000組數(shù)據(jù)m={(n11j,n10j,n01j,n00j,xj,yj):j=1,2},在顯著性水平α=0.05下,對(duì)于每個(gè)檢驗(yàn)統(tǒng)計(jì)量Ti(i=w1,w2,sc,l,tan1,tan2),其犯第一類錯(cuò)誤的概率可通過以下公式計(jì)算:基于統(tǒng)計(jì)量Ti拒絕原假設(shè)的次數(shù)/5 000(δ=0),經(jīng)驗(yàn)功效通過以下公式計(jì)算基于統(tǒng)計(jì)量Ti拒絕原假設(shè)的次數(shù)/5 000(δ≠0)。犯第一類錯(cuò)誤概率的模擬結(jié)果見表4~7。由于篇幅的限制,只列出了小樣本和中等樣本下的部分功效的模擬結(jié)果,見表8~9。

    表4 小樣本(n1,n2,N1,N2)=(20,20,30,30)下犯第一類錯(cuò)誤的概率(顯著性水平α=0.05)

    表5 小樣本(n1,n2,N1,N2)=(30,30,50,50)下犯第一類錯(cuò)誤的概率(顯著性水平α=0.05)

    續(xù)表(表5)

    表6 中等樣本(n1,n2,N1,N2)=(50,50,100,100)下犯第一類錯(cuò)誤的概率(顯著性水平α=0.05)

    續(xù)表(表6)

    表7 大樣本(n1,n2,N1,N2)=(200,200,500,500)下犯第一類錯(cuò)誤的概率(α=0.05)

    表8 小樣本(n1,n2,N1,N2)=(20,20,30,30)下的經(jīng)驗(yàn)功效(顯著性水平α=0.05)

    表9 中等樣本(n1,n2,N1,N2)=(50,50,100,100)下的經(jīng)驗(yàn)功效(顯著性水平α=0.05)

    續(xù)表(表9)

    通過模擬結(jié)果可以得到以下結(jié)論:① 對(duì)于2個(gè)模型,即使在極小樣本(如(n1,n2,N1,N2)=(20,20,30,30))下,Score檢驗(yàn)統(tǒng)計(jì)量(Tsc)犯第一類錯(cuò)誤概率接近給定的顯著性水平(α=0.05),相應(yīng)的功效也比較大,隨著樣本量的增大,犯第一類錯(cuò)誤的概率越來越接近顯著性水平0.05且功效越來越大;② 對(duì)于模型1,除了小樣本情形,考慮原假設(shè)下方差的Wald檢驗(yàn)統(tǒng)計(jì)量(Tw2)、Score檢驗(yàn)(Tsc)以及反雙曲正切變換統(tǒng)計(jì)量(Ttan2)都表現(xiàn)良好,它們犯第一類錯(cuò)誤的概率都比較接近于給定的顯著性水平0.05,且有較大的功效;③ 對(duì)于模型2,在中等樣本和大樣本且犯病率不是極小時(shí),所有檢驗(yàn)統(tǒng)計(jì)量犯第一類錯(cuò)誤的概率都基本接近名義水平且具有較高的功效。綜上所述,在模型1下,推薦使用Tw2、Tsc和Ttan1;在模型2下,當(dāng)樣本量極小時(shí)(如(n1,n2,N1,N2)=(20,20,30,30)),推薦使用Score統(tǒng)計(jì)量,在樣本量不是極小時(shí),所有統(tǒng)計(jì)量都可以推薦使用。

    4 實(shí)例分析

    為了驗(yàn)證本文提出方法的有效性,利用Nedelman[11]研究中的瘧疾數(shù)據(jù)進(jìn)行分析??紤]了未成年組(即9~18歲)的2個(gè)不同組的瘧疾數(shù)據(jù),數(shù)據(jù)結(jié)構(gòu)見表10。

    表10 實(shí)例數(shù)據(jù)

    對(duì)于假設(shè)檢驗(yàn):H0:δ=0?H1:δ≠0,在模型1下,基于本文考慮的檢驗(yàn)統(tǒng)計(jì)量Tw1、Tw2、Tl、Tsc、Ttan1、Ttan2的檢驗(yàn)p值分別為:0.113、0.130、0.129、0.184、0.115、0.132。由此,沒有充分理由拒絕原假設(shè),認(rèn)為瘧疾的患病率與不同的調(diào)查組無關(guān)。在模型2下,本文所提出的檢驗(yàn)統(tǒng)計(jì)量的檢驗(yàn)p值分別為:0.127、0.141、0.143、0.145、0.128、0.142,同樣認(rèn)為沒有充分理由拒絕原假設(shè),即認(rèn)為瘧疾的患病率與不同的調(diào)查組無關(guān)。

    5 結(jié)束語

    本文基于2個(gè)有誤判分類器下的二重抽樣數(shù)據(jù)對(duì)2組患病率進(jìn)行了比較研究,分別考慮2個(gè)二重抽樣模型的6種檢驗(yàn)統(tǒng)計(jì)量,并通過不同樣本量和不同參數(shù)設(shè)置考察了它們犯第一類錯(cuò)誤的概率和檢驗(yàn)功效。對(duì)于模型1,即使在極小樣本下,Score檢驗(yàn)以及反雙曲正切統(tǒng)計(jì)量都表現(xiàn)良好,在樣本量不太小時(shí),似然比統(tǒng)計(jì)量也能很好地控制犯第一類錯(cuò)誤的概率;對(duì)于模型2,本文所考慮的所有檢驗(yàn)統(tǒng)計(jì)量(在中等樣本和大樣本下且犯病率不是極小時(shí))都表現(xiàn)很好。在2個(gè)模型中,考慮初級(jí)分類器和高級(jí)分類器不存在假陽性,在實(shí)際的模型中該條件可能不存在,但如果2種分類器都存在假陽性分類時(shí),本文考慮的模型是不可識(shí)別的。此時(shí),可以考慮采用貝葉斯統(tǒng)計(jì)方法進(jìn)行研究,這是一個(gè)將來值得研究的問題。

    猜你喜歡
    假設(shè)檢驗(yàn)樣本量分類器
    醫(yī)學(xué)研究中樣本量的選擇
    航空裝備測試性試驗(yàn)樣本量確定方法
    BP-GA光照分類器在車道線識(shí)別中的應(yīng)用
    電子測試(2018年1期)2018-04-18 11:52:35
    Sample Size Calculations for Comparing Groups with Binary Outcomes
    統(tǒng)計(jì)推斷的研究
    雙冪變換下正態(tài)線性回歸模型參數(shù)的假設(shè)檢驗(yàn)
    Primary Question and Hypothesis Testing in Randomized Controlled Clinical Trials
    統(tǒng)計(jì)學(xué)教學(xué)中關(guān)于假設(shè)檢驗(yàn)問題探討
    加權(quán)空-譜與最近鄰分類器相結(jié)合的高光譜圖像分類
    結(jié)合模糊(C+P)均值聚類和SP-V-支持向量機(jī)的TSK分類器
    天美传媒精品一区二区| 免费高清在线观看视频在线观看| 天堂中文最新版在线下载 | 非洲黑人性xxxx精品又粗又长| 中文乱码字字幕精品一区二区三区 | 久久久久性生活片| 欧美3d第一页| 午夜久久久久精精品| 亚洲三级黄色毛片| 色5月婷婷丁香| 国产伦一二天堂av在线观看| 十八禁网站网址无遮挡 | 乱码一卡2卡4卡精品| 国产精品女同一区二区软件| 九九久久精品国产亚洲av麻豆| 哪个播放器可以免费观看大片| 欧美日韩精品成人综合77777| 国产精品无大码| 在线a可以看的网站| 国产伦在线观看视频一区| 国产成人a区在线观看| 久久精品夜色国产| 国产精品女同一区二区软件| 日韩三级伦理在线观看| 只有这里有精品99| 蜜臀久久99精品久久宅男| 男女下面进入的视频免费午夜| 国产成年人精品一区二区| 国产精品爽爽va在线观看网站| 女人十人毛片免费观看3o分钟| 亚洲欧美清纯卡通| 最近的中文字幕免费完整| 18禁动态无遮挡网站| 国产色爽女视频免费观看| 三级国产精品欧美在线观看| 久久精品国产亚洲av天美| 国产伦精品一区二区三区四那| 国产精品99久久久久久久久| 黄色日韩在线| 国产三级在线视频| 99热网站在线观看| 青春草亚洲视频在线观看| 69人妻影院| 一区二区三区免费毛片| 91精品伊人久久大香线蕉| 99久久中文字幕三级久久日本| 亚洲欧美成人综合另类久久久| 久久久久久久久大av| 天堂√8在线中文| 国产精品久久久久久av不卡| 少妇裸体淫交视频免费看高清| 欧美性感艳星| 免费观看的影片在线观看| 国产精品国产三级国产专区5o| 午夜激情福利司机影院| 国产在视频线在精品| 日韩成人伦理影院| www.av在线官网国产| 免费无遮挡裸体视频| 黄色日韩在线| 亚洲三级黄色毛片| 欧美高清性xxxxhd video| 亚洲精品成人av观看孕妇| 日本欧美国产在线视频| 精品久久久久久久末码| 能在线免费看毛片的网站| 国产69精品久久久久777片| kizo精华| 久久99蜜桃精品久久| 舔av片在线| 人妻系列 视频| 国产伦理片在线播放av一区| 国产一区二区三区av在线| 偷拍熟女少妇极品色| 亚洲成人一二三区av| 国产精品99久久久久久久久| 亚洲av中文av极速乱| 欧美bdsm另类| 国产色婷婷99| 亚洲丝袜综合中文字幕| 99久久九九国产精品国产免费| 肉色欧美久久久久久久蜜桃 | 亚洲精品一区蜜桃| 精品国内亚洲2022精品成人| 成年女人看的毛片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 日本免费a在线| 色视频www国产| 国产日韩欧美在线精品| 人人妻人人澡人人爽人人夜夜 | 久久久精品94久久精品| 亚洲精品久久午夜乱码| 久久久久网色| 久久久久性生活片| 免费看av在线观看网站| 欧美成人精品欧美一级黄| 欧美丝袜亚洲另类| 成人鲁丝片一二三区免费| 青春草视频在线免费观看| 国产精品一区二区三区四区免费观看| 亚洲经典国产精华液单| 亚洲经典国产精华液单| 我的老师免费观看完整版| 只有这里有精品99| 久久久久久久久大av| 综合色av麻豆| 一级毛片我不卡| 国产精品蜜桃在线观看| 成人高潮视频无遮挡免费网站| 26uuu在线亚洲综合色| 18禁在线无遮挡免费观看视频| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 国产高清国产精品国产三级 | 亚洲国产精品成人久久小说| 国产成人一区二区在线| 国产精品国产三级国产av玫瑰| 精品欧美国产一区二区三| 亚洲丝袜综合中文字幕| 91av网一区二区| 日韩av在线大香蕉| 成人欧美大片| 男人舔女人下体高潮全视频| 精品久久久精品久久久| 午夜免费观看性视频| 免费高清在线观看视频在线观看| 美女内射精品一级片tv| 欧美变态另类bdsm刘玥| 亚洲欧美成人综合另类久久久| 亚洲欧美清纯卡通| 麻豆成人av视频| 久久精品久久久久久噜噜老黄| 午夜老司机福利剧场| 午夜精品一区二区三区免费看| 国产永久视频网站| 国精品久久久久久国模美| videos熟女内射| 两个人的视频大全免费| 69人妻影院| 欧美成人午夜免费资源| 国产精品99久久久久久久久| 国产免费一级a男人的天堂| 美女主播在线视频| 丰满乱子伦码专区| 毛片女人毛片| 三级毛片av免费| 国产精品蜜桃在线观看| 亚洲在线自拍视频| 国产伦一二天堂av在线观看| 两个人的视频大全免费| 毛片一级片免费看久久久久| 久久精品国产亚洲网站| 中文字幕人妻熟人妻熟丝袜美| 中国美白少妇内射xxxbb| 大话2 男鬼变身卡| 少妇被粗大猛烈的视频| 日日撸夜夜添| 亚洲国产精品成人综合色| 中文字幕人妻熟人妻熟丝袜美| 久久久久久九九精品二区国产| 综合色丁香网| 中文字幕久久专区| 又大又黄又爽视频免费| 好男人在线观看高清免费视频| 欧美高清成人免费视频www| 在线观看一区二区三区| 午夜激情欧美在线| 看免费成人av毛片| 国产成人福利小说| 午夜福利成人在线免费观看| 中文字幕av成人在线电影| 观看美女的网站| 精品人妻一区二区三区麻豆| 国产av在哪里看| 男女视频在线观看网站免费| 久久久久久久久中文| 五月伊人婷婷丁香| 五月天丁香电影| 91狼人影院| 菩萨蛮人人尽说江南好唐韦庄| 好男人在线观看高清免费视频| 欧美人与善性xxx| 美女国产视频在线观看| 成年免费大片在线观看| 成人毛片60女人毛片免费| av在线天堂中文字幕| 尤物成人国产欧美一区二区三区| 日韩欧美精品v在线| 一级毛片我不卡| 免费看不卡的av| 亚洲美女视频黄频| 成年免费大片在线观看| 日韩欧美 国产精品| 日韩成人av中文字幕在线观看| 尤物成人国产欧美一区二区三区| 国产一区二区三区av在线| 亚洲综合精品二区| av在线亚洲专区| 日韩制服骚丝袜av| 免费看日本二区| 日本黄大片高清| 综合色丁香网| 欧美变态另类bdsm刘玥| 亚洲自偷自拍三级| 亚洲精品乱码久久久v下载方式| 永久免费av网站大全| 青青草视频在线视频观看| 伊人久久国产一区二区| 国产女主播在线喷水免费视频网站 | 欧美 日韩 精品 国产| 永久网站在线| 亚洲婷婷狠狠爱综合网| 亚洲在久久综合| xxx大片免费视频| 又爽又黄无遮挡网站| 大陆偷拍与自拍| 国产亚洲最大av| 直男gayav资源| 日韩精品青青久久久久久| 大片免费播放器 马上看| 日韩欧美一区视频在线观看 | 精品一区二区三卡| 亚洲精品成人久久久久久| 一本一本综合久久| 国产一区二区三区综合在线观看 | 亚洲一级一片aⅴ在线观看| 久久久久精品性色| 国产视频内射| av黄色大香蕉| 国产伦一二天堂av在线观看| 国产极品天堂在线| 国产成人91sexporn| 国产免费视频播放在线视频 | 亚洲精品久久久久久婷婷小说| 亚洲欧洲日产国产| 五月天丁香电影| 大又大粗又爽又黄少妇毛片口| 欧美三级亚洲精品| 国产精品一及| 午夜福利视频1000在线观看| 天堂影院成人在线观看| 欧美区成人在线视频| 亚洲av男天堂| 能在线免费看毛片的网站| 中文在线观看免费www的网站| 全区人妻精品视频| 亚洲国产欧美在线一区| 国产精品嫩草影院av在线观看| 一个人看视频在线观看www免费| 午夜激情久久久久久久| 亚洲久久久久久中文字幕| 婷婷色麻豆天堂久久| 日本一二三区视频观看| 黄色日韩在线| 亚洲aⅴ乱码一区二区在线播放| av专区在线播放| 亚洲国产av新网站| 国产v大片淫在线免费观看| 在线观看av片永久免费下载| 毛片女人毛片| 99视频精品全部免费 在线| 国产精品女同一区二区软件| 18+在线观看网站| 日韩三级伦理在线观看| 极品少妇高潮喷水抽搐| 搡女人真爽免费视频火全软件| 国产精品国产三级国产专区5o| 身体一侧抽搐| 午夜免费激情av| 免费看光身美女| 精品一区在线观看国产| 亚洲丝袜综合中文字幕| 午夜激情久久久久久久| 乱码一卡2卡4卡精品| 美女主播在线视频| 91精品国产九色| 亚洲精品视频女| 最近手机中文字幕大全| av天堂中文字幕网| 丝袜美腿在线中文| 国产成人freesex在线| 国产一区二区亚洲精品在线观看| 午夜精品一区二区三区免费看| 亚洲国产av新网站| 欧美高清性xxxxhd video| 中文欧美无线码| 国产人妻一区二区三区在| 在线天堂最新版资源| 国产高清有码在线观看视频| 伊人久久精品亚洲午夜| 乱人视频在线观看| 赤兔流量卡办理| 国产一区亚洲一区在线观看| 又黄又爽又刺激的免费视频.| 26uuu在线亚洲综合色| 人人妻人人澡人人爽人人夜夜 | 中文字幕人妻熟人妻熟丝袜美| 久久久久久久亚洲中文字幕| 亚洲成人中文字幕在线播放| 国产视频内射| 只有这里有精品99| 韩国av在线不卡| 成人毛片a级毛片在线播放| 丰满少妇做爰视频| 色网站视频免费| 亚洲国产精品成人综合色| 国产成人一区二区在线| 免费大片18禁| 熟妇人妻久久中文字幕3abv| 综合色丁香网| 日韩国内少妇激情av| 青春草亚洲视频在线观看| 午夜免费观看性视频| xxx大片免费视频| 极品少妇高潮喷水抽搐| 亚洲综合色惰| 又爽又黄a免费视频| 久久精品夜夜夜夜夜久久蜜豆| 免费看日本二区| 人人妻人人澡欧美一区二区| 久久久久性生活片| 人人妻人人澡欧美一区二区| 麻豆av噜噜一区二区三区| 久久热精品热| 亚洲国产色片| 欧美精品一区二区大全| 国产成人aa在线观看| 人妻夜夜爽99麻豆av| 国产亚洲av嫩草精品影院| 欧美日韩综合久久久久久| 久久韩国三级中文字幕| 久久久久久久久久久免费av| 精品久久久久久久末码| 人体艺术视频欧美日本| 春色校园在线视频观看| 夜夜看夜夜爽夜夜摸| 能在线免费看毛片的网站| 黄色日韩在线| 午夜老司机福利剧场| av福利片在线观看| 一区二区三区免费毛片| 国产有黄有色有爽视频| 最近手机中文字幕大全| 大话2 男鬼变身卡| 99久国产av精品国产电影| 99视频精品全部免费 在线| 亚洲高清免费不卡视频| 波野结衣二区三区在线| 国产在视频线精品| 精品一区二区三区人妻视频| 亚洲国产色片| 亚洲精品乱码久久久久久按摩| 精品人妻熟女av久视频| 国产黄色免费在线视频| 韩国高清视频一区二区三区| 日韩欧美精品v在线| 欧美丝袜亚洲另类| 亚洲欧美日韩无卡精品| 久久草成人影院| 美女大奶头视频| 成人鲁丝片一二三区免费| 黄片无遮挡物在线观看| 超碰97精品在线观看| 国产黄频视频在线观看| 如何舔出高潮| 草草在线视频免费看| 中文字幕免费在线视频6| 3wmmmm亚洲av在线观看| 高清毛片免费看| 最新中文字幕久久久久| 精品久久久久久久久久久久久| 国产精品久久久久久精品电影| 高清在线视频一区二区三区| 亚洲av免费在线观看| 免费黄频网站在线观看国产| 亚洲av福利一区| 精品国内亚洲2022精品成人| av一本久久久久| 久久久欧美国产精品| 久久久精品欧美日韩精品| 国产精品不卡视频一区二区| 国产综合精华液| 国产探花极品一区二区| 欧美zozozo另类| 大香蕉久久网| 一区二区三区乱码不卡18| 熟女电影av网| 成人av在线播放网站| 午夜福利高清视频| 韩国高清视频一区二区三区| 国产真实伦视频高清在线观看| 免费看av在线观看网站| 免费播放大片免费观看视频在线观看| 午夜久久久久精精品| 久久久久久久久久黄片| 国产伦精品一区二区三区四那| 激情五月婷婷亚洲| 亚洲在久久综合| 久久6这里有精品| 免费无遮挡裸体视频| 一区二区三区免费毛片| 国内精品美女久久久久久| 精品一区二区三区视频在线| 国产极品天堂在线| 亚洲av不卡在线观看| 国产精品av视频在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲av免费在线观看| 啦啦啦韩国在线观看视频| 成人特级av手机在线观看| 国产亚洲精品av在线| 一级二级三级毛片免费看| 亚洲在久久综合| 91久久精品国产一区二区成人| .国产精品久久| 午夜福利视频1000在线观看| 久久久久久九九精品二区国产| 亚洲精品成人av观看孕妇| 日本与韩国留学比较| 国产精品蜜桃在线观看| 日日撸夜夜添| 日韩一区二区视频免费看| 淫秽高清视频在线观看| 成人高潮视频无遮挡免费网站| 欧美成人a在线观看| 九九在线视频观看精品| 麻豆成人午夜福利视频| 久久99蜜桃精品久久| 午夜福利成人在线免费观看| 在线免费观看的www视频| 啦啦啦韩国在线观看视频| 噜噜噜噜噜久久久久久91| 十八禁网站网址无遮挡 | 日韩av不卡免费在线播放| 青春草国产在线视频| 日韩成人av中文字幕在线观看| 人妻制服诱惑在线中文字幕| av国产久精品久网站免费入址| 夜夜看夜夜爽夜夜摸| 午夜爱爱视频在线播放| 国产 一区精品| 蜜桃久久精品国产亚洲av| 丰满人妻一区二区三区视频av| 啦啦啦韩国在线观看视频| 我要看日韩黄色一级片| 看免费成人av毛片| 99久久精品国产国产毛片| 国产乱人视频| 一级av片app| 热99在线观看视频| 亚洲在线自拍视频| 青春草亚洲视频在线观看| 岛国毛片在线播放| 高清欧美精品videossex| 91午夜精品亚洲一区二区三区| 麻豆av噜噜一区二区三区| 国内揄拍国产精品人妻在线| 精品久久久噜噜| av一本久久久久| 成人高潮视频无遮挡免费网站| 能在线免费看毛片的网站| 成人毛片a级毛片在线播放| 久久亚洲国产成人精品v| 深爱激情五月婷婷| 久久久久久久久久久免费av| 好男人在线观看高清免费视频| 蜜臀久久99精品久久宅男| 欧美极品一区二区三区四区| xxx大片免费视频| 国产精品不卡视频一区二区| 嫩草影院入口| 国产永久视频网站| 中文字幕制服av| 爱豆传媒免费全集在线观看| 国产精品人妻久久久久久| 免费观看无遮挡的男女| 午夜福利在线观看免费完整高清在| 在线免费观看不下载黄p国产| 精品人妻视频免费看| 亚洲av在线观看美女高潮| 2018国产大陆天天弄谢| 麻豆成人av视频| 欧美xxxx黑人xx丫x性爽| 国产精品一区二区性色av| 建设人人有责人人尽责人人享有的 | 又黄又爽又刺激的免费视频.| a级毛片免费高清观看在线播放| 成人高潮视频无遮挡免费网站| 美女cb高潮喷水在线观看| 国产乱人视频| 99热这里只有精品一区| 国产黄a三级三级三级人| 亚洲在久久综合| 国产在视频线精品| 一级爰片在线观看| or卡值多少钱| 久久久久久久久久人人人人人人| 国产人妻一区二区三区在| 国产三级在线视频| 老女人水多毛片| 嫩草影院精品99| av福利片在线观看| 亚洲精品第二区| 26uuu在线亚洲综合色| 午夜精品一区二区三区免费看| 偷拍熟女少妇极品色| 亚洲成人av在线免费| 97精品久久久久久久久久精品| 久久精品久久久久久噜噜老黄| 黄色欧美视频在线观看| 日本爱情动作片www.在线观看| av.在线天堂| 国产精品国产三级国产专区5o| 99久久精品热视频| 国产在线男女| 亚洲婷婷狠狠爱综合网| 日韩精品青青久久久久久| 欧美成人午夜免费资源| 成人特级av手机在线观看| h日本视频在线播放| 水蜜桃什么品种好| 日韩不卡一区二区三区视频在线| 国产激情偷乱视频一区二区| 在线观看免费高清a一片| 亚洲精品aⅴ在线观看| 国产v大片淫在线免费观看| 精品午夜福利在线看| 搞女人的毛片| 亚洲最大成人手机在线| 免费观看a级毛片全部| 亚洲人成网站高清观看| 一级毛片久久久久久久久女| 一级毛片黄色毛片免费观看视频| 欧美丝袜亚洲另类| 日韩三级伦理在线观看| 不卡视频在线观看欧美| 国产精品嫩草影院av在线观看| 久久国产乱子免费精品| 国产探花在线观看一区二区| 一二三四中文在线观看免费高清| 国内少妇人妻偷人精品xxx网站| 亚洲久久久久久中文字幕| 男女边吃奶边做爰视频| 欧美不卡视频在线免费观看| 亚洲欧美中文字幕日韩二区| av免费观看日本| 日韩视频在线欧美| 久久97久久精品| 午夜福利高清视频| 久久久久久久久久久丰满| 亚洲最大成人av| av福利片在线观看| 天堂俺去俺来也www色官网 | 国产在线一区二区三区精| 淫秽高清视频在线观看| 国产av码专区亚洲av| 国产精品国产三级国产专区5o| 免费电影在线观看免费观看| 亚洲性久久影院| 国产精品嫩草影院av在线观看| 欧美日韩精品成人综合77777| av专区在线播放| 18禁在线无遮挡免费观看视频| 汤姆久久久久久久影院中文字幕 | 亚洲欧洲国产日韩| 日韩在线高清观看一区二区三区| 午夜免费激情av| 哪个播放器可以免费观看大片| 免费人成在线观看视频色| 午夜福利视频精品| 国产极品天堂在线| 国产日韩欧美在线精品| 欧美激情在线99| 男人舔奶头视频| 又爽又黄a免费视频| 在线观看美女被高潮喷水网站| 国产一级毛片七仙女欲春2| 久久久久久久大尺度免费视频| 床上黄色一级片| 成人亚洲精品一区在线观看 | 天美传媒精品一区二区| 91久久精品电影网| 亚洲不卡免费看| 国产高清国产精品国产三级 | 久久久久久久久久人人人人人人| 麻豆av噜噜一区二区三区| 亚洲av成人av| 人人妻人人澡人人爽人人夜夜 | 国产永久视频网站| 午夜久久久久精精品| 免费观看性生交大片5| 国产永久视频网站| 国产亚洲av片在线观看秒播厂 | 啦啦啦啦在线视频资源| 日韩一本色道免费dvd| 国产亚洲91精品色在线| 国产午夜精品一二区理论片| 亚洲国产欧美人成| 十八禁网站网址无遮挡 | 精品欧美国产一区二区三| 久久久亚洲精品成人影院| a级毛色黄片| 精品久久久久久电影网| 欧美激情国产日韩精品一区| 简卡轻食公司| 99久久中文字幕三级久久日本| 内射极品少妇av片p| 国产麻豆成人av免费视频| 在线免费观看的www视频| 少妇高潮的动态图| 久久99热6这里只有精品| 天堂√8在线中文| 岛国毛片在线播放| 在现免费观看毛片| 亚洲真实伦在线观看|