• <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
    美女主播在线视频| 久久精品久久久久久久性| 纵有疾风起免费观看全集完整版| 一本一本综合久久| 免费观看在线日韩| 水蜜桃什么品种好| 亚洲精品自拍成人| 久久毛片免费看一区二区三区| 最近2019中文字幕mv第一页| 在线观看国产h片| 日本色播在线视频| 三级国产精品欧美在线观看| av专区在线播放| 免费看av在线观看网站| 国产精品三级大全| 国产av精品麻豆| 国产在线免费精品| 国产成人精品福利久久| 亚洲欧美精品自产自拍| 欧美日本中文国产一区发布| 蜜桃在线观看..| 欧美激情极品国产一区二区三区 | 三级国产精品欧美在线观看| 国产 一区精品| 欧美精品一区二区免费开放| 色吧在线观看| 国国产精品蜜臀av免费| 国产综合精华液| 亚洲国产精品成人久久小说| 国产深夜福利视频在线观看| 欧美人与善性xxx| 国产免费一级a男人的天堂| 日韩中字成人| 中国国产av一级| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美日韩另类电影网站| 日韩视频在线欧美| 最近中文字幕2019免费版| 日产精品乱码卡一卡2卡三| 一区在线观看完整版| 国产成人91sexporn| 青青草视频在线视频观看| 亚洲国产色片| 99热这里只有精品一区| 免费人成在线观看视频色| 伦精品一区二区三区| 国产白丝娇喘喷水9色精品| 久久国产精品大桥未久av | 丰满迷人的少妇在线观看| 91久久精品国产一区二区三区| 欧美另类一区| 精品少妇内射三级| 亚洲色图综合在线观看| 亚洲一级一片aⅴ在线观看| 女人久久www免费人成看片| 国产高清三级在线| 51国产日韩欧美| 国产成人免费无遮挡视频| 美女国产视频在线观看| 国产一区二区在线观看日韩| 成人毛片a级毛片在线播放| 国产伦精品一区二区三区四那| 免费黄网站久久成人精品| 亚洲精品456在线播放app| 一区二区三区乱码不卡18| 精品久久久久久久久av| 国产亚洲欧美精品永久| 国产成人91sexporn| 亚洲三级黄色毛片| a级一级毛片免费在线观看| 久久国内精品自在自线图片| 91精品伊人久久大香线蕉| 国产精品欧美亚洲77777| av天堂久久9| 亚洲va在线va天堂va国产| 女性被躁到高潮视频| 国产无遮挡羞羞视频在线观看| 99热这里只有精品一区| 人妻 亚洲 视频| 久久毛片免费看一区二区三区| 国产又色又爽无遮挡免| 美女脱内裤让男人舔精品视频| 免费人妻精品一区二区三区视频| 国内少妇人妻偷人精品xxx网站| 久久久亚洲精品成人影院| 日产精品乱码卡一卡2卡三| 国产精品成人在线| 菩萨蛮人人尽说江南好唐韦庄| 亚洲,一卡二卡三卡| 一本一本综合久久| 三级国产精品片| 亚洲第一av免费看| 性高湖久久久久久久久免费观看| 波野结衣二区三区在线| 国产极品粉嫩免费观看在线 | 国产精品一区二区在线观看99| 亚洲精品自拍成人| 桃花免费在线播放| 精品午夜福利在线看| 国产精品一区二区在线不卡| av国产久精品久网站免费入址| 女性被躁到高潮视频| 午夜精品国产一区二区电影| 另类亚洲欧美激情| 精品熟女少妇av免费看| 日韩一本色道免费dvd| 精品国产乱码久久久久久小说| av免费观看日本| 成人漫画全彩无遮挡| 国模一区二区三区四区视频| 深夜a级毛片| 男女边吃奶边做爰视频| 99久久精品国产国产毛片| 最近2019中文字幕mv第一页| 国产成人aa在线观看| 国产色婷婷99| 亚洲成色77777| 久久热精品热| 色婷婷av一区二区三区视频| 中文字幕精品免费在线观看视频 | 男女边吃奶边做爰视频| 久久 成人 亚洲| 国产精品一区www在线观看| 一级av片app| 久久久a久久爽久久v久久| 在线免费观看不下载黄p国产| 国产美女午夜福利| 国产精品一二三区在线看| 极品教师在线视频| 中国美白少妇内射xxxbb| 人人妻人人澡人人爽人人夜夜| 美女主播在线视频| 搡女人真爽免费视频火全软件| 夫妻午夜视频| 18禁在线播放成人免费| 亚洲欧美日韩卡通动漫| 毛片一级片免费看久久久久| 99久久精品国产国产毛片| 久久国产亚洲av麻豆专区| 中国三级夫妇交换| 最新中文字幕久久久久| 夫妻午夜视频| 国产日韩一区二区三区精品不卡 | 亚洲精品aⅴ在线观看| 亚洲欧美一区二区三区黑人 | 交换朋友夫妻互换小说| 免费观看a级毛片全部| 日本av免费视频播放| 搡女人真爽免费视频火全软件| 校园人妻丝袜中文字幕| 成年av动漫网址| 久久久欧美国产精品| 亚洲av日韩在线播放| 免费播放大片免费观看视频在线观看| 国产精品一区www在线观看| 丁香六月天网| 亚洲第一区二区三区不卡| 国产黄色免费在线视频| 欧美性感艳星| 国产精品伦人一区二区| 高清视频免费观看一区二区| 精品久久久久久久久亚洲| 久久久精品94久久精品| 简卡轻食公司| 少妇被粗大的猛进出69影院 | 99热6这里只有精品| 久久久精品免费免费高清| 亚洲电影在线观看av| 免费观看a级毛片全部| 精品卡一卡二卡四卡免费| 青春草国产在线视频| 少妇人妻一区二区三区视频| av天堂久久9| 国产乱来视频区| 又黄又爽又刺激的免费视频.| 亚洲精品aⅴ在线观看| 国产亚洲91精品色在线| 国产 精品1| 秋霞伦理黄片| 日本av免费视频播放| 偷拍熟女少妇极品色| 精品99又大又爽又粗少妇毛片| 麻豆成人av视频| 国产精品偷伦视频观看了| 高清av免费在线| 在线观看www视频免费| 国产免费又黄又爽又色| 亚洲精品自拍成人| 国产黄片视频在线免费观看| 黑人巨大精品欧美一区二区蜜桃 | 熟妇人妻不卡中文字幕| 黄色一级大片看看| 国产精品成人在线| 久久午夜福利片| 日本黄色片子视频| 美女脱内裤让男人舔精品视频| 亚洲av成人精品一二三区| 黄色毛片三级朝国网站 | 又大又黄又爽视频免费| videossex国产| 免费看日本二区| 免费av不卡在线播放| 久久久久精品久久久久真实原创| 嫩草影院新地址| 观看免费一级毛片| 99精国产麻豆久久婷婷| 中国三级夫妇交换| 精品久久久精品久久久| av专区在线播放| 22中文网久久字幕| 激情五月婷婷亚洲| 中文字幕制服av| 欧美激情极品国产一区二区三区 | 中国美白少妇内射xxxbb| 亚洲欧美精品专区久久| 精品久久久精品久久久| 99国产精品免费福利视频| 伦精品一区二区三区| 91午夜精品亚洲一区二区三区| 桃花免费在线播放| 日韩欧美精品免费久久| 国产色婷婷99| 国产一区二区在线观看日韩| 国模一区二区三区四区视频| 久久久久久久亚洲中文字幕| 性高湖久久久久久久久免费观看| 一边亲一边摸免费视频| 男人和女人高潮做爰伦理| 久久久a久久爽久久v久久| 精品熟女少妇av免费看| 久久 成人 亚洲| 18禁动态无遮挡网站| 制服丝袜香蕉在线| 五月玫瑰六月丁香| 纯流量卡能插随身wifi吗| 亚洲经典国产精华液单| 精品久久国产蜜桃| 水蜜桃什么品种好| 日韩欧美一区视频在线观看 | 99九九线精品视频在线观看视频| 日韩免费高清中文字幕av| 一级黄片播放器| 国产免费又黄又爽又色| 日韩一本色道免费dvd| 黄色毛片三级朝国网站 | 日产精品乱码卡一卡2卡三| 色哟哟·www| 丝袜在线中文字幕| 亚洲性久久影院| 国产精品成人在线| 美女福利国产在线| 精品一品国产午夜福利视频| 精华霜和精华液先用哪个| 中国国产av一级| 晚上一个人看的免费电影| 看十八女毛片水多多多| 亚洲精品乱久久久久久| 在线免费观看不下载黄p国产| 少妇的逼水好多| 这个男人来自地球电影免费观看 | 男的添女的下面高潮视频| 国产亚洲午夜精品一区二区久久| 少妇熟女欧美另类| 成人亚洲精品一区在线观看| 国产精品国产三级专区第一集| 免费大片18禁| 亚洲不卡免费看| 肉色欧美久久久久久久蜜桃| 日韩中文字幕视频在线看片| 国产淫片久久久久久久久| 美女中出高潮动态图| 在线观看av片永久免费下载| 成人影院久久| 午夜精品国产一区二区电影| 亚洲电影在线观看av| 国国产精品蜜臀av免费| 国产爽快片一区二区三区| 久久久久国产网址| 街头女战士在线观看网站| 色视频在线一区二区三区| av黄色大香蕉| 三级国产精品欧美在线观看| 狂野欧美激情性xxxx在线观看| 啦啦啦啦在线视频资源| 免费大片黄手机在线观看| 一区二区av电影网| 熟女人妻精品中文字幕| 国产精品国产av在线观看| 免费人妻精品一区二区三区视频| 毛片一级片免费看久久久久| 免费人成在线观看视频色| 国产午夜精品久久久久久一区二区三区| 99九九线精品视频在线观看视频| 国产一级毛片在线| 免费av不卡在线播放| 大码成人一级视频| 超碰97精品在线观看| 亚洲欧美日韩东京热| 日本vs欧美在线观看视频 | 精品久久久久久久久亚洲| 午夜福利影视在线免费观看| 成年人免费黄色播放视频 | 成人18禁高潮啪啪吃奶动态图 | 搡女人真爽免费视频火全软件| 久久韩国三级中文字幕| 99久久精品一区二区三区| 少妇人妻久久综合中文| 狂野欧美激情性bbbbbb| 黄色视频在线播放观看不卡| 亚洲精品中文字幕在线视频 | 丰满饥渴人妻一区二区三| 久久精品国产鲁丝片午夜精品| 亚洲无线观看免费| 91成人精品电影| 久久久a久久爽久久v久久| 国语对白做爰xxxⅹ性视频网站| 春色校园在线视频观看| 乱系列少妇在线播放| 国产精品久久久久久久久免| 欧美日本中文国产一区发布| 亚洲av日韩在线播放| 亚洲自偷自拍三级| 高清不卡的av网站| 亚洲国产成人一精品久久久| 少妇人妻精品综合一区二区| 亚洲精品乱码久久久v下载方式| 日日摸夜夜添夜夜添av毛片| 18+在线观看网站| 色网站视频免费| 欧美精品一区二区免费开放| 国产亚洲午夜精品一区二区久久| 在线观看美女被高潮喷水网站| 嫩草影院新地址| 精品久久久久久久久av| 亚洲,一卡二卡三卡| 成年人午夜在线观看视频| 乱人伦中国视频| 精品人妻熟女毛片av久久网站| 亚洲电影在线观看av| 人人妻人人澡人人看| 国产一区二区在线观看av| 一级a做视频免费观看| av国产久精品久网站免费入址| 涩涩av久久男人的天堂| 在线观看一区二区三区激情| 在线观看av片永久免费下载| 国产免费一区二区三区四区乱码| 久久av网站| 亚洲伊人久久精品综合| 男人狂女人下面高潮的视频| 日韩av不卡免费在线播放| 少妇猛男粗大的猛烈进出视频| 看免费成人av毛片| 精品一区在线观看国产| 又黄又爽又刺激的免费视频.| 亚洲真实伦在线观看| 国产亚洲精品久久久com| 久久鲁丝午夜福利片| 97超视频在线观看视频| 国产精品伦人一区二区| 亚洲一区二区三区欧美精品| 在现免费观看毛片| 成人亚洲精品一区在线观看| 亚洲色图综合在线观看| 91精品一卡2卡3卡4卡| 国产亚洲最大av| 精品亚洲成a人片在线观看| 国产黄频视频在线观看| 国产亚洲精品久久久com| av.在线天堂| 久久久久久久精品精品| 秋霞伦理黄片| 精品国产一区二区久久| 国产亚洲一区二区精品| 亚洲激情五月婷婷啪啪| 亚洲四区av| 尾随美女入室| 天美传媒精品一区二区| 国产在线男女| 777米奇影视久久| 中文字幕免费在线视频6| 日产精品乱码卡一卡2卡三| 欧美 亚洲 国产 日韩一| 看免费成人av毛片| 777米奇影视久久| 亚洲成人一二三区av| 国产精品久久久久久久电影| 国产成人a∨麻豆精品| av专区在线播放| 桃花免费在线播放| 欧美97在线视频| 一区二区av电影网| 建设人人有责人人尽责人人享有的| 欧美三级亚洲精品| 在线 av 中文字幕| av福利片在线观看| 免费观看性生交大片5| 高清毛片免费看| 大话2 男鬼变身卡| 中文字幕精品免费在线观看视频 | 久久久精品94久久精品| 观看免费一级毛片| 国产一区二区在线观看日韩| 久久精品久久久久久久性| 亚洲av欧美aⅴ国产| 国产成人免费无遮挡视频| 久久女婷五月综合色啪小说| 久久精品久久精品一区二区三区| 成人午夜精彩视频在线观看| 99久久中文字幕三级久久日本| 亚洲人成网站在线播| 婷婷色麻豆天堂久久| 国产淫片久久久久久久久| 高清午夜精品一区二区三区| 久久国产亚洲av麻豆专区| 久久国产精品大桥未久av | www.av在线官网国产| 国产亚洲最大av| 精品久久久噜噜| 少妇人妻 视频| 男人狂女人下面高潮的视频| 男女边摸边吃奶| 麻豆乱淫一区二区| 国产成人精品无人区| 国产在视频线精品| 黄色一级大片看看| 亚洲欧美日韩卡通动漫| 久久久欧美国产精品| 777米奇影视久久| 亚洲精品一区蜜桃| 国产一区二区三区av在线| 青春草国产在线视频| 日韩一本色道免费dvd| 99久国产av精品国产电影| 免费观看在线日韩| 美女cb高潮喷水在线观看| 人人妻人人澡人人爽人人夜夜| 久久精品久久精品一区二区三区| 国产精品99久久久久久久久| 久久久精品94久久精品| 中文字幕久久专区| 国产高清国产精品国产三级| 欧美变态另类bdsm刘玥| 波野结衣二区三区在线| 91久久精品国产一区二区成人| 亚洲无线观看免费| 成人漫画全彩无遮挡| 另类亚洲欧美激情| 一区二区av电影网| 大又大粗又爽又黄少妇毛片口| 国产精品不卡视频一区二区| 18禁裸乳无遮挡动漫免费视频| 成年av动漫网址| 亚洲欧美精品专区久久| a级毛色黄片| 国产成人精品无人区| av不卡在线播放| 日本免费在线观看一区| 五月玫瑰六月丁香| 一本久久精品| 亚洲在久久综合| 免费大片18禁| 大又大粗又爽又黄少妇毛片口| 美女国产视频在线观看| 亚洲av综合色区一区| 成年美女黄网站色视频大全免费 | 亚洲精品日本国产第一区| 日韩熟女老妇一区二区性免费视频| av国产精品久久久久影院| 欧美变态另类bdsm刘玥| 亚洲精品国产成人久久av| 成人午夜精彩视频在线观看| 一级毛片 在线播放| 亚洲图色成人| 九九久久精品国产亚洲av麻豆| 久久久久久久国产电影| 亚洲美女黄色视频免费看| 精品人妻熟女av久视频| 人妻少妇偷人精品九色| 九九在线视频观看精品| 午夜91福利影院| 成人免费观看视频高清| av卡一久久| 亚洲电影在线观看av| 国产精品久久久久成人av| 亚洲成色77777| 亚洲高清免费不卡视频| 最近中文字幕2019免费版| 一区二区三区乱码不卡18| 伦理电影大哥的女人| 日本爱情动作片www.在线观看| 国产av精品麻豆| 我要看黄色一级片免费的| 亚洲第一区二区三区不卡| 婷婷色综合www| 成年美女黄网站色视频大全免费 | 国产亚洲精品久久久com| 肉色欧美久久久久久久蜜桃| 免费观看的影片在线观看| 亚洲精品自拍成人| 99热6这里只有精品| 成人国产麻豆网| 久久精品国产亚洲av天美| 少妇被粗大猛烈的视频| 亚洲国产毛片av蜜桃av| 爱豆传媒免费全集在线观看| 成人美女网站在线观看视频| 有码 亚洲区| 日韩欧美一区视频在线观看 | 一个人免费看片子| 国产成人免费观看mmmm| 亚洲欧美日韩另类电影网站| 亚洲美女黄色视频免费看| 丝袜在线中文字幕| 少妇被粗大猛烈的视频| 亚洲不卡免费看| 亚洲av.av天堂| 97在线人人人人妻| 日韩三级伦理在线观看| 免费黄频网站在线观看国产| 亚洲自偷自拍三级| 成人美女网站在线观看视频| 日本猛色少妇xxxxx猛交久久| 在线观看免费视频网站a站| av国产精品久久久久影院| 伦理电影大哥的女人| 午夜福利视频精品| 夫妻性生交免费视频一级片| 好男人视频免费观看在线| 日本免费在线观看一区| 我要看黄色一级片免费的| 亚洲精品久久久久久婷婷小说| 欧美丝袜亚洲另类| 免费观看a级毛片全部| 久久久久久人妻| 亚洲av.av天堂| 精品久久久噜噜| av又黄又爽大尺度在线免费看| 欧美精品高潮呻吟av久久| 免费av中文字幕在线| 国产精品不卡视频一区二区| 精品少妇久久久久久888优播| 性高湖久久久久久久久免费观看| av国产久精品久网站免费入址| 久久人妻熟女aⅴ| 男人和女人高潮做爰伦理| 黄片无遮挡物在线观看| 女的被弄到高潮叫床怎么办| 亚洲欧美成人综合另类久久久| 国产精品蜜桃在线观看| 亚洲欧美精品专区久久| 午夜精品国产一区二区电影| 亚洲精品,欧美精品| 性色av一级| 男人舔奶头视频| 老司机亚洲免费影院| 日韩在线高清观看一区二区三区| 在线观看国产h片| 久久毛片免费看一区二区三区| 亚州av有码| 美女视频免费永久观看网站| 亚洲真实伦在线观看| 国产无遮挡羞羞视频在线观看| 国产极品粉嫩免费观看在线 | 在线观看www视频免费| av福利片在线| 亚洲图色成人| 爱豆传媒免费全集在线观看| 久久精品熟女亚洲av麻豆精品| a级毛片免费高清观看在线播放| 日韩在线高清观看一区二区三区| 国产精品女同一区二区软件| 国产av精品麻豆| 能在线免费看毛片的网站| 大话2 男鬼变身卡| 国产亚洲一区二区精品| 国产精品国产三级专区第一集| 在线免费观看不下载黄p国产| 国产欧美日韩综合在线一区二区 | 夫妻性生交免费视频一级片| 亚洲精品第二区| 久久99精品国语久久久| 麻豆成人午夜福利视频| 欧美bdsm另类| av在线app专区| 噜噜噜噜噜久久久久久91| 性高湖久久久久久久久免费观看| 大香蕉97超碰在线| 男人狂女人下面高潮的视频| 亚洲av男天堂| 欧美高清成人免费视频www| 在线观看免费高清a一片| 国产精品一区二区三区四区免费观看| 久久久久国产精品人妻一区二区| 亚洲欧美精品专区久久| 久久精品夜色国产| 爱豆传媒免费全集在线观看| 国产又色又爽无遮挡免| 日韩欧美精品免费久久| 亚洲av国产av综合av卡| 亚洲精品成人av观看孕妇| 男女啪啪激烈高潮av片| 99热这里只有是精品在线观看| 26uuu在线亚洲综合色| 亚洲第一av免费看| 观看美女的网站| 精品一区二区三卡| 亚洲婷婷狠狠爱综合网| 日韩人妻高清精品专区| 丰满饥渴人妻一区二区三| 国产精品秋霞免费鲁丝片| a级一级毛片免费在线观看| 一区二区三区四区激情视频| 波野结衣二区三区在线| 中文乱码字字幕精品一区二区三区|