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

    微觀調(diào)查數(shù)據(jù)抽樣權(quán)數(shù)的可忽略性檢驗(yàn)及實(shí)證研究

    2020-06-03 06:57:04米子川
    統(tǒng)計(jì)與信息論壇 2020年5期
    關(guān)鍵詞:權(quán)數(shù)估計(jì)量顯著性

    王 峰,米子川

    (山西財(cái)經(jīng)大學(xué) 統(tǒng)計(jì)學(xué)院,山西 太原 030006)

    一、引 言

    在大數(shù)據(jù)日益成為歸納性統(tǒng)計(jì)研究的主流時(shí),以小數(shù)據(jù)為特征的微觀調(diào)查數(shù)據(jù)也再次成為演繹性抽樣估計(jì)關(guān)注的焦點(diǎn)。一般認(rèn)為,大數(shù)據(jù)告訴我們“是什么”,小數(shù)據(jù)則可以告訴我們“為什么”。2018年中國國家統(tǒng)計(jì)局與清華大學(xué)共同建立了國家統(tǒng)計(jì)局-清華大學(xué)數(shù)據(jù)開發(fā)中心,對(duì)“規(guī)模以上工業(yè)企業(yè)財(cái)務(wù)狀況年度調(diào)查”“住戶收支調(diào)查”等多個(gè)微觀調(diào)查數(shù)據(jù)和普查數(shù)據(jù)進(jìn)行開發(fā)應(yīng)用研究。此外還有諸如中國家庭追蹤調(diào)查(CFPS)、中國健康與養(yǎng)老追蹤調(diào)查(CHARLS)、中國綜合社會(huì)調(diào)查(CGSS)以及中國家庭金融調(diào)查(CHFS)等等一大批有影響力的微觀調(diào)查數(shù)據(jù)。這些眾多的微觀調(diào)查數(shù)據(jù),幾乎全部是采用復(fù)雜抽樣方法得到的。在實(shí)際的調(diào)查中,由于受到調(diào)查目的、估計(jì)精度、調(diào)查費(fèi)用和可操作性等調(diào)查要求的限制,需要采用分層、整群、不等概率和多階段等多種抽樣方法相結(jié)合的抽樣設(shè)計(jì),來提高抽樣效率和樣本的代表性。我們把除單純的簡單隨機(jī)抽樣方法外,由其他多種抽樣方法組合而成的抽樣方式稱為復(fù)雜抽樣,由此得到的樣本為復(fù)雜樣本。這里的“復(fù)雜”指的是抽樣設(shè)計(jì)的特征偏離了放回簡單隨機(jī)抽樣的設(shè)計(jì)特征,而在有限總體的抽樣框中放回簡單隨機(jī)抽樣所獲得的數(shù)據(jù)就是獨(dú)立同分布的數(shù)據(jù)[1]。顯然復(fù)雜抽樣數(shù)據(jù)不是獨(dú)立同分布數(shù)據(jù),其入樣概率多為不等概率,因此需要引入抽樣權(quán)數(shù)對(duì)復(fù)雜抽樣數(shù)據(jù)做系列調(diào)整。抽樣權(quán)數(shù)中的基礎(chǔ)權(quán)數(shù),也就是單元入樣概率的倒數(shù),可用于調(diào)整由于復(fù)雜設(shè)計(jì)帶來的不等概率的偏差。更進(jìn)一步,統(tǒng)計(jì)學(xué)家在考慮無回答、抽樣框覆蓋偏差等因素時(shí),需要對(duì)抽樣權(quán)數(shù)做進(jìn)一步調(diào)整,形成最終的抽樣權(quán)數(shù)。

    很顯然,在估計(jì)總體均值或比例等一些描述性分析時(shí),抽樣權(quán)數(shù)是避免偏差必不可少的重要因素[2-4]。但是,在研究解釋變量和被解釋變量關(guān)系的模型中是否應(yīng)該考慮抽樣權(quán)數(shù),一直存有爭議[5-7]。近年來,對(duì)于抽樣權(quán)數(shù)和回歸模型也有集中性的討論[8-15]。應(yīng)用抽樣權(quán)數(shù)的優(yōu)勢是明顯的。首先能使樣本較好地代表總體,至少在一些重要特征上樣本的分布與總體的分布趨于一致,且在一般情況下,加權(quán)后的估計(jì)量是其總體參數(shù)的無偏估計(jì);其次,引入抽樣權(quán)數(shù)可以方便調(diào)整樣本信息,使樣本數(shù)據(jù)的應(yīng)用更加靈活。當(dāng)然,應(yīng)用抽樣權(quán)數(shù)進(jìn)行調(diào)整也存在明顯缺陷:首先增加了估計(jì)量的標(biāo)準(zhǔn)誤,從而導(dǎo)致估計(jì)的不穩(wěn)定性與精度的降低;其次是抽樣權(quán)數(shù)變化越大,其設(shè)計(jì)效應(yīng)(Kishs deff)也越大。因此,在沒有必要應(yīng)用抽樣權(quán)數(shù)的情況下,引入抽樣權(quán)數(shù)可能不僅不會(huì)使偏差減少,反而會(huì)造成估計(jì)量有效性的降低。另一方面,如果在應(yīng)該加入抽樣權(quán)數(shù)的情況下,沒有加入抽樣權(quán)數(shù),可能會(huì)使估計(jì)量產(chǎn)生更大偏差。

    那么,在什么情況下應(yīng)該使用抽樣權(quán)數(shù)呢?更進(jìn)一步,使用復(fù)雜調(diào)查數(shù)據(jù)建立解釋變量和被解釋變量關(guān)系模型中是否應(yīng)該使用抽樣權(quán)數(shù)?Bollen認(rèn)為大體分三類:一類主要來自生物統(tǒng)計(jì)、公共衛(wèi)生和抽樣調(diào)查領(lǐng)域的學(xué)者,一般使用抽樣權(quán)數(shù);另一類來自經(jīng)濟(jì)社會(huì)領(lǐng)域(包括計(jì)量經(jīng)濟(jì))的學(xué)者,一般不考慮抽樣權(quán)數(shù);還有一小部分學(xué)者可能用也可能不用權(quán)數(shù)[2]。Chambers和Skinner認(rèn)為目前還沒有一個(gè)明確的答案[16]。本文將從抽樣權(quán)數(shù)的低效性和抽樣權(quán)數(shù)影響的顯著性入手,在平衡兩者得失的基礎(chǔ)上,給出一個(gè)判斷是否使用權(quán)數(shù)的思路,作為解決上述問題的一個(gè)參考,這也是本文的一個(gè)創(chuàng)新之處。

    二、抽樣權(quán)數(shù)影響的檢驗(yàn)

    一般地,引入抽樣權(quán)數(shù)的估計(jì)量會(huì)滿足無偏性,但在通常情況下,其方差也會(huì)變大。Korn給出了引入抽樣權(quán)數(shù)所帶來的方差增大情況的測度指標(biāo),即抽樣權(quán)數(shù)低效性的測量,并給出了近似計(jì)算公式[17]。但是多數(shù)統(tǒng)計(jì)學(xué)家仍傾向于通過抽樣權(quán)數(shù)對(duì)模型造成的影響進(jìn)行測度,以此來決定抽樣權(quán)數(shù)的取舍。從該測度的方法角度可以分為直接測度(系數(shù)差異性的檢驗(yàn))和間接測度(抽樣權(quán)數(shù)輔助檢驗(yàn))。接下來,通過分析抽樣權(quán)數(shù)低效性和抽樣權(quán)數(shù)影響的顯著性,在其各自適用性和優(yōu)缺點(diǎn)的基礎(chǔ)上,得到抽樣權(quán)數(shù)可忽略性的檢驗(yàn)思路。

    (一)抽樣權(quán)數(shù)低效性的測量

    (1)

    (2)

    該公式來源于Korn,這里的CV是樣本權(quán)數(shù)的變異系數(shù)。低效性的檢驗(yàn)來自如下邏輯:我們當(dāng)然希望有一個(gè)無偏估計(jì)量,但是如果因?yàn)闊o偏而造成的方差(均方誤差)過大,我們還是選擇一個(gè)有偏估計(jì)量。這個(gè)低效性的計(jì)算就是測算這個(gè)方差(均方誤差)大多少,如果沒大多少,我們應(yīng)該選擇無偏估計(jì)量,即應(yīng)該使用抽樣權(quán)數(shù)。如果低效性很大了,說明加入權(quán)數(shù)后造成的方差(均方誤差)太大了,我們放棄無偏性,選擇有效性,即不應(yīng)該加入權(quán)數(shù)。

    那么到底方差(均方誤差)大多少就是太大了呢?目前還沒有一個(gè)客觀的標(biāo)準(zhǔn),但是可以從以下幾個(gè)方面去考慮:

    從低效性的絕對(duì)程度看:Korn認(rèn)為,5%到10%的低效性在實(shí)際應(yīng)用中不是很大[17]。本文認(rèn)為,考慮到當(dāng)前抽樣設(shè)計(jì)的復(fù)雜性,尤其是為降低抽樣成本以及無回答因素所帶來的抽樣設(shè)計(jì)上的復(fù)雜化,不超過30%的低效性,在實(shí)際中都會(huì)認(rèn)為不是很大。這一點(diǎn)在本文后續(xù)的國內(nèi)外數(shù)據(jù)測算中得到驗(yàn)證。

    從低效性對(duì)估計(jì)量期望精度的影響看:例如,在75%的低效性下,使估計(jì)2%的標(biāo)準(zhǔn)誤從0.01%增加一倍到了0.02%,不會(huì)認(rèn)為是不可接受的大。同樣75%的低效性下,使估計(jì)2%的標(biāo)準(zhǔn)誤從0.5%增加一倍到了1%,這一般會(huì)認(rèn)為大到了無法接受的程度[17]。

    從未加權(quán)估計(jì)量偏差的期望大小來看:如果能夠準(zhǔn)確估計(jì)未加權(quán)估計(jì)量的偏差,則可以考慮用偏差幫助選擇是否加權(quán)。即使不能從數(shù)據(jù)中足夠準(zhǔn)確地估計(jì)出偏差,若能從理論上說明對(duì)于某些類型的參數(shù),偏差可能會(huì)很小。例如,在估計(jì)兩個(gè)均值或比例的差異時(shí),偏差會(huì)比估計(jì)單個(gè)均值或比例時(shí)小。在這些情況下,即期望偏差很小時(shí),一般不會(huì)特別反對(duì)使用未加權(quán)估計(jì)量。

    總的來看,從抽樣權(quán)數(shù)低效性測量的角度來判斷是否應(yīng)該使用抽樣權(quán)數(shù),還是需要研究者的主觀判斷,不同的使用者可能會(huì)有不同的結(jié)果;另外,在實(shí)際的應(yīng)用中,出于計(jì)算方便,一般采用近似計(jì)算,在某些情況下可能存在誤導(dǎo)(參見例子[17])。但是,低效性測量的優(yōu)勢也是明顯的,首先低效性測量不依賴于所采用的推斷方法或模型,這樣適用性就很強(qiáng);另外低效性的計(jì)算只依賴于抽樣權(quán)數(shù)的分布且計(jì)算簡單方便。

    以上是根據(jù)抽樣權(quán)數(shù)對(duì)估計(jì)量方差(均方誤差)的影響,來判斷是否使用抽樣權(quán)數(shù)。另外一個(gè)思路是通過比較回歸模型中加權(quán)和未加權(quán)的系數(shù)是否有顯著差異來判斷是否需要加權(quán),或者將權(quán)數(shù)置于模型中判斷其是否顯著。

    (二)抽樣權(quán)數(shù)影響的顯著性測量

    2.抽樣權(quán)數(shù)輔助檢驗(yàn)(Weight Association Tests)。Hausman曾建議評(píng)估其他形式的方程Y=Xβ+XMβM+ε中βM的顯著性來判斷設(shè)定誤差[18]。這里的XM是X的一個(gè)合適變形。用F檢驗(yàn)H0:βM=0來檢驗(yàn)設(shè)定誤差。另外,與一般多元回歸的假定類似,使用F檢驗(yàn)需要我們假定ε來自于正態(tài)分布。盡管Hausman建議用這樣一個(gè)形式去檢驗(yàn)各種設(shè)定誤差,但他并沒有考慮用它做加權(quán)檢驗(yàn)。

    Dumouchel等利用Hausman的回歸方法并將其用于確定是否加權(quán)[20-21]。再次考慮方程Y=Xβu+Xwβw+ε,這里Y是被解釋變量的向量值,X是未加權(quán)的解釋變量陣對(duì)應(yīng)的系數(shù)為βu,Xw是相同解釋變量陣對(duì)應(yīng)的加權(quán)系數(shù)為βw,ε是誤差向量。Dumouchel and Duncan建議用OLS估計(jì)這個(gè)回歸模型然后用F檢驗(yàn)H0:βM=0去檢驗(yàn)是否需要權(quán)數(shù),拒絕原假設(shè)意味著需要加權(quán)。可以看出盡管Hausman并沒有這么做,作為權(quán)數(shù)輔助檢驗(yàn)(WA檢驗(yàn))的Dumouchel和Duncan和Fuller的F檢驗(yàn)回歸法都是遵循Hausman替代基于回歸的設(shè)定誤差檢驗(yàn)。Fuller將這一方法做了一個(gè)改變,推薦用回歸Y=Xβu+Wα+ε,這里W是權(quán)數(shù)變量,檢驗(yàn)其系數(shù)α是否顯著非零。

    三、抽樣權(quán)數(shù)可忽略性檢驗(yàn)的操作路徑

    在得到調(diào)查數(shù)據(jù)的抽樣權(quán)數(shù)后,其分布就容易直觀得到。加之低效性的判斷不依賴于所采用的模型方法,因此在判斷抽樣權(quán)數(shù)是否可以忽略時(shí),首先計(jì)算抽樣權(quán)數(shù)的低效性,并根據(jù)前述內(nèi)容做出低效性是否可以接受的判斷。

    為給出低效性的一個(gè)一般經(jīng)驗(yàn)水平,這里給出了一些常見分布下的抽樣權(quán)數(shù)頻數(shù)分布圖。樣本量取1 000,見圖1,從上往下依次為:均勻分布 [U(a=1,b=1 000)] 、正態(tài)分布 [N(μ=1 000,σ2=10)] 、二項(xiàng)分布 [B(n=50,p=0.25)] 、指數(shù)分布 [e(λ=0.1)] 、卡方分布 [χ2(3)] 、Possion分布 [P(λ=4)]抽樣權(quán)數(shù)的頻數(shù)分布圖。因?yàn)檫@里都是權(quán)數(shù)所以不是整數(shù)的全部取最大整數(shù)。

    圖1常見分布下的權(quán)數(shù)分布圖

    考慮到常數(shù)分布很簡單,為一條水平直線,這里沒有畫出。對(duì)于常數(shù)分布,樣本權(quán)數(shù)的低效性顯然為零。因?yàn)榇藭r(shí)的樣本就是等概率抽選,也可以認(rèn)為樣本是自加權(quán)的,不存在因?yàn)闄?quán)數(shù)的引入而造成估計(jì)量標(biāo)準(zhǔn)誤的增大。因此,抽樣權(quán)數(shù)的低效性為零,公式的計(jì)算結(jié)果也為零。可以認(rèn)為,在不考慮抽樣權(quán)數(shù)的調(diào)查數(shù)據(jù)分析時(shí),都是默認(rèn)抽樣權(quán)數(shù)為常數(shù),這顯然可能會(huì)低估估計(jì)量的標(biāo)準(zhǔn)誤,更主要的是通常會(huì)失去估計(jì)量無偏性。一般情況下,抽樣權(quán)數(shù)的分布不會(huì)是常數(shù),除樣本是自加權(quán)樣本外,均會(huì)產(chǎn)生抽樣權(quán)數(shù)低效性的發(fā)生,這里測算了以上幾種常見參數(shù)分布的抽樣權(quán)數(shù)的低效性,見表1。

    表1 抽樣權(quán)數(shù)在常見分布下的低效性

    ① 這里需要說明的是“再抽樣權(quán)數(shù)”與抽樣權(quán)數(shù)一樣僅僅是一個(gè)公布的權(quán)數(shù)。因?yàn)橹袊彝プ粉櫿{(diào)查(CFPS)抽樣設(shè)計(jì)的復(fù)雜性,調(diào)查數(shù)據(jù)清理和加權(quán)調(diào)整也異常復(fù)雜,耗時(shí)較長,無法滿足研究者對(duì)調(diào)查數(shù)據(jù)進(jìn)行快速分析的需求。因此,在中國家庭追蹤調(diào)查的數(shù)據(jù)庫中包含了一個(gè)再整合數(shù)據(jù)庫,即對(duì)5個(gè)“大省”進(jìn)行再抽樣調(diào)整樣本,使5個(gè)“大省”的抽樣比與“小省”的抽樣比近似相同,以便在沒有及時(shí)獲得抽樣權(quán)數(shù)的情況下利用再抽樣調(diào)查權(quán)數(shù)去推斷總體。

    以上圖形可以從便捷的角度給出實(shí)踐中抽樣權(quán)數(shù)分布的直觀認(rèn)知。即實(shí)際中可以很方便的畫出抽樣權(quán)數(shù)的頻數(shù)分布圖,如果這個(gè)分布圖與上述的某個(gè)圖形類似,就可以相應(yīng)的得到抽樣權(quán)數(shù)低效性的一般認(rèn)知。

    由于抽樣權(quán)數(shù)的分布圖形只是對(duì)抽樣權(quán)數(shù)低效性判斷的一般認(rèn)知,具有很強(qiáng)的主觀性,多數(shù)情況下甚至是無法做出判斷,因此這里綜合抽樣權(quán)數(shù)的低效性和抽樣權(quán)數(shù)影響的顯著性,提出一個(gè)檢驗(yàn)抽樣權(quán)數(shù)可忽略性的判斷路徑,見圖2。

    當(dāng)抽樣權(quán)數(shù)的分布為正態(tài)分布[N(μ=1 000,σ2= 10) ]、二項(xiàng)分布[B(n= 50,p= 0.25) ]和Possion分布[P(λ=4)]時(shí),抽樣權(quán)數(shù)的低效性比較小,也就是說考慮抽樣權(quán)數(shù)對(duì)估計(jì)量標(biāo)準(zhǔn)誤的影響很小,甚至可以忽略不計(jì),這時(shí)候應(yīng)該使用抽樣權(quán)數(shù),即在保證估計(jì)量無偏的情況下,不會(huì)對(duì)標(biāo)準(zhǔn)誤產(chǎn)生很大影響;如果抽樣權(quán)數(shù)的分布為指數(shù)分布[e(λ=0.1)]時(shí),此時(shí)的抽樣權(quán)數(shù)無效性達(dá)到了44.6%在這種情況下就要格外謹(jǐn)慎,雖然考慮抽樣權(quán)數(shù)能使估計(jì)量無偏,但此時(shí)估計(jì)量的穩(wěn)健性會(huì)較差。通常會(huì)選擇穩(wěn)健性,而放棄考慮抽樣權(quán)數(shù)。當(dāng)然也可以選擇重新構(gòu)建研究的子總體,獲取新的抽樣權(quán)數(shù)重新判斷; 如果抽樣分布類似均勻分布[U(a=1,b=1 000) ]和卡方分布[χ2(3)]時(shí),低效性達(dá)到了24.5%和25%,此時(shí)考慮抽樣權(quán)數(shù)可能會(huì)影響到估計(jì)量的標(biāo)準(zhǔn)誤,至于影響的大小是否能夠接受,可以用上述假設(shè)檢驗(yàn)的方法去檢驗(yàn)。但為什么不直接用該假設(shè)檢驗(yàn)來判斷呢?原因有二:一是低效性的判斷方便快捷不依賴于所采用的方法;二是假設(shè)檢驗(yàn)的判斷,目前還局限于回歸模型的分析中,雖然對(duì)其他模型的構(gòu)建可以提供參考,但當(dāng)前還沒有證明該假設(shè)檢驗(yàn)可以應(yīng)用于其他一切模型。

    圖2 抽樣權(quán)數(shù)可忽略性的判斷路徑圖

    四、實(shí)證研究

    CFPS是北京大學(xué)中國社會(huì)科學(xué)調(diào)查中心主持的追蹤調(diào)查項(xiàng)目之一,調(diào)查對(duì)象是中國的25個(gè)省市自治區(qū)的家庭戶和家庭戶中的所有滿足調(diào)查條件的家庭成員。在抽樣設(shè)計(jì)上采用三階段、不等概率的整群抽樣設(shè)計(jì)[22]。本文以CFPS2010的數(shù)據(jù)為例,來說明抽樣權(quán)數(shù)可忽略性檢驗(yàn)。

    (一)低效性檢驗(yàn)

    首先利用抽樣權(quán)數(shù)的分布圖做一個(gè)直觀的分析,圖3中左側(cè)圖形為中國家庭抽樣權(quán)數(shù)分布圖右側(cè)為中國家庭再抽樣權(quán)數(shù)①分布圖。這兩個(gè)圖形都與常數(shù)分布差的很遠(yuǎn),與圖1中正態(tài)分布[N(μ=1 000,σ2=10)]、二項(xiàng)分布[B(n=50,p=0.25)]相差也比較大,與卡方分布[χ2(3)]類似,因此其低效性都不會(huì)很小。對(duì)比兩個(gè)圖形,再抽樣權(quán)數(shù)的低效性要大于抽樣權(quán)數(shù)低效性。

    圖3 全國家庭抽樣和再抽樣權(quán)數(shù)分布圖

    這里根據(jù)式(2),對(duì)CFPS2010的抽樣調(diào)查數(shù)據(jù)中的家庭權(quán)數(shù)和家庭再抽樣權(quán)數(shù)的低效性做進(jìn)一步分類測算,結(jié)果見表2。可以看出中國家庭的抽樣權(quán)數(shù)的低效性為26.208%;其中城鎮(zhèn)家庭抽樣權(quán)數(shù)的低效性低于全國水平為25.422%;農(nóng)村家庭抽樣權(quán)數(shù)的低效性高于全國水平為26.627%。同時(shí)也可以看出,家庭再抽樣的權(quán)數(shù)低效性較高,全國、城鎮(zhèn)和農(nóng)村家庭的再抽樣權(quán)數(shù)低效性分別為34.537%、34.676%和33.702%。無論是抽樣權(quán)數(shù)還是再抽樣的權(quán)數(shù),低效性雖然不小,但也沒有大到無法接受的程度。根據(jù)筆者掌握的國外抽樣權(quán)數(shù)低效性的測算結(jié)果:美國國家健康和營養(yǎng)調(diào)查(NHANES)中,NHANESII中2~18歲男孩抽樣權(quán)數(shù)的低效性是34%;NHANESI中25~74歲女性抽樣權(quán)數(shù)的低效性是48%;美國國民健康訪問調(diào)查(NHIS1987)18歲以上男性的抽樣權(quán)數(shù)的低效性是22%。為便于與國外比較,我們也計(jì)算了CFPS2010中25~74歲女性抽樣權(quán)數(shù)的低效性是26.451%,18歲以上男性的抽樣權(quán)數(shù)的低效性是26.131%,前者比美國同一指標(biāo)低,后者略高。

    表2 農(nóng)村、城鎮(zhèn)及全國家庭抽樣權(quán)數(shù)的低效性單位:%

    ① 對(duì)于缺失數(shù)據(jù)可以考慮多種方法去插補(bǔ),由于這些缺失值對(duì)權(quán)數(shù)的分布影響不大,因此本文直接刪除處理;對(duì)于家庭消費(fèi)小于100元的家庭,CFPS用戶手冊(cè)中說明對(duì)其做過插補(bǔ)處理,但數(shù)據(jù)中仍然有6個(gè)家庭的消費(fèi)小于100元,這里也刪除處理。

    為了做進(jìn)一步檢驗(yàn),我們計(jì)算了家庭消費(fèi)和家庭收入調(diào)整后的抽樣權(quán)數(shù)的低效性。這里面的調(diào)整主要包括刪除了家庭消費(fèi)和家庭收入中的缺失數(shù)據(jù),也刪除了家庭消費(fèi)在100元以下的數(shù)據(jù)①。經(jīng)過這樣的調(diào)整,抽樣權(quán)數(shù)略微有些變化,見表2。調(diào)整后城鎮(zhèn)家庭抽樣權(quán)數(shù)的低效性均高于農(nóng)村家庭抽樣權(quán)數(shù)的低效性,但低效性的變化都不是很大。

    接下來我們測算了不同規(guī)模家庭的抽樣權(quán)數(shù)的低效性。根據(jù)家庭成員數(shù)將家庭規(guī)模劃分為:家1為1位家庭成員、家2為2位家庭成員,以此得到家6及以上為家庭成員數(shù)為6位及以上。從表3可以看出,雖然不同家庭規(guī)模的樣本量差異較大,一位家庭成員的家庭樣本數(shù)為728,三位家庭成員的家庭樣本數(shù)為33 556,但它們的抽樣權(quán)數(shù)的低效性都不是非常大,都沒有超過30%。

    表3 按家庭人口數(shù)劃分的抽樣權(quán)數(shù)的低效性

    因此,無論是家庭抽樣權(quán)數(shù)還是成人抽樣權(quán)數(shù),以及不同規(guī)模家庭的抽樣權(quán)數(shù),低效性都不是特別的大。因此,我們?cè)谑褂谜{(diào)查數(shù)據(jù)時(shí),應(yīng)該考慮抽樣權(quán)數(shù),這樣既可以得到一個(gè)無偏估計(jì)量,同時(shí)其有效性也不會(huì)降低很多。

    (二)抽樣權(quán)數(shù)影響的顯著性水平檢驗(yàn)

    由前面的分析可知,在考慮家庭收入和消費(fèi)時(shí),由于刪除了缺失數(shù)據(jù)和家庭消費(fèi)小于100元的數(shù)據(jù)后,抽樣權(quán)數(shù)的分布有了略微的變化。但抽樣權(quán)數(shù)的低效性依然不是很大,因此從低效性的角度看,我們應(yīng)該考慮抽樣權(quán)數(shù)。為對(duì)上述思路與路徑做完整的分析,接下來從抽樣權(quán)數(shù)影響的顯著性水平來檢驗(yàn)抽樣權(quán)數(shù)是否可以忽略。即假設(shè)認(rèn)為抽樣權(quán)數(shù)低效性的大小無法判斷,則通過假設(shè)檢驗(yàn)來判斷。

    這里利用權(quán)數(shù)輔助檢驗(yàn)的方法,選取Dumouchel和Duncan所述方法對(duì)上述數(shù)據(jù)的抽樣權(quán)數(shù)是否應(yīng)該引入模型進(jìn)行檢驗(yàn)。Dumouchel和Duncan是將權(quán)數(shù)及權(quán)數(shù)與解釋變量的交互作用引入模型,利用F檢驗(yàn)其系數(shù)是否全部為零來判斷抽樣權(quán)數(shù)是否可以忽略。因?yàn)闄?quán)數(shù)輔助檢驗(yàn)需要依據(jù)模型來檢驗(yàn),由前述數(shù)據(jù),這里對(duì)家庭消費(fèi)和支出建立線性模型來檢驗(yàn)抽樣權(quán)數(shù)的可忽略性。

    表4 家庭支出關(guān)于家庭收入的線性模型比較

    注:***表示在1%水平上顯著。

    為全面展示抽樣權(quán)數(shù)可忽略性的判斷路徑圖,這里給出了另外一種情況,即抽樣權(quán)數(shù)是可忽略的。依然以CFPS2010的數(shù)據(jù)為例,選取cyear指標(biāo)為抽樣權(quán)數(shù)。這其實(shí)相當(dāng)于認(rèn)為CFPS2010的數(shù)據(jù)是等概率抽選的,這也是在實(shí)際中,微觀調(diào)查數(shù)據(jù)建模時(shí)經(jīng)常默認(rèn)的,即數(shù)據(jù)是等概率抽樣獲取的,因此抽樣權(quán)數(shù)可以忽略。利用Dumouchel和Duncan所建議的檢驗(yàn)方法,計(jì)算得到F(2,8 188) = 1.411,P=0.244,在0.1的顯著性水平下,我們不能拒絕,由此得出cyear作為權(quán)數(shù)是可忽略的。

    然而事實(shí)是這些數(shù)據(jù)并不是等概率獲取的,所以說,在給微觀調(diào)查數(shù)據(jù)建模的時(shí)候,如果對(duì)抽樣權(quán)數(shù)視而不見,其所建立的模型就值得懷疑,因?yàn)槎鄶?shù)大型、規(guī)范的調(diào)查數(shù)據(jù)很少是等概率抽取樣本,同時(shí)還要考慮因?yàn)闊o回答和抽樣框覆蓋偏差等因素所帶來的“過抽樣”和“欠抽樣”等問題。當(dāng)然不是說對(duì)于這類調(diào)查數(shù)據(jù)就必須在模型中考慮抽樣權(quán)數(shù),如果因?yàn)橐霗?quán)數(shù)帶來方差的增大,足以影響到數(shù)據(jù)分析,或者假設(shè)檢驗(yàn)沒有被拒絕,這時(shí)候可以如圖2所示,返回重新選擇研究的數(shù)據(jù)或者選擇不考慮抽樣權(quán)數(shù)處理。但無論怎么樣,對(duì)抽樣權(quán)數(shù)總歸要有一個(gè)處理,用要有用的理由,不用也要有個(gè)說明,這是我們的寫作動(dòng)機(jī)之一,本文就此給出了解決該問題的一個(gè)路徑。

    五、結(jié)論與建議

    權(quán)數(shù)是把雙刃劍,在提高精度的同時(shí),有可能提高估計(jì)量的誤差[15,23]。因此,在使用微觀調(diào)查數(shù)據(jù)時(shí),是否使用抽樣權(quán)數(shù),本文給出了一個(gè)判斷思路。即從抽樣權(quán)數(shù)的低效性和抽樣權(quán)數(shù)影響的顯著性兩方面入手,給出了一個(gè)判斷路徑圖。第一,通過分析抽樣權(quán)數(shù)低效性的計(jì)算過程,得出了抽樣權(quán)數(shù)低效性判斷的優(yōu)缺點(diǎn),并給出了抽樣權(quán)數(shù)在幾種常見分布下,低效性的一個(gè)一般經(jīng)驗(yàn)判斷。第二,通過對(duì)抽樣權(quán)數(shù)影響的顯著性檢驗(yàn)方法的綜述,根據(jù)方法的不同,分為系數(shù)差異性檢驗(yàn)和權(quán)數(shù)輔助檢驗(yàn),并得出這些檢驗(yàn)與Hausman的模型設(shè)定誤差檢驗(yàn)的關(guān)系。最后,通過CFPS2010驗(yàn)證了本文提出的檢驗(yàn)路徑,并給出了CFPS2010家庭抽樣調(diào)查數(shù)據(jù)抽樣權(quán)數(shù)的低效性,與國外類似調(diào)查數(shù)據(jù)的低效性進(jìn)行了比較。通過低效性判斷和Dumouchel-Duncan檢驗(yàn)說明其抽樣權(quán)數(shù)的不可忽略性,即得到在使用CFPS2010家庭抽樣調(diào)查數(shù)據(jù)做統(tǒng)計(jì)推斷或數(shù)據(jù)建模時(shí),應(yīng)當(dāng)考慮抽樣權(quán)數(shù)。當(dāng)然,本文所運(yùn)用的抽樣權(quán)數(shù)影響的顯著性檢驗(yàn)還局限于回歸分析,如何將其推廣到其他模型,比如結(jié)構(gòu)方程模型中,還需要做進(jìn)一步研究;文中抽樣權(quán)數(shù)低效性大小的判斷,還缺乏客觀統(tǒng)一的標(biāo)準(zhǔn),需要更多經(jīng)驗(yàn)數(shù)據(jù)的積累和方法的創(chuàng)新。

    猜你喜歡
    權(quán)數(shù)估計(jì)量顯著性
    豬肉在CPI中的權(quán)數(shù)被調(diào)低了嗎?
    豬肉在CPI中的權(quán)數(shù)被調(diào)低了嗎?
    我國CPI權(quán)數(shù)調(diào)整中存在的問題及對(duì)策思考
    基于顯著性權(quán)重融合的圖像拼接算法
    電子制作(2019年24期)2019-02-23 13:22:26
    權(quán)數(shù)可靠性的假設(shè)檢驗(yàn)探討
    基于視覺顯著性的視頻差錯(cuò)掩蓋算法
    一種基于顯著性邊緣的運(yùn)動(dòng)模糊圖像復(fù)原方法
    淺談估計(jì)量的優(yōu)良性標(biāo)準(zhǔn)
    論商標(biāo)固有顯著性的認(rèn)定
    基于配網(wǎng)先驗(yàn)信息的諧波狀態(tài)估計(jì)量測點(diǎn)最優(yōu)配置
    電測與儀表(2015年6期)2015-04-09 12:00:50
    久久精品成人免费网站| 狂野欧美激情性xxxx| 91麻豆精品激情在线观看国产| 美国免费a级毛片| 久久午夜综合久久蜜桃| 欧美丝袜亚洲另类 | 国产亚洲欧美在线一区二区| 婷婷六月久久综合丁香| 好男人电影高清在线观看| 国产乱人伦免费视频| 亚洲精品久久国产高清桃花| 国产蜜桃级精品一区二区三区| 777久久人妻少妇嫩草av网站| 欧美久久黑人一区二区| 亚洲色图 男人天堂 中文字幕| 99久久无色码亚洲精品果冻| 黑丝袜美女国产一区| 90打野战视频偷拍视频| 女生性感内裤真人,穿戴方法视频| 国产私拍福利视频在线观看| 国产精品久久久av美女十八| 自线自在国产av| 色播在线永久视频| 中文字幕人妻熟女乱码| 亚洲中文av在线| 97人妻精品一区二区三区麻豆 | 久久久久亚洲av毛片大全| 国产亚洲精品久久久久久毛片| 丁香六月欧美| 50天的宝宝边吃奶边哭怎么回事| 成人国语在线视频| 两人在一起打扑克的视频| bbb黄色大片| or卡值多少钱| 每晚都被弄得嗷嗷叫到高潮| 一本综合久久免费| 在线av久久热| 国产单亲对白刺激| 在线看三级毛片| 欧美性猛交╳xxx乱大交人| 草草在线视频免费看| 久久国产精品男人的天堂亚洲| 99re在线观看精品视频| 亚洲 国产 在线| 精品人妻1区二区| 嫩草影视91久久| 大型黄色视频在线免费观看| 免费在线观看影片大全网站| 国产99白浆流出| 国产精品野战在线观看| 亚洲成人久久爱视频| 男女之事视频高清在线观看| 国产成年人精品一区二区| 97人妻精品一区二区三区麻豆 | 午夜亚洲福利在线播放| 久久精品影院6| 长腿黑丝高跟| 宅男免费午夜| 黄片小视频在线播放| 免费av毛片视频| 亚洲成av片中文字幕在线观看| 在线观看免费视频日本深夜| 黑人欧美特级aaaaaa片| 两人在一起打扑克的视频| 两性夫妻黄色片| 欧美日韩亚洲综合一区二区三区_| 欧美日本亚洲视频在线播放| 99国产极品粉嫩在线观看| 久久人人精品亚洲av| 国产三级在线视频| 日韩精品青青久久久久久| 999久久久国产精品视频| 老汉色∧v一级毛片| 日本免费一区二区三区高清不卡| 久久 成人 亚洲| 亚洲欧美精品综合久久99| 一级作爱视频免费观看| 中亚洲国语对白在线视频| 日韩视频一区二区在线观看| 一个人免费在线观看的高清视频| 国产av不卡久久| 亚洲中文av在线| 露出奶头的视频| 90打野战视频偷拍视频| 亚洲九九香蕉| 琪琪午夜伦伦电影理论片6080| 亚洲av成人av| 久热爱精品视频在线9| 脱女人内裤的视频| 精品久久久久久久久久久久久 | 99热这里只有精品一区 | 欧美精品亚洲一区二区| 日韩精品中文字幕看吧| 日本熟妇午夜| 成年版毛片免费区| 日韩精品中文字幕看吧| 白带黄色成豆腐渣| 女同久久另类99精品国产91| 757午夜福利合集在线观看| 国产aⅴ精品一区二区三区波| 国产区一区二久久| 日韩有码中文字幕| 大香蕉久久成人网| 亚洲精品久久国产高清桃花| www.www免费av| 欧美人与性动交α欧美精品济南到| 欧美日韩福利视频一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| bbb黄色大片| 国产97色在线日韩免费| 黄色丝袜av网址大全| 丁香六月欧美| 在线国产一区二区在线| 91麻豆av在线| 久久久久国内视频| 亚洲熟妇熟女久久| 国产蜜桃级精品一区二区三区| 久久精品人妻少妇| 此物有八面人人有两片| 日韩 欧美 亚洲 中文字幕| 宅男免费午夜| 色综合站精品国产| 亚洲久久久国产精品| 女生性感内裤真人,穿戴方法视频| 国产精品自产拍在线观看55亚洲| 日日干狠狠操夜夜爽| 99久久无色码亚洲精品果冻| 他把我摸到了高潮在线观看| 精品国产亚洲在线| 国产乱人伦免费视频| 日日夜夜操网爽| 精品一区二区三区av网在线观看| 在线av久久热| 女人被狂操c到高潮| 桃色一区二区三区在线观看| 欧美丝袜亚洲另类 | 日韩欧美在线二视频| 欧美一级a爱片免费观看看 | 9191精品国产免费久久| 两人在一起打扑克的视频| 麻豆国产av国片精品| 久久久国产成人免费| 国产欧美日韩一区二区精品| 久久中文字幕人妻熟女| 亚洲av五月六月丁香网| 99久久久亚洲精品蜜臀av| 美女免费视频网站| 男女床上黄色一级片免费看| 99在线人妻在线中文字幕| 国产成人精品无人区| 两性夫妻黄色片| av片东京热男人的天堂| 中文在线观看免费www的网站 | 黄片大片在线免费观看| 黑丝袜美女国产一区| 久久中文字幕一级| 啪啪无遮挡十八禁网站| 久久天躁狠狠躁夜夜2o2o| 男人操女人黄网站| 久久久久国产一级毛片高清牌| 正在播放国产对白刺激| 法律面前人人平等表现在哪些方面| 91老司机精品| www日本在线高清视频| 在线观看免费视频日本深夜| 亚洲av成人不卡在线观看播放网| 啦啦啦免费观看视频1| 欧美日韩黄片免| 日韩中文字幕欧美一区二区| 精品国产亚洲在线| 日韩欧美 国产精品| 午夜视频精品福利| 国产精品久久久av美女十八| 99热只有精品国产| 久久中文字幕人妻熟女| 一区二区三区精品91| aaaaa片日本免费| 国内久久婷婷六月综合欲色啪| 久久国产亚洲av麻豆专区| 美女国产高潮福利片在线看| 天天躁夜夜躁狠狠躁躁| 国产蜜桃级精品一区二区三区| 国产野战对白在线观看| 午夜影院日韩av| 啦啦啦免费观看视频1| 人人妻人人澡欧美一区二区| 国产一区在线观看成人免费| 国产免费av片在线观看野外av| 亚洲无线在线观看| 成人午夜高清在线视频 | 美女国产高潮福利片在线看| 男人舔女人的私密视频| 国产蜜桃级精品一区二区三区| 亚洲狠狠婷婷综合久久图片| 久久久久国内视频| 美女 人体艺术 gogo| 欧美成人一区二区免费高清观看 | 麻豆一二三区av精品| 超碰成人久久| 亚洲五月色婷婷综合| 99久久久亚洲精品蜜臀av| 亚洲片人在线观看| 1024香蕉在线观看| 免费观看人在逋| 可以在线观看的亚洲视频| svipshipincom国产片| 91九色精品人成在线观看| 欧美三级亚洲精品| 国产又爽黄色视频| 哪里可以看免费的av片| 日韩欧美一区二区三区在线观看| 亚洲国产精品sss在线观看| 国产又爽黄色视频| 亚洲av美国av| 91麻豆精品激情在线观看国产| 国产精品九九99| 亚洲精华国产精华精| aaaaa片日本免费| 精品一区二区三区av网在线观看| 女人爽到高潮嗷嗷叫在线视频| 12—13女人毛片做爰片一| 中文字幕人成人乱码亚洲影| 19禁男女啪啪无遮挡网站| 99riav亚洲国产免费| 成人三级做爰电影| 亚洲国产日韩欧美精品在线观看 | 久久久久精品国产欧美久久久| 99国产精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 国产乱人伦免费视频| www日本在线高清视频| 好看av亚洲va欧美ⅴa在| 免费一级毛片在线播放高清视频| 精品久久久久久久久久免费视频| 岛国在线观看网站| 午夜福利高清视频| 久久精品国产亚洲av香蕉五月| 国产精品自产拍在线观看55亚洲| a在线观看视频网站| 亚洲成国产人片在线观看| 麻豆一二三区av精品| 亚洲一区二区三区色噜噜| 国产精品久久久人人做人人爽| 亚洲精品美女久久av网站| 91麻豆精品激情在线观看国产| 日韩欧美国产一区二区入口| 国产成人欧美| 亚洲国产看品久久| 久久中文字幕一级| 精品熟女少妇八av免费久了| 一级毛片女人18水好多| 老汉色∧v一级毛片| 黄片小视频在线播放| 亚洲精品中文字幕一二三四区| 国产成人一区二区三区免费视频网站| 久久 成人 亚洲| 日日摸夜夜添夜夜添小说| 美国免费a级毛片| 亚洲国产高清在线一区二区三 | 黄色女人牲交| 19禁男女啪啪无遮挡网站| 久久精品91无色码中文字幕| 99久久综合精品五月天人人| 国产成+人综合+亚洲专区| 色哟哟哟哟哟哟| 久久久久久久久久黄片| 又大又爽又粗| 人妻久久中文字幕网| 熟妇人妻久久中文字幕3abv| 国产男靠女视频免费网站| 黑丝袜美女国产一区| 国产真实乱freesex| 欧美色欧美亚洲另类二区| www.精华液| 在线国产一区二区在线| 老司机在亚洲福利影院| 精品电影一区二区在线| 国产真人三级小视频在线观看| 在线天堂中文资源库| 男人操女人黄网站| 黄片播放在线免费| 欧美日韩黄片免| 国产欧美日韩精品亚洲av| 俺也久久电影网| 91大片在线观看| 一区二区三区激情视频| 久久天堂一区二区三区四区| 两性夫妻黄色片| 最新在线观看一区二区三区| 日本一本二区三区精品| 老司机福利观看| 成人欧美大片| 亚洲一区中文字幕在线| 亚洲全国av大片| 美女高潮喷水抽搐中文字幕| av超薄肉色丝袜交足视频| 高清在线国产一区| 女同久久另类99精品国产91| 国产私拍福利视频在线观看| 久久久久国内视频| 啦啦啦韩国在线观看视频| 麻豆久久精品国产亚洲av| 免费高清视频大片| 午夜精品在线福利| 亚洲人成网站高清观看| 日韩欧美 国产精品| 国产黄a三级三级三级人| 亚洲成人久久爱视频| 国产单亲对白刺激| 亚洲精品av麻豆狂野| 国产视频一区二区在线看| 美女午夜性视频免费| 国产精品精品国产色婷婷| 国产精品精品国产色婷婷| 香蕉国产在线看| 免费一级毛片在线播放高清视频| 久久精品人妻少妇| av在线天堂中文字幕| xxxwww97欧美| 国产精品乱码一区二三区的特点| 丁香欧美五月| 大型av网站在线播放| 欧美日韩亚洲国产一区二区在线观看| 欧美乱色亚洲激情| 91成人精品电影| 亚洲精品粉嫩美女一区| 男女下面进入的视频免费午夜 | 亚洲aⅴ乱码一区二区在线播放 | 99riav亚洲国产免费| 成年女人毛片免费观看观看9| 欧美 亚洲 国产 日韩一| 亚洲精品一区av在线观看| 制服丝袜大香蕉在线| 成熟少妇高潮喷水视频| 久热爱精品视频在线9| 草草在线视频免费看| 欧美黑人巨大hd| 国产久久久一区二区三区| 老司机深夜福利视频在线观看| 久久久国产欧美日韩av| 满18在线观看网站| www.精华液| 亚洲电影在线观看av| 亚洲成av片中文字幕在线观看| 久久中文看片网| 69av精品久久久久久| 精品欧美一区二区三区在线| 美国免费a级毛片| 精品人妻1区二区| 国产精品自产拍在线观看55亚洲| 热99re8久久精品国产| 中文资源天堂在线| 国产野战对白在线观看| 精品国产乱码久久久久久男人| 成熟少妇高潮喷水视频| 国产视频一区二区在线看| av欧美777| 18禁裸乳无遮挡免费网站照片 | 国产亚洲精品一区二区www| 在线视频色国产色| 中文字幕另类日韩欧美亚洲嫩草| www国产在线视频色| 久久午夜综合久久蜜桃| 中文字幕av电影在线播放| 人人妻人人澡欧美一区二区| 久久久久九九精品影院| 精品久久久久久久人妻蜜臀av| 国产av一区二区精品久久| 免费在线观看黄色视频的| 观看免费一级毛片| 亚洲av五月六月丁香网| 久久亚洲精品不卡| 成人午夜高清在线视频 | 色哟哟哟哟哟哟| av在线天堂中文字幕| 色尼玛亚洲综合影院| 每晚都被弄得嗷嗷叫到高潮| 亚洲狠狠婷婷综合久久图片| 亚洲第一欧美日韩一区二区三区| 少妇 在线观看| 18禁国产床啪视频网站| 桃红色精品国产亚洲av| 午夜免费鲁丝| 人人妻人人澡欧美一区二区| 久久亚洲真实| 免费观看精品视频网站| 国产亚洲精品久久久久久毛片| 国产精品二区激情视频| 香蕉丝袜av| 亚洲成人久久爱视频| 一级毛片精品| 亚洲全国av大片| 亚洲性夜色夜夜综合| 成人三级做爰电影| 一二三四在线观看免费中文在| 亚洲av成人av| 国产精品亚洲av一区麻豆| 午夜两性在线视频| 午夜福利免费观看在线| 欧美一级毛片孕妇| 亚洲一区二区三区色噜噜| 听说在线观看完整版免费高清| 国产av一区在线观看免费| 在线观看免费视频日本深夜| 淫秽高清视频在线观看| 亚洲第一电影网av| 人人妻,人人澡人人爽秒播| 欧美日韩亚洲综合一区二区三区_| 中文在线观看免费www的网站 | 精品高清国产在线一区| 色哟哟哟哟哟哟| 90打野战视频偷拍视频| 欧美中文日本在线观看视频| 欧美激情久久久久久爽电影| 国产精品免费一区二区三区在线| 99热6这里只有精品| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品美女久久av网站| 成人亚洲精品av一区二区| 999精品在线视频| 亚洲av第一区精品v没综合| 50天的宝宝边吃奶边哭怎么回事| 97人妻精品一区二区三区麻豆 | 99久久精品国产亚洲精品| 国产精品一区二区精品视频观看| 美女高潮到喷水免费观看| 午夜福利成人在线免费观看| 麻豆久久精品国产亚洲av| 亚洲国产欧美一区二区综合| 精品久久久久久成人av| 亚洲欧美精品综合一区二区三区| 宅男免费午夜| 午夜老司机福利片| 757午夜福利合集在线观看| 国产私拍福利视频在线观看| 黄色片一级片一级黄色片| 亚洲一区中文字幕在线| 一二三四社区在线视频社区8| 中出人妻视频一区二区| 亚洲国产欧美日韩在线播放| 满18在线观看网站| 日日摸夜夜添夜夜添小说| 亚洲国产看品久久| 在线观看日韩欧美| 女性被躁到高潮视频| 波多野结衣高清作品| 女人被狂操c到高潮| 熟女电影av网| 亚洲av第一区精品v没综合| 看免费av毛片| 久久婷婷人人爽人人干人人爱| 69av精品久久久久久| 国产欧美日韩精品亚洲av| x7x7x7水蜜桃| 国产亚洲精品久久久久久毛片| 国产精品精品国产色婷婷| 免费观看精品视频网站| 精品国产一区二区三区四区第35| 国产精品 国内视频| 国产激情欧美一区二区| 久久国产精品男人的天堂亚洲| 亚洲五月婷婷丁香| 国语自产精品视频在线第100页| 精品久久久久久久毛片微露脸| 午夜免费鲁丝| 亚洲国产毛片av蜜桃av| 91国产中文字幕| 亚洲最大成人中文| 午夜福利免费观看在线| 久久久久久久久免费视频了| 啦啦啦免费观看视频1| 在线观看免费午夜福利视频| 国产精品一区二区精品视频观看| 国内毛片毛片毛片毛片毛片| 免费在线观看日本一区| 国产高清激情床上av| 在线看三级毛片| 午夜视频精品福利| 亚洲精华国产精华精| 韩国av一区二区三区四区| 长腿黑丝高跟| 国产aⅴ精品一区二区三区波| 国产一区在线观看成人免费| 中文字幕精品亚洲无线码一区 | 999精品在线视频| 日本a在线网址| 啦啦啦韩国在线观看视频| 中文字幕人成人乱码亚洲影| 久久香蕉激情| 精品福利观看| 一个人免费在线观看的高清视频| 久久人妻av系列| 午夜福利欧美成人| 视频在线观看一区二区三区| 最新美女视频免费是黄的| 男女做爰动态图高潮gif福利片| 久久精品成人免费网站| 欧美久久黑人一区二区| 18禁黄网站禁片免费观看直播| 国产国语露脸激情在线看| 精品少妇一区二区三区视频日本电影| 免费女性裸体啪啪无遮挡网站| 真人做人爱边吃奶动态| 深夜精品福利| 美女扒开内裤让男人捅视频| 国产三级在线视频| 操出白浆在线播放| 国产99白浆流出| 欧美成人一区二区免费高清观看 | 18禁国产床啪视频网站| 免费看日本二区| 午夜福利高清视频| 一边摸一边抽搐一进一小说| 亚洲va日本ⅴa欧美va伊人久久| 99在线人妻在线中文字幕| 国产av一区在线观看免费| 久久国产乱子伦精品免费另类| 国产欧美日韩一区二区精品| 成人特级黄色片久久久久久久| av超薄肉色丝袜交足视频| 日本熟妇午夜| 亚洲九九香蕉| 国产1区2区3区精品| 波多野结衣高清作品| 国产片内射在线| 巨乳人妻的诱惑在线观看| 国产视频内射| 少妇粗大呻吟视频| 后天国语完整版免费观看| 熟女电影av网| 久久亚洲真实| 免费无遮挡裸体视频| 免费看美女性在线毛片视频| 人妻丰满熟妇av一区二区三区| 欧美黑人精品巨大| e午夜精品久久久久久久| 精品一区二区三区视频在线观看免费| 亚洲aⅴ乱码一区二区在线播放 | 高潮久久久久久久久久久不卡| 丰满的人妻完整版| 国产高清videossex| 色综合亚洲欧美另类图片| 亚洲avbb在线观看| 啦啦啦免费观看视频1| 亚洲免费av在线视频| 日韩精品青青久久久久久| 最近在线观看免费完整版| 妹子高潮喷水视频| 99riav亚洲国产免费| 亚洲五月色婷婷综合| 成人午夜高清在线视频 | 嫩草影院精品99| 在线国产一区二区在线| 国产精品九九99| 怎么达到女性高潮| 深夜精品福利| 亚洲欧美日韩高清在线视频| 久久中文看片网| 成人亚洲精品av一区二区| 久久久国产欧美日韩av| 非洲黑人性xxxx精品又粗又长| 一区二区三区高清视频在线| 日韩大码丰满熟妇| 婷婷精品国产亚洲av| 亚洲精品中文字幕一二三四区| 成人永久免费在线观看视频| av在线天堂中文字幕| 91九色精品人成在线观看| 亚洲九九香蕉| 欧美乱码精品一区二区三区| 国产精品爽爽va在线观看网站 | 午夜视频精品福利| 1024手机看黄色片| 亚洲精品国产区一区二| 国产精品美女特级片免费视频播放器 | 黄色成人免费大全| 欧美黑人巨大hd| 亚洲精品一卡2卡三卡4卡5卡| 国产成人影院久久av| 午夜两性在线视频| 在线观看午夜福利视频| 黄色成人免费大全| 男女视频在线观看网站免费 | 波多野结衣巨乳人妻| 十八禁网站免费在线| 久久人妻av系列| 亚洲第一青青草原| xxx96com| 国产av不卡久久| 制服人妻中文乱码| 亚洲精品色激情综合| 久久久久久久久免费视频了| 88av欧美| 日日爽夜夜爽网站| 久久久精品欧美日韩精品| 亚洲av成人一区二区三| 日韩欧美免费精品| 欧美精品啪啪一区二区三区| 少妇裸体淫交视频免费看高清 | 特大巨黑吊av在线直播 | 美女 人体艺术 gogo| 免费看a级黄色片| 精品久久蜜臀av无| avwww免费| 日韩欧美国产在线观看| 宅男免费午夜| 精品久久久久久久久久免费视频| 欧美黑人欧美精品刺激| 啦啦啦韩国在线观看视频| 久久香蕉国产精品| 久久久久久久久中文| 黄色视频,在线免费观看| 无限看片的www在线观看| 搡老妇女老女人老熟妇| 亚洲一区中文字幕在线|