• <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é)模擬試題(二)
    亚洲国产精品国产精品| 国产黄频视频在线观看| 久久久精品94久久精品| 在线 av 中文字幕| 亚洲精品一区蜜桃| 久久久久久久久久久免费av| 一区二区av电影网| 中文字幕av电影在线播放| 秋霞在线观看毛片| 成年女人在线观看亚洲视频| 精品国产一区二区三区久久久樱花| 波多野结衣一区麻豆| 成人影院久久| 国产成人aa在线观看| 一个人免费看片子| 亚洲第一区二区三区不卡| 看非洲黑人一级黄片| 黑人欧美特级aaaaaa片| 99热6这里只有精品| 97精品久久久久久久久久精品| 老熟女久久久| 97在线人人人人妻| 欧美激情极品国产一区二区三区 | 妹子高潮喷水视频| 久久久久精品人妻al黑| 看非洲黑人一级黄片| 午夜免费鲁丝| 香蕉精品网在线| 亚洲,欧美精品.| 亚洲一码二码三码区别大吗| 国产男女内射视频| 狂野欧美激情性bbbbbb| 国产午夜精品一二区理论片| 国产极品天堂在线| 亚洲精品国产色婷婷电影| 亚洲国产精品成人久久小说| 一区二区日韩欧美中文字幕 | 国产av国产精品国产| 韩国高清视频一区二区三区| 不卡视频在线观看欧美| 日韩伦理黄色片| 欧美性感艳星| 亚洲成人一二三区av| 国产精品不卡视频一区二区| 91国产中文字幕| 国产片内射在线| 亚洲欧美中文字幕日韩二区| 免费看不卡的av| 久久久久久久久久久久大奶| 狂野欧美激情性xxxx在线观看| 下体分泌物呈黄色| 97超碰精品成人国产| av网站免费在线观看视频| 少妇的丰满在线观看| 久久精品国产综合久久久 | 国产男女内射视频| 国产 精品1| 精品国产一区二区三区四区第35| 男女边摸边吃奶| 国产精品欧美亚洲77777| 国产精品嫩草影院av在线观看| 国产亚洲av片在线观看秒播厂| 国产淫语在线视频| 成人毛片a级毛片在线播放| 国产精品人妻久久久影院| 最后的刺客免费高清国语| 亚洲国产精品国产精品| 国产福利在线免费观看视频| 激情五月婷婷亚洲| 欧美成人精品欧美一级黄| 欧美人与善性xxx| 91国产中文字幕| 综合色丁香网| 久久精品国产亚洲av天美| 久久精品国产亚洲av涩爱| 日本av免费视频播放| 男女高潮啪啪啪动态图| 国产成人欧美| av视频免费观看在线观看| 欧美+日韩+精品| 老司机影院毛片| 国产免费福利视频在线观看| 欧美bdsm另类| 亚洲内射少妇av| 国产高清国产精品国产三级| 9191精品国产免费久久| 欧美日韩av久久| 青春草国产在线视频| 美女视频免费永久观看网站| 日日啪夜夜爽| 亚洲精品美女久久久久99蜜臀 | 免费观看a级毛片全部| 一区二区av电影网| 激情五月婷婷亚洲| av电影中文网址| 成年人午夜在线观看视频| 久久人人爽人人片av| 国产一级毛片在线| 精品一品国产午夜福利视频| 色94色欧美一区二区| 欧美日韩一区二区视频在线观看视频在线| 国产黄频视频在线观看| 欧美精品国产亚洲| 久久 成人 亚洲| 另类精品久久| 老司机影院毛片| 成人国语在线视频| 日韩 亚洲 欧美在线| 精品人妻在线不人妻| 亚洲精品自拍成人| 人人澡人人妻人| 蜜桃国产av成人99| 插逼视频在线观看| 又粗又硬又长又爽又黄的视频| 中文乱码字字幕精品一区二区三区| 秋霞伦理黄片| 国产 一区精品| 啦啦啦中文免费视频观看日本| 欧美 日韩 精品 国产| 国产欧美亚洲国产| 亚洲精品乱久久久久久| 免费日韩欧美在线观看| 日韩成人伦理影院| 亚洲av成人精品一二三区| 亚洲国产最新在线播放| 少妇的逼水好多| 日产精品乱码卡一卡2卡三| 美女xxoo啪啪120秒动态图| 热re99久久精品国产66热6| 国产在线视频一区二区| 亚洲av.av天堂| 99久国产av精品国产电影| 青春草视频在线免费观看| 亚洲精品视频女| 两个人看的免费小视频| 国产乱人偷精品视频| 永久网站在线| 视频中文字幕在线观看| 精品亚洲成a人片在线观看| 欧美日韩一区二区视频在线观看视频在线| av天堂久久9| 久热久热在线精品观看| 中文字幕人妻丝袜制服| 欧美性感艳星| 老女人水多毛片| 交换朋友夫妻互换小说| 国产综合精华液| 在线亚洲精品国产二区图片欧美| 在线天堂最新版资源| 久久人人爽人人爽人人片va| 亚洲av综合色区一区| 看免费成人av毛片| 国产成人91sexporn| 欧美人与性动交α欧美软件 | 狂野欧美激情性xxxx在线观看| 蜜桃国产av成人99| 日本免费在线观看一区| 欧美性感艳星| 国产成人精品一,二区| 男女啪啪激烈高潮av片| 日本av手机在线免费观看| 日韩制服丝袜自拍偷拍| 精品人妻偷拍中文字幕| 日本av手机在线免费观看| 99久久人妻综合| 人体艺术视频欧美日本| 在线观看三级黄色| 女的被弄到高潮叫床怎么办| 午夜福利网站1000一区二区三区| 人人澡人人妻人| 亚洲图色成人| 国产成人精品无人区| 中文字幕最新亚洲高清| 亚洲av在线观看美女高潮| 天堂中文最新版在线下载| √禁漫天堂资源中文www| 少妇熟女欧美另类| 成人亚洲欧美一区二区av| 中文精品一卡2卡3卡4更新| 欧美成人午夜精品| 久久国产精品大桥未久av| 中文精品一卡2卡3卡4更新| 久久人人爽人人片av| 久久鲁丝午夜福利片| 亚洲av.av天堂| 丝袜脚勾引网站| 久久av网站| 国产麻豆69| 色婷婷久久久亚洲欧美| 久久av网站| 精品国产一区二区久久| 少妇 在线观看| 免费看光身美女| 99热网站在线观看| 麻豆精品久久久久久蜜桃| 国产色爽女视频免费观看| 日本黄色日本黄色录像| 又黄又爽又刺激的免费视频.| 亚洲精品美女久久av网站| 亚洲av欧美aⅴ国产| 亚洲欧美成人综合另类久久久| 久久久久人妻精品一区果冻| 少妇 在线观看| 国产免费又黄又爽又色| 男女免费视频国产| 日韩av免费高清视频| 9191精品国产免费久久| 男女午夜视频在线观看 | 久久99热这里只频精品6学生| 黄色视频在线播放观看不卡| 精品人妻在线不人妻| 亚洲美女黄色视频免费看| 毛片一级片免费看久久久久| 午夜福利,免费看| 亚洲国产av影院在线观看| 人妻少妇偷人精品九色| 久久女婷五月综合色啪小说| 少妇的逼好多水| 久久久久久久国产电影| 日韩一本色道免费dvd| 亚洲精品av麻豆狂野| 国产精品久久久av美女十八| 26uuu在线亚洲综合色| 精品一区二区三区视频在线| 看十八女毛片水多多多| 观看美女的网站| 精品亚洲成国产av| 观看av在线不卡| 五月天丁香电影| 欧美人与善性xxx| 日韩中文字幕视频在线看片| 国产69精品久久久久777片| 成人亚洲欧美一区二区av| 国产日韩一区二区三区精品不卡| 亚洲av在线观看美女高潮| 2022亚洲国产成人精品| 国内精品宾馆在线| 成人亚洲精品一区在线观看| 人人澡人人妻人| 又黄又爽又刺激的免费视频.| 香蕉国产在线看| 亚洲四区av| 岛国毛片在线播放| 国产一区亚洲一区在线观看| 久久久精品区二区三区| 免费观看av网站的网址| 亚洲精品av麻豆狂野| 80岁老熟妇乱子伦牲交| www.色视频.com| 麻豆乱淫一区二区| 满18在线观看网站| 一区二区三区四区激情视频| 亚洲欧美成人精品一区二区| 777米奇影视久久| 天天躁夜夜躁狠狠躁躁| 国产极品天堂在线| 9热在线视频观看99| 久久久久久久久久久免费av| 观看美女的网站| 亚洲丝袜综合中文字幕| 日本vs欧美在线观看视频| 水蜜桃什么品种好| 制服丝袜香蕉在线| 一级毛片我不卡| 国产伦理片在线播放av一区| 多毛熟女@视频| 精品亚洲成国产av| 少妇人妻 视频| 一边摸一边做爽爽视频免费| 国产成人精品久久久久久| 亚洲 欧美一区二区三区| 国产免费又黄又爽又色| av片东京热男人的天堂| 中文字幕免费在线视频6| 日韩av在线免费看完整版不卡| 欧美日韩国产mv在线观看视频| 亚洲,欧美精品.| 国产国拍精品亚洲av在线观看| 18禁动态无遮挡网站| 伊人亚洲综合成人网| 精品卡一卡二卡四卡免费| 久久亚洲国产成人精品v| 国产成人精品婷婷| 丝袜喷水一区| 精品一区二区免费观看| 三上悠亚av全集在线观看| 国产精品国产三级专区第一集| 一区二区三区乱码不卡18| 午夜福利,免费看| 一级片'在线观看视频| 中文字幕免费在线视频6| 欧美日韩成人在线一区二区| 曰老女人黄片| 永久免费av网站大全| 国产1区2区3区精品| 国产伦理片在线播放av一区| 91久久精品国产一区二区三区| 亚洲少妇的诱惑av| 亚洲综合色惰| 亚洲精品国产av蜜桃| 亚洲精品久久久久久婷婷小说| 夜夜骑夜夜射夜夜干| 国产一区二区三区综合在线观看 | 久久久久久久久久成人| 男女国产视频网站| 一区二区日韩欧美中文字幕 | 久热久热在线精品观看| 丝袜人妻中文字幕| 99热全是精品| 亚洲精品自拍成人| 少妇人妻精品综合一区二区| 精品熟女少妇av免费看| 少妇猛男粗大的猛烈进出视频| 国产精品国产av在线观看| 色吧在线观看| 亚洲图色成人| 看非洲黑人一级黄片| 热re99久久精品国产66热6| 成人毛片60女人毛片免费| 国产成人免费无遮挡视频| 国产精品不卡视频一区二区| 国产亚洲午夜精品一区二区久久| 日韩一区二区视频免费看| 亚洲国产精品专区欧美| 欧美丝袜亚洲另类| 欧美亚洲 丝袜 人妻 在线| www日本在线高清视频| 免费高清在线观看日韩| 久久国产亚洲av麻豆专区| 深夜精品福利| 中国国产av一级| 80岁老熟妇乱子伦牲交| 国产精品三级大全| 人妻 亚洲 视频| 各种免费的搞黄视频| 97在线人人人人妻| 国产成人91sexporn| 九九爱精品视频在线观看| 91aial.com中文字幕在线观看| 韩国av在线不卡| 免费观看在线日韩| 国产精品久久久久成人av| 久久久国产精品麻豆| 免费观看无遮挡的男女| 国产在线免费精品| 毛片一级片免费看久久久久| 国产精品.久久久| 亚洲,一卡二卡三卡| 午夜福利,免费看| 青春草亚洲视频在线观看| 日韩一区二区视频免费看| 国产国拍精品亚洲av在线观看| 国产色爽女视频免费观看| av片东京热男人的天堂| videos熟女内射| 日日撸夜夜添| 老女人水多毛片| 一级毛片 在线播放| 久久久久久人妻| 国产欧美日韩综合在线一区二区| 中文字幕免费在线视频6| 美女主播在线视频| 纯流量卡能插随身wifi吗| 最黄视频免费看| 999精品在线视频| 免费观看a级毛片全部| 天堂俺去俺来也www色官网| 亚洲av国产av综合av卡| 国产乱来视频区| 亚洲激情五月婷婷啪啪| 亚洲精品成人av观看孕妇| 亚洲av中文av极速乱| 日韩 亚洲 欧美在线| 久久免费观看电影| 九色亚洲精品在线播放| 男女下面插进去视频免费观看 | 满18在线观看网站| 亚洲欧美中文字幕日韩二区| 视频在线观看一区二区三区| 国产乱来视频区| 免费人妻精品一区二区三区视频| 欧美成人午夜精品| 哪个播放器可以免费观看大片| 免费看av在线观看网站| 美女主播在线视频| 免费黄网站久久成人精品| 日韩不卡一区二区三区视频在线| 99久久人妻综合| 精品一区二区三区四区五区乱码 | 欧美亚洲 丝袜 人妻 在线| 侵犯人妻中文字幕一二三四区| 亚洲av欧美aⅴ国产| 色婷婷久久久亚洲欧美| av网站免费在线观看视频| 九色亚洲精品在线播放| freevideosex欧美| 性色av一级| 蜜臀久久99精品久久宅男| 欧美精品亚洲一区二区| 丝袜人妻中文字幕| 亚洲成人一二三区av| 国产一区二区三区av在线| 天堂俺去俺来也www色官网| 日产精品乱码卡一卡2卡三| 人妻系列 视频| 亚洲一区二区三区欧美精品| 韩国av在线不卡| 免费人妻精品一区二区三区视频| 性色av一级| 制服诱惑二区| 熟妇人妻不卡中文字幕| 99热网站在线观看| 在线观看一区二区三区激情| 亚洲综合色网址| 日韩电影二区| 欧美精品av麻豆av| 亚洲,欧美精品.| 综合色丁香网| 成人免费观看视频高清| 成人无遮挡网站| 黑丝袜美女国产一区| 人妻人人澡人人爽人人| 99视频精品全部免费 在线| 婷婷色综合www| 欧美xxⅹ黑人| 欧美xxxx性猛交bbbb| 丝袜人妻中文字幕| 激情视频va一区二区三区| 欧美日韩视频精品一区| 免费黄网站久久成人精品| 99热网站在线观看| 国产xxxxx性猛交| 男女高潮啪啪啪动态图| 欧美+日韩+精品| 国产1区2区3区精品| 五月开心婷婷网| av卡一久久| 精品国产露脸久久av麻豆| av电影中文网址| 肉色欧美久久久久久久蜜桃| 免费日韩欧美在线观看| 欧美精品人与动牲交sv欧美| 国产成人aa在线观看| 18在线观看网站| 妹子高潮喷水视频| 国产精品无大码| 青春草国产在线视频| 欧美日韩视频精品一区| av线在线观看网站| 日韩三级伦理在线观看| av视频免费观看在线观看| 成人漫画全彩无遮挡| 在线精品无人区一区二区三| 亚洲,一卡二卡三卡| 免费av不卡在线播放| 一区二区三区乱码不卡18| 寂寞人妻少妇视频99o| 老司机亚洲免费影院| 亚洲av欧美aⅴ国产| 麻豆乱淫一区二区| 精品国产乱码久久久久久小说| 九九爱精品视频在线观看| 久久久久视频综合| 午夜免费观看性视频| 久久人人爽av亚洲精品天堂| 久久久久精品人妻al黑| 国产1区2区3区精品| 美女脱内裤让男人舔精品视频| 美女国产高潮福利片在线看| 美女大奶头黄色视频| 99热这里只有是精品在线观看| 男的添女的下面高潮视频| 久久久久国产精品人妻一区二区| 久久久精品94久久精品| 国产一区有黄有色的免费视频| 欧美国产精品一级二级三级| 亚洲精品美女久久久久99蜜臀 | 寂寞人妻少妇视频99o| 51国产日韩欧美| 如日韩欧美国产精品一区二区三区| 免费黄网站久久成人精品| 国产精品国产三级国产av玫瑰| 午夜福利乱码中文字幕| av国产精品久久久久影院| 99精国产麻豆久久婷婷| av一本久久久久| 国产日韩欧美在线精品| 黑人欧美特级aaaaaa片| 亚洲久久久国产精品| 嫩草影院入口| 寂寞人妻少妇视频99o| 在线免费观看不下载黄p国产| 国产一区二区在线观看日韩| 国产色婷婷99| 在线观看国产h片| 国产免费现黄频在线看| av女优亚洲男人天堂| 免费人妻精品一区二区三区视频| 丝瓜视频免费看黄片| 九色成人免费人妻av| a级片在线免费高清观看视频| 亚洲精品久久午夜乱码| 日本免费在线观看一区| 久久人人爽人人片av| 成人国产麻豆网| 成人综合一区亚洲| 校园人妻丝袜中文字幕| 国产精品久久久久久久电影| 亚洲精品成人av观看孕妇| 中文字幕人妻熟女乱码| 男人添女人高潮全过程视频| 国产成人av激情在线播放| 大陆偷拍与自拍| 人妻少妇偷人精品九色| 婷婷色综合大香蕉| 青春草视频在线免费观看| 久久人人爽av亚洲精品天堂| 中文字幕另类日韩欧美亚洲嫩草| 侵犯人妻中文字幕一二三四区| 亚洲欧洲国产日韩| 精品一品国产午夜福利视频| 一二三四中文在线观看免费高清| 天堂中文最新版在线下载| 在线观看www视频免费| 日韩不卡一区二区三区视频在线| 精品人妻在线不人妻| 国精品久久久久久国模美| 成年美女黄网站色视频大全免费| 免费看av在线观看网站| 中国三级夫妇交换| 人妻一区二区av| 午夜老司机福利剧场| 汤姆久久久久久久影院中文字幕| 乱人伦中国视频| 伦理电影免费视频| 久久女婷五月综合色啪小说| 国产成人aa在线观看| 成人国产麻豆网| 日韩精品免费视频一区二区三区 | 免费在线观看完整版高清| 美女内射精品一级片tv| 一级,二级,三级黄色视频| 日韩制服丝袜自拍偷拍| 色哟哟·www| 久久久国产精品麻豆| 丝袜脚勾引网站| 成人国产麻豆网| 秋霞伦理黄片| 99久久综合免费| 日韩av不卡免费在线播放| 男人舔女人的私密视频| 精品午夜福利在线看| 中国美白少妇内射xxxbb| 一级毛片我不卡| 日产精品乱码卡一卡2卡三| 91成人精品电影| av天堂久久9| 国产精品久久久久久久电影| 亚洲精品中文字幕在线视频| 中国三级夫妇交换| 久久精品国产亚洲av涩爱| 2021少妇久久久久久久久久久| 亚洲美女黄色视频免费看| 一个人免费看片子| 日本爱情动作片www.在线观看| 久久99一区二区三区| 飞空精品影院首页| 国产精品三级大全| 久久久国产精品麻豆| 久久人人爽人人片av| 少妇人妻久久综合中文| 日韩在线高清观看一区二区三区| 人妻人人澡人人爽人人| 街头女战士在线观看网站| 国产成人精品一,二区| 久久久久久久精品精品| 中文天堂在线官网| av线在线观看网站| 搡老乐熟女国产| 免费女性裸体啪啪无遮挡网站| 亚洲av综合色区一区| 一边亲一边摸免费视频| 日韩,欧美,国产一区二区三区| 在线 av 中文字幕| 精品国产一区二区久久| 97超碰精品成人国产| 国产亚洲精品久久久com| 亚洲天堂av无毛| 日韩电影二区| 午夜福利视频精品| 一区二区三区精品91| 精品亚洲成国产av| 免费观看性生交大片5| 母亲3免费完整高清在线观看 | 亚洲第一av免费看| 亚洲 欧美一区二区三区| 母亲3免费完整高清在线观看 | av视频免费观看在线观看| 国产免费一区二区三区四区乱码| 成人手机av| 精品一区二区三卡| 18禁在线无遮挡免费观看视频| 这个男人来自地球电影免费观看 | 久久久久久人妻| 亚洲精品国产色婷婷电影| 免费观看在线日韩| 亚洲国产成人一精品久久久| 伦理电影免费视频| 亚洲欧美一区二区三区黑人 | 女性生殖器流出的白浆| 少妇人妻精品综合一区二区| 亚洲精品乱久久久久久| 久久鲁丝午夜福利片| 亚洲欧美日韩另类电影网站| 亚洲一级一片aⅴ在线观看|