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

    基于Monte Carlo模擬的非參數(shù)多重比較方法評價

    2011-07-27 09:29:40濱州醫(yī)學(xué)院264003孫紅衛(wèi)韓春蕾
    中國衛(wèi)生統(tǒng)計 2011年5期
    關(guān)鍵詞:組數(shù)均數(shù)樣本量

    濱州醫(yī)學(xué)院(264003) 孫紅衛(wèi) 王 玖 韓春蕾

    在實(shí)際工作中,經(jīng)常要用Kruskal-Wallis檢驗(yàn)對多于2組的數(shù)據(jù)進(jìn)行比較,如果有統(tǒng)計學(xué)意義,需要采用非參數(shù)的多重比較方法(multiple comparison procedures,MCP)來鑒別到底哪兩組有差異。相比參數(shù)檢驗(yàn)的MCP,非參數(shù)MCP研究較少,提供的方法不多,且統(tǒng)計軟件沒有直接調(diào)用的過程可供實(shí)現(xiàn)。本文收集了在實(shí)際中常常采用的5種非參數(shù)MCP,用Monte-Carlo模擬來考察這些方法的第一、二類錯誤和判對率,并對實(shí)際中應(yīng)用MCP提出建議。

    非參數(shù)多重比較方法

    1.Dunn_z

    這是Dunn〔1〕在1964年提出的一種用于非參數(shù)多重比較的方法。第i組和第j組進(jìn)行比較的統(tǒng)計量為:

    z服從正態(tài)分布。要對檢驗(yàn)水準(zhǔn)進(jìn)行調(diào)整,α'=2α/k(k-1),z與za'比較。為兩對比組中第i組、第 j組的平均秩和;ni、nj為 i組、j組的樣本含量;N為總例數(shù)。

    2.擴(kuò)展的t檢驗(yàn)法

    擴(kuò)展的t檢驗(yàn)法是衛(wèi)生統(tǒng)計學(xué)第6版教材〔2〕中介紹的方法。第i組和第j組進(jìn)行比較的統(tǒng)計量為:

    t服從自由度為N-k的t分布。其中,Ri、Rj、ni、nj、N的意義同上,k為組數(shù),H為Kruskal-Wallis檢驗(yàn)所得的統(tǒng)計量。這種方法也需要編程實(shí)現(xiàn)。

    3.秩次轉(zhuǎn)換之后的LSD(R_LSD)、SNK(R_SNK)和Bonferroni(R_BON)

    Conover和 Iman〔3〕認(rèn)為,對于非參數(shù)多重比較的問題,可以對數(shù)據(jù)取秩次之后,再對秩次用參數(shù)多重比較的方法來解決。這也是實(shí)際工作中經(jīng)常采用的方法。本文考察了R_LSD,R_SNK和R_BON三種方法的性能。

    在以上五種方法中,經(jīng)過證明,R_LSD與擴(kuò)展的t檢驗(yàn)法是等價的。所以,本文主要考察了四種非參數(shù)多重比較方法,即Dunn_z,R_LSD,R_SNK和R_BON。

    評價指標(biāo)

    理想的多重比較方法不但控制第一類錯誤在一定的檢驗(yàn)水準(zhǔn)下,而且能夠盡量減少第二類錯誤,提高檢驗(yàn)效能。

    1.第一類錯誤〔4〕

    (1)CER(comparisonwise error rate),指在單個的假設(shè)檢驗(yàn)中犯第一類錯誤的概率大小

    (2)FWE(familywise error rate),為實(shí)驗(yàn)的第一類錯誤。在檢驗(yàn)集合的所有假設(shè)檢驗(yàn)中出現(xiàn)第一類錯誤的概率,即檢驗(yàn)集合中至少發(fā)生了一個CER錯誤的概率。這是多重比較方法需要控制的指標(biāo)。

    2.第二類錯誤〔5〕

    (1)CP(complete power),即完全效能,指在檢驗(yàn)集合的比較中發(fā)現(xiàn)所有的差異的能力。

    (2)PP(proportional power),即比例效能,指在一次兩兩比較中發(fā)現(xiàn)差異的能力。

    3.TM(true model)〔6〕

    判對率指對于檢驗(yàn)集合的統(tǒng)計推斷完全正確的概率,其同時考察了犯第一、二類錯誤的概率。

    模擬實(shí)驗(yàn)設(shè)計

    由于非參數(shù)多重比較主要用于不滿足參數(shù)檢驗(yàn)的條件時,比如不符從正態(tài)分布,或者雖然正態(tài)分布但是方差不齊時,所以考慮了兩種分布,一是對數(shù)正態(tài)分布,二是正態(tài)分布但是方差不齊。

    考慮到多重比較方法的性能受到多個因素的影響,比如樣本含量、組數(shù)等,還有就是均數(shù)間的差異、方差大小等。為了盡量模擬實(shí)際應(yīng)用情況,實(shí)驗(yàn)設(shè)計見表1、表2。

    表1 正態(tài)分布方差不齊時的實(shí)驗(yàn)參數(shù)設(shè)定

    表2 對數(shù)正態(tài)分布時的實(shí)驗(yàn)參數(shù)設(shè)定

    在不同組合下根據(jù)對應(yīng)的分布產(chǎn)生隨機(jī)數(shù),并用4種多重比較方法進(jìn)行檢驗(yàn),考察4種方法的一類錯誤,二類錯誤以及判對率。隨機(jī)數(shù)的產(chǎn)生在SAS8.2上實(shí)現(xiàn),4種多重比較方法及各種指標(biāo)的計算均利用SAS的宏功能。每種組合下均模擬1000次。結(jié)果的分析使用SAS8.2和SPSS13.0實(shí)現(xiàn)。

    結(jié) 果

    本文兩種分布皆考慮了總體均數(shù)全部相等,總體均數(shù)部分相等和總體均數(shù)全部不等三種情況。經(jīng)過分析發(fā)現(xiàn),正態(tài)分布方差不齊和對數(shù)正態(tài)分布的結(jié)果相似,所以兩種分布的結(jié)果合在一起分析。

    1.總體均數(shù)全部相等

    由于所有總體均數(shù)全部相等,不可能犯第二類錯誤。這里要考察的指標(biāo)只有FWE。

    從表3可以看出,Dunn_z和R_BON對FWE的控制在0.1以下,而R_SNK的FWE超過了0.1的部分占17.4%,且當(dāng)樣本量小于20時,其檢驗(yàn)水準(zhǔn)都控制在0.1以內(nèi)。而R_LSD的FWE失去了控制。

    表3 總體均數(shù)全部相等時FWE的模擬結(jié)果

    2.總體均數(shù)部分相等

    由于總體均數(shù)部分相等,這里要考察的指標(biāo)有FWE,CP,PP及 TM(結(jié)果見表4~6,圖1~4)。

    表4 總體均數(shù)部分相等時FWE的模擬結(jié)果

    表4和表3的結(jié)果相近,但是在部分相等時,四種方法的FWE都減小了。其中Dunn_z控制在FWE在0.1以下,而R_BON超過0.1的部分只有2.8%,R_SNK超過0.1的部分有23.5%,但當(dāng)樣本量小于20時,其FWE超過0.1的部分只有7.4%,而這時R_LSD的FWE超過0.1的部分有19.8%。

    從圖1可以看出,Dunn_z法FWE中位數(shù)隨樣本量的增長變化不大,R_BON次之,而R_SNK和R_LSD增長較大。

    圖1 總體均數(shù)部分相等時FWE的中位數(shù)隨樣本量的變化情況

    圖2 總體均數(shù)部分相等時FWE的中位數(shù)隨組數(shù)的變化情況

    從圖2可以看出,Dunn_z和R_BON的FWE在組數(shù)增大時變化不大,這由于兩種方法都根據(jù)組數(shù)來調(diào)整了檢驗(yàn)水準(zhǔn)。R_SNK的FWE變化不大,而R_LSD的FWE隨組數(shù)的增長而增大。

    表5 總體均數(shù)部分相等時CP和PP的模擬結(jié)果

    從表5可以看出,R_LSD和R_SNK的完全效能和比例效能接近,都較高,R_BON次之,Dunn_z的檢驗(yàn)效能最低。

    從圖3可以看出,四種方法的檢驗(yàn)效能都隨著樣本量的增長而增大,且之間的差距越來越小。當(dāng)樣本量大于40時,四種方法的檢驗(yàn)效能相近。

    圖4顯示四種方法的完全效能都隨著組數(shù)的增長而減小,但是減小的程度不一樣,Dunn_z的完全效能下降幅度最大 ,R_BON的下降幅度也較大,而R_SNK和R_LSD的完全效能下降幅度較小。

    圖3 總體均數(shù)部分相等時CP的中位數(shù)隨樣本量的變化情況

    圖4 總體均數(shù)部分相等時CP的中位數(shù)隨組數(shù)的變化情況

    表6 總體均數(shù)部分相等時TM的模擬結(jié)果

    判對率綜合了第一類錯誤和第二類錯誤,從判對率的角度看,其中最高的是R_BON和R_SNK。

    3.總體均數(shù)全部不等

    由于所有總體均數(shù)全部不等,不可能犯第一類錯誤。這里要考察的指標(biāo)只有CP和PP。

    表7 總體均數(shù)全部不等時CP和PP的模擬結(jié)果

    R_LSD和R_SNK的完全效能和比例效能接近,都較高,R_BON次之,Dunn_z最低。

    結(jié) 論

    選擇合適的檢驗(yàn)方法的原則是能將FWE控制在一定的水準(zhǔn)下,選擇檢驗(yàn)效能最高的方法。

    從上面的結(jié)果中可以看出,R_LSD、R_SNK、R_BON以及Dunn_z犯第一類錯誤的概率依次減小,但檢驗(yàn)效能也依次減小。

    R_LSD以控制CER為目的,所以它的FWE不受控制,會遠(yuǎn)大于檢驗(yàn)水準(zhǔn),而R_SNK在樣本量比較小時,基本控制在0.1內(nèi)。而兩者的檢驗(yàn)效能相差不大。所以綜合來看,R_SNK要優(yōu)于R_LSD。

    R_BON與Dunn_z都將FWE控制在0.1內(nèi),但是R_BON的檢驗(yàn)效能要比Dunn_z高出很多,所以,在兩者中R_BON較優(yōu)。

    R_BON與R_SNK比較,R_BON控制FWE的能力較高,但是檢驗(yàn)效能較低。從前面的結(jié)果可以看到,R_BON的FWE隨組數(shù)和樣本量變化比較平穩(wěn),但是檢驗(yàn)效能卻隨著組數(shù)的增大和樣本量的減小而減小,在組數(shù)較小或樣本量很大時,其檢驗(yàn)效能與R_SNK接近。所以當(dāng)組數(shù)較小時,可以選用R_BON,而在組數(shù)較大時,如果樣本量也較大,可以考慮用R_BON,而如果組數(shù)較大但樣本量不大時用R_BON,檢驗(yàn)效能會較低。這時可以考慮用R_SNK,其檢驗(yàn)效能較高,而在樣本量較小時,其FWE基本控制在0.1以內(nèi)。

    綜上,五種方法的性能如下:擴(kuò)展的t檢驗(yàn)法與R_LSD 等價,R_LSD、R_SNK、R_BON 以及 Dunn_z犯第一類錯誤的概率依次減小,但檢驗(yàn)效能也依次減小。四種方法受樣本量和組數(shù)的影響不同。R_BON和R_SNK是兩種較優(yōu)的方法。在實(shí)際中,當(dāng)組數(shù)較小時,或組數(shù)較大,且樣本量也較大時,可以選用R_BON;而如果組數(shù)較大,樣本量較小時,可以選用用R_SNK。

    1.Dunn O.Multiple comparisons using rank sums.Technometrics,1964,6:241-252.

    2.方積乾主編.衛(wèi)生統(tǒng)計學(xué).第6版.北京:人民衛(wèi)生出版社,2008:183-187.

    3.Conover W,Ronald L.Rank Transformations as a Bridge Between Parametric and Nonparametric Statistics.The American Statistician,1981,35(3):124-129.

    4.Westfall P.Multiple comparisons and multiple tests:using the SAS system.Cary,NC:SAS Institute,Inc,1999:16-21.

    5.Seaman M.New developments in pairwise multiple comparisons:Some powerful and practicable procedures.Psychological Bulletin,1991,110(3):577-586.

    6.Cribbie R.Pairwise multiple comparisons:A model comparison approach versus stepwise procedures.British Journal of Mathematical and Statistical Psychology,2003,55:27-39.

    猜你喜歡
    組數(shù)均數(shù)樣本量
    組數(shù)
    醫(yī)學(xué)研究中樣本量的選擇
    航空裝備測試性試驗(yàn)樣本量確定方法
    一類求不定方程正整數(shù)解的組數(shù)問題的解法及推廣
    Sample Size Calculations for Comparing Groups with Binary Outcomes
    關(guān)于均數(shù)與偏差
    關(guān)于均數(shù)與偏差
    關(guān)于均數(shù)與偏差
    關(guān)于均數(shù)與偏差
    自適應(yīng)樣本量調(diào)整中Fisher合并P值法和傳統(tǒng)檢驗(yàn)法的模擬比較
    国产白丝娇喘喷水9色精品| 99re6热这里在线精品视频| 国产在线一区二区三区精| 亚洲图色成人| 国产一级毛片在线| 97精品久久久久久久久久精品| 人妻少妇偷人精品九色| 亚洲av免费在线观看| 18+在线观看网站| 国产成人a区在线观看| 成人综合一区亚洲| 免费看不卡的av| 人人妻人人看人人澡| 精品人妻一区二区三区麻豆| 日韩欧美精品免费久久| 精品久久久久久电影网| 免费看光身美女| 国产91av在线免费观看| 丰满少妇做爰视频| 国产一级毛片在线| 日韩欧美精品免费久久| 亚洲四区av| 国产片特级美女逼逼视频| 日本欧美国产在线视频| 午夜激情欧美在线| 国产久久久一区二区三区| 国产在视频线精品| 免费观看的影片在线观看| 99久久精品国产国产毛片| 日本欧美国产在线视频| 婷婷色av中文字幕| av在线老鸭窝| 国产91av在线免费观看| 日本爱情动作片www.在线观看| 欧美3d第一页| 男人和女人高潮做爰伦理| 国产精品久久久久久久久免| 黄色一级大片看看| 久久久成人免费电影| 午夜亚洲福利在线播放| 国模一区二区三区四区视频| 青青草视频在线视频观看| 高清日韩中文字幕在线| 国产一区二区在线观看日韩| 91av网一区二区| 国产成人精品久久久久久| 国产成人精品福利久久| 麻豆成人av视频| 亚洲国产日韩欧美精品在线观看| 搞女人的毛片| 久久久午夜欧美精品| 日本wwww免费看| 国产一区二区在线观看日韩| 日本黄大片高清| 成人高潮视频无遮挡免费网站| 中文在线观看免费www的网站| 男女那种视频在线观看| 国产精品久久久久久久久免| 国产一区有黄有色的免费视频 | 精品久久久久久久人妻蜜臀av| 色5月婷婷丁香| 十八禁国产超污无遮挡网站| 一级av片app| 欧美一区二区亚洲| 99久久中文字幕三级久久日本| 久久久久久久久久黄片| 在线观看一区二区三区| 国产成年人精品一区二区| 性色avwww在线观看| 午夜福利成人在线免费观看| 日韩视频在线欧美| 国产欧美日韩精品一区二区| 国产大屁股一区二区在线视频| 亚洲一区高清亚洲精品| 国产片特级美女逼逼视频| 国产亚洲一区二区精品| 亚洲不卡免费看| 99久久精品一区二区三区| 午夜免费观看性视频| videossex国产| 国语对白做爰xxxⅹ性视频网站| 99热这里只有是精品在线观看| 成人亚洲精品av一区二区| 九九爱精品视频在线观看| 精品久久久久久久久av| 久久久久久久久久久免费av| 久久亚洲国产成人精品v| 看黄色毛片网站| 国产日韩欧美在线精品| 免费av观看视频| 午夜精品一区二区三区免费看| 久久这里有精品视频免费| 亚洲人成网站在线观看播放| 97人妻精品一区二区三区麻豆| xxx大片免费视频| 一夜夜www| av在线播放精品| 综合色av麻豆| 夫妻午夜视频| 亚洲av成人av| 国产欧美另类精品又又久久亚洲欧美| 国产精品人妻久久久影院| 亚洲av一区综合| 国产单亲对白刺激| 赤兔流量卡办理| 三级经典国产精品| 国产精品99久久久久久久久| 国产成人精品一,二区| 精品久久久久久久久久久久久| 最近最新中文字幕大全电影3| 丝袜喷水一区| 尾随美女入室| 免费观看a级毛片全部| 一区二区三区四区激情视频| 九九爱精品视频在线观看| 看黄色毛片网站| 国产精品一区二区在线观看99 | 特大巨黑吊av在线直播| 男人爽女人下面视频在线观看| 97精品久久久久久久久久精品| 国产黄色视频一区二区在线观看| 欧美日韩亚洲高清精品| 国产精品一区二区在线观看99 | 少妇熟女欧美另类| 国产精品一区二区在线观看99 | 乱系列少妇在线播放| 欧美成人午夜免费资源| 成年av动漫网址| 三级毛片av免费| 国模一区二区三区四区视频| 亚洲国产最新在线播放| 国产在视频线精品| 久久久久久久午夜电影| 亚洲欧洲日产国产| 少妇人妻精品综合一区二区| 国产大屁股一区二区在线视频| 亚洲丝袜综合中文字幕| 国产一区有黄有色的免费视频 | 天堂√8在线中文| 亚洲国产欧美在线一区| 精品一区在线观看国产| 好男人视频免费观看在线| 亚洲精品成人av观看孕妇| 干丝袜人妻中文字幕| 国产精品熟女久久久久浪| 午夜亚洲福利在线播放| 久久久国产一区二区| 天堂√8在线中文| 视频中文字幕在线观看| 亚洲婷婷狠狠爱综合网| 视频中文字幕在线观看| 亚洲精品日韩在线中文字幕| 欧美性感艳星| 亚洲综合色惰| 3wmmmm亚洲av在线观看| 免费看a级黄色片| 大陆偷拍与自拍| 欧美成人a在线观看| 成人午夜精彩视频在线观看| 国产伦理片在线播放av一区| 禁无遮挡网站| 久久久国产一区二区| 国产黄色免费在线视频| 天天一区二区日本电影三级| 黄色一级大片看看| 亚洲成人av在线免费| 国产女主播在线喷水免费视频网站 | 久久99热6这里只有精品| 欧美极品一区二区三区四区| 男人狂女人下面高潮的视频| a级毛色黄片| 狂野欧美白嫩少妇大欣赏| 亚洲av不卡在线观看| 午夜激情福利司机影院| 99久国产av精品| 99久久九九国产精品国产免费| 九九久久精品国产亚洲av麻豆| 国产精品久久视频播放| 国产欧美日韩精品一区二区| 亚洲av成人精品一二三区| 综合色丁香网| 2021少妇久久久久久久久久久| 毛片女人毛片| 亚洲婷婷狠狠爱综合网| 国精品久久久久久国模美| 嫩草影院新地址| 天美传媒精品一区二区| 久久精品国产亚洲网站| 久久鲁丝午夜福利片| 亚洲美女视频黄频| 亚洲aⅴ乱码一区二区在线播放| 亚洲成人久久爱视频| 亚洲图色成人| 国产91av在线免费观看| 国产淫语在线视频| 高清视频免费观看一区二区 | 我要看日韩黄色一级片| 三级毛片av免费| 久久久久国产网址| 色综合站精品国产| 国产亚洲av片在线观看秒播厂 | 禁无遮挡网站| 熟妇人妻久久中文字幕3abv| 欧美bdsm另类| 69av精品久久久久久| 国产视频内射| 国产麻豆成人av免费视频| 男女视频在线观看网站免费| 亚洲欧洲日产国产| 激情五月婷婷亚洲| 男女国产视频网站| 男女国产视频网站| 最近中文字幕高清免费大全6| 少妇熟女欧美另类| 十八禁网站网址无遮挡 | 街头女战士在线观看网站| 黄片wwwwww| 2022亚洲国产成人精品| 免费大片黄手机在线观看| 97超碰精品成人国产| 天堂网av新在线| 亚洲人成网站高清观看| 国产精品麻豆人妻色哟哟久久 | 少妇人妻一区二区三区视频| 国产一区二区亚洲精品在线观看| 欧美日韩视频高清一区二区三区二| 深爱激情五月婷婷| 日韩强制内射视频| 国产精品爽爽va在线观看网站| 美女大奶头视频| 国产成人免费观看mmmm| 超碰av人人做人人爽久久| 国语对白做爰xxxⅹ性视频网站| 中文字幕免费在线视频6| 三级经典国产精品| 黄片wwwwww| 日韩成人av中文字幕在线观看| 亚洲成人一二三区av| 日本-黄色视频高清免费观看| 亚洲国产欧美人成| 国产精品一二三区在线看| 国产伦精品一区二区三区四那| 人妻夜夜爽99麻豆av| 男女边摸边吃奶| 婷婷色麻豆天堂久久| av免费在线看不卡| 一二三四中文在线观看免费高清| 午夜福利成人在线免费观看| 特级一级黄色大片| 欧美日韩在线观看h| 日日摸夜夜添夜夜爱| 99久久精品一区二区三区| 亚洲国产av新网站| 亚洲四区av| 亚洲在线观看片| 成人亚洲精品一区在线观看 | 中文字幕av成人在线电影| 国产精品一区二区三区四区久久| 一本久久精品| av在线老鸭窝| 22中文网久久字幕| 毛片女人毛片| 亚洲最大成人中文| 国产成人91sexporn| 最后的刺客免费高清国语| 亚洲精品国产av蜜桃| 夜夜爽夜夜爽视频| 国产伦精品一区二区三区视频9| 色播亚洲综合网| 青春草国产在线视频| 国产伦精品一区二区三区四那| 夜夜看夜夜爽夜夜摸| 亚洲av免费在线观看| 最近2019中文字幕mv第一页| 深夜a级毛片| 美女国产视频在线观看| 国产av国产精品国产| 国产精品国产三级国产专区5o| 99热这里只有是精品50| 熟女电影av网| 国产一区亚洲一区在线观看| 春色校园在线视频观看| 尤物成人国产欧美一区二区三区| 色网站视频免费| 最近视频中文字幕2019在线8| 午夜精品国产一区二区电影 | 精品久久国产蜜桃| 久久韩国三级中文字幕| 欧美xxxx性猛交bbbb| 欧美成人午夜免费资源| 99re6热这里在线精品视频| 神马国产精品三级电影在线观看| 免费看美女性在线毛片视频| 精品国内亚洲2022精品成人| 国产精品人妻久久久影院| 18禁动态无遮挡网站| 国产熟女欧美一区二区| 中文精品一卡2卡3卡4更新| 日本一二三区视频观看| 国产视频内射| 欧美变态另类bdsm刘玥| 亚洲电影在线观看av| 亚洲av成人av| 日韩 亚洲 欧美在线| 久久久久久久久中文| 午夜日本视频在线| 五月天丁香电影| 欧美区成人在线视频| 搡老妇女老女人老熟妇| 中文字幕久久专区| 自拍偷自拍亚洲精品老妇| 一夜夜www| 一本久久精品| 韩国av在线不卡| 可以在线观看毛片的网站| 成年女人看的毛片在线观看| 日韩欧美三级三区| 91久久精品国产一区二区成人| 中文字幕制服av| 肉色欧美久久久久久久蜜桃 | 国产精品久久久久久精品电影小说 | 久久久久久久久中文| 九草在线视频观看| 欧美精品一区二区大全| 日本猛色少妇xxxxx猛交久久| 亚洲精品第二区| 久久这里有精品视频免费| 国产av在哪里看| 亚洲精品日本国产第一区| 最近中文字幕2019免费版| 久久6这里有精品| 国产成人午夜福利电影在线观看| a级毛片免费高清观看在线播放| 高清在线视频一区二区三区| ponron亚洲| 超碰av人人做人人爽久久| 免费看av在线观看网站| 大又大粗又爽又黄少妇毛片口| 精品不卡国产一区二区三区| 99热网站在线观看| 久久国内精品自在自线图片| 国产又色又爽无遮挡免| 只有这里有精品99| 日日干狠狠操夜夜爽| 又黄又爽又刺激的免费视频.| 在线观看一区二区三区| 五月天丁香电影| 亚洲欧美日韩东京热| 亚洲乱码一区二区免费版| 男的添女的下面高潮视频| 国精品久久久久久国模美| 女的被弄到高潮叫床怎么办| 国产白丝娇喘喷水9色精品| 日本猛色少妇xxxxx猛交久久| 婷婷色av中文字幕| 99视频精品全部免费 在线| 亚洲最大成人av| 久久久a久久爽久久v久久| 成人亚洲欧美一区二区av| av在线播放精品| 亚洲精品一二三| 免费观看a级毛片全部| 久久久久精品性色| 国产高清不卡午夜福利| 九色成人免费人妻av| 亚洲美女搞黄在线观看| 日本与韩国留学比较| 亚洲综合精品二区| 麻豆成人av视频| 国产精品一区二区三区四区久久| 91精品国产九色| 纵有疾风起免费观看全集完整版 | 国产精品嫩草影院av在线观看| 不卡视频在线观看欧美| 久久人人爽人人爽人人片va| 中文字幕免费在线视频6| 中文字幕av成人在线电影| 亚洲高清免费不卡视频| 国产精品不卡视频一区二区| 亚州av有码| 久久久久久久大尺度免费视频| 麻豆久久精品国产亚洲av| 免费观看性生交大片5| 联通29元200g的流量卡| 久久久久精品久久久久真实原创| 国产精品久久久久久av不卡| 内地一区二区视频在线| 黄色日韩在线| 免费电影在线观看免费观看| 国产淫片久久久久久久久| 看黄色毛片网站| 亚洲精品自拍成人| 国产精品一区二区三区四区久久| 麻豆成人午夜福利视频| 午夜久久久久精精品| 熟女电影av网| 最近最新中文字幕免费大全7| av卡一久久| 禁无遮挡网站| 亚洲av男天堂| 亚洲高清免费不卡视频| 亚洲av在线观看美女高潮| 久久精品人妻少妇| 精品久久久久久久久亚洲| 久久99热这里只频精品6学生| a级毛色黄片| 少妇熟女aⅴ在线视频| 国产精品人妻久久久影院| 日本黄色片子视频| 麻豆久久精品国产亚洲av| 国产黄片视频在线免费观看| 国产成人a区在线观看| 亚洲精品久久久久久婷婷小说| 综合色丁香网| 国产精品麻豆人妻色哟哟久久 | 国产亚洲精品久久久com| 欧美+日韩+精品| 五月玫瑰六月丁香| 中文乱码字字幕精品一区二区三区 | 亚洲精品第二区| 真实男女啪啪啪动态图| 免费观看av网站的网址| 日韩成人av中文字幕在线观看| 国产在线一区二区三区精| 国产精品国产三级专区第一集| 午夜精品在线福利| 亚洲国产高清在线一区二区三| 精品一区二区三区视频在线| 黄色一级大片看看| 欧美一区二区亚洲| 联通29元200g的流量卡| 精华霜和精华液先用哪个| 亚洲高清免费不卡视频| 国产高清三级在线| 91久久精品电影网| 国产黄片视频在线免费观看| 免费观看精品视频网站| 精品国产三级普通话版| 国产人妻一区二区三区在| 免费无遮挡裸体视频| xxx大片免费视频| 高清毛片免费看| 亚洲真实伦在线观看| 国精品久久久久久国模美| 精品99又大又爽又粗少妇毛片| 国产白丝娇喘喷水9色精品| xxx大片免费视频| 国产一级毛片在线| 国产午夜福利久久久久久| 久久久精品欧美日韩精品| 国产成人精品一,二区| 国产成人免费观看mmmm| 女的被弄到高潮叫床怎么办| 亚洲真实伦在线观看| 丰满少妇做爰视频| 成年女人在线观看亚洲视频 | 精品少妇黑人巨大在线播放| 亚洲av男天堂| 欧美成人一区二区免费高清观看| 亚洲国产精品国产精品| 黄色日韩在线| 久久99热这里只有精品18| 亚洲四区av| 高清在线视频一区二区三区| 午夜视频国产福利| 国产男女超爽视频在线观看| 婷婷色综合大香蕉| 久久99热这里只频精品6学生| 午夜爱爱视频在线播放| 国产久久久一区二区三区| 亚洲久久久久久中文字幕| 超碰av人人做人人爽久久| 欧美人与善性xxx| 免费看美女性在线毛片视频| 人人妻人人看人人澡| 国内揄拍国产精品人妻在线| 亚洲精品中文字幕在线视频 | 最近的中文字幕免费完整| 亚洲欧美日韩东京热| 99九九线精品视频在线观看视频| 大又大粗又爽又黄少妇毛片口| 亚洲欧美成人综合另类久久久| 亚洲成人一二三区av| 日本与韩国留学比较| 超碰97精品在线观看| 久久99热这里只有精品18| 欧美97在线视频| 国产精品人妻久久久影院| 国产精品麻豆人妻色哟哟久久 | 亚洲国产成人一精品久久久| 神马国产精品三级电影在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲人成网站在线观看播放| 高清日韩中文字幕在线| 国产一级毛片在线| ponron亚洲| 一本一本综合久久| 欧美区成人在线视频| 校园人妻丝袜中文字幕| 亚洲性久久影院| 最近视频中文字幕2019在线8| kizo精华| 久久久精品欧美日韩精品| 亚洲精品一二三| 国产亚洲最大av| 少妇的逼水好多| 久久久久久久久久黄片| 成人亚洲精品av一区二区| 2022亚洲国产成人精品| 成人av在线播放网站| 久久久久久久久久人人人人人人| av在线天堂中文字幕| 亚洲丝袜综合中文字幕| 亚洲av免费在线观看| 亚洲av在线观看美女高潮| 久久久久精品久久久久真实原创| 国产精品蜜桃在线观看| 国产成人免费观看mmmm| 97超碰精品成人国产| 亚洲高清免费不卡视频| 日韩av在线免费看完整版不卡| 日产精品乱码卡一卡2卡三| 熟女人妻精品中文字幕| 纵有疾风起免费观看全集完整版 | 又爽又黄a免费视频| 久久这里只有精品中国| 国产精品国产三级国产专区5o| 日本-黄色视频高清免费观看| 六月丁香七月| 精品酒店卫生间| 搡老妇女老女人老熟妇| 欧美日本视频| 我要看日韩黄色一级片| 看黄色毛片网站| 久久久久久久午夜电影| 国产又色又爽无遮挡免| 搡老乐熟女国产| 国产成人一区二区在线| 国产在视频线精品| 2021天堂中文幕一二区在线观| 又爽又黄a免费视频| 淫秽高清视频在线观看| 久久亚洲国产成人精品v| 欧美性感艳星| 国产午夜精品论理片| 欧美日韩视频高清一区二区三区二| 久久久精品免费免费高清| 丰满人妻一区二区三区视频av| 中文字幕免费在线视频6| 日韩欧美精品v在线| 美女高潮的动态| 亚洲精品久久久久久婷婷小说| 97精品久久久久久久久久精品| 男女边摸边吃奶| 乱系列少妇在线播放| 日韩精品青青久久久久久| 建设人人有责人人尽责人人享有的 | 国产精品久久久久久精品电影| 亚洲精品国产av蜜桃| 亚洲乱码一区二区免费版| 国产乱来视频区| 男人舔奶头视频| 麻豆乱淫一区二区| 日产精品乱码卡一卡2卡三| 午夜激情福利司机影院| 免费看a级黄色片| 狂野欧美白嫩少妇大欣赏| 久久精品人妻少妇| 午夜激情福利司机影院| 国产精品福利在线免费观看| 国产单亲对白刺激| 国产一区二区在线观看日韩| 嘟嘟电影网在线观看| 26uuu在线亚洲综合色| 精品久久久久久成人av| av免费观看日本| 特大巨黑吊av在线直播| 五月伊人婷婷丁香| 特大巨黑吊av在线直播| 少妇丰满av| 中文字幕制服av| 99久久精品热视频| 极品少妇高潮喷水抽搐| 欧美精品国产亚洲| 人人妻人人看人人澡| 校园人妻丝袜中文字幕| 白带黄色成豆腐渣| 成人午夜高清在线视频| 一个人看的www免费观看视频| 亚洲高清免费不卡视频| 欧美最新免费一区二区三区| 一区二区三区免费毛片| 成人高潮视频无遮挡免费网站| 少妇猛男粗大的猛烈进出视频 | 免费观看精品视频网站| 久久久久久久午夜电影| av在线亚洲专区| 久久久色成人| 国产亚洲5aaaaa淫片| 亚洲av不卡在线观看| 国产成人免费观看mmmm| 精品一区二区三卡| 少妇的逼好多水| 日韩在线高清观看一区二区三区| 99九九线精品视频在线观看视频| 尾随美女入室| 女人十人毛片免费观看3o分钟| 三级男女做爰猛烈吃奶摸视频| 18禁裸乳无遮挡免费网站照片| 麻豆久久精品国产亚洲av| 韩国av在线不卡| 3wmmmm亚洲av在线观看| 又爽又黄a免费视频| 国产亚洲91精品色在线| 中文资源天堂在线|