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

    如何正確運用χ2檢驗
    ——秩和檢驗與SAS實現(xiàn)

    2021-11-04 08:31:34胡純嚴(yán)胡良平
    四川精神衛(wèi)生 2021年5期
    關(guān)鍵詞:見式中位數(shù)線性

    胡純嚴(yán) ,胡良平 ,2*

    (1.軍事科學(xué)院研究生院,北京 100850;2.世界中醫(yī)藥學(xué)會聯(lián)合會臨床科研統(tǒng)計學(xué)專業(yè)委員會,北京 100029*通信作者:胡良平,E-mail:lphu927@163.com)

    在對單因素資料進(jìn)行差異性分析時,若發(fā)現(xiàn)資料不滿足參數(shù)檢驗的前提條件,宜選用適用面寬的非參數(shù)檢驗。其中,秩和檢驗可能是常被選用的方法之一。本文將介紹兩樣本資料秩和檢驗[1-6]、多樣本資料秩和檢驗[1-6]以及前述兩種情形下都可能會用到的10種評分方法[1]。

    1 兩樣本資料秩和檢驗

    1.1 簡單線性秩檢驗

    1.1.1 未分層資料的簡單線性秩檢驗

    未分層資料的簡單線性秩檢驗的檢驗統(tǒng)計量[1]見式(1):

    在式(1)中,z是一個服從標(biāo)準(zhǔn)正態(tài)分布的檢驗統(tǒng)計量(即隨機變量),S、E0(S)和Var0(S)分別代表未分層資料的“簡單線性秩統(tǒng)計量”“簡單線性秩統(tǒng)計量的期望值”和“簡單線性秩統(tǒng)計量的方差”,定義分別見式(2)、式(3)、式(4):

    在上述三式中,各變量的含義如下:Rj是第j個觀測或個體的秩;a(Rj)是基于第j個觀測的秩的得分(注意:評分的具體方法有多種,將在后文介紹);Cj是指示變量代表第j個觀測所在的組;n代表總觀測數(shù)(即總樣本含量);n1代表第1組(較小樣本含量)樣本含量;n2代表第2組樣本含量;是平均得分。

    基于標(biāo)準(zhǔn)正態(tài)分布理論和式(1)計算所得到的z值,可計算出與z對應(yīng)的正態(tài)分布曲線下尾端的概率P值,右單側(cè)概率、左單側(cè)概率和雙側(cè)概率分別見式(5)、式(6)、式(7):

    連續(xù)性校正:當(dāng)基于Wilcoxon和Siegel-Tukey評分且進(jìn)行漸近雙側(cè)檢驗時,SAS/STAT的NPAR1WAY過程將默認(rèn)進(jìn)行連續(xù)性校正,即當(dāng)式(1)分子的計算結(jié)果大于0時,分子減掉0.5;當(dāng)式(1)分子的計算結(jié)果小于0時,分子加上0.5。若想取消連續(xù)性校正,需要在“PROC NPAR1WAY”語句中增加選項“CORRECT=NO”。

    1.1.2 分層資料的簡單線性秩檢驗

    若資料中有一個分層因素(通常稱其為重要非試驗因素),設(shè)其有K(K>2)個水平,在分層因素的每個水平下,試驗因素均有兩個水平(即兩個對比組)。于是,分層資料的簡單線性秩檢驗的檢驗統(tǒng)計量[1]見式(8):

    式(8)中,z是一個服從標(biāo)準(zhǔn)正態(tài)分布的檢驗統(tǒng)計量(即隨機變量),T、E0(T)和Var0(T)分別代表分層資料的“簡單線性秩統(tǒng)計量”“簡單線性秩統(tǒng)計量的期望值”和“簡單線性秩統(tǒng)計量的方差”,其定義分別見式(9)、式(10)、式(11):

    在以上三式中,Sk為第k層的“簡單線性秩統(tǒng)計量”,wk為第k層的權(quán)重,其定義見式(12):

    在式(12)中,nk為第k層的樣本含量。如果在“STRATA”語句中,指定“WEIGHTS=STRATUM”,則wk=1/(nk+1 );如果指定“WEIGHTS=EQUAL”,則wk=1。

    1.2 Fligner-Policello檢驗

    1.2.1 概述

    Fligner和Policello于1981年提出的比較兩組定量資料中位數(shù)的檢驗方法[1],被稱為“Fligner-Policello檢驗法”。該法假定每組定量資料服從對稱分布,但不要求兩組定量資料具有相同的分布,也不要求兩組定量資料的方差相等。該法是基于Orban和Wolfe于1979年提出的“配置得分”而構(gòu)建。設(shè)有X與Y兩個組,來自X組的觀測Xi的配置得分記為P(Xi),其取值定義如下:P(Xi)=Y組中取值小于Xi的數(shù)據(jù)個數(shù);如果遇到相等的數(shù)值,需要對P(Xi)進(jìn)行校正,即在已知P(Xi)的基礎(chǔ)上增加Y組中取值等于Xi的數(shù)據(jù)個數(shù)的一半。對來自Y組的觀測Yj的配置得分記為P(Yj),其取值定義與P(Xi)相同。

    1.2.2 配置得分的定義

    在式(13)、式(14)中,nX和nY分別代表X組與Y組的樣本含量;I(·)是指示函數(shù)。于是,兩組各自的平均配置得分的計算公式分別見式(15)、式(16):

    1.2.3 Fligner-Policello檢驗統(tǒng)計量

    Fligner-Policello檢驗統(tǒng)計量見式(17):

    在式(17)中,z是一個服從標(biāo)準(zhǔn)正態(tài)分布的檢驗統(tǒng)計量(即隨機變量);VX和VY的計算分別見式(18)、式(19):

    X和Y與兩個組的配置得分的標(biāo)準(zhǔn)差分別見式(20)、式(21):

    【說明】P值的定義與式(5)、式(6)、式(7)相同,此處從略。

    2 多樣本資料秩和檢驗

    2.1 概述

    對多組定量資料進(jìn)行比較的秩和檢驗法常有下面兩個名稱,第一個叫做“單因素ANOVA檢驗”;第二個叫做“Kruskal-Wallis檢驗(采取Wilcoxon評分法)”。其實,它們本質(zhì)上都屬于“χ2檢驗”。當(dāng)對多組定量資料進(jìn)行整體比較時,其檢驗假設(shè)為“H0:各組之間沒有差別”。

    2.2 多組之間的整體比較

    設(shè)有一個具有r個水平的試驗因素,對定量資料進(jìn)行r組之間的整體比較時,所需要的檢驗統(tǒng)計量[1]見式(22):

    在式(22)中,C是一個服從自由度為df=r-1的χ2分布的檢驗統(tǒng)計量;ni是第i個水平組的樣本含量;Ti是第i個水平組的得分之和;E0(Ti)是在H0成立的條件下第i個水平組的期望秩和;S2是得分的樣本方差。Ti、E0(Ti)和S2的計算公式分別見式(23)、式(24)、式(25):

    2.3 多組之間的兩兩比較

    2.3.1 概述

    由Dwass、Steel、Critchlow和Fligner提出的多重比較方法,簡稱為“DSCF檢驗法”。此法從r(r>2)個組中每次抽取兩組進(jìn)行比較,故總共需要比較r×(r-1)/2次。每次比較都基于標(biāo)準(zhǔn)化的威爾科克森檢驗統(tǒng)計量,即采取威爾科克森法評分,并采用式(1)計算z統(tǒng)計量。

    2.3.2 DSCF檢驗統(tǒng)計量

    基于標(biāo)準(zhǔn)化的威爾科克森z檢驗統(tǒng)計量構(gòu)造出DSCF檢驗統(tǒng)計量見式(26):

    在式(26)中,z是采取威爾科克森法評分,并采用式(1)計算的結(jié)果(注意:每次比較只涉及兩組定量資料);而DSCF是一個近似服從于“r個標(biāo)準(zhǔn)正態(tài)變量的學(xué)生化極差分布”的檢驗統(tǒng)計量。兩樣本DSCF比較的P值可以通過下面的方法求出,即將DSCF統(tǒng)計量的值視為學(xué)生化極差分布的百分位數(shù),從而,基于學(xué)生化極差分布下特定百分位數(shù)計算出分布曲線下尾端的概率,即為所求的P值。

    3 秩和檢驗中的評分方法

    3.1 概述

    秩和檢驗的一個特點就是不直接利用原始數(shù)據(jù),而是先根據(jù)原始數(shù)據(jù)的大小給它們編秩。所謂編秩,就是給每個原始數(shù)據(jù)賦予一個自然數(shù),代表每個原始數(shù)據(jù)在一組和整個資料中的“相對位置”。然后再依據(jù)不同的數(shù)學(xué)原理,對每個“秩”進(jìn)行“評分”或“賦值”。SAS/STAT的NPAR1WAY過程[1]中介紹了十多種評分方法,現(xiàn)呈現(xiàn)其主要內(nèi)容。

    3.2 用于位置比較的評分方法

    3.2.1 威爾科克森(Wilcoxon)評分法

    威爾科克森評分是觀測的秩,可用公式表示如下:

    在式(27)中,Rj是第j個觀測的秩,而a(Rj)是第j個觀測的評分。

    【說明】在兩樣本資料的線性秩統(tǒng)計量中采用威爾科克森評分產(chǎn)生Mann-Whitney-Wilcoxon檢驗的秩和統(tǒng)計量;在多樣本資料的單因素ANOVA統(tǒng)計量中采用威爾科克森評分產(chǎn)生Kruskal-Wallis檢驗的秩和統(tǒng)計量;對于logistic分布的位置改變來說,威爾科克森評分是局部最有效能的。

    3.2.2 中位數(shù)(Median)評分法

    當(dāng)資料中的觀測值大于中位數(shù)時,則該觀測的中位數(shù)評分等于1;否則,中位數(shù)評分等于0。依據(jù)觀測的秩,中位數(shù)評分的定義見下式:

    【說明】在兩樣本資料的線性秩統(tǒng)計量中采用中位數(shù)評分產(chǎn)生兩樣本中位數(shù)檢驗的秩和統(tǒng)計量;在多樣本資料的單因素ANOVA統(tǒng)計量中采用中位數(shù)評分產(chǎn)生Brown-Mood檢驗的秩和統(tǒng)計量;中位數(shù)評分用于尾部抬高且對稱分布時,效能特別高。

    3.2.3 Van der Waerden(正態(tài))評分

    Van der Waerden評分是標(biāo)準(zhǔn)正態(tài)分布的分位數(shù),也被稱為分位數(shù)正態(tài)評分。該評分的計算公式見式(29):

    在式(29)中,Φ是標(biāo)準(zhǔn)正態(tài)分布的累計分布函數(shù)。對于正態(tài)分布而言,這些評分的效能極高。

    3.2.4 Savage評分

    Savage評分是來自指數(shù)分布的順序統(tǒng)計量的期望值,通過減掉1使評分的中心位于0附近。該評分的計算公式見式(30):

    Savage評分在以下兩種情形中具有高效能,其一,在指數(shù)分布中比較尺度差異;其二,在極值分布中比較位置變化。

    3.3 用于尺度比較的評分方法

    3.3.1 Siegel-Tukey評分

    Siegel-Tukey評分的定義如下:

    a(1)=1,a(n)=2,a(n-1)=3,a(2)=4

    a(3)=5,a(n-2)=6,a(n-3)=7,a(4)=8,…

    這里得分值按此模式朝著中間連續(xù)增加,直到全部觀測中的每個觀測都被賦予一個得分值為止。

    【說明】當(dāng)進(jìn)行Siegel-Tukey兩樣本檢驗的計算時,SAS/STAT中NPAR1WAY過程默認(rèn)需要進(jìn)行校正;如果用戶不想進(jìn)行校正,需要在“PROC NPAR1WAY”語句中增加選項“CORRECT=NO”。

    3.3.2 Ansari-Bradley評分

    Ansari-Bradley評分為對應(yīng)的極端秩賦予相同的得分,其定義如下:

    a(1)=1,a(n)=1,a(2)=2,a(n-1)=2

    a(3)=3,a(n-2)=3,a(4)=4,a(n-3)=4,…

    等價地,Ansari-Bradley評分可用如下通式表示:

    3.3.3 Klotz評分

    Klotz評分是Van der Waerden評分的平方,其定義如下:

    在式(32)中,Φ是標(biāo)準(zhǔn)正態(tài)分布的累計分布函數(shù)。

    3.3.4 Mood評分

    Mood評分按照觀測的秩與平均秩之差量的平方進(jìn)行計算,其定義如下:

    3.4 用于位置和尺度比較的評分方法

    Conove評分是基于觀測值與其樣本算術(shù)平均值之離差絕對值的秩的平方,對于第j個觀測而言,其定義如下:

    在式(34)中,Uj的計算見式(35):

    在式(35)中,i代表第i個樣本(組);j代表第i個樣本中第j個觀測;Xj(i)代表第i樣本中第j個觀測的觀測值;i代表第i個樣本的算術(shù)平均值;Uj代表第i個樣本中第j個觀測的秩。

    【說明】Conove于1999年提出,若在第i個樣本的全部Uj中出現(xiàn)了相同的數(shù)值(稱為“ties”),則先按無相同數(shù)據(jù)編秩(即給予編號),再求那幾個相同數(shù)據(jù)所對應(yīng)秩的算術(shù)平均值,并以此平均值作為它們的“秩”。Conove評分檢驗也被稱為“方差的平方秩檢驗”。

    4 實例與SAS實現(xiàn)

    4.1 問題與數(shù)據(jù)

    【例1】某地59例女性類風(fēng)濕性關(guān)節(jié)炎患者參加了一項臨床試驗[1],她們被隨機分配進(jìn)入試驗組(n=27)與安慰劑對照組(n=32)。結(jié)果變量有5種不同的取值,即“療效特好=5”“療效尚好=4”“療效中等=3”“療效一般=2”和“療效差=1”。記錄每位患者所接受的處理和療效的具體取值,臨床試驗結(jié)果以頻數(shù)表形式呈現(xiàn),詳見后面的SAS程序,此處從略。試對兩組有序資料進(jìn)行秩和檢驗,以評價兩種治療方法的效果差異是否有統(tǒng)計學(xué)意義。

    4.2 SAS實現(xiàn)

    4.2.1 對例1的SAS實現(xiàn)

    【分析與解答】設(shè)所需要的SAS程序如下:

    【程序說明】“從ab到fp”這10個選項是秩和檢驗中的10種評分方法;其中,“fp”是前文介紹的Fligner-Policello檢驗法。

    【SAS輸出結(jié)果及解釋】因篇幅所限,以下僅呈現(xiàn)“威爾科克森檢驗結(jié)果”,其他檢驗方法輸出的結(jié)果從略。

    以上是威爾科克森檢驗輸出的第1部分結(jié)果,即輸出兩組描述性統(tǒng)計量的計算結(jié)果,治療組的平均秩為37.00分,對照組的平均秩為24.09分。

    以上是威爾科克森檢驗輸出的第2部分結(jié)果,即以標(biāo)準(zhǔn)正態(tài)分布為理論根據(jù)計算得到的結(jié)果,雙側(cè)檢驗的P=0.0046<α=0.05。

    以上是威爾科克森檢驗輸出的第3部分結(jié)果,即以χ2分布為理論根據(jù)計算得到的結(jié)果,P=0.0031<α=0.05。

    【統(tǒng)計結(jié)論與專業(yè)結(jié)論】由輸出結(jié)果可知,前4種檢驗方法和第10種檢驗方法給出的檢驗結(jié)果均為P<0.05,說明治療組與對照組療效的“平均秩”或“中位數(shù)”之間差異有統(tǒng)計學(xué)意義;由于治療組的“平均秩”或“中位數(shù)”大于對照組的“平均秩”或“中位數(shù)”,又因評分值越大標(biāo)志著療效越好,故可以認(rèn)為“active”治療方法的效果優(yōu)于安慰劑。

    第5~8種檢驗方法的檢驗結(jié)果均為P>0.05,說明治療組與對照組療效的“離散度(即尺度參數(shù))”之間差異無統(tǒng)計學(xué)意義,即兩組有序資料的變化范圍接近一致。

    第9種檢驗方法的檢驗結(jié)果為P>0.05,說明治療組與對照組療效的“位置參數(shù)”和“離散度(即尺度參數(shù))”綜合的指標(biāo)之間差異無統(tǒng)計學(xué)意義。

    5 討論與小結(jié)

    5.1 討論

    秩和檢驗有兩個優(yōu)點:其一,對資料的要求不高;其二,選擇不同的評分方法可以分別實現(xiàn)“位置參數(shù)(如平均值、平均秩、中位數(shù))”“尺度參數(shù)(如標(biāo)準(zhǔn)差、分位數(shù)間距)”和“位置參數(shù)以及尺度參數(shù)”的比較。其缺點在于:不適合分析多因素資料。為了保留對資料要求不高的優(yōu)點,又能夠處理多因素資料,需要選擇復(fù)雜的非參數(shù)統(tǒng)計分析方法[6-8]。

    5.2 小結(jié)

    本文介紹了適用于分析單因素資料的秩和檢驗方法,包括分析單因素兩水平設(shè)計資料的“簡單線性秩檢驗”和單因素多水平設(shè)計資料的“單因素多水平ANOVA檢驗”。詳細(xì)介紹了在前述兩類檢驗中都不可缺少的10種評分方法。通過一個實例并借助SAS軟件,實現(xiàn)了單因素兩水平設(shè)計資料的簡單線性秩檢驗,呈現(xiàn)了10種評分方法計算所得到的結(jié)果,對輸出結(jié)果作出了解釋,并給出了統(tǒng)計結(jié)論和專業(yè)結(jié)論。

    猜你喜歡
    見式中位數(shù)線性
    高速公路下穿既有鐵路橋橋墩基底承載力驗算*
    公路與汽運(2024年1期)2024-03-07 03:02:06
    低溫下船用鋼材彈塑性曲線研究
    河南科技(2023年1期)2023-02-11 12:17:04
    漸近線性Klein-Gordon-Maxwell系統(tǒng)正解的存在性
    Effects of Landau damping and collision on stimulated Raman scattering with various phase-space distributions
    線性回歸方程的求解與應(yīng)用
    橋(門)式起重機起升機構(gòu)高速浮動軸設(shè)計
    二階線性微分方程的解法
    中位數(shù)計算公式及數(shù)學(xué)性質(zhì)的新認(rèn)識
    2015年中考數(shù)學(xué)模擬試題(五)
    2015年中考數(shù)學(xué)模擬試題(二)
    亚洲欧美日韩高清在线视频| 91精品三级在线观看| 操出白浆在线播放| 国产97色在线日韩免费| 精品久久久精品久久久| 桃红色精品国产亚洲av| 两个人看的免费小视频| 久久人妻福利社区极品人妻图片| 久久久国产精品麻豆| 国产成人精品无人区| 美女 人体艺术 gogo| 国产av又大| 夫妻午夜视频| 在线十欧美十亚洲十日本专区| 不卡一级毛片| 国产91精品成人一区二区三区| 精品免费久久久久久久清纯 | 美国免费a级毛片| 搡老熟女国产l中国老女人| 成人特级黄色片久久久久久久| 97人妻天天添夜夜摸| 老司机靠b影院| 91成人精品电影| 成年人免费黄色播放视频| 久久亚洲真实| 中文字幕人妻熟女乱码| 亚洲 欧美一区二区三区| 在线观看66精品国产| 新久久久久国产一级毛片| 欧美日本中文国产一区发布| 后天国语完整版免费观看| 欧美中文综合在线视频| 国产免费男女视频| 精品亚洲成国产av| 欧美激情高清一区二区三区| 欧美黄色淫秽网站| 成人永久免费在线观看视频| 午夜福利视频在线观看免费| 精品国产乱码久久久久久男人| 免费观看a级毛片全部| 国产人伦9x9x在线观看| 亚洲熟女毛片儿| 久久狼人影院| 久久亚洲真实| 国产高清videossex| 高清欧美精品videossex| 成人国产一区最新在线观看| 免费高清在线观看日韩| 很黄的视频免费| 99精品在免费线老司机午夜| 热re99久久精品国产66热6| 日韩 欧美 亚洲 中文字幕| 丝袜美腿诱惑在线| 久久久久视频综合| 99国产精品一区二区三区| 亚洲专区中文字幕在线| 精品久久久久久久毛片微露脸| 在线观看日韩欧美| 91精品国产国语对白视频| 超碰97精品在线观看| 黄色女人牲交| 中国美女看黄片| 亚洲综合色网址| netflix在线观看网站| 国产高清视频在线播放一区| 中文字幕人妻丝袜一区二区| 桃红色精品国产亚洲av| 精品久久久久久久久久免费视频 | 成年人午夜在线观看视频| 村上凉子中文字幕在线| 黄频高清免费视频| 在线av久久热| 岛国在线观看网站| 在线永久观看黄色视频| 国产精品偷伦视频观看了| 视频区欧美日本亚洲| 超碰成人久久| 亚洲av欧美aⅴ国产| 亚洲精品在线美女| 天天躁日日躁夜夜躁夜夜| 久久精品国产99精品国产亚洲性色 | svipshipincom国产片| 欧美色视频一区免费| 国产成人精品无人区| videosex国产| 日韩欧美一区二区三区在线观看 | 黄网站色视频无遮挡免费观看| 国产激情久久老熟女| 狠狠狠狠99中文字幕| 男人舔女人的私密视频| 免费在线观看日本一区| av片东京热男人的天堂| 丝袜人妻中文字幕| 91麻豆av在线| 又紧又爽又黄一区二区| 国产主播在线观看一区二区| 久久久久国内视频| 人妻丰满熟妇av一区二区三区 | 中文字幕另类日韩欧美亚洲嫩草| 精品午夜福利视频在线观看一区| 午夜激情av网站| av欧美777| 精品福利永久在线观看| 丰满迷人的少妇在线观看| 国产精品1区2区在线观看. | 国产亚洲欧美精品永久| 在线观看66精品国产| 国产蜜桃级精品一区二区三区 | 飞空精品影院首页| videosex国产| 777米奇影视久久| 99国产精品一区二区三区| 搡老岳熟女国产| 最近最新中文字幕大全免费视频| 伦理电影免费视频| 亚洲国产欧美日韩在线播放| 国产av精品麻豆| 99热网站在线观看| 激情在线观看视频在线高清 | 脱女人内裤的视频| 国产日韩欧美亚洲二区| 久久精品人人爽人人爽视色| 亚洲成av片中文字幕在线观看| 亚洲精品一卡2卡三卡4卡5卡| 两人在一起打扑克的视频| 成人亚洲精品一区在线观看| 亚洲精品一二三| 亚洲熟女毛片儿| 国产男女内射视频| 久久香蕉激情| 一级,二级,三级黄色视频| 国产一区二区三区在线臀色熟女 | 精品久久久久久电影网| 不卡av一区二区三区| 亚洲成人免费av在线播放| 狠狠狠狠99中文字幕| 欧美日韩成人在线一区二区| 欧美一级毛片孕妇| 女人爽到高潮嗷嗷叫在线视频| 久久香蕉精品热| 99国产精品免费福利视频| 99精品久久久久人妻精品| 亚洲中文日韩欧美视频| 国产淫语在线视频| 欧美日韩一级在线毛片| 国产精品国产高清国产av | 中文字幕高清在线视频| 黄色女人牲交| 老司机午夜十八禁免费视频| 久久精品亚洲熟妇少妇任你| 99精品欧美一区二区三区四区| 日本一区二区免费在线视频| www.自偷自拍.com| 久久久精品免费免费高清| 一级,二级,三级黄色视频| xxxhd国产人妻xxx| 国产一区在线观看成人免费| 亚洲人成电影免费在线| 亚洲专区字幕在线| 中文字幕人妻丝袜一区二区| 大码成人一级视频| 每晚都被弄得嗷嗷叫到高潮| 久久中文字幕人妻熟女| 黑人欧美特级aaaaaa片| 亚洲第一av免费看| 99精品久久久久人妻精品| 欧美人与性动交α欧美软件| 国产一区二区三区综合在线观看| 美女扒开内裤让男人捅视频| 伊人久久大香线蕉亚洲五| 国产精品av久久久久免费| 亚洲国产中文字幕在线视频| www.自偷自拍.com| 高清毛片免费观看视频网站 | 老司机亚洲免费影院| 欧美午夜高清在线| 久久精品国产99精品国产亚洲性色 | 久久久久国内视频| 美国免费a级毛片| 一进一出好大好爽视频| 久久久久久亚洲精品国产蜜桃av| 国产精品亚洲一级av第二区| 国内毛片毛片毛片毛片毛片| 午夜免费鲁丝| 亚洲中文字幕日韩| 99国产综合亚洲精品| 大码成人一级视频| 老熟妇仑乱视频hdxx| 久久精品人人爽人人爽视色| 亚洲精华国产精华精| 69精品国产乱码久久久| 黄色片一级片一级黄色片| 大码成人一级视频| 午夜久久久在线观看| 人妻丰满熟妇av一区二区三区 | 他把我摸到了高潮在线观看| 夜夜爽天天搞| 大码成人一级视频| 高清欧美精品videossex| 精品久久久久久,| 国产精品自产拍在线观看55亚洲 | 精品电影一区二区在线| a级片在线免费高清观看视频| 视频在线观看一区二区三区| 国产精品欧美亚洲77777| 手机成人av网站| 久久精品熟女亚洲av麻豆精品| 亚洲七黄色美女视频| 色尼玛亚洲综合影院| 中文亚洲av片在线观看爽 | 国产色视频综合| 黄网站色视频无遮挡免费观看| 99国产极品粉嫩在线观看| 国产一区二区三区在线臀色熟女 | 国产成人欧美在线观看 | 麻豆av在线久日| 久久影院123| 国产亚洲精品久久久久久毛片 | 午夜成年电影在线免费观看| 国产精品98久久久久久宅男小说| 老汉色av国产亚洲站长工具| 欧美日韩一级在线毛片| 在线十欧美十亚洲十日本专区| 久久中文看片网| 久久久国产成人免费| 天天躁日日躁夜夜躁夜夜| 日韩欧美国产一区二区入口| 国产精品自产拍在线观看55亚洲 | 欧美亚洲日本最大视频资源| 亚洲av美国av| 欧美精品一区二区免费开放| 久久热在线av| 建设人人有责人人尽责人人享有的| 午夜免费观看网址| 欧美日韩成人在线一区二区| xxx96com| 悠悠久久av| 女人久久www免费人成看片| 国产欧美日韩综合在线一区二区| 久久国产亚洲av麻豆专区| 国产在线一区二区三区精| 大陆偷拍与自拍| 亚洲欧美激情综合另类| 国产精品1区2区在线观看. | 啦啦啦在线免费观看视频4| 久久精品熟女亚洲av麻豆精品| bbb黄色大片| 变态另类成人亚洲欧美熟女 | 捣出白浆h1v1| 久久久久精品国产欧美久久久| 日本五十路高清| 香蕉国产在线看| 美女国产高潮福利片在线看| 国产有黄有色有爽视频| 怎么达到女性高潮| 在线免费观看的www视频| 久久精品国产99精品国产亚洲性色 | 中文字幕人妻丝袜制服| 精品久久久久久,| 国产欧美日韩精品亚洲av| 18在线观看网站| 久久人妻av系列| 亚洲午夜精品一区,二区,三区| 国产又爽黄色视频| 十分钟在线观看高清视频www| 国产99白浆流出| 亚洲国产精品一区二区三区在线| 欧美日韩国产mv在线观看视频| 亚洲熟妇中文字幕五十中出 | 久久久久久久久免费视频了| 国产欧美日韩一区二区三| 女人被狂操c到高潮| 午夜视频精品福利| 亚洲熟女精品中文字幕| 女同久久另类99精品国产91| 亚洲aⅴ乱码一区二区在线播放 | 日韩成人在线观看一区二区三区| 很黄的视频免费| 国产一区二区三区视频了| 母亲3免费完整高清在线观看| 午夜日韩欧美国产| 99国产精品免费福利视频| 免费av中文字幕在线| 精品亚洲成国产av| 免费在线观看亚洲国产| 丝袜在线中文字幕| tocl精华| 丰满饥渴人妻一区二区三| 天天操日日干夜夜撸| 午夜91福利影院| 日本黄色日本黄色录像| 人人澡人人妻人| 精品人妻在线不人妻| 久久久国产成人精品二区 | 久久久国产精品麻豆| 亚洲情色 制服丝袜| 久久精品熟女亚洲av麻豆精品| 国产亚洲欧美98| 少妇粗大呻吟视频| 国产熟女午夜一区二区三区| 美女视频免费永久观看网站| 香蕉国产在线看| 在线观看66精品国产| 亚洲七黄色美女视频| 午夜91福利影院| 19禁男女啪啪无遮挡网站| 黑人欧美特级aaaaaa片| 80岁老熟妇乱子伦牲交| 国产精品 欧美亚洲| 国产成人系列免费观看| 久久久国产成人精品二区 | 亚洲九九香蕉| 亚洲午夜理论影院| 午夜福利在线免费观看网站| 99久久99久久久精品蜜桃| 亚洲熟妇熟女久久| 国精品久久久久久国模美| 亚洲,欧美精品.| 在线天堂中文资源库| 久久国产精品人妻蜜桃| 久久精品国产清高在天天线| 亚洲欧洲精品一区二区精品久久久| 久久国产精品男人的天堂亚洲| 日本vs欧美在线观看视频| 少妇粗大呻吟视频| 操出白浆在线播放| 伊人久久大香线蕉亚洲五| 午夜精品久久久久久毛片777| 亚洲熟女毛片儿| 国产精品久久视频播放| 一级a爱视频在线免费观看| 极品教师在线免费播放| 欧美久久黑人一区二区| 精品国产超薄肉色丝袜足j| 岛国毛片在线播放| 搡老乐熟女国产| 欧美黑人精品巨大| 成人18禁在线播放| 久久精品国产清高在天天线| 久久人妻av系列| 国产欧美日韩一区二区三| 欧美日韩亚洲国产一区二区在线观看 | 黄色成人免费大全| 国产精品一区二区精品视频观看| 美女高潮喷水抽搐中文字幕| 色婷婷av一区二区三区视频| 亚洲精品成人av观看孕妇| 捣出白浆h1v1| 久久精品国产99精品国产亚洲性色 | 国产成人精品久久二区二区免费| 国产又爽黄色视频| 亚洲 国产 在线| 18禁美女被吸乳视频| 精品第一国产精品| 亚洲精品在线观看二区| 亚洲av日韩在线播放| 国产精品一区二区免费欧美| 黄网站色视频无遮挡免费观看| 啦啦啦 在线观看视频| 可以免费在线观看a视频的电影网站| 女性被躁到高潮视频| www.999成人在线观看| 日韩有码中文字幕| 欧美日韩av久久| 亚洲久久久国产精品| 国产亚洲欧美精品永久| 热re99久久精品国产66热6| 午夜福利在线免费观看网站| 色尼玛亚洲综合影院| 男人舔女人的私密视频| 亚洲av日韩在线播放| 亚洲色图综合在线观看| 天天影视国产精品| 香蕉国产在线看| 无限看片的www在线观看| 成年动漫av网址| 如日韩欧美国产精品一区二区三区| 欧美日韩国产mv在线观看视频| 少妇被粗大的猛进出69影院| 国产极品粉嫩免费观看在线| 色婷婷久久久亚洲欧美| 黄片小视频在线播放| 99热网站在线观看| 亚洲欧美日韩另类电影网站| 国产欧美亚洲国产| 美女扒开内裤让男人捅视频| 国产精品久久久久久精品古装| 一级a爱片免费观看的视频| tocl精华| 亚洲国产精品合色在线| 亚洲中文av在线| 黄色丝袜av网址大全| 国产精品久久久久久精品古装| 啪啪无遮挡十八禁网站| 久久久久久久久免费视频了| 国产又色又爽无遮挡免费看| 一进一出抽搐动态| 亚洲欧美日韩另类电影网站| 久久99一区二区三区| 午夜免费鲁丝| 国产一区在线观看成人免费| 久久久久久久久久久久大奶| 国产成人免费无遮挡视频| 午夜免费成人在线视频| 国产精品偷伦视频观看了| 男男h啪啪无遮挡| 久久天堂一区二区三区四区| 精品国产一区二区久久| 色老头精品视频在线观看| 99精品欧美一区二区三区四区| 久久久国产成人精品二区 | 国产精品乱码一区二三区的特点 | 国产精品98久久久久久宅男小说| 99精国产麻豆久久婷婷| 中文字幕高清在线视频| 黄频高清免费视频| 变态另类成人亚洲欧美熟女 | 人成视频在线观看免费观看| 久久香蕉精品热| 午夜精品在线福利| 乱人伦中国视频| 久久久国产成人精品二区 | av超薄肉色丝袜交足视频| 黑人巨大精品欧美一区二区mp4| 18禁裸乳无遮挡免费网站照片 | 亚洲avbb在线观看| 亚洲专区字幕在线| 天天添夜夜摸| 亚洲情色 制服丝袜| 日本五十路高清| 欧美日韩一级在线毛片| 飞空精品影院首页| 超色免费av| 美女国产高潮福利片在线看| 午夜亚洲福利在线播放| 看片在线看免费视频| 一a级毛片在线观看| 操出白浆在线播放| 自线自在国产av| 国产麻豆69| 欧美+亚洲+日韩+国产| 人人妻人人澡人人看| 女人被躁到高潮嗷嗷叫费观| 久99久视频精品免费| 亚洲三区欧美一区| 在线观看舔阴道视频| 国产亚洲av高清不卡| 性色av乱码一区二区三区2| 91在线观看av| 久久久久精品人妻al黑| 国产aⅴ精品一区二区三区波| 久久精品熟女亚洲av麻豆精品| 亚洲综合色网址| 国产在线观看jvid| 亚洲欧美激情综合另类| 国产高清视频在线播放一区| 18禁裸乳无遮挡免费网站照片 | 免费在线观看日本一区| 在线永久观看黄色视频| 国产亚洲欧美精品永久| 精品亚洲成a人片在线观看| 人人妻人人澡人人看| 国产高清视频在线播放一区| 99热网站在线观看| 午夜福利乱码中文字幕| 精品一区二区三区av网在线观看| 女人久久www免费人成看片| 中文字幕精品免费在线观看视频| 国产精品久久视频播放| 每晚都被弄得嗷嗷叫到高潮| 亚洲免费av在线视频| 国产精品香港三级国产av潘金莲| 国产精品一区二区免费欧美| 午夜福利在线免费观看网站| 一边摸一边抽搐一进一出视频| 欧美亚洲 丝袜 人妻 在线| 久久久久久久国产电影| 亚洲人成伊人成综合网2020| 国产真人三级小视频在线观看| 国产高清国产精品国产三级| 女性生殖器流出的白浆| 757午夜福利合集在线观看| 嫁个100分男人电影在线观看| 国产欧美日韩一区二区三区在线| 久久久久视频综合| 黄片小视频在线播放| 成人精品一区二区免费| 国产成人精品无人区| 999久久久精品免费观看国产| 99久久国产精品久久久| av在线播放免费不卡| 国产日韩一区二区三区精品不卡| 国产成人精品在线电影| 午夜激情av网站| 久久香蕉激情| 天天影视国产精品| 亚洲在线自拍视频| 极品教师在线免费播放| 麻豆成人av在线观看| bbb黄色大片| 国产又色又爽无遮挡免费看| 我的亚洲天堂| 精品国产乱码久久久久久男人| 少妇 在线观看| 久99久视频精品免费| 久久国产精品男人的天堂亚洲| 久久热在线av| 国内毛片毛片毛片毛片毛片| 亚洲欧美激情在线| 久久久精品国产亚洲av高清涩受| 欧美黑人精品巨大| 91精品国产国语对白视频| 亚洲欧美一区二区三区黑人| 建设人人有责人人尽责人人享有的| 精品久久久久久电影网| 国产精品国产av在线观看| 一级毛片精品| 欧美在线一区亚洲| 欧美日韩亚洲高清精品| 欧美一级毛片孕妇| 久久久久精品人妻al黑| 欧美精品啪啪一区二区三区| 久久九九热精品免费| 三上悠亚av全集在线观看| 成熟少妇高潮喷水视频| 国产欧美亚洲国产| 9热在线视频观看99| 久99久视频精品免费| 精品视频人人做人人爽| 亚洲在线自拍视频| 久久久水蜜桃国产精品网| 久久青草综合色| 国产在视频线精品| 男女床上黄色一级片免费看| 亚洲情色 制服丝袜| 欧美激情 高清一区二区三区| 法律面前人人平等表现在哪些方面| 在线视频色国产色| 久久精品国产99精品国产亚洲性色 | 少妇粗大呻吟视频| 久99久视频精品免费| 日韩免费av在线播放| 视频区图区小说| 国产精品亚洲一级av第二区| 久久性视频一级片| 国产成人欧美| 女警被强在线播放| 国产野战对白在线观看| 亚洲精品乱久久久久久| av欧美777| 成人影院久久| 亚洲av成人av| 老司机亚洲免费影院| 一级作爱视频免费观看| 久久精品熟女亚洲av麻豆精品| 日日摸夜夜添夜夜添小说| 母亲3免费完整高清在线观看| 91字幕亚洲| 黄色丝袜av网址大全| 黄色视频,在线免费观看| 欧美成人午夜精品| 国产亚洲欧美98| 男女床上黄色一级片免费看| 岛国在线观看网站| 视频区图区小说| 欧美午夜高清在线| 一本大道久久a久久精品| 如日韩欧美国产精品一区二区三区| 国产xxxxx性猛交| 国产成人免费无遮挡视频| 精品福利永久在线观看| av线在线观看网站| xxx96com| 最近最新中文字幕大全免费视频| 飞空精品影院首页| 亚洲成a人片在线一区二区| 日韩 欧美 亚洲 中文字幕| 久久热在线av| 亚洲av日韩在线播放| 飞空精品影院首页| 最新在线观看一区二区三区| 久久精品熟女亚洲av麻豆精品| 免费观看精品视频网站| a级毛片在线看网站| 国产男靠女视频免费网站| 我的亚洲天堂| 电影成人av| 国产精品综合久久久久久久免费 | 久久人妻福利社区极品人妻图片| 黄色毛片三级朝国网站| 丝袜美腿诱惑在线| 亚洲色图av天堂| 香蕉国产在线看| 婷婷精品国产亚洲av在线 | 国产不卡一卡二| 50天的宝宝边吃奶边哭怎么回事| 90打野战视频偷拍视频| 欧美 亚洲 国产 日韩一| 欧美在线一区亚洲| 精品一区二区三区视频在线观看免费 | 涩涩av久久男人的天堂| 熟女少妇亚洲综合色aaa.| 精品乱码久久久久久99久播| 亚洲成人手机| а√天堂www在线а√下载 | 一边摸一边抽搐一进一小说 | 人妻久久中文字幕网| 91成年电影在线观看| 国产99白浆流出| 亚洲国产中文字幕在线视频| 俄罗斯特黄特色一大片| 女人被躁到高潮嗷嗷叫费观| 欧美精品av麻豆av| 亚洲av日韩在线播放| 婷婷精品国产亚洲av在线 |