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

    研究設(shè)計中樣本量的確定

    2014-12-02 04:22:54
    外國語文 2014年5期
    關(guān)鍵詞:錯誤率樣本量平均數(shù)

    鮑 貴

    (南京工業(yè)大學(xué) 英語系,江蘇 南京 211816)

    1.引言

    任何推理統(tǒng)計都涉及到兩類錯誤率,即第一類錯誤率(Type I error rate,記作α)和第二類錯誤率(Type II error rate,記作β)。第一類錯誤率是零假設(shè)(H0)為真卻被錯誤拒絕的概率,因而不拒絕零假設(shè)的置信度為1-α。第二類錯誤率是零假設(shè)為誤卻被錯誤接受的概率,因而拒絕零假設(shè)的置信度為1-β。以往的研究重視第一類錯誤率,將α值設(shè)得很低,譬如.05或.01等,以期獲得科學(xué)的新發(fā)現(xiàn)。傳統(tǒng)上對第二類錯誤率的重視程度不足,主要是因為研究者通常希望證實與零假設(shè)相對立的備擇假設(shè)(Ha,又稱研究假設(shè))而非零假設(shè)本身,零假設(shè)的提出與拒絕只是為備擇假設(shè)的成立提供反證。

    實際上,第一類和第二類錯誤率是緊密聯(lián)系的,重視一類錯誤、忽略另一類錯誤不是很好的統(tǒng)計思維方法。近40年來,尤其在 Cohen(1969;1977;1988)的力作問世之后,研究者們逐漸意識到第二類錯誤率和第一類錯誤率一樣值得重視,因為它們均對能否得到科學(xué)的新發(fā)現(xiàn)起著至關(guān)重要的作用。不過,研究者們在研究第二類錯誤率時卻通常換用其補(bǔ)數(shù)(1-β),將之稱作統(tǒng)計效力(statistical power,簡稱效力),以此反映零假設(shè)為誤時研究拒絕零假設(shè)的能力。在科學(xué)的實證研究中,由于并不真正知道會犯哪類錯誤,因而研究者既需要控制第一類錯誤率,以避免在零假設(shè)為真的情況下得到虛假的研究發(fā)現(xiàn),又需要控制第二類錯誤率(即提高統(tǒng)計效力),以便在研究假設(shè)為真的情況下能夠得到研究發(fā)現(xiàn)。關(guān)于實證研究中效力的重要性,Hallahan&Rosenthal(1996:491)做了很好的概括:(1)如果計劃研究時不考慮效力,研究能夠發(fā)現(xiàn)實際存在的效應(yīng)的可能性或許很小。其結(jié)果是,由于不可能拒絕零假設(shè)而造成時間和資源的浪費。(2)由于解釋結(jié)果時沒有考慮效力,研究者可能將無顯著性結(jié)果解釋為零假設(shè)為真,從而過早地放棄很有前景的研究取向。Shadish等.(2002)和 Heppner等(2008)等將低統(tǒng)計效力視作統(tǒng)計結(jié)論效度的重要威脅。美國心理學(xué)會(American Psychological As-sociation,2010:30)建議,使用推理統(tǒng)計時,研究者應(yīng)考慮與假設(shè)檢驗相關(guān)的統(tǒng)計效力。這關(guān)系到在α水平、一定效應(yīng)量(effect size,簡稱ES)和樣本量(sample size)條件下正確拒絕所要檢驗的假設(shè)的可能性,因此,研究者應(yīng)常規(guī)性地提供證據(jù),表明研究有足夠效力發(fā)現(xiàn)有實質(zhì)性意義的效應(yīng)。

    在研究設(shè)計中,要保證統(tǒng)計效力,就要考慮選用合適的樣本量。外語教學(xué)研究者在研究設(shè)計中通常對選擇多大的樣本量沒有把握,樣本量選擇過大或過小的情況時常出現(xiàn)。樣本量選擇過大雖然不會削弱統(tǒng)計效力,但是往往會造成人力和物力資源的浪費,使研究設(shè)計顯得不經(jīng)濟(jì)。而且,樣本量過大也更易產(chǎn)生統(tǒng)計顯著性結(jié)果(即概率p<α),可能導(dǎo)致研究發(fā)現(xiàn)沒有實際意義的效應(yīng)量。當(dāng)然,習(xí)慣于直接解釋效應(yīng)量的研究者不會落入根據(jù)p值判斷效應(yīng)是否重要這一陷阱(Ellis,2010:52)。另一方面,樣本量過小則會使統(tǒng)計效力降低,減少了發(fā)現(xiàn)總體效應(yīng)的可能性,不僅造成人力和物力資源的浪費,還會導(dǎo)致錯誤的結(jié)論,譬如將統(tǒng)計不顯著的結(jié)果錯誤地解釋為接受零假設(shè)。因此,樣本量過大和過小都會降低研究的效率。優(yōu)化研究設(shè)計必須考慮效力、效應(yīng)量和樣本量之間的關(guān)系,選擇適當(dāng)?shù)臉颖玖俊1狙芯恳苑讲罘治?含t檢驗)為例介紹統(tǒng)計效力和效應(yīng)量的基本概念,分析樣本量與其他影響參數(shù)之間的關(guān)系,為樣本量的確定提供必要的方法。

    2.決定研究所需樣本量的參數(shù)及其相互關(guān)系

    2.1 樣本量、第一類錯誤率、效力和非中心參數(shù)

    就方差分析而言,傳統(tǒng)上影響效力的參數(shù)為第一類錯誤率、樣本量和非中心參數(shù)(noncentrality parameter,簡稱NCP,記作λ)。非中心參數(shù)表示研究假設(shè)偏離零假設(shè)的程度,計算公式為:

    公式中,SSHa是研究假設(shè)中期望平均數(shù)的平方和,n是每個比較總體平均數(shù)μi的樣本量,μ是各個比較總體平均數(shù)μi的平均數(shù),k是比較總體數(shù)。σ2是總體方差,常用方差分析中的誤差均方(MSE)來估計。零假設(shè)為真時,λ =0,F(xiàn)分布為中心分布。研究假設(shè)為真時,λ≠0,F(xiàn)分布為非中心分布。自由度確定時,零假設(shè)條件下的F分布只有一個,即中心F分布,但是同樣自由度時的非中心F分布卻有多個,分布的位置取決于λ。非中心參數(shù)λ值越大,非中心F分布就越遠(yuǎn)離中心分布,統(tǒng)計效力就越高。圖1顯示單因素方差分析時中心和非中心F分布中各個參數(shù)之間的關(guān)系。

    圖1 中心與非中心F分布比較

    左圖中心F分布F(k-1,k(n-1))臨界值右邊的面積等于α值。右圖中的β值為非中心F分布中小于或等于α值所對應(yīng)的F臨界值的F統(tǒng)計量的概率,效力為1-β。中心F分布臨界值是F分布的分子、分母的自由度和α的函數(shù),其中,分子、分母的自由度完全由比較總體數(shù)k和每組樣本量n決定。中心F分布的臨界值越大,臨界線就越向右移,β值就會越大,效力也就越低,反之亦然。α值不變時,對于比較組數(shù)(k)相同的方差分析,每組樣本量(n)越大,F(xiàn)分布臨界值就越小(圖中的臨界線左移),β值隨之減小,效力也就增加。如果α值也增加,F(xiàn)分布臨界值就會更小,β值也會更小,效力也就會再增加。

    下面仍以單因素方差分析為例探討不同參數(shù)之間的關(guān)系。假設(shè)三個比較總體的平均數(shù)分別為μ1=85、μ2=90 和 μ3=80,標(biāo)準(zhǔn)差 σ =10。我們利用R軟件(鮑貴,2012)編寫程序考察在α =.05、每組樣本量n取10-30區(qū)間21個不同值時,樣本量、非中心參數(shù)和效力之間的關(guān)系,統(tǒng)計結(jié)果見表1。

    表1 顯示,k=3、n=10 時,效力為 0.46,即在總體效應(yīng)存在的情況下,基于該樣本量發(fā)現(xiàn)它的概率卻不到50%。要提高效力,就要增加樣本量。當(dāng)每組樣本量增至20時,效力水平基本達(dá)到0.8。在此之前,效力隨樣本量增加而增加的速度較快,而在此之后,效力隨樣本量增加的速度有減緩的趨勢,樣本量的影響力減小。

    總體上,α值不變時,增大任何兩個參數(shù)值,都會增大第三個參數(shù)值。在一個參數(shù)值減小的情況下,要使另一個參數(shù)值保持不變,就要增大第三個參數(shù)值。譬如,如果非中心參數(shù)值小,要使效力保持在一個較高的水平,那就需要增加樣本量。實際研究中,由于α值通常設(shè)定為.05,研究者往往根據(jù)n和λ確定效力或根據(jù)λ和效力確定研究所需的樣本量。

    表1 不同參數(shù)之間的關(guān)系

    2.2 樣本量、第一類錯誤率、效力和效應(yīng)量

    從λ的計算公式可以看到,λ不獨立于樣本量n,即λ不是反映總體特征的一個固定值。這兩個參數(shù)的相互依賴性往往使得研究者不便依據(jù)λ來確定研究所需的樣本量。當(dāng)今研究中,往往用反映總體特征的一個固定參數(shù)來替代λ。這一固定參數(shù)就是Cohen(1969;1977;1988)提出的獨立于樣本量的效應(yīng)量。效應(yīng)量是指某個現(xiàn)象存在于總體中的程度或零假設(shè)錯誤的程度(Cohen,1988:9-10),測量上表現(xiàn)為標(biāo)準(zhǔn)化平均數(shù)的標(biāo)準(zhǔn)差。方差分析中,效應(yīng)量(用f表示)的計算公式為:

    零假設(shè)為真時,f=0;零假設(shè)為誤時,f≠ 0。在n和α值不變時,f值越大,效力就越高,反之亦然。表1的21個樣本中,雖然λ值各不相同,但是效應(yīng)量均相等(f≈0.41)。在f和α值不變時,如果要使效力保持在一個較高的水平,切實可行的辦法是增加樣本量。

    由于f是反映總體平均數(shù)差異大小的一個穩(wěn)定的特征,所以用它來估計研究所需的樣本量比使用λ更為方便。以下仍以三組單因素方差分析為例探討不同參數(shù)之間的關(guān)系。圖2顯示單因素方差分析(k=3,α =.05)中效應(yīng)量與樣本量之間的關(guān)系隨效力變化的趨勢。

    圖2 不同效力中效應(yīng)量與樣本量的關(guān)系

    圖2中的9條曲線自下而上反映效力水平為0.1-0.9時研究需要的每組樣本量??傮w上,效力恒定時,隨著效應(yīng)量的增加,研究所需的樣本量呈下降趨勢。效應(yīng)量介于0.1-0.25之間時,圖中曲線下降較為陡峭,說明樣本量受效力水平的影響較大。效應(yīng)量大于0.25之后,曲線變化較為平緩,說明樣本量受效力水平的影響減弱。效應(yīng)量大于0.4之后,9條線幾乎重合,效力水平對樣本量的影響大大減弱。效應(yīng)量小時,要使研究的效力保持在較高的水平,則必須擴(kuò)大樣本量。如果效應(yīng)量大,即便要保持較高的效力水平,樣本量也不必很大。

    表2以具體數(shù)值的形式反映單因素方差分析(k=3,α =.05)中每組樣本量、效力和效應(yīng)量之間的關(guān)系。

    表2 不同效力和效應(yīng)量條件下的每組樣本量

    當(dāng)f值為最小值0.1、每組樣本量為22時,效力只有0.1,因而雖然總體效應(yīng)存在,基于該樣本量發(fā)現(xiàn)它的概率卻只有10%。隨著樣本量的增大,效力也在增大。如果要達(dá)到0.8這一效力水平,每組樣本量則應(yīng)為323。隨著效應(yīng)量的提高,達(dá)到同等效力所需的每組樣本量則減少。譬如,在效應(yīng)量為0.3時,要達(dá)到效力水平0.8,每組所需樣本量降至37。如果效應(yīng)量高達(dá)0.5,每組樣本量為14時便可達(dá)到效力水平0.8。

    3.確定研究所需樣本量的方法

    3.1 確定α、效力水平和效應(yīng)量

    效力分析涉及四個參數(shù),即樣本量(n)、α、效力和效應(yīng)量(或非中心參數(shù))。知道其中的任何三個參數(shù)便可求得另一個參數(shù)。效力分析的主要目的之一是在研究設(shè)計階段根據(jù)參數(shù)之間的關(guān)系確定研究所需的樣本量。要確定研究所需的樣本量,就要確定其他三個參數(shù)。為了同時控制第一類和第二類錯誤率,α和β通常設(shè)定為很小的水平。按常規(guī),α =.05、.01或.001,其中 α 最常用值為.05。對于β值設(shè)定的規(guī)約性沒有α那么強(qiáng),但是當(dāng)今研究中β的通常值設(shè)定為0.1,即效力取值為0.9,0.8常被視作可接受的最小效力值(即 β =0.2)(Batterham & Atkinson,2005:158)。傳統(tǒng)上,β 取值大于α之值,反映出研究者對錯誤接受零假設(shè)不那么保守。正如Cohen(1988:56)所說,通常情況下,行為科學(xué)家判定第一類錯誤比第二類錯誤嚴(yán)重,因此需要更加嚴(yán)防。沒有發(fā)現(xiàn)比虛假的發(fā)現(xiàn)危害程度小這一認(rèn)識與傳統(tǒng)的科學(xué)觀是一致的。研究者通常固定α和β值,其目的是為了控制基于樣本的統(tǒng)計推理可能犯的不同類型錯誤。因此,研究者最重要的工作就是確定研究總體效應(yīng)量或非中心參數(shù)。由于非中心參數(shù)與樣本量不獨立,因而,比較切合實際的做法是估計總體效應(yīng)量。一旦明確了總體效應(yīng)量估計,研究所需的樣本量便可確定??傮w效應(yīng)量通常是未知的,準(zhǔn)確估計它至關(guān)重要,因為它不僅直接影響對研究所需樣本量的估計,而且也是科學(xué)研究追求的目標(biāo)。

    研究者通??梢杂貌煌椒ü烙嬁傮w效應(yīng)量。最好的方法是回顧相關(guān)實證研究,計算平均效應(yīng)量作為總體效應(yīng)量的最接近的估計。如果前期研究沒有提供基于樣本的效應(yīng)量,研究者可以根據(jù)它們提供的基本數(shù)據(jù)(比如平均數(shù)和標(biāo)準(zhǔn)差)計算估計的效應(yīng)量。計算平均效應(yīng)量可采用加權(quán)(weight)方法,最簡便的方法是以樣本量作為權(quán)重(Hunter&Schmidt,2004;Ellis,2010)。利用前期各項研究的效應(yīng)量估計平均效應(yīng)量采用加權(quán)的方式而不采用簡單地求它們的算術(shù)平均數(shù)(即各個研究效應(yīng)量之和除以研究的數(shù)量)的理由是,來自大樣本的效應(yīng)量估計因為取樣誤差小,所以比來自小樣本的估計更準(zhǔn)確,在平均效應(yīng)量估計中應(yīng)占有更大的權(quán)重。現(xiàn)以實驗研究中常用的獨立樣本t檢驗為例。該檢驗用于比較實驗組和對照組平均數(shù)差異,是單因素方差分析的簡化形式。效應(yīng)量(d)的計算公式為:

    公式中,ME和Mc分別代表實驗組和控制組總體平均數(shù),σ為每個總體的標(biāo)準(zhǔn)差(假設(shè)每個總體的標(biāo)準(zhǔn)差相等)。該公式也適用于非實驗條件下兩個獨立組比較時的效應(yīng)量計算。如果前期研究對兩個獨立樣本平均數(shù)差異采用單因素方差分析,則d=2f。實際研究中,用樣本平均數(shù)(和)和合并標(biāo)準(zhǔn)差(SD)估計對應(yīng)總體參數(shù)。兩組樣本量相等時,合并標(biāo)準(zhǔn)差為兩個樣本方差(VE和VC)平均數(shù)的平方根。假定某研究者回顧前期五項研究得到表3數(shù)據(jù)。

    表3 五項研究的結(jié)果和效應(yīng)量

    根據(jù)d的計算公式,求得各項研究的效應(yīng)量,如表3最后一列所示。從效應(yīng)量估計值來看,各研究之間有很大差異。總體效應(yīng)量的最優(yōu)估計為各研究效應(yīng)量加權(quán)平均數(shù)(),計算公式為:

    公式中,wi是權(quán)重,di是第i個研究的效應(yīng)量,k是研究的數(shù)量。本研究使用的權(quán)重是樣本量Ni(第i個研究兩個樣本量之和),經(jīng)計算得到≈ 0.35。

    由于d可以有方向性地比較實驗組和控制組標(biāo)準(zhǔn)化平均數(shù)差異的大小,因而在元分析(metaanalysis)中常采用加權(quán)的方法估計總體效應(yīng)。方差分析是無方向性的。在比較多組(三組或三組以上)平均數(shù)差異時,該分析只能回答各組之間有無顯著性差異存在,不能明確差異的具體位置,除非進(jìn)行多重比較(multiple comparisons)。正因為如此,統(tǒng)計學(xué)家不對多組方差分析中的效應(yīng)量f進(jìn)行元分析。不過,在研究設(shè)計階段,為了確定研究所需的樣本量,研究者也可以用以上對d的加權(quán)方法估計平均效應(yīng)量f,作為對總體效應(yīng)量的最近似估計。如果研究者不采用這一估算方法,也可以只計算出與自身研究相似的前期研究的效應(yīng)量f,再決定合理的效應(yīng)量應(yīng)該有多大。

    另一種方法是在沒有前期實證研究結(jié)果做參考的情況下開展先導(dǎo)研究,對研究總體的效應(yīng)量進(jìn)行恰當(dāng)?shù)墓烙?。第三種方法是結(jié)合研究實際合理地主觀估計效應(yīng)量。該方法適用于既沒有前期相關(guān)實證研究結(jié)果作為支撐、條件又不允許開展先導(dǎo)研究的情形。由此得到的總體效應(yīng)量估計可能不是很精確。最后,如果研究者認(rèn)為合理的估計很難,不妨采用Cohen(1988)的建議。Cohen(1988)針對各種統(tǒng)計分析方法提出的效應(yīng)量大小的標(biāo)準(zhǔn)可以作為合理的參照。譬如,Cohen(1988:25-27,284-288)將t檢驗中小、中、大效應(yīng)量分別操作定義為d=0.2、d=0.5 和d=0.8;方差分析中的小、中、大效應(yīng)量分別操作定義為f=0.1、f=0.25和f=0.4。在缺乏判斷效應(yīng)量大小的依據(jù)時,建議研究者采用小效應(yīng)量計算研究所需的樣本量,以免效應(yīng)量估計過高導(dǎo)致實際統(tǒng)計效力降低。不管研究的實際效應(yīng)量有多大,如果研究有足夠的效力確保發(fā)現(xiàn)小效應(yīng),那么犯第一類或第二類嚴(yán)重錯誤的風(fēng)險就會很小(Murphy&Myors,2004:59)。

    3.2 計算研究所需樣本量的方法

    在設(shè)定第一類錯誤率α和效力水平,并且估計出總體效應(yīng)量之后,便可以確定研究所需的樣本量了。研究所需樣本量的確定方法主要有以下幾種。第一種是利用R和SAS等軟件的外置或附加效力分析程序計算出不同統(tǒng)計分析需要的樣本量。譬如,表3中的平均效應(yīng)量珔d≈0.35,設(shè)α=.05,效力水平為0.8。使用R的外置程序pwr,輸入命令:pwr.t.test(d=0.35,sig.level=.05,power=.8,type= “two.sample”,alternative= “greater”),則單側(cè)t檢驗所需要的每組樣本量為102。這意味著,要使研究在80%的情況下能夠發(fā)現(xiàn)總體效應(yīng)(0.35),每組的樣本量應(yīng)為102。如果研究者使用雙側(cè)t檢驗,則只需將命令中的“greater”改為“two.sided”,由此得到n=130,比單側(cè)t檢驗需要的樣本量要大一些。

    PASS和G*Power等效力分析專用軟件也是很方便的選擇。這些軟件能夠計算各種統(tǒng)計分析的效力或需要的樣本量,并能繪制出各種效力分析圖形。G*Power是免費使用軟件(網(wǎng)址:http://www.psycho.uni-duesseldorf.de/abteilungen/aap/gpower3/download-and-register),且界面操作簡便,因此建議研究者使用該軟件。關(guān)于該軟件各種統(tǒng)計分析功能的詳細(xì)介紹,研究者可參考Faulet al.(2007)、Mayret al.(2007)和 Lan & Lian(2010)。

    第三種方法是查閱效力分析樣本量表。以Cohen(1988)提供的不同統(tǒng)計分析(包括t檢驗、方差分析和相關(guān)分析等)樣本量表為例。如果研究的某個效應(yīng)量在表中有顯示,則查表比較方便:在某個α水平上,某個效應(yīng)量所在的列與設(shè)定的某個效力水平所在的排交叉得到的數(shù)值即為樣本量。譬如,d=0.30、α =.05、效力水平為 0.8 時,查表得到單側(cè)獨立樣本t檢驗需要的每組樣本量n=138。但是,如果研究的效應(yīng)量d不在表中,則需要采用近似計算的方法得到樣本量,計算公式為:

    公式中,分子n.10是在某個α和效力水平上d=0.10時所需的樣本量,分母中的d為表中沒有顯示的效應(yīng)量,常數(shù)1是調(diào)整所需樣本量的經(jīng)驗平均值。根據(jù)表3求出的加權(quán)平均效應(yīng)量0.35不在效力分析表中,但是表中提供了單側(cè)獨立樣本t檢驗在 α =.05、效力水平為 0.8 時的n.10值(1237)。在α =.05(單側(cè)t檢驗)時,樣本略微更精確的估計可采用常數(shù)0.7代替1,經(jīng)計算得到本例所需的樣本量為n≈102。如果采用雙側(cè)t檢驗,α值和效力水平不變,則n.10=1571(表中值),經(jīng)計算得到n≈ 130。

    關(guān)于方差分析樣本量表需要做點說明。Cohen(1988)提供的表為單因素方差分析樣本量表。表中列出在分子自由度不同,α和效應(yīng)量f取不同值時不同效力水平對應(yīng)的樣本量。如果研究的效應(yīng)量f不在表中,則也需要采用近似計算的方法得到樣本量,計算公式為:

    公式中,分子n.05是在某個α和效力水平上f=0.05時所需的樣本量,分母中的f為表中沒有顯示的效應(yīng)量,1為調(diào)整樣本量的常數(shù)。前面提到,獨立樣本t檢驗是單因素方差分析的簡化形式(比較組限于2個)。用方差分析同樣可以得到研究所需的樣本量。兩組比較時,表3求出的效應(yīng)量d≈0.35等同于f≈0.175。在 α =.05,效力為0.8時,n.05=1571(表中值),經(jīng)計算得到n≈ 130,與上面得到的樣本量相同。

    Cohen(1988:396-403)談到因素方差分析中在計算某個效力水平上發(fā)現(xiàn)主效應(yīng)和交互效應(yīng)所需的樣本量時如何利用單因素方差分析樣本量表的問題。各個效應(yīng)量單元格樣本量(nc)的計算公式為:

    公式中,n'數(shù)值上等于表中的n,number of cells為單元格總數(shù),u是效應(yīng)自由度(等于因素水平數(shù)-1),末尾的1是校正單元格樣本〗uB=2,fB=0.4;uA×B=2,fA×B=0.25。為了確定各個效應(yīng)需要的樣本量,研究者決定將α和效力水平分別設(shè)定為.05和0.8。對于A因素效應(yīng),查表得到n=n'=45,利用上面的公式得到nc≈16,那么研究所需的總樣本量N應(yīng)為96。同樣方法可以得到:因素B的單元格nc=11,總樣本量N為66;交互作用A×B的單元格nc≈27,總樣本量N為162。如果研究者要在效力為0.8的水平上發(fā)現(xiàn)中等效應(yīng)的交互作用,則總樣本量應(yīng)為162。由于162是三個效應(yīng)量所需樣本量的最大值,所以使用此樣本量能夠保證兩個主效應(yīng)的效力水平在0.8以上。用G*Power計算得到的各效應(yīng)所需總樣本量分別為90(A因素)、64(B因素)和 158(A×B交互作用)。由于舍入和計算方法的不同,這些值略低于利用表格計算得到的對應(yīng)總樣本量。在算法上,Cohen(1988)表中的樣本量采用近似算法,G*Power采用精確算法??傮w上,這兩種方法得到的結(jié)果基本一致。對于因素設(shè)計,Cohen(1988)算法得到的樣本量估計偏高(Erdfelderet al.,1996)。

    4.結(jié)語

    統(tǒng)計學(xué)是發(fā)展迅猛的一門學(xué)科。傳統(tǒng)的統(tǒng)計分析方法和統(tǒng)計觀念不斷地被調(diào)整,新的統(tǒng)計分析方法不斷涌現(xiàn)。既然外語教學(xué)研究是以統(tǒng)計分析為主導(dǎo)的,那么外語教學(xué)研究者就有必要跟得上統(tǒng)計學(xué)發(fā)展的步伐,更新研究方法和統(tǒng)計觀念。效力分析是外語教學(xué)研究必不可少的研究設(shè)計與分析環(huán)節(jié)。其主要目的之一就是為了在規(guī)定第一類和第二類錯誤率的前提下選擇適當(dāng)?shù)臉颖玖浚_??傮w效應(yīng)量存在時一項研究在很大程度上能夠發(fā)現(xiàn)它。

    外語教學(xué)研究者對研究所需樣本量問題的認(rèn)識不夠深刻,對效力分析還比較陌生。鑒于此,本文以單因素方差分析為例,介紹了決定研究所需樣本量的三個重要參數(shù),即第一類錯誤率、效力和效應(yīng)量(或非中心參數(shù)),探討了它們之間的關(guān)系,詳細(xì)分析了確定效應(yīng)量的方法以及計算研究所需樣本量的多種途徑,以期引起研究者對確定樣本量問題的重視,并在未來研究設(shè)計中能夠確定恰當(dāng)?shù)臉颖玖浚岣呓y(tǒng)計分析的水平。

    [1] American Psychological Association.Publication Manual of the American Psychological Association(6th ed.)[M].Washington,DC:Author,2010.

    [2]Batterham,A.M.& G.Atkinson.How Big does My Sample Need to Be?A Primer on the Murky World of Sample Size Estimation [J].Physical Therapy in Sport,2005(6)3:153-163.

    [3]Cohen,J.Statistical Power Analysis for the Behavioral Sciences[M].New York:Academic Press,1969.

    [4]Cohen,J.Statistical Power Analysis for the Behavioral Sciences[M].New York:Academic Press,1977.

    [5]Cohen,J.Statistical Power Analysis for the Behavioral Sciences[M].Hillsdale,NJ:Lawrence Erlbaum Associates,1988.

    [6]Ellis,P.D.The Essential Guide to Effect Sizes:Statistical Power,Meta-Analysis,and the Interpretation of Research Results[M]. Cambridge: Cambridge University Press,2010.

    [7]Erdfelder,E.,F(xiàn).Faul& A.Buchner.G*POWER:A General Power Analysis Program [J].Behavior Research Methods,Instruments,&Computers,1996,(28)1:1–11.

    [8]Faul,F(xiàn).,E.Erdfelder,A.G.Lang& A.Buchner.G*Power 3:A Flexible Statistical Power Analysis Program for the Social,Behavioral,and Biomedical Sciences[J].Be-havior Research Methods,2007(39)2:175-191.

    [9]Hallahan,M.& R.Rosenthal.Statistical Power:Concepts,Procedures,and Applications[J].Behaviour Research and Therapy,1996,(34)5/6:489-499.

    [10]Heppner,P.P.,B.E.Wampold& D.M.Kivlighan,Jr.Research Design in Counseling(3th ed.)[M].Belmont,CA:Thomson Brooks/Cole,2008.

    [11]Hunter,J.E.& F.L.Schmidt.Methods of Meta-Analysis:Correcting Error and Bias in Research Findings(2nd ed.) [M]. London:Sage Publications Ltd.,2004.

    [12]Lan,L.& ZW.Lian.Application of Statistical Power A-nalysis–How to Determine the Right Sample Size in Human Health,Comfort and Productivity Research[J].Building and Environment,2010,(45)5:1202-1213.

    [13] Mayr,S.,E.Erdfelder,A.Buchner & F.Faul.A Short Tutorial of G*Power[J].Tutorials in Quantitative Methods for Psychology,2007,(3)2:51-59.

    [14]Murphy,K.R.& B.Myors.Statistical Power Analysis:A Simple and General Model for Traditional and Modern Hypothesis Tests(2nd ed.)[M].Mahwah,NJ:Lawrence Erlbaum Associates,2004.

    [15]Shadish,W.R.,T.D.Cook& D.T.Campbell.Experimental and Quasi-Experimental Designs for Generalized Causal Inference[M].Boston:Houghton Mifflin Company,2002.

    [16]鮑貴.多元回歸分析中的交互作用問題——以語言閾限假設(shè)檢驗為例[J].外國語文,2012(4):63-68.

    猜你喜歡
    錯誤率樣本量平均數(shù)
    限制性隨機(jī)試驗中選擇偏倚導(dǎo)致的一類錯誤率膨脹*
    加權(quán)平均數(shù)的應(yīng)用
    醫(yī)學(xué)研究中樣本量的選擇
    航空裝備測試性試驗樣本量確定方法
    Sample Size Calculations for Comparing Groups with Binary Outcomes
    正視錯誤,尋求策略
    教師·中(2017年3期)2017-04-20 21:49:49
    關(guān)注加權(quán)平均數(shù)中的“權(quán)”
    平均數(shù)應(yīng)用舉隅
    說說加權(quán)平均數(shù)
    解析小學(xué)高段學(xué)生英語單詞抄寫作業(yè)錯誤原因
    最近2019中文字幕mv第一页| 香蕉国产在线看| 精品亚洲乱码少妇综合久久| 日韩精品免费视频一区二区三区| 国产精品av久久久久免费| 在线观看免费高清a一片| 国产精品一区二区精品视频观看| 国产精品国产av在线观看| 亚洲国产中文字幕在线视频| 捣出白浆h1v1| 久久久久久久久久久久大奶| 国产无遮挡羞羞视频在线观看| 国产免费又黄又爽又色| 2018国产大陆天天弄谢| 午夜91福利影院| 亚洲国产欧美一区二区综合| 中国三级夫妇交换| avwww免费| 90打野战视频偷拍视频| 七月丁香在线播放| 日韩电影二区| 欧美最新免费一区二区三区| 操美女的视频在线观看| 91成人精品电影| 亚洲五月色婷婷综合| 亚洲第一av免费看| av女优亚洲男人天堂| 中文字幕人妻丝袜制服| 日日爽夜夜爽网站| 99热国产这里只有精品6| 久久青草综合色| 深夜精品福利| 亚洲美女搞黄在线观看| 日韩中文字幕欧美一区二区 | 亚洲精品日韩在线中文字幕| 90打野战视频偷拍视频| 免费在线观看黄色视频的| 亚洲国产精品国产精品| 国产亚洲欧美精品永久| 丰满少妇做爰视频| 精品人妻在线不人妻| 国产99久久九九免费精品| 制服人妻中文乱码| 国产精品二区激情视频| 综合色丁香网| 在线观看一区二区三区激情| 成年人午夜在线观看视频| 午夜福利网站1000一区二区三区| av又黄又爽大尺度在线免费看| 国产色婷婷99| 男女免费视频国产| 久久久精品94久久精品| 亚洲第一区二区三区不卡| 国产一区二区在线观看av| 成年动漫av网址| 丝袜美腿诱惑在线| 另类精品久久| 精品视频人人做人人爽| 丝袜美足系列| 韩国精品一区二区三区| 国产男女内射视频| 久久久久精品性色| 久久久久精品久久久久真实原创| 日韩精品免费视频一区二区三区| 在线看a的网站| 男人爽女人下面视频在线观看| 免费人妻精品一区二区三区视频| 国产高清国产精品国产三级| 精品国产一区二区三区久久久樱花| 久久精品亚洲熟妇少妇任你| 一级爰片在线观看| 精品一区二区三区av网在线观看 | 亚洲,一卡二卡三卡| 99热全是精品| 日本色播在线视频| 日本爱情动作片www.在线观看| 性色av一级| 亚洲熟女精品中文字幕| 综合色丁香网| 久久性视频一级片| 一级毛片电影观看| 在线精品无人区一区二区三| 国产精品偷伦视频观看了| 国产亚洲精品第一综合不卡| 女性生殖器流出的白浆| 777米奇影视久久| 国产又色又爽无遮挡免| 亚洲欧美一区二区三区黑人| 免费黄网站久久成人精品| 日本黄色日本黄色录像| 亚洲一区中文字幕在线| 另类亚洲欧美激情| 大码成人一级视频| 一级毛片黄色毛片免费观看视频| 日韩欧美一区视频在线观看| 婷婷成人精品国产| 天天躁夜夜躁狠狠躁躁| 国产又爽黄色视频| 十八禁网站网址无遮挡| 永久免费av网站大全| 99精品久久久久人妻精品| 美女扒开内裤让男人捅视频| 午夜福利视频精品| 亚洲在久久综合| 天天添夜夜摸| 九九爱精品视频在线观看| 国产精品秋霞免费鲁丝片| 日韩一本色道免费dvd| 99精品久久久久人妻精品| 免费在线观看黄色视频的| 久久久精品区二区三区| 欧美人与性动交α欧美精品济南到| 国产伦人伦偷精品视频| 亚洲熟女毛片儿| 99久国产av精品国产电影| 超碰成人久久| 午夜福利网站1000一区二区三区| 满18在线观看网站| 热re99久久精品国产66热6| 亚洲精品成人av观看孕妇| 久久精品熟女亚洲av麻豆精品| 丝袜人妻中文字幕| 9191精品国产免费久久| 亚洲精华国产精华液的使用体验| 日本色播在线视频| 中文乱码字字幕精品一区二区三区| 日韩一区二区视频免费看| 又大又爽又粗| 国产免费视频播放在线视频| 日韩精品免费视频一区二区三区| 中文字幕人妻熟女乱码| 9热在线视频观看99| 国产精品香港三级国产av潘金莲 | 国产成人免费观看mmmm| 免费观看a级毛片全部| 国产男人的电影天堂91| 免费观看性生交大片5| 韩国精品一区二区三区| 麻豆av在线久日| 日本爱情动作片www.在线观看| 国产乱人偷精品视频| 国产福利在线免费观看视频| 人人澡人人妻人| 老司机靠b影院| 夫妻性生交免费视频一级片| 欧美另类一区| 三上悠亚av全集在线观看| 狠狠婷婷综合久久久久久88av| svipshipincom国产片| 国产一级毛片在线| 国产精品久久久久久人妻精品电影 | 成人三级做爰电影| 国产亚洲一区二区精品| 色网站视频免费| 成人三级做爰电影| 悠悠久久av| 成人亚洲欧美一区二区av| 99热网站在线观看| 亚洲成人一二三区av| 欧美黑人欧美精品刺激| 国产探花极品一区二区| 欧美人与性动交α欧美软件| 亚洲精品国产av蜜桃| 看免费成人av毛片| 免费在线观看黄色视频的| www.熟女人妻精品国产| 国产精品女同一区二区软件| 一二三四中文在线观看免费高清| 久久久久精品性色| 久久久久久久国产电影| 亚洲在久久综合| 最新的欧美精品一区二区| 欧美另类一区| 老司机在亚洲福利影院| 国产免费又黄又爽又色| 一级毛片黄色毛片免费观看视频| 人人妻人人澡人人看| 成人18禁高潮啪啪吃奶动态图| 国产黄色视频一区二区在线观看| 天天影视国产精品| 精品国产乱码久久久久久男人| 十分钟在线观看高清视频www| 天天躁日日躁夜夜躁夜夜| 91成人精品电影| 女人久久www免费人成看片| 丝袜美腿诱惑在线| 赤兔流量卡办理| 如日韩欧美国产精品一区二区三区| 国产视频首页在线观看| 久热爱精品视频在线9| 日韩视频在线欧美| www.av在线官网国产| av免费观看日本| 国精品久久久久久国模美| 国产野战对白在线观看| 久久99一区二区三区| 久久午夜综合久久蜜桃| 美女视频免费永久观看网站| 国产一区二区三区综合在线观看| 成人黄色视频免费在线看| 观看美女的网站| 啦啦啦啦在线视频资源| 久热这里只有精品99| 久久免费观看电影| 成人漫画全彩无遮挡| 黄网站色视频无遮挡免费观看| 80岁老熟妇乱子伦牲交| av网站免费在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 99久国产av精品国产电影| 一级爰片在线观看| 免费高清在线观看日韩| 最黄视频免费看| 久久久久久久精品精品| 在线免费观看不下载黄p国产| 涩涩av久久男人的天堂| 999久久久国产精品视频| 丝袜美腿诱惑在线| 大片电影免费在线观看免费| 天堂俺去俺来也www色官网| a级片在线免费高清观看视频| 婷婷色综合大香蕉| 欧美人与性动交α欧美软件| 国产日韩欧美在线精品| 欧美人与性动交α欧美精品济南到| 久久综合国产亚洲精品| 国产精品蜜桃在线观看| 亚洲成人av在线免费| 久久久久视频综合| 少妇 在线观看| 日韩,欧美,国产一区二区三区| 十分钟在线观看高清视频www| 久热爱精品视频在线9| 欧美人与善性xxx| 亚洲国产av影院在线观看| 女性被躁到高潮视频| 久久毛片免费看一区二区三区| 中文字幕人妻熟女乱码| 亚洲成人手机| 伊人久久国产一区二区| 制服丝袜香蕉在线| 久久ye,这里只有精品| 丝袜在线中文字幕| 最近中文字幕2019免费版| 日本猛色少妇xxxxx猛交久久| 性色av一级| 精品国产一区二区三区久久久樱花| 国产视频首页在线观看| 亚洲伊人色综图| 天天躁夜夜躁狠狠久久av| 9色porny在线观看| 伊人久久大香线蕉亚洲五| 亚洲国产精品一区二区三区在线| 国产片内射在线| 男人舔女人的私密视频| 最近最新中文字幕大全免费视频 | 亚洲国产欧美在线一区| 人妻人人澡人人爽人人| 日韩精品免费视频一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 在线观看免费午夜福利视频| 咕卡用的链子| 国产亚洲精品第一综合不卡| 香蕉丝袜av| 中国三级夫妇交换| 下体分泌物呈黄色| 蜜桃在线观看..| 国产精品成人在线| 日本爱情动作片www.在线观看| 18禁观看日本| 亚洲中文av在线| 国产一区二区在线观看av| 日日爽夜夜爽网站| 十八禁高潮呻吟视频| 狠狠婷婷综合久久久久久88av| 国产成人啪精品午夜网站| 亚洲三区欧美一区| 久久久久精品人妻al黑| 亚洲伊人久久精品综合| 亚洲自偷自拍图片 自拍| avwww免费| 久久99精品国语久久久| 久久人人爽av亚洲精品天堂| 久久天堂一区二区三区四区| 欧美另类一区| 人妻 亚洲 视频| 日韩中文字幕视频在线看片| av在线老鸭窝| 亚洲精品乱久久久久久| www.自偷自拍.com| 欧美日韩福利视频一区二区| 悠悠久久av| 亚洲人成电影观看| 国产日韩欧美视频二区| 汤姆久久久久久久影院中文字幕| 亚洲欧美日韩另类电影网站| 伊人亚洲综合成人网| 丰满少妇做爰视频| 老司机影院毛片| 国产av国产精品国产| 免费在线观看完整版高清| 亚洲国产中文字幕在线视频| 在线观看免费视频网站a站| 操美女的视频在线观看| 精品亚洲成a人片在线观看| 国产精品一国产av| 日本欧美国产在线视频| 欧美人与善性xxx| 免费黄网站久久成人精品| 欧美激情 高清一区二区三区| 日韩人妻精品一区2区三区| 国产1区2区3区精品| 日韩一区二区视频免费看| 欧美精品一区二区免费开放| 黑人欧美特级aaaaaa片| 亚洲国产欧美日韩在线播放| 黑丝袜美女国产一区| 亚洲色图 男人天堂 中文字幕| 综合色丁香网| 麻豆精品久久久久久蜜桃| 人妻 亚洲 视频| 国产在线免费精品| 在线观看免费午夜福利视频| 熟女少妇亚洲综合色aaa.| 一区二区三区乱码不卡18| 欧美精品高潮呻吟av久久| 纵有疾风起免费观看全集完整版| 成人国语在线视频| 午夜激情av网站| 人人妻人人爽人人添夜夜欢视频| 日韩制服骚丝袜av| 亚洲四区av| 国产欧美亚洲国产| 久久青草综合色| 亚洲第一青青草原| 国产日韩一区二区三区精品不卡| 80岁老熟妇乱子伦牲交| 熟女少妇亚洲综合色aaa.| 久久久国产欧美日韩av| 国产成人免费观看mmmm| 少妇精品久久久久久久| 国产男人的电影天堂91| 日韩一区二区视频免费看| 久久久国产一区二区| videos熟女内射| 一区二区av电影网| 久久久久久人妻| 国产一区二区三区av在线| 一边摸一边做爽爽视频免费| 久久久亚洲精品成人影院| 日本欧美国产在线视频| 亚洲欧美精品自产自拍| 亚洲综合精品二区| 欧美激情高清一区二区三区 | 国产伦理片在线播放av一区| 久久热在线av| 国产不卡av网站在线观看| 成人三级做爰电影| 日韩 欧美 亚洲 中文字幕| 中文字幕色久视频| 国产日韩欧美在线精品| av又黄又爽大尺度在线免费看| 高清在线视频一区二区三区| 亚洲av男天堂| 久久久久久人妻| 免费在线观看黄色视频的| 三上悠亚av全集在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 自线自在国产av| 精品国产一区二区久久| 亚洲成人免费av在线播放| 色婷婷av一区二区三区视频| 女人高潮潮喷娇喘18禁视频| 十分钟在线观看高清视频www| 日本91视频免费播放| 亚洲av男天堂| 欧美 日韩 精品 国产| 看十八女毛片水多多多| 亚洲av欧美aⅴ国产| 亚洲一级一片aⅴ在线观看| 久热爱精品视频在线9| 亚洲av电影在线进入| 大码成人一级视频| 亚洲成人一二三区av| 99国产综合亚洲精品| 一级片'在线观看视频| 美女视频免费永久观看网站| 狠狠婷婷综合久久久久久88av| 一级毛片电影观看| 亚洲图色成人| 电影成人av| 丝袜美足系列| 亚洲精华国产精华液的使用体验| 操美女的视频在线观看| 性色av一级| 最近最新中文字幕大全免费视频 | 欧美中文综合在线视频| 久久久久久人妻| 汤姆久久久久久久影院中文字幕| 亚洲综合精品二区| 18禁国产床啪视频网站| 桃花免费在线播放| 18禁动态无遮挡网站| 日韩av不卡免费在线播放| 亚洲国产欧美网| 久久精品国产a三级三级三级| 中文字幕人妻丝袜制服| 人人澡人人妻人| 日本欧美视频一区| 亚洲一卡2卡3卡4卡5卡精品中文| a级毛片黄视频| 波野结衣二区三区在线| 日韩av免费高清视频| 男人舔女人的私密视频| 国产精品国产三级专区第一集| av不卡在线播放| a级片在线免费高清观看视频| 观看av在线不卡| 亚洲欧美一区二区三区国产| 人人妻人人澡人人爽人人夜夜| 国产精品一区二区在线不卡| 久久精品国产a三级三级三级| 国产黄频视频在线观看| 五月天丁香电影| 午夜福利在线免费观看网站| 黄片播放在线免费| 国产淫语在线视频| 一级毛片 在线播放| 久久久久久久大尺度免费视频| 99热国产这里只有精品6| 高清在线视频一区二区三区| 亚洲精品国产一区二区精华液| 精品国产一区二区三区四区第35| 久久久精品94久久精品| 亚洲精品美女久久av网站| 高清不卡的av网站| 免费高清在线观看日韩| 下体分泌物呈黄色| 如日韩欧美国产精品一区二区三区| 18禁观看日本| 亚洲国产欧美在线一区| 熟女av电影| 日日摸夜夜添夜夜爱| 久久久国产精品麻豆| 欧美 亚洲 国产 日韩一| 亚洲精品aⅴ在线观看| 国产精品国产三级国产专区5o| 69精品国产乱码久久久| videos熟女内射| 国产激情久久老熟女| 国产成人一区二区在线| av国产久精品久网站免费入址| 在线看a的网站| 亚洲欧美激情在线| 在线观看免费午夜福利视频| 国产精品久久久av美女十八| 国产片内射在线| 日韩,欧美,国产一区二区三区| 亚洲国产最新在线播放| 日韩一本色道免费dvd| 999精品在线视频| 色婷婷av一区二区三区视频| 人妻 亚洲 视频| 精品少妇久久久久久888优播| 国产男女超爽视频在线观看| 亚洲精品日本国产第一区| 久久99精品国语久久久| 夫妻午夜视频| 欧美黄色片欧美黄色片| 亚洲国产欧美网| 纵有疾风起免费观看全集完整版| 在线观看www视频免费| 成人国产麻豆网| 国产色婷婷99| 久久人人爽av亚洲精品天堂| 大陆偷拍与自拍| 久久久久网色| 亚洲天堂av无毛| 日本91视频免费播放| www.熟女人妻精品国产| 日韩熟女老妇一区二区性免费视频| 亚洲精品国产av成人精品| 伊人亚洲综合成人网| 老汉色∧v一级毛片| 国产97色在线日韩免费| 1024香蕉在线观看| 国产免费一区二区三区四区乱码| 日韩人妻精品一区2区三区| 美女视频免费永久观看网站| 性少妇av在线| 热99久久久久精品小说推荐| 嫩草影院入口| 日韩av不卡免费在线播放| 18禁动态无遮挡网站| 大码成人一级视频| 99国产精品免费福利视频| 美国免费a级毛片| 国产黄色免费在线视频| 一区二区三区激情视频| 亚洲人成电影观看| 亚洲天堂av无毛| 国产精品99久久99久久久不卡 | 大码成人一级视频| 18禁观看日本| 成人国产麻豆网| 亚洲精品国产av蜜桃| 欧美人与性动交α欧美精品济南到| 亚洲成人一二三区av| 波多野结衣一区麻豆| 欧美日韩亚洲高清精品| 亚洲熟女精品中文字幕| 国产成人欧美| 少妇被粗大猛烈的视频| 久久久久久人人人人人| 精品国产乱码久久久久久小说| 欧美 日韩 精品 国产| av福利片在线| 女人被躁到高潮嗷嗷叫费观| 操出白浆在线播放| 伊人久久大香线蕉亚洲五| 亚洲国产av新网站| 精品人妻一区二区三区麻豆| 成人国产av品久久久| 欧美激情 高清一区二区三区| 成人18禁高潮啪啪吃奶动态图| 欧美激情 高清一区二区三区| 久久久久人妻精品一区果冻| 亚洲,欧美精品.| 亚洲av电影在线观看一区二区三区| 久久久欧美国产精品| 亚洲精品国产色婷婷电影| 1024香蕉在线观看| 黑人欧美特级aaaaaa片| 久久精品国产亚洲av涩爱| 日日摸夜夜添夜夜爱| 免费观看人在逋| 美女国产高潮福利片在线看| 美女大奶头黄色视频| 国产亚洲av高清不卡| 青草久久国产| 免费看av在线观看网站| 69精品国产乱码久久久| 成年av动漫网址| av片东京热男人的天堂| 午夜影院在线不卡| 欧美中文综合在线视频| 国产又色又爽无遮挡免| 黑人欧美特级aaaaaa片| 九草在线视频观看| 亚洲国产欧美网| 国产免费福利视频在线观看| 久久久久精品性色| 亚洲精品久久久久久婷婷小说| 亚洲精品中文字幕在线视频| 99热国产这里只有精品6| 亚洲美女搞黄在线观看| 伦理电影免费视频| 欧美日韩亚洲综合一区二区三区_| 久久久精品区二区三区| 中文精品一卡2卡3卡4更新| 九草在线视频观看| 亚洲欧美精品综合一区二区三区| 国产成人91sexporn| 人妻一区二区av| 国产精品嫩草影院av在线观看| 国产 一区精品| 免费少妇av软件| 中文字幕av电影在线播放| 日韩av不卡免费在线播放| 精品人妻一区二区三区麻豆| 国产成人av激情在线播放| 国产福利在线免费观看视频| 久久久久久人妻| 搡老乐熟女国产| 成年美女黄网站色视频大全免费| 国产精品一二三区在线看| 精品午夜福利在线看| 亚洲中文av在线| 亚洲精品成人av观看孕妇| 亚洲精品乱久久久久久| www.av在线官网国产| 侵犯人妻中文字幕一二三四区| 亚洲免费av在线视频| 亚洲久久久国产精品| 老司机影院成人| 亚洲欧美成人综合另类久久久| 亚洲欧美成人精品一区二区| 久久精品亚洲熟妇少妇任你| 日韩精品有码人妻一区| 宅男免费午夜| 国产精品嫩草影院av在线观看| bbb黄色大片| 黄片小视频在线播放| 国产精品嫩草影院av在线观看| 视频区图区小说| 国产伦理片在线播放av一区| 国产老妇伦熟女老妇高清| 国产麻豆69| 中文字幕人妻丝袜一区二区 | 日韩av在线免费看完整版不卡| 丁香六月天网| 亚洲,一卡二卡三卡| 这个男人来自地球电影免费观看 | 黑人欧美特级aaaaaa片| 国产女主播在线喷水免费视频网站| 久久久久久人妻| 国产精品女同一区二区软件| 欧美在线一区亚洲| 精品人妻一区二区三区麻豆| 久久青草综合色| 婷婷色麻豆天堂久久| 欧美在线黄色| 久久精品国产综合久久久| 亚洲熟女精品中文字幕| 亚洲人成网站在线观看播放|