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

    方差不齊時(shí)兩組及以上均數(shù)比較時(shí)不同分析方法的穩(wěn)健性和把握度比較

    2016-12-27 08:49:10復(fù)旦大學(xué)公共衛(wèi)生學(xué)院生物統(tǒng)計(jì)教研室衛(wèi)生部衛(wèi)生技術(shù)評(píng)估重點(diǎn)實(shí)驗(yàn)室
    關(guān)鍵詞:穩(wěn)健性均數(shù)樣本量

    復(fù)旦大學(xué)公共衛(wèi)生學(xué)院生物統(tǒng)計(jì)教研室,衛(wèi)生部衛(wèi)生技術(shù)評(píng)估重點(diǎn)實(shí)驗(yàn)室

    蔣 丹 趙耐青

    方差不齊時(shí)兩組及以上均數(shù)比較時(shí)不同分析方法的穩(wěn)健性和把握度比較

    復(fù)旦大學(xué)公共衛(wèi)生學(xué)院生物統(tǒng)計(jì)教研室,衛(wèi)生部衛(wèi)生技術(shù)評(píng)估重點(diǎn)實(shí)驗(yàn)室

    蔣 丹 趙耐青

    目的 本文著重比較秩和檢驗(yàn)、調(diào)整自由度的t′檢驗(yàn)、混合效應(yīng)模型(mixed model)以及方差加權(quán)最小二乘法(VWLS)等方法在方差不齊時(shí),用于兩組/多組獨(dú)立樣本均數(shù)比較時(shí)的穩(wěn)健性和把握度。方法 本文通過(guò)模擬分析方法,分別設(shè)計(jì)總體均數(shù)相等或不等時(shí),在不同標(biāo)準(zhǔn)差和樣本量的條件下,用幾種統(tǒng)計(jì)方法比較2組及3組樣本均數(shù)的Ⅰ類(lèi)錯(cuò)誤和Power度。結(jié)果 (1)證實(shí)樣本量相等時(shí),t檢驗(yàn)對(duì)于方差不齊的2組樣本均數(shù)比較具有穩(wěn)健性,但是樣本量相等方差不齊的3組獨(dú)立樣本均數(shù)比較時(shí),方差分析方法卻不具有穩(wěn)健性。(2)不論是2組還是3組樣本均數(shù)比較,秩和檢驗(yàn)在特定條件下對(duì)于方差不齊具有穩(wěn)健性。(3)兩組方差不齊樣本均數(shù)比較時(shí),t′檢驗(yàn)和mixed model因?yàn)棰耦?lèi)錯(cuò)誤更穩(wěn)健,比VWLS方法更穩(wěn)定,且這三種方法的Power值相互比較接近。(4)三組方差不齊樣本均數(shù)比較,mixed model方法在樣本量較少時(shí)比VWLS方法Ⅰ類(lèi)錯(cuò)誤更穩(wěn)健,但是隨著樣本量增加,這一優(yōu)勢(shì)消失,而VWLS的Power值明顯高于mixed model統(tǒng)計(jì)方法。結(jié)論 2組方差不齊樣本均數(shù)比較時(shí),可以使用t′檢驗(yàn)、mixed model及VWLS等方法,其中首選更為穩(wěn)健的t′檢驗(yàn)、mixed model,3組方差不齊樣本均數(shù)比較時(shí)可以使用mixed model及VWLS等方法,當(dāng)樣本量較小時(shí)首選mixed model方法,樣本量增大時(shí),以VWLS方法更優(yōu)。

    方差不齊 均數(shù)比較 第Ⅰ類(lèi)錯(cuò)誤 把握度

    比較兩組總體均數(shù)是否相同,通常使用t檢驗(yàn)進(jìn)行統(tǒng)計(jì)推斷,比較三組及以上總體均數(shù)是否相同,則可使用ANOVA方差分析。然而,不論t檢驗(yàn)還是ANOVA,均要求數(shù)據(jù)滿(mǎn)足獨(dú)立性、正態(tài)性、方差齊性等前提假設(shè)。而數(shù)據(jù)是否具有獨(dú)立性在研究設(shè)計(jì)時(shí)已經(jīng)明確;大數(shù)定律則告訴我們:在樣本量足夠的情況下不論數(shù)據(jù)本身是否服從正態(tài)分布,樣本均數(shù)的抽樣分布都可以逼近正態(tài)分布[1]。因此,在樣本量較大情況下,偏態(tài)分布的資料也可以用t檢驗(yàn)或方差分析。所以,對(duì)于t檢驗(yàn)或方差分析,樣本資料是否滿(mǎn)足方差齊性就顯得尤為重要。

    倘若方差齊性檢驗(yàn)的統(tǒng)計(jì)學(xué)結(jié)論為方差不齊,則傳統(tǒng)醫(yī)學(xué)統(tǒng)計(jì)教材往往會(huì)給出如下建議:2組總體均數(shù)比較使用調(diào)整自由度的t′檢驗(yàn)或者非參數(shù)統(tǒng)計(jì)分析方法如Wilcoxon秩和檢驗(yàn),3組及以上總體均數(shù)比較使用非參數(shù)統(tǒng)計(jì)分析方法如Kruskal-Wallis檢驗(yàn)。然而,Wilcoxon秩和檢驗(yàn)和Kruskal-Wallis檢驗(yàn)的無(wú)效假設(shè)為:兩個(gè)或多個(gè)樣本來(lái)自同一總體[2]。當(dāng)出現(xiàn)方差不齊情況時(shí),也就是對(duì)應(yīng)的已經(jīng)不是同一總體而是兩個(gè)總體或多個(gè)總體了。此時(shí),如果還用上述秩和檢驗(yàn)就存在明顯的邏輯問(wèn)題,當(dāng)然不排除在特殊情況下,方差不齊對(duì)于秩和檢驗(yàn)也可能是穩(wěn)健的(robust)。正如2組總體均數(shù)比較的情況下,2組樣本量相同,t檢驗(yàn)對(duì)于方差不齊具有穩(wěn)健性[3](即2組樣本量相同情況下,可以忽略方差不齊情況,t檢驗(yàn)仍然有效)。從線性模型的角度考慮,混合效應(yīng)模型(mixed model)和方差加權(quán)最小二乘法(variance-weighted least squares,VWLS)也可以用于方差不齊的2組或多組均數(shù)比較。綜上所述,本文研究目的有二:一是探究方差不齊時(shí),使用以下4種方法比較2組獨(dú)立樣本均數(shù)的穩(wěn)健性和把握度:調(diào)整自由度的t′檢驗(yàn)、Wilcoxon秩和檢驗(yàn)、混合效應(yīng)模型(mixed model)以及方差加權(quán)最小二乘法(VWLS);二是探究方差不齊時(shí),3組及以上獨(dú)立樣本均數(shù)比較用以下3種分析方法的穩(wěn)健性和把握度:混合效應(yīng)模型(mixed model)、Kruskal-Wallis檢驗(yàn)、方差加權(quán)最小二乘法(VWLS)。

    研究方法

    根據(jù)研究目的,有必要對(duì)方差不齊程度、樣本量之比歸類(lèi)來(lái)觀察不同條件下Ⅰ類(lèi)錯(cuò)誤和把握度(power)的大小變化情況,從而對(duì)各種方法的穩(wěn)健性和把握度進(jìn)行比較分析,因而需要在不同的情況下進(jìn)行模擬分析,下面對(duì)各種模擬情境作簡(jiǎn)單介紹。

    A、成組獨(dú)立樣本均數(shù)比較的模擬情境

    A1:設(shè)兩組總體均數(shù)相等為25,總體標(biāo)準(zhǔn)差之比為 1、1.5、2、2.5、3、4,使樣本量之比為 1、2、3、4,在不同情況下產(chǎn)生兩組服從正態(tài)分布的獨(dú)立樣本,觀察Ⅰ類(lèi)錯(cuò)誤大小。

    A2:設(shè)兩組總體均數(shù)不等,第1組為24,第2組為25,總體標(biāo)準(zhǔn)差之比和樣本量之比的考慮同上,在不同情況下產(chǎn)生均服從正態(tài)分布的兩個(gè)獨(dú)立樣本,比較power值大小。

    B、三組獨(dú)立樣本均數(shù)比較的模擬情境

    B1:設(shè)三組總體均數(shù)相等均為25,考慮σ1=σ2<σ3和σ1<σ2<σ32種情形,并考慮3組樣本量相等、3組中有2組樣本量相等、3組樣本量互不相等的三種情況,在不同情形下產(chǎn)生均服從正態(tài)分布的3個(gè)樣本,觀察Ⅰ類(lèi)錯(cuò)誤大小的變化。

    B2:設(shè)三組總體均數(shù)互不相等,分別為24.5、25、25.5,總體標(biāo)準(zhǔn)差和樣本量的考慮同上,得到不同情況下均服從正態(tài)分布的3個(gè)樣本,比較power值大小。

    上述4種模擬情境,在每種情形下獨(dú)立模擬2000次,將Ⅰ類(lèi)錯(cuò)誤α設(shè)為0.05;使用SAS9.3和stata11.0統(tǒng)計(jì)軟件進(jìn)行模擬分析。

    結(jié) 果

    1.兩組獨(dú)立樣本連續(xù)型定量資料比較

    ——使用 t′檢驗(yàn)、VWLS、mixed model和秩和檢驗(yàn)統(tǒng)計(jì)方法

    (1)觀察Ⅰ類(lèi)錯(cuò)誤的變化

    由表1可知 t′檢驗(yàn)、VWLS、mixed model這三種方法始終表現(xiàn)出穩(wěn)健性,Ⅰ類(lèi)錯(cuò)誤大小隨機(jī)波動(dòng),與標(biāo)準(zhǔn)差之比、樣本量之比無(wú)任何關(guān)聯(lián);但是 VWLS統(tǒng)計(jì)方法分析時(shí)Ⅰ類(lèi)錯(cuò)誤的發(fā)生比例高于其他2種方法;此外,表1中t檢驗(yàn)Ⅰ類(lèi)錯(cuò)誤大小的穩(wěn)定性,證實(shí)了文獻(xiàn)中提到的當(dāng)樣本量相等時(shí),t檢驗(yàn)對(duì)方差不齊的穩(wěn)健性;因?yàn)闃颖玖肯嗟葧r(shí)t檢驗(yàn)是mixed model的特例,所以由模擬結(jié)果可知兩種方法的Ⅰ類(lèi)錯(cuò)誤大小完全一致。

    表1 兩組樣本量和總體標(biāo)準(zhǔn)差不同比例下幾種方法Ⅰ類(lèi)錯(cuò)誤比較

    雖然秩和檢驗(yàn)在樣本量相等時(shí)Ⅰ類(lèi)錯(cuò)誤會(huì)隨著標(biāo)準(zhǔn)差之比的擴(kuò)大而增加,但是表1結(jié)果顯示秩和檢驗(yàn)的Ⅰ類(lèi)錯(cuò)誤大小在特定條件下增幅不十分顯著,比如總體標(biāo)準(zhǔn)差之比為2時(shí),Ⅰ類(lèi)錯(cuò)誤發(fā)生的比例增大5%;如果樣本量之比與標(biāo)準(zhǔn)差之比大小變化方向一致,那么隨著兩組樣本量或總體標(biāo)準(zhǔn)差之間差距拉大,Ⅰ類(lèi)錯(cuò)誤呈遞減趨勢(shì)(見(jiàn)表1上半部分),相反地,樣本量之比與總體標(biāo)準(zhǔn)差之比大小變化方向相反,那么隨著差距拉大,Ⅰ類(lèi)錯(cuò)誤呈遞增趨勢(shì)(見(jiàn)表1),其變化幅度與樣本量之比成正比。所以當(dāng)樣本量相等且標(biāo)準(zhǔn)差之比不超過(guò)2時(shí),秩和檢驗(yàn)具有穩(wěn)健性。

    (2)觀察power值的變化

    不論樣本量之比和總體標(biāo)準(zhǔn)差之比大小如何變化,t′檢驗(yàn)、VWLS、mixed model統(tǒng)計(jì)方法的 power值相互都比較接近;秩和檢驗(yàn)當(dāng)樣本量相等和標(biāo)準(zhǔn)差之比不超過(guò)2時(shí),power值同樣靠近其他方法。因?yàn)闃?biāo)準(zhǔn)差之比大于2后秩和檢驗(yàn)不具有穩(wěn)健性,所以power值不被納入比較范疇。

    表2 兩組樣本量和總體標(biāo)準(zhǔn)差不同比例下幾種方法power值比較

    2.三組獨(dú)立樣本連續(xù)型定量資料比較

    ——使用秩和檢驗(yàn)\mixed model\VWLS統(tǒng)計(jì)方法

    (1)觀察Ⅰ類(lèi)錯(cuò)誤變化(3組標(biāo)準(zhǔn)差存在形式如下:σ1=σ2<σ3\σ1<σ2<σ3)

    當(dāng)樣本量相等時(shí),表3和表4的結(jié)果說(shuō)明發(fā)現(xiàn)方差分析、秩和檢驗(yàn)的Ⅰ類(lèi)錯(cuò)誤隨著標(biāo)準(zhǔn)差之比的增加而增加;mixed model方法和VWLS方法的Ⅰ類(lèi)錯(cuò)誤大小穩(wěn)定,且其變化隨機(jī)波動(dòng),與標(biāo)準(zhǔn)差之比無(wú)任何關(guān)聯(lián);但是在表3和表4中,當(dāng)每組樣本量為30時(shí),VWLS方法Ⅰ類(lèi)錯(cuò)誤的發(fā)生比例明顯高于mixed model方法,尤其當(dāng)標(biāo)準(zhǔn)差形式是 σ1<σ2<σ3,由此說(shuō)明當(dāng)樣本量較小時(shí),mixed model方法比VWLS方法對(duì)于方差不齊更為穩(wěn)健。

    遇到樣本量不等時(shí),秩和檢驗(yàn)Ⅰ類(lèi)錯(cuò)誤大體的變化趨勢(shì)如下:當(dāng)樣本量之比和標(biāo)準(zhǔn)差之比大小變化方向一致時(shí),Ⅰ類(lèi)錯(cuò)誤隨著樣本量或者標(biāo)準(zhǔn)差之間差距增大呈現(xiàn)降低趨勢(shì);當(dāng)樣本量之比與標(biāo)準(zhǔn)差之比大小變化方向相反時(shí),Ⅰ類(lèi)錯(cuò)誤隨著樣本量或者標(biāo)準(zhǔn)差之間差距拉大呈現(xiàn)遞增趨勢(shì)(見(jiàn)表3、4后半部分),變化幅度與樣本量之比成正比;表中結(jié)果說(shuō)明了mixed model和VWLS的Ⅰ類(lèi)錯(cuò)誤發(fā)生比例的穩(wěn)定性,與樣本量和標(biāo)準(zhǔn)差之比無(wú)關(guān)聯(lián)的趨勢(shì)性改變。(見(jiàn)表3、4后半部分)。

    表3 3組獨(dú)立樣本在不同情況下幾種方法Ⅰ類(lèi)錯(cuò)誤比較(1)

    *:每種情況重復(fù)獨(dú)立模擬2000次,標(biāo)準(zhǔn)差存在形式為σ1=σ2<σ3

    (2)比較power值大?。?組標(biāo)準(zhǔn)差存在形式如下:σ1=σ2<σ3\σ1<σ2<σ3)

    由表5、表6可知,不論標(biāo)準(zhǔn)差之比和樣本量之比大小如何變化,VWLS方法的 power值始終高于mixed model方法。

    討論與結(jié)論

    從統(tǒng)計(jì)學(xué)背景看,方差不齊的多組獨(dú)立樣本比較已經(jīng)不符合非參數(shù)秩和檢驗(yàn)的原假設(shè)H0,即總體分布相同。當(dāng)方差不齊時(shí),總體分布已經(jīng)為偏離H0假設(shè),為符合H1的情況。所以,方差不齊的程度越大,即越偏離H0假設(shè),Ⅰ類(lèi)錯(cuò)誤越不穩(wěn)健。如本研究發(fā)現(xiàn),當(dāng)樣本量相等時(shí),Ⅰ類(lèi)錯(cuò)誤大小隨方差之比的增加而增加,樣本量不等時(shí),Ⅰ類(lèi)錯(cuò)誤大小隨著樣本量之比的增大而拒絕比例漸漸遠(yuǎn)離檢驗(yàn)水準(zhǔn) α(本研究 α=0.05);然而縱觀 t′檢驗(yàn)、mixed model、VWLS,原假設(shè)是總體均數(shù)相同,方差不齊不屬于這些統(tǒng)計(jì)檢驗(yàn)方法的目標(biāo),且這3種方法容許存在方差不齊的情況而進(jìn)行分析計(jì)算,所以才使得Ⅰ類(lèi)錯(cuò)誤保持穩(wěn)健性,未隨著方差之比的變化而變化。

    從模擬結(jié)果來(lái)看,在表1中樣本量之比為1∶2,標(biāo)準(zhǔn)差之比 3時(shí),t′檢驗(yàn)、mixed model、VWLS的Ⅰ類(lèi)錯(cuò)誤全部增至0.06左右,因?yàn)樯鲜霈F(xiàn)象是獨(dú)立出現(xiàn)的,樣本量之比改變使其恢復(fù)正常,又因?yàn)樵摤F(xiàn)象涉及了研究所用的所有參數(shù)檢驗(yàn)方法,且這些方法的數(shù)學(xué)計(jì)算容許方差不齊的情況存在,因此考慮種子數(shù)設(shè)置引起的誤差所致,調(diào)整模擬所需種子數(shù)后,上述現(xiàn)象消失,Ⅰ類(lèi)錯(cuò)誤回歸至0.05附近(結(jié)果未顯示)。

    方差不齊的2組獨(dú)立樣本均數(shù)比較證實(shí)了樣本量相等的t檢驗(yàn)對(duì)于方差不齊具有穩(wěn)健性,同時(shí)說(shuō)明秩和檢驗(yàn)在樣本量相等、標(biāo)準(zhǔn)差之比不超過(guò)2時(shí)是穩(wěn)健的。因?yàn)?t′檢驗(yàn)、mixed model、VWLS這 3種參數(shù)檢驗(yàn)中,VWLS方法的Ⅰ類(lèi)錯(cuò)誤發(fā)生比例最大,所以mixed model和t′檢驗(yàn)比VWLS方法更穩(wěn)健,并且這三種方法的把握度又相互接近,所以t′檢驗(yàn)和mixed model方法更優(yōu)。

    當(dāng)方差不齊的3組獨(dú)立樣本均數(shù)比較時(shí),證實(shí)當(dāng)樣本量相等時(shí),方差分析對(duì)于方差不齊是不穩(wěn)健的,此結(jié)果與ROGAN文中的結(jié)果一致[4]。而秩和檢驗(yàn)在樣本量相等、方差之比較小如不超過(guò)2時(shí)才算穩(wěn)健,Ⅰ類(lèi)錯(cuò)誤偏頗不至太大。

    3組獨(dú)立樣本均數(shù)比較的模擬結(jié)果顯示,當(dāng)樣本量偏小時(shí),VWLS方法Ⅰ類(lèi)錯(cuò)誤發(fā)生比例易過(guò)大增至0.06左右,主要是因?yàn)閂WLS的加權(quán)方差由樣本的方差決定,所以樣本量較小時(shí)方差偏差較大,導(dǎo)致Ⅰ類(lèi)錯(cuò)誤發(fā)生比例增大;因?yàn)閙ixed model方法的Ⅰ類(lèi)錯(cuò)誤仍在0.05(本研究設(shè)置的α值)附近,所以樣本量較小時(shí)mixed model更穩(wěn)??;但是樣本量增加使VWLS的Ⅰ類(lèi)錯(cuò)誤大小恢復(fù)至設(shè)定值附近,又因?yàn)閂WLS的把握度明顯高于mixed model,這時(shí)以VWLS方法更優(yōu)。

    表4 3組獨(dú)立樣本在不同情況下幾種方法Ⅰ類(lèi)錯(cuò)誤比較(2)

    表5 3組獨(dú)立樣本在不同情況下幾種方法的power值大?。?)

    表6 3組獨(dú)立樣本在不同情況下幾種方法的Power值大小(2)

    [1]趙耐青,陳鋒主編.衛(wèi)生統(tǒng)計(jì)學(xué).北京:高等教育出版社,2008:49-61.

    [2]方積乾主編.衛(wèi)生統(tǒng)計(jì)學(xué).北京:人民衛(wèi)生出版社,2008:177-192.

    [3]Posten HO,ChengHY,Owen DB,et al.Robustness of the two-sample t-test under violations of the homogeneity of variance assumption.Communications in Statistics Theory and Methods,1982,11(2):109-126.

    [4]Rogan JC,Keselman HJ.Is the ANOVA F-Test Robust to Variance Heterogeneity When Sample Sizes Are Equal?:An Investigation via a Coefficient of Variation.American Educational Research Journal,1977,14(4):493-498.

    A Comparison of Different Statistical Methods of Their Robustness and Power When Analyzing the Population Mean Difference of Unequal-variance Data

    Jiang Dan,Zhao Naiqing
    (Biostatistical Department of Public Health School of Fudan University;Key Lab of Health Technology Assessment,Ministry of Health(Fudan University),(200032)Shanghai)

    Objective To compare the Type I error and power of different statistical methods,such as Rank Sum Test/t′test/Mixed Model and VWLS,for analyzing independent unequal variance data.Methods Under the hypothesis of equal population means and unequal population means,two/three independent groups with various ratios of standard deviation and sample size are simulated and compared using above statistical ways for thousands of times to acquire their each Type I Error and Power.Results (1)t test is robust to unequal variance data with equal sample size.However,ANOVA lacks its robustness to unequal variance data with equal sample size.(2)The method of Rank Sum Test is robust to unequal variance data only when special conditions are met.(3)t′test and Mixed Model are better than VWLSbecause of their model stable Type I Error when analyzing independent samples with unequal variance.Plus,the Power of these method is very close to each other.(4)When analyzing three independent unequal variance groups,Mixed Model is more robust than VWLS if the sample size is small;with the size getting bigger,VWLS gets robust and is better than Mixed Model for higher power.Conclusion When it comes to the analysis of two independent quantitative samples with unequal variance,t′test/Mixed Model/VWLS all can be used but the first two methods are better.When it comes to the analysis of three or more unequal variance groups,Mixed Model is better than VWLSunder the situation of small sample size,but with the size getting bigger,VWLSis getting better than Mixed Model.

    Unequal variance;Means compare;Type I error;Power

    郭海強(qiáng))

    猜你喜歡
    穩(wěn)健性均數(shù)樣本量
    醫(yī)學(xué)研究中樣本量的選擇
    會(huì)計(jì)穩(wěn)健性的定義和計(jì)量
    商情(2019年3期)2019-03-29 12:04:52
    航空裝備測(cè)試性試驗(yàn)樣本量確定方法
    會(huì)計(jì)穩(wěn)健性的文獻(xiàn)綜述
    財(cái)訊(2018年22期)2018-05-14 08:55:57
    Sample Size Calculations for Comparing Groups with Binary Outcomes
    貨幣政策、會(huì)計(jì)穩(wěn)健性與銀行信貸關(guān)系探析
    關(guān)于均數(shù)與偏差
    Effect of Ammonia on the Performance of Catalysts for Selective Hydrogenation of 1-Methylnaphthalene
    關(guān)于均數(shù)與偏差
    關(guān)于均數(shù)與偏差
    黄色一级大片看看| 精品人妻熟女av久视频| 七月丁香在线播放| 三级国产精品欧美在线观看| 精品久久久噜噜| 国产69精品久久久久777片| 国产美女午夜福利| 婷婷色综合www| 69人妻影院| 18禁在线播放成人免费| 亚洲怡红院男人天堂| 亚洲欧美日韩东京热| 欧美潮喷喷水| 九九爱精品视频在线观看| 建设人人有责人人尽责人人享有的 | 亚洲欧美日韩无卡精品| 成人免费观看视频高清| av黄色大香蕉| 香蕉精品网在线| 毛片女人毛片| 禁无遮挡网站| 插阴视频在线观看视频| 黄色怎么调成土黄色| a级毛片免费高清观看在线播放| 一级片'在线观看视频| 日韩精品有码人妻一区| 黄色日韩在线| 人人妻人人看人人澡| 乱系列少妇在线播放| 亚洲综合精品二区| 精品久久久久久久末码| 国产成人免费无遮挡视频| 久久99热这里只频精品6学生| 激情五月婷婷亚洲| 国产精品国产三级国产av玫瑰| 26uuu在线亚洲综合色| 免费看a级黄色片| 男女边摸边吃奶| 九九爱精品视频在线观看| 麻豆乱淫一区二区| av在线天堂中文字幕| a级毛片免费高清观看在线播放| 国产av不卡久久| 国产成人91sexporn| 免费播放大片免费观看视频在线观看| 婷婷色麻豆天堂久久| 日韩av不卡免费在线播放| 麻豆久久精品国产亚洲av| 人人妻人人澡人人爽人人夜夜| 亚洲在久久综合| av在线亚洲专区| 日韩成人伦理影院| 国产精品人妻久久久久久| 少妇的逼水好多| 亚洲在线观看片| 成人国产麻豆网| 一区二区三区精品91| 久久久午夜欧美精品| 亚洲精品,欧美精品| 日韩欧美 国产精品| 亚洲电影在线观看av| 精品国产一区二区三区久久久樱花 | 99久久精品热视频| 亚洲一级一片aⅴ在线观看| 国产一区亚洲一区在线观看| 精品一区二区免费观看| 日本色播在线视频| eeuss影院久久| 婷婷色综合大香蕉| 99热国产这里只有精品6| av国产免费在线观看| 精品久久久久久久久亚洲| 婷婷色综合大香蕉| 亚洲欧洲国产日韩| 少妇丰满av| 亚洲欧洲国产日韩| 精品一区二区免费观看| 国产精品国产三级国产av玫瑰| 色吧在线观看| 久久久国产一区二区| 丰满乱子伦码专区| 国产综合精华液| 久久99热这里只频精品6学生| 亚洲欧美日韩卡通动漫| 最近中文字幕高清免费大全6| 免费观看性生交大片5| 能在线免费看毛片的网站| 一本久久精品| 2021少妇久久久久久久久久久| 好男人在线观看高清免费视频| 国产精品99久久99久久久不卡 | 自拍偷自拍亚洲精品老妇| 边亲边吃奶的免费视频| 亚洲av日韩在线播放| 两个人的视频大全免费| 日韩欧美 国产精品| 国产黄频视频在线观看| 街头女战士在线观看网站| 欧美日韩视频精品一区| 亚洲精品成人久久久久久| 26uuu在线亚洲综合色| 丰满人妻一区二区三区视频av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产视频首页在线观看| 综合色丁香网| 六月丁香七月| 午夜老司机福利剧场| 日日摸夜夜添夜夜爱| 少妇猛男粗大的猛烈进出视频 | 国产成人精品福利久久| 三级经典国产精品| 在线看a的网站| av在线观看视频网站免费| 精品国产乱码久久久久久小说| 草草在线视频免费看| 日韩不卡一区二区三区视频在线| 2021天堂中文幕一二区在线观| 国产黄色免费在线视频| av.在线天堂| 精品熟女少妇av免费看| 午夜福利高清视频| 婷婷色麻豆天堂久久| 精品国产乱码久久久久久小说| 18禁裸乳无遮挡动漫免费视频 | 丝瓜视频免费看黄片| 丝袜喷水一区| 亚洲精品一二三| 欧美最新免费一区二区三区| 欧美区成人在线视频| 1000部很黄的大片| 麻豆久久精品国产亚洲av| 日产精品乱码卡一卡2卡三| 少妇人妻一区二区三区视频| 亚洲欧美成人综合另类久久久| 99久久精品热视频| 久久人人爽人人爽人人片va| a级毛片免费高清观看在线播放| 少妇裸体淫交视频免费看高清| 在线 av 中文字幕| 国产精品国产三级国产专区5o| 久久久久久久大尺度免费视频| 中文在线观看免费www的网站| 亚洲最大成人中文| 成人免费观看视频高清| 麻豆久久精品国产亚洲av| 久久99精品国语久久久| 免费播放大片免费观看视频在线观看| www.色视频.com| 中文精品一卡2卡3卡4更新| 成人漫画全彩无遮挡| 亚洲,一卡二卡三卡| 一级二级三级毛片免费看| 男女边吃奶边做爰视频| 国产高清有码在线观看视频| 嫩草影院入口| 在线免费十八禁| 下体分泌物呈黄色| 免费看日本二区| 国产乱人偷精品视频| 黄色欧美视频在线观看| 国产精品99久久久久久久久| 人体艺术视频欧美日本| 神马国产精品三级电影在线观看| 日韩伦理黄色片| 久久人人爽人人爽人人片va| 国产精品女同一区二区软件| 一二三四中文在线观看免费高清| 婷婷色综合www| 麻豆精品久久久久久蜜桃| 久久精品国产亚洲av天美| 蜜桃久久精品国产亚洲av| 久久这里有精品视频免费| 在线a可以看的网站| 成人毛片60女人毛片免费| 亚洲成人中文字幕在线播放| 午夜激情久久久久久久| 午夜亚洲福利在线播放| 美女cb高潮喷水在线观看| 久久人人爽人人片av| 中文字幕人妻熟人妻熟丝袜美| 美女视频免费永久观看网站| 久久精品国产a三级三级三级| 色5月婷婷丁香| 成人免费观看视频高清| 纵有疾风起免费观看全集完整版| 我的老师免费观看完整版| 最新中文字幕久久久久| 精品久久久久久久久av| av国产精品久久久久影院| 深爱激情五月婷婷| 亚洲婷婷狠狠爱综合网| 乱系列少妇在线播放| 日本三级黄在线观看| kizo精华| 全区人妻精品视频| 日韩欧美精品v在线| 国产av不卡久久| 欧美成人a在线观看| 精品人妻熟女av久视频| 别揉我奶头 嗯啊视频| 国产美女午夜福利| 乱码一卡2卡4卡精品| 联通29元200g的流量卡| 99久久九九国产精品国产免费| 91久久精品电影网| 久久热精品热| 少妇人妻久久综合中文| 国产av不卡久久| 国语对白做爰xxxⅹ性视频网站| 日韩 亚洲 欧美在线| 国产69精品久久久久777片| 亚洲三级黄色毛片| 黄色一级大片看看| 亚洲av国产av综合av卡| 精品久久久久久电影网| 丝袜美腿在线中文| 亚洲精品成人av观看孕妇| av黄色大香蕉| 搞女人的毛片| 国产黄片美女视频| av播播在线观看一区| 日日撸夜夜添| 久久97久久精品| 色视频www国产| 精品久久久噜噜| 黄色日韩在线| 亚洲熟女精品中文字幕| 听说在线观看完整版免费高清| 日韩成人av中文字幕在线观看| 亚洲欧美清纯卡通| 极品教师在线视频| 久久这里有精品视频免费| 成人亚洲欧美一区二区av| 亚洲欧美精品专区久久| 色哟哟·www| 国产精品熟女久久久久浪| 18禁动态无遮挡网站| 欧美老熟妇乱子伦牲交| 麻豆久久精品国产亚洲av| 久久99热这里只有精品18| 久久精品久久精品一区二区三区| 秋霞在线观看毛片| 观看免费一级毛片| 久久久久国产网址| 男女边摸边吃奶| 自拍欧美九色日韩亚洲蝌蚪91 | 天天一区二区日本电影三级| 国产免费一级a男人的天堂| 亚洲精品久久久久久婷婷小说| 国产av不卡久久| 国模一区二区三区四区视频| 晚上一个人看的免费电影| 国产色婷婷99| 欧美日韩一区二区视频在线观看视频在线 | 国产成人精品久久久久久| 少妇人妻久久综合中文| 欧美性感艳星| 最近2019中文字幕mv第一页| 精品一区二区三卡| av.在线天堂| 少妇人妻久久综合中文| 久久97久久精品| 91精品国产九色| 中文资源天堂在线| 午夜福利网站1000一区二区三区| 国产伦精品一区二区三区视频9| 国内揄拍国产精品人妻在线| 欧美一区二区亚洲| 国产黄色视频一区二区在线观看| 免费少妇av软件| 免费观看av网站的网址| 黑人高潮一二区| 嫩草影院新地址| 少妇人妻精品综合一区二区| 又黄又爽又刺激的免费视频.| 免费观看av网站的网址| 精品久久久久久电影网| 69av精品久久久久久| 少妇 在线观看| 免费黄频网站在线观看国产| 国产成人精品福利久久| 亚洲欧洲国产日韩| 91精品伊人久久大香线蕉| 中文字幕免费在线视频6| 欧美老熟妇乱子伦牲交| 久久久午夜欧美精品| 搡老乐熟女国产| 亚洲欧美日韩另类电影网站 | 国产免费一区二区三区四区乱码| 亚洲人成网站在线播| 色综合色国产| 久久精品久久久久久噜噜老黄| 久久99热这里只频精品6学生| 亚洲精品中文字幕在线视频 | 亚洲精品亚洲一区二区| 色哟哟·www| 国产成人aa在线观看| 日本一本二区三区精品| 99精国产麻豆久久婷婷| 色吧在线观看| 丰满少妇做爰视频| 97在线人人人人妻| av天堂中文字幕网| 日韩不卡一区二区三区视频在线| 丝袜美腿在线中文| 欧美成人精品欧美一级黄| 亚洲,欧美,日韩| 免费高清在线观看视频在线观看| 大码成人一级视频| 热99国产精品久久久久久7| 成人综合一区亚洲| 我的老师免费观看完整版| 一边亲一边摸免费视频| 亚洲人与动物交配视频| 人人妻人人看人人澡| 婷婷色麻豆天堂久久| 国产亚洲最大av| 乱系列少妇在线播放| 看黄色毛片网站| av播播在线观看一区| 视频中文字幕在线观看| 女人被狂操c到高潮| 十八禁网站网址无遮挡 | 老司机影院成人| 久久这里有精品视频免费| 亚洲精品自拍成人| 一级毛片电影观看| 欧美丝袜亚洲另类| 另类亚洲欧美激情| 日韩欧美精品v在线| 精品少妇久久久久久888优播| 蜜桃亚洲精品一区二区三区| 国产欧美亚洲国产| 最近最新中文字幕免费大全7| 日本一二三区视频观看| 精品一区二区三卡| 精品人妻偷拍中文字幕| 色视频在线一区二区三区| 如何舔出高潮| 网址你懂的国产日韩在线| 精品一区二区免费观看| 18禁裸乳无遮挡免费网站照片| 成人美女网站在线观看视频| 国内精品美女久久久久久| 777米奇影视久久| 国产精品国产三级国产专区5o| 大香蕉久久网| 国产精品av视频在线免费观看| 成人无遮挡网站| 成人综合一区亚洲| 天天躁夜夜躁狠狠久久av| 国产亚洲最大av| 国产69精品久久久久777片| 女人久久www免费人成看片| 亚洲国产av新网站| 高清欧美精品videossex| 欧美高清成人免费视频www| 国产真实伦视频高清在线观看| 一个人观看的视频www高清免费观看| 精品亚洲乱码少妇综合久久| 美女脱内裤让男人舔精品视频| av又黄又爽大尺度在线免费看| 下体分泌物呈黄色| 亚洲精华国产精华液的使用体验| 亚洲av国产av综合av卡| 久久综合国产亚洲精品| 国产一级毛片在线| av福利片在线观看| 国产成人91sexporn| 搡老乐熟女国产| 美女视频免费永久观看网站| 成人欧美大片| 欧美xxⅹ黑人| 成人特级av手机在线观看| 热re99久久精品国产66热6| 亚洲在久久综合| 草草在线视频免费看| 国产亚洲av嫩草精品影院| 80岁老熟妇乱子伦牲交| 黄片无遮挡物在线观看| 91精品伊人久久大香线蕉| 97在线人人人人妻| 91久久精品电影网| 久久人人爽av亚洲精品天堂 | 男人狂女人下面高潮的视频| 人人妻人人爽人人添夜夜欢视频 | 最近中文字幕高清免费大全6| 高清毛片免费看| 日本免费在线观看一区| 干丝袜人妻中文字幕| 成年版毛片免费区| 久久99热这里只有精品18| 2021少妇久久久久久久久久久| 国产成人午夜福利电影在线观看| 五月伊人婷婷丁香| 午夜老司机福利剧场| 日韩一区二区视频免费看| av播播在线观看一区| av国产久精品久网站免费入址| 欧美97在线视频| 国产精品成人在线| 国产免费一区二区三区四区乱码| 黄色欧美视频在线观看| 国产午夜精品久久久久久一区二区三区| 久久精品久久精品一区二区三区| av黄色大香蕉| 午夜日本视频在线| 69av精品久久久久久| 男女那种视频在线观看| 男的添女的下面高潮视频| 一区二区三区四区激情视频| 99热这里只有是精品50| 汤姆久久久久久久影院中文字幕| 久久6这里有精品| 观看免费一级毛片| 亚洲激情五月婷婷啪啪| 丝袜美腿在线中文| 国产探花在线观看一区二区| 国产精品成人在线| 最后的刺客免费高清国语| 国产伦精品一区二区三区四那| 少妇被粗大猛烈的视频| 在线 av 中文字幕| 人妻制服诱惑在线中文字幕| 各种免费的搞黄视频| 青春草国产在线视频| 人体艺术视频欧美日本| 狂野欧美激情性xxxx在线观看| 中文字幕免费在线视频6| 永久网站在线| 黄色怎么调成土黄色| 夜夜爽夜夜爽视频| 黄片无遮挡物在线观看| 色婷婷久久久亚洲欧美| 国产成人一区二区在线| 一区二区av电影网| 久久精品夜色国产| 国产高清不卡午夜福利| 国产黄a三级三级三级人| 蜜桃亚洲精品一区二区三区| 大话2 男鬼变身卡| 精品少妇久久久久久888优播| 精品人妻一区二区三区麻豆| 午夜福利在线在线| 亚洲va在线va天堂va国产| 欧美变态另类bdsm刘玥| 少妇丰满av| 狂野欧美白嫩少妇大欣赏| 欧美最新免费一区二区三区| 日本猛色少妇xxxxx猛交久久| 亚洲伊人久久精品综合| 久久久久久九九精品二区国产| 一个人观看的视频www高清免费观看| 联通29元200g的流量卡| 又大又黄又爽视频免费| 少妇人妻精品综合一区二区| 国产精品一区二区三区四区免费观看| 日韩三级伦理在线观看| 又爽又黄无遮挡网站| 你懂的网址亚洲精品在线观看| 有码 亚洲区| 狠狠精品人妻久久久久久综合| 国产女主播在线喷水免费视频网站| 我要看日韩黄色一级片| 亚洲人成网站高清观看| 国内揄拍国产精品人妻在线| 欧美丝袜亚洲另类| 插逼视频在线观看| 午夜免费鲁丝| 在线观看一区二区三区激情| 22中文网久久字幕| 一二三四中文在线观看免费高清| 99久久精品一区二区三区| 精品国产乱码久久久久久小说| 亚洲精品第二区| 欧美日韩综合久久久久久| 国产综合精华液| 成人欧美大片| 亚洲国产av新网站| 亚洲欧美日韩另类电影网站 | 青春草国产在线视频| 免费av不卡在线播放| 久久ye,这里只有精品| 涩涩av久久男人的天堂| 2021天堂中文幕一二区在线观| 2021少妇久久久久久久久久久| 日韩国内少妇激情av| 在线观看三级黄色| 成人亚洲精品av一区二区| 久久热精品热| 91久久精品国产一区二区成人| 欧美变态另类bdsm刘玥| 高清日韩中文字幕在线| 久久久久久久国产电影| 永久免费av网站大全| 亚洲va在线va天堂va国产| 成人高潮视频无遮挡免费网站| 一边亲一边摸免费视频| 涩涩av久久男人的天堂| 成人黄色视频免费在线看| 国产成人aa在线观看| 亚洲电影在线观看av| 性插视频无遮挡在线免费观看| 国产色爽女视频免费观看| av网站免费在线观看视频| 91在线精品国自产拍蜜月| 亚洲人成网站高清观看| 身体一侧抽搐| 国产免费一区二区三区四区乱码| 国产精品99久久久久久久久| 麻豆国产97在线/欧美| 国产高清三级在线| 亚洲国产av新网站| 又爽又黄无遮挡网站| 午夜日本视频在线| 麻豆成人午夜福利视频| 成人亚洲欧美一区二区av| 亚洲国产精品成人综合色| 日韩在线高清观看一区二区三区| 日日撸夜夜添| 国产久久久一区二区三区| 亚洲精品成人久久久久久| 国产视频内射| 天天躁夜夜躁狠狠久久av| 99久国产av精品国产电影| 看十八女毛片水多多多| 精品一区在线观看国产| 99热这里只有是精品50| 精品一区二区免费观看| 成人鲁丝片一二三区免费| 午夜日本视频在线| 亚洲国产精品成人久久小说| 肉色欧美久久久久久久蜜桃 | 三级经典国产精品| 一级a做视频免费观看| 丝袜喷水一区| 国产精品秋霞免费鲁丝片| 亚洲精品中文字幕在线视频 | 黑人高潮一二区| av在线app专区| 九九久久精品国产亚洲av麻豆| 超碰97精品在线观看| 久久99精品国语久久久| 下体分泌物呈黄色| 大陆偷拍与自拍| 国产精品久久久久久av不卡| 亚洲最大成人中文| 久久人人爽av亚洲精品天堂 | 黄色欧美视频在线观看| 亚洲国产精品专区欧美| 免费黄网站久久成人精品| 少妇人妻一区二区三区视频| 免费黄色在线免费观看| 亚洲精品自拍成人| 麻豆成人午夜福利视频| 午夜福利在线观看免费完整高清在| 97热精品久久久久久| 久久久久久久久久人人人人人人| 少妇猛男粗大的猛烈进出视频 | 亚洲av电影在线观看一区二区三区 | 特级一级黄色大片| 观看美女的网站| 欧美激情国产日韩精品一区| 我要看日韩黄色一级片| 国产午夜精品一二区理论片| 777米奇影视久久| 大又大粗又爽又黄少妇毛片口| 亚洲一级一片aⅴ在线观看| 欧美成人精品欧美一级黄| 神马国产精品三级电影在线观看| 亚洲人与动物交配视频| 22中文网久久字幕| 国产精品99久久99久久久不卡 | 国产亚洲精品久久久com| 黄色一级大片看看| 国产免费视频播放在线视频| 精品酒店卫生间| 大片免费播放器 马上看| 欧美最新免费一区二区三区| 好男人视频免费观看在线| 亚洲色图av天堂| 亚洲高清免费不卡视频| 看十八女毛片水多多多| 天天躁日日操中文字幕| 国产黄片美女视频| av卡一久久| .国产精品久久| 亚洲精品国产av成人精品| 寂寞人妻少妇视频99o| 成年人午夜在线观看视频| 亚洲av日韩在线播放| 日韩亚洲欧美综合| 一级毛片黄色毛片免费观看视频| 久久国内精品自在自线图片| 在线看a的网站| 熟女人妻精品中文字幕| 亚洲精品日韩在线中文字幕| 99热全是精品| 免费观看a级毛片全部| 亚洲人成网站在线播| 久久99热这里只有精品18| 又爽又黄a免费视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美变态另类bdsm刘玥| 国产精品精品国产色婷婷| 麻豆成人av视频| 亚洲欧美日韩卡通动漫| 五月开心婷婷网| 午夜老司机福利剧场| 亚洲av日韩在线播放| 亚洲欧美日韩另类电影网站 | a级一级毛片免费在线观看| 国产视频首页在线观看| 日韩国内少妇激情av| 国产成人a∨麻豆精品|