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

    人口普查遺漏的組合式估計方法

    2024-02-04 02:49:42胡桂華LPEZCARRDavid黃艷華
    統(tǒng)計與信息論壇 2024年2期
    關(guān)鍵詞:估計量人口數(shù)遺漏

    胡桂華,LPEZ-CARR David,黃艷華,吳 笛

    (1.重慶工商大學(xué) a.數(shù)學(xué)與統(tǒng)計學(xué)院,b.成渝地區(qū)雙城經(jīng)濟(jì)圈建設(shè)研究院,重慶 400067;2.加州大學(xué)圣芭芭拉分校 地理系,加利福尼亞 圣芭芭拉,93106)

    一、引 言

    人口普查不可避免發(fā)生遺漏,導(dǎo)致普查登記人口數(shù)低于實際人口數(shù)。人口普查遺漏指應(yīng)該在普查中登記卻未能登記在人口普查表中,要估計普查遺漏人口數(shù),需要使用其他資料(如覆蓋調(diào)查資料等),通過比對其他資料與人口普查資料來發(fā)現(xiàn)遺漏人口,或者在某假設(shè)條件下利用已有資料構(gòu)造遺漏估計量。

    人口普查質(zhì)量評估實踐表明,遺漏在整個人群中的分布并非均勻。不同人口群體或地區(qū)之間的遺漏率差異較大,美國2010年人口普查遺漏估計結(jié)果顯示[1],西班牙裔少數(shù)民族的遺漏率高于非西班牙族白人;租房者遺漏率(8.5%)比有房者(3.7%)遺漏率高;在各州中,遺漏率從愛荷華州的2.6%到密西西比州的8.9%不等;大城市的遺漏率往往高于其他地區(qū)的普查遺漏率。

    每次人口普查后,政府統(tǒng)計部門需要估計總?cè)丝诘倪z漏及各類別人口的遺漏。與遺漏相關(guān)的一個概念是凈遺漏。雖然遺漏和凈遺漏都反映了人口普查數(shù)據(jù)準(zhǔn)確性的各個方面,但揭示的內(nèi)容不同,遺漏反映的是本應(yīng)在本次人口普查中登記卻未登記的人數(shù);凈遺漏率反映的是遺漏人數(shù)的百分比與在普查中重復(fù)登記人數(shù)的百分比之差。比如,美國2010年全國的遺漏率為5.3%,而凈遺漏率為0.01%。

    造成人口普查遺漏的原因較多,有些受訪者認(rèn)為不應(yīng)該將幼兒填寫在人口普查表中[2]。美國2010年的一項專項調(diào)查結(jié)果表明,在低收入家庭的人口普查登記中,只有82%的家庭將幼兒填寫在人口普查表中,另外18%的家庭認(rèn)為幼兒不應(yīng)該納入人口普查登記范圍。雖然人口普查操作指南手冊明確規(guī)定應(yīng)該登記普查日前出生的幼兒和收養(yǎng)的棄嬰,但實際上并未嚴(yán)格遵守,從而導(dǎo)致遺漏幼兒。在人口普查登記期間,有些住戶成員外出經(jīng)商、旅游或?qū)W習(xí),當(dāng)普查員多次上門登記時不在家,從而造成遺漏;有些居民住宅修建在偏遠(yuǎn)地區(qū),造成普查員登記困難,在普查表中未登記這樣的住宅及其家庭成員;有些普查員不適當(dāng)?shù)貙⒕幼∮米≌?dāng)作商業(yè)用房,從而未登記居住在其中的家庭和個人造成遺漏;多個家庭共用同一個地址和同一個門牌號碼,普查員只登記了這個地址上的一個家庭,而未登記該地址上其他應(yīng)該登記的家庭等。這些情況和其他尚未列舉的情況表明,人口普查遺漏客觀存在,有必要開展相關(guān)研究。

    筆者研究過普查人口名單及覆蓋調(diào)查人口名單情形下的完整遺漏估計量,以及這兩項調(diào)查人口名單和行政記錄人口名單情形下的合成遺漏估計量。完整遺漏估計量由一個單重遺漏估計量和兩個雙重遺漏估計量構(gòu)成[2];合成遺漏估計量包括兩個單重遺漏估計量、一個雙重遺漏估計量和一個三重遺漏估計量[3]。這兩種遺漏估計量適合于人口普查質(zhì)量評估水平較高的發(fā)達(dá)國家使用。本文研究的組合式遺漏估計量包括一個單重遺漏估計量和一個雙重遺漏估計量,適合于人口普查質(zhì)量評估水平相對較低的發(fā)展中國家使用。

    本文研究目標(biāo)是,使用所提出的組合式遺漏估計量替代目前聯(lián)合國統(tǒng)計司建議各國使用的單重遺漏估計量,以解決后者低估總體人口普查遺漏人口數(shù)的問題。

    研究具有一定的理論意義與現(xiàn)實意義。理論意義表現(xiàn)在兩個方面:一方面,組合式遺漏估計量由單重遺漏估計量和統(tǒng)計獨立的雙重遺漏估計量構(gòu)成,雙重遺漏估計量的理論基礎(chǔ)是捕獲—再捕獲模型。為滿足該模型等概率要求,需要使用體現(xiàn)人口在普查中登記概率大小的人口統(tǒng)計特征變量和居住位置變量對總體人口分層。在等概率層構(gòu)造雙重遺漏估計量,否則產(chǎn)生異質(zhì)性偏差。在分層二重抽樣下,雙重遺漏估計量的構(gòu)成要素使用再加權(quán)擴(kuò)張估計量構(gòu)造。另一方面,組合式遺漏估計量較為復(fù)雜,其抽樣方差使用分層刀切抽樣方差估計量近似估計,組合式遺漏估計量有偏,因此需估計其偏差。

    現(xiàn)實意義表現(xiàn)在三個方面。第一,組合式遺漏估計量除了包括未登記在普查人口名單但登記在覆蓋調(diào)查人口名單的單重遺漏人口外,還包括了雙重遺漏人口,即同時遺漏于這兩份人口名單的人口,因而覆蓋了總體絕大部分遺漏人口(未包括覆蓋調(diào)查與普查非獨立情況下的雙重遺漏估計量)。未包括這類雙重遺漏估計量基于兩點考慮,一是筆者已經(jīng)研究過這類雙重遺漏估計量,二是估計這類雙重遺漏人口數(shù)需要雙系統(tǒng)估計量和人口統(tǒng)計分析模型估計的總體實際人口數(shù),以及覆蓋調(diào)查數(shù)據(jù)處理誤差數(shù)據(jù)。中國目前缺少這些數(shù)據(jù),而本文為未來中國使用組合式遺漏估計量提供參考;第二,為構(gòu)造組合式遺漏估計量,需要比對兩份人口名單,在比對過程中,可查明遺漏者特征、遺漏原因、遺漏程度及遺漏在總體中的分布情況;第三,組合式遺漏估計量利用了普查人口名單輔助信息,使用高質(zhì)量的、相關(guān)程度高的輔助信息是統(tǒng)計推斷的一個重要原則。這三個方面的情況表明,組合式遺漏估計量有望提高人口普查遺漏估計精度,應(yīng)用于政府統(tǒng)計部門人口普查遺漏估計。

    本文的創(chuàng)新體現(xiàn)在兩個方面,一是使用樣本普查小區(qū)的普查人口名單和覆蓋調(diào)查人口名單的原始數(shù)據(jù),全面演示了組合式遺漏估計量的計算過程;二是討論了組合式遺漏估計量的統(tǒng)計性質(zhì)。尚未發(fā)現(xiàn)政府統(tǒng)計部門和國內(nèi)外其他學(xué)者研究人口普查遺漏估計量的統(tǒng)計性質(zhì)。

    二、文獻(xiàn)綜述

    估計人口普查遺漏的方法較多,如剩余法、行政記錄比較法和覆蓋調(diào)查法等[2]?;谑S喾ǖ钠詹檫z漏估計量為雙系統(tǒng)估計量與基于比率估計量的普查正確登記人口數(shù)估計量之差[4-7]。普查正確登記人口數(shù)估計量為最終普查登記人口數(shù)與普查正確登記率的乘積。普查正確登記率為普查正確登記人口數(shù)估計量與至少填寫了姓名和兩個調(diào)查項目的普查登記人口數(shù)估計量之比。剩余法的優(yōu)勢是能夠利用已經(jīng)算出的結(jié)果快速取得普查遺漏估計值,劣勢是無法分析遺漏者的人口統(tǒng)計特征。行政記錄比較法的基本思想是,將總體中的相關(guān)人口行政記錄合并為一份名單,以家庭為抽樣單位,從這份名單中抽取若干家庭,并與人口普查微觀記錄進(jìn)行比較。如果某人被登記在這份名單,而未登記在普查記錄中,收判斷這個人在普查中遺漏。這種方法的優(yōu)勢是確保了行政記錄人口名單獨立于人口普查,而且無需花費數(shù)據(jù)采集成本,劣勢是獲取行政記錄人口名單難度大,另外還需要剔除行政記錄人口名單中的重復(fù)人口,以及定期更新行政記錄人口名單。覆蓋調(diào)查是為了估計普查遺漏和其他指標(biāo)在人口普查之后進(jìn)行的一項抽樣調(diào)查,其抽樣單位不是人或住房單元,而是小范圍的地理區(qū)域,例如中國的某一普查小區(qū)。在覆蓋調(diào)查中,一方面可以得到樣本普查小區(qū)的普查人口名單,另一方面還可以獲得覆蓋調(diào)查人口名單。普查人口名單由樣本普查小區(qū)每個住房單元填寫的普查表匯編而成,覆蓋調(diào)查人口名單由樣本普查小區(qū)每個住房單元填寫的覆蓋調(diào)查表編制而成,依據(jù)這兩份調(diào)查人口名單可以構(gòu)造單重遺漏估計量、雙重遺漏估計量和組合式遺漏估計量。覆蓋調(diào)查的優(yōu)勢是,可以構(gòu)造多種形式的遺漏估計量,劣勢是需要采取措施確保普查與覆蓋調(diào)查之間的獨立性,否則據(jù)此構(gòu)造的遺漏估計量存在交互作用偏差。

    單重遺漏估計量是目前許多國家估計普查遺漏的主要方法[8-10]。除美國和加拿大等少數(shù)國家和地區(qū)外,包括中國、南非和盧旺達(dá)在內(nèi)的許多國家使用此方法[11-12]。單重遺漏估計量為樣本普查小區(qū)的未匹配人口與其抽樣權(quán)數(shù)的線性估計量。未匹配人口是指登記在覆蓋調(diào)查人口名單而未登記在普查人口名單的遺漏人口[13],這里有一個假設(shè),即樣本普查小區(qū)的覆蓋調(diào)查人口名單本身未遺漏人口。如果這一假設(shè)不成立,那就意味著,有些人同時遺漏于這兩份調(diào)查人口名單(稱為雙重遺漏人口)。單重遺漏估計量未包括雙重遺漏人口,因而低估遺漏人口數(shù)。雖然覆蓋調(diào)查規(guī)模比普查小許多,而且調(diào)查員比普查員專業(yè)性更強(qiáng),但覆蓋調(diào)查依然可能遺漏人口。美國人口普查局設(shè)計的2000年、2010年和2020年人口普查質(zhì)量評估方案中均有一章專門論述覆蓋調(diào)查遺漏人口的處理方法。中國2000年、2010年和2020年覆蓋調(diào)查樣本規(guī)模分別為602個、402個和406個普查小區(qū)。如此小的樣本規(guī)模,原則上應(yīng)該做到無人口遺漏,然而實際上覆蓋調(diào)查難以100%登記人口,這意味著,當(dāng)構(gòu)造普查遺漏估計量時,不能想當(dāng)然認(rèn)為覆蓋調(diào)查不遺漏人口,而應(yīng)該包括雙重遺漏人口。與單重遺漏人口所不同的是,雙重遺漏人口既未登記在普查人口名單,也未登記在覆蓋調(diào)查人口名單,因此找到雙重遺漏人口難度更大。要找到雙重遺漏人口,首先,要估計研究區(qū)域內(nèi)的雙重遺漏人口數(shù),如果估計結(jié)果為零,就放棄尋找雙重遺漏人口;其次,如果估計結(jié)果不為零,就通過社區(qū)負(fù)責(zé)人篩選可能的雙重遺漏人口,一般來說,雙重遺漏人口主要是本社區(qū)的無固定住所者、犯罪人員、在逃人員、獨住者和不關(guān)心國家大事者;最后,在確定雙重遺漏人口后,通過面訪調(diào)查或者人口行政資料(如戶籍資料)獲得其姓名、性別、年齡、文化程度、婚姻狀況、職業(yè)等人口特征。

    三、人口普查遺漏估計理論

    (一)組合式遺漏估計量及其方差估計

    為便于計算,將組合式遺漏估計量(Combined Omission Estimator,COE)的單重遺漏估計量和統(tǒng)計獨立情況下的雙重遺漏估計量放在同一等概率人口層(用v表示)建立,盡管單重遺漏估計量無需在等概率人口層構(gòu)造,等概率人口層是通過對總體按照體現(xiàn)人口在普查中登記概率大小的變量進(jìn)行分層得到的。變量值相同或大致相同的人在同一層,分層在覆蓋調(diào)查樣本抽取后實施,以確保樣本中的每一個人有一個對應(yīng)的層。分層目標(biāo)是減少組合式遺漏估計量的異質(zhì)性偏差[14]。分層變量越多,層內(nèi)的同質(zhì)性越強(qiáng),異質(zhì)性偏差越小,每一層的覆蓋調(diào)查樣本量也越少,估計遺漏人口數(shù)的抽樣方差越大,因此,在確定對總體人口等概率分層變量及其變量值時,要綜合考慮層內(nèi)的異質(zhì)性偏差、抽樣方差及覆蓋調(diào)查的樣本規(guī)模。

    為構(gòu)造組合式遺漏估計量,引入不完整二維列聯(lián)表,將同一樣本普查小區(qū)的普查人口名單和覆蓋調(diào)查人口名單的比對結(jié)果填入該表,見表1。

    表1 等概率人口層v的不完整二維列聯(lián)表

    比對通常在同一樣本普查小區(qū)內(nèi)進(jìn)行,其目標(biāo)之一是查找登記在覆蓋調(diào)查人口名單的人是否也登記在普查人口名單。如果未登記在普查人口名單,就判斷這個人在普查中遺漏。為避免虛增遺漏,兩份名單的比對范圍應(yīng)由此樣本普查小區(qū)擴(kuò)大到周圍區(qū)域,即搜索區(qū)域(2)最理想的比對范圍是全國普查微觀記錄數(shù)據(jù)庫。在這個數(shù)據(jù)庫搜索覆蓋調(diào)查人口名單的人口,確定其是否在普查人口名單中登記。如果未登記,則作為普查遺漏人口。,判斷是否能夠在搜索區(qū)域的普查人口名單找到與覆蓋調(diào)查人口名單一致的人,如果未找到則判作覆蓋調(diào)查人口名單的這個人在普查中遺漏,如果找到了則不能夠判作遺漏,有些人在普查中并未登記在所屬的樣本普查小區(qū),而是登記在其搜索區(qū)域。在這種情況下,如果只是將比對范圍局限在此樣本普查小區(qū),勢必虛增遺漏人口。

    (1)

    (2)

    (3)

    (4)

    表明,在普查中登記的人認(rèn)為已經(jīng)答復(fù)了普查問題,因而相比未參加普查的人,更加不愿意參加覆蓋調(diào)查,使Y11v的人口數(shù)減少,從而高估等概率人口層v的實際人口數(shù)。

    (5)

    意味著,在普查中登記的人相比未登記的人更加意識到普查的重要性,從而更加愿意參加覆蓋調(diào)查,使Y11v的人口數(shù)增加,從而低估等概率人口層v的實際人口數(shù)。本文構(gòu)造的雙重遺漏估計量為式(3)。

    如果采取第二種方法,在同時滿足兩個假設(shè)條件的情況下(覆蓋調(diào)查與普查獨立,通過在這兩項調(diào)查中使用不同的調(diào)查方法、不同的調(diào)查員和不同組織機(jī)構(gòu)實現(xiàn);總體中的每一個人有同樣的概率登記在普查或覆蓋調(diào)查中,通過對總體人口等概率分層來實現(xiàn)),(Yv,p1+v,p+1v)的最大似然函數(shù)為:

    (6)

    式(6)中的p1+v,p+1v分別為等概率人口層的人在普查和覆蓋調(diào)查中登記的邊際概率。

    依據(jù)最大似然估計方法,得到式(6)中的Yv,P1+v,P+1v的最大似然估計量為:

    (7)

    考慮到覆蓋調(diào)查實際為抽樣調(diào)查,式(1)可寫作式(8):

    (8)

    (9)

    式(9)中,h表示第一重抽樣層的任意層,h=1,2,…,H;g是對第一重樣本進(jìn)一步分的任意層,g=1,2,…,Gh。如果第一重樣本小區(qū)i進(jìn)入層g,那么xhgi=1,否則為0;如果繼續(xù)進(jìn)入第二重樣本,那么shgi=1,否則為0;yhgiv為第二重樣本普查小區(qū)hgi在等概率人口層v的觀察值;αhgi為經(jīng)過兩重抽樣后樣本普查小區(qū)hgi的抽樣權(quán)數(shù),其功能是將樣本指標(biāo)擴(kuò)張到總體指標(biāo)。

    式(8)為復(fù)雜估計量,采用分層刀切抽樣方差估計量近似計算其抽樣方差,表示為:

    (10)

    (11)

    (12)

    (13)

    式(13)中,Chg,chg分別表示交叉層hg的普查小區(qū)集合和樣本普查小區(qū)集合;Nh,nh分別是層h的普查小區(qū)總數(shù)目和樣本普查小區(qū)數(shù)目;Mhg,mhg分別是層hg的普查小區(qū)總數(shù)目和樣本普查小區(qū)數(shù)目。

    (14)

    (15)

    式(15)中的協(xié)方差計算公式如下:

    (16)

    (二)組合式遺漏估計量的統(tǒng)計性質(zhì)

    1.無偏性

    關(guān)于估計量的無偏性,應(yīng)該從理論上證明[21]。如果估計量相對簡單,那么完成其無偏性理論證明較容易,例如,簡單隨機(jī)抽樣下樣本均值是總體均值的無偏估計量,但當(dāng)估計量復(fù)雜時,完成其無偏性的理論證明則會遇到較大困難。從式(8)和式(9)可以看出,要從理論上證明組合式遺漏估計量的無偏性是一項較難完成的工作。首先,這兩個估計量的概率分布無從知曉;其次,在分層二重抽樣下,采用雙重擴(kuò)張估計量構(gòu)造,尤其是組合式遺漏估計量中的雙重遺漏估計量是一個分?jǐn)?shù)估計量,而且分母與分子都是依據(jù)雙重擴(kuò)張估計量構(gòu)造[22]。鑒于這兩個遺漏估計量難以從理論上進(jìn)行無偏性證明,于是本文采用模擬的方法予以討論。首先,根據(jù)抽樣方法確定可能的樣本個數(shù),計算每個樣本的組合式遺漏估計值;其次,計算所有可能樣本的組合式遺漏估計值的平均值,將其作為組合式遺漏估計量的數(shù)學(xué)期望;最后,將組合式遺漏估計量的估計值作為其真實值,依據(jù)偏差公式(估計量的均值與真實值之差)計算這兩個估計量的偏差,如果偏差為0,則為無偏估計量,反之為有偏估計量。組合式遺漏估計量可能為有偏估計量,需計算其均方誤差。

    2.有效性

    有效性是指估計量與總體參數(shù)的離散程度。如果兩個估計量都是無偏的,那么離散程度較小的估計量相對來說是有效的,離散程度用方差來衡量。從數(shù)理統(tǒng)計理論來看,討論估計量有效性的前提條件是這兩個估計量均為無偏估計量。單重遺漏估計量和組合式遺漏估計量是同一總體人口普查遺漏參數(shù)的兩個估計量。在后面的模擬研究中,發(fā)現(xiàn)這兩個遺漏估計量的偏差均不為零,即它們?yōu)橛衅烙嬃?。然?并不意味著對這兩個有偏估計量就不能進(jìn)行有效性比較,仍可以使用均方誤差來比較它們的有效性。

    3.一致性

    一致性是指隨著樣本規(guī)模的增大,估計量越來越接近總體參數(shù)的真值。在人口普查質(zhì)量評估中,覆蓋調(diào)查的樣本規(guī)模受到嚴(yán)格限制,達(dá)不到一致性所要求的樣本規(guī)模。因此,單重遺漏估計量和組合式遺漏估計量不具備一致性。

    4.充分性

    如果一個估計量利用了需要估計的總體參數(shù)的全部信息,則稱其具有充分性。從理論角度論證估計量的充分性有兩種方法,一是從充分估計量的定義出發(fā),確定既定估計量取值后樣本的條件分布;二是使用因子分解定理證明,使用這兩種方法的前提條件是總體的概率函數(shù)已知,然而,要從理論上論證單重遺漏估計量和組合式遺漏估計量是否為充分估計量十分困難。困難之處在于,人口普查標(biāo)準(zhǔn)時點上的人口總體的概率函數(shù)未知,既無法確定條件分布,也無法使用因子分解定理,因此,在這里只從直觀上討論單重遺漏估計量和組合式遺漏估計量的充分性。在兩份名單獨立的情況下,為了估計總體遺漏,應(yīng)該用到的全部信息包括登記在覆蓋調(diào)查人口名單但未登記在普查人口名單的人,同時未登記在這兩份調(diào)查人口名單的人。單重遺漏估計量利用了第一種信息,雙重遺漏估計量利用了第二種信息,組合式遺漏估計量則同時利用了第一種和第二種信息,可見,單重遺漏估計量并沒有利用全部信息,不具備充分性,而組合式遺漏估計量利用了全部信息,具有充分性。因此,從充分性看,組合式遺漏估計量優(yōu)于單重遺漏估計量和雙重遺漏估計量。

    相較于對估計量統(tǒng)計性質(zhì)的重視,政府統(tǒng)計部門更關(guān)注估計量是否覆蓋了研究總體。在人口普查凈誤差估計中,盡管用來估計總體實際人口數(shù)的雙系統(tǒng)估計量是一個有偏估計量,但1980年以來,部分國家一直使用雙系統(tǒng)估計量[23-25]。美國在人口普查質(zhì)量評估報告均未討論雙系統(tǒng)估計量的統(tǒng)計性質(zhì),主要原因在于不僅工作難度大,而且對實際工作并無多大益處。

    四、實證與模擬

    第一部分為實證分析,為政府統(tǒng)計部門提供組合式遺漏估計量及其抽樣方差、偏差和均方誤差估計量完整的計算程序,推廣組合式遺漏估計量;第二部分為模擬分析,討論單重遺漏估計量和組合式遺漏估計量的無偏性。

    (一)實證分析

    1.基本情況及數(shù)據(jù)來源

    以XXX省級單位XXX市XXX區(qū)XXX街道為實證范圍,估計該街道2010年普查遺漏人口數(shù)。采取分層二重抽樣,抽樣單位為普查小區(qū),從街道的100個普查小區(qū)中抽取8個。在第一重抽樣中,按城鄉(xiāng)屬性,將普查小區(qū)分為兩層,即城市層(h=1)和鄉(xiāng)村層(h=2),使用Nh表示層h的普查小區(qū)總數(shù),nh為從層h抽取的第一重樣本普查小區(qū)數(shù);在第二重抽樣中,對抽取的第一重樣本普查小區(qū),按照住房單元數(shù)目進(jìn)一步分為兩層,即70~90個層(g=1),以及70個以下和90個以上層(g=2),Mhg和mhg分別表示層hg的普查小區(qū)總數(shù)和樣本普查小區(qū)數(shù)。

    采用較為簡單的2X2分層模式有兩個原因。一是實證對象為郊區(qū),既有城市普查小區(qū),也有鄉(xiāng)村普查小區(qū),因此選擇城鄉(xiāng)屬性對普查小區(qū)分層具有合理性,中國國家統(tǒng)計局一直按照城鄉(xiāng)對普查小區(qū)分層;二是中國一個普查小區(qū)平均包括80個住房單元,根據(jù)普查小區(qū)實際規(guī)模對第一重樣本普查小區(qū)分為上面的兩個新g層,也具有一定合理性。實際中,采取何種分層模式與數(shù)據(jù)的可得性有直接關(guān)系。

    在覆蓋調(diào)查樣本普查小區(qū)抽取之后和使用組合式遺漏估計量估計普查遺漏人口數(shù)之前,需要做的一項工作是對總體人口使用性別、年齡、戶籍所在地、文化程度、民族等變量進(jìn)行等概率分層。不難看出,分層變量越多,等概率人口層內(nèi)部的同質(zhì)性就越強(qiáng),但是,在覆蓋調(diào)查樣本規(guī)模一定的情況下,分層變量過多會導(dǎo)致有些等概率人口層的抽樣方差過大,為避免抽樣方差過大,本文只選擇性別對總體人口分層。

    有關(guān)抽樣及其結(jié)果和樣本數(shù)據(jù),見表2和表3。

    表2 抽樣層、樣本和抽樣權(quán)數(shù)

    表3 層及樣本數(shù)據(jù)

    表2中的數(shù)據(jù)有兩個用途,一是用來計算表4中的單元的加權(quán)人數(shù),以及表5的遺漏估計值;二是用來計算表6進(jìn)入第二重樣本的每個樣本普查小區(qū)的復(fù)制權(quán)數(shù),以及表7的單重遺漏估計值和組合式遺漏估計值的抽樣方差。

    表4 等概率人口層的單元加權(quán)人數(shù)

    表5 等概率人口層及總體的遺漏估計值

    表6 第二重樣本普查小區(qū)復(fù)制權(quán)數(shù)

    表7 基于抽樣方差的變異系數(shù)估計值

    表3中的y10v,y01v,y11v是在比對同一樣本普查小區(qū)的普查人口名單和覆蓋調(diào)查人口名單,以及將比對結(jié)果劃分到男性層和女性層的結(jié)果。為提高匹配人口數(shù)y11v,比對范圍至少應(yīng)該擴(kuò)大到樣本普查小區(qū)的鄰近小區(qū)。由于沒有獲得鄰近普查小區(qū)的這兩份人口名單,比對只在樣本普查小區(qū)內(nèi)進(jìn)行,y11v可能偏低。

    2.估計結(jié)果

    根據(jù)式(9)和表2~3樣本數(shù)據(jù),計算式(8)等號右邊每個單元的加權(quán)人口數(shù),計算結(jié)果見表4:

    依據(jù)式(8)和表4計算男性層和女性層及總體的普查遺漏人口數(shù),計算結(jié)果見表5:

    從表5可以看出:(1)存在雙重遺漏人口。這便是提出組合式遺漏估計量的原因,雙重遺漏人口大多在普查和覆蓋調(diào)查期間外出經(jīng)商或旅游,或者故意躲避調(diào)查員。(2)如果采取組合式遺漏估計量,估計的總體遺漏人口數(shù)為173人;如果采取單重遺漏估計量,估計的總體遺漏人口數(shù)為170人??梢?單重遺漏估計量低估總體遺漏人口數(shù)3人。(3)無論是采取單重遺漏估計量還是組合式遺漏估計量,男性遺漏人口數(shù)均多于女性遺漏人口數(shù),例如,如果采取組合式遺漏估計量,男性遺漏人口數(shù)為107人,而女性只有66人,可見,相比女性,男性更容易在普查中遺漏,這源于兩方面的原因:一方面,男性外出經(jīng)商、務(wù)工的人數(shù)多于女性,男性是流動人口的主要人群,普查員上門登記時,難以遇到他們,因此男性比女性更容易在普查中遺漏;另一方面,男性對人口普查的重視程度不及女性,認(rèn)為人口普查與自身利益關(guān)系不大,因而不愿意抽出專門時間接受普查員調(diào)查,有意或無意躲避普查,造成遺漏,而女性比起男性更愿意接受調(diào)查。以上情況表明,在制定人口普查方案時,對男性應(yīng)該給予足夠的重視,采取有效措施防止男性普查遺漏。

    計算遺漏估計值的抽樣方差。使用表2和式(13)計算復(fù)制權(quán)數(shù),結(jié)果見表6:

    基于式(10)和式(15)以及表4~6數(shù)據(jù),得到等概率人口層及總體的單重和組合式遺漏估計值的抽樣方差等,結(jié)果見表7:

    從表7可以看出,使用單重遺漏估計量估計的男性層、女性層及總體的變異系數(shù)分別為0.18、0.34、0.20;使用組合式遺漏估計量估計的變異系數(shù)分別為0.19、0.33、0.19。除男性層外,使用組合式遺漏估計量得到的女性和總體的遺漏人口數(shù)估計值的變異系數(shù)均小于相應(yīng)的單重遺漏估計量,表明組合式遺漏估計量的估計精度高于單重遺漏估計量。因此,應(yīng)該選擇組合式遺漏估計量。

    (二)模擬分析

    為了進(jìn)一步比較單重遺漏估計量和組合式遺漏估計量的抽樣估計精度,對這兩個有偏估計量,要考慮其偏差,不能只依據(jù)其抽樣方差來判斷其抽樣估計精度,而要使用均方誤差表示其抽樣估計精度,采用模擬方法實現(xiàn)這個目標(biāo)。根據(jù)前文理論中模擬分析的步驟,計算單重遺漏估計量和組合式遺漏估計量的偏差,討論單重遺漏估計量和組合式遺漏估計量的無偏性問題。

    由于單重漏登估計量及雙重漏登估計量的偏差均不為零,需進(jìn)一步計算其均方誤差及變異系數(shù)。結(jié)果見表8。單重遺漏估計量男性層、女性層、總體的均方誤差分別為546.03、579.63、1 629.7;組合式遺漏估計量男性層、女性層、總體的均方誤差分別為496.89、551.81、1 456.82。單重遺漏估計量男性層、女性層、總體的變異系數(shù)分別為0.198、0.321、0.209;組合式遺漏估計量男性層、女性層、總體的變異系數(shù)分別為0.191、0.317、0.200。

    表8 基于均方誤差的變異系數(shù)估計值

    表8表明三點重要信息,一是單重遺漏估計量及組合式遺漏估計量,其偏差均不為零,因而它們都是有偏估計量,應(yīng)使用均方誤差表示其抽樣估計精度;二是無論是男性層、女性層以及總體,使用單重遺漏估計量的偏差均大于組合式遺漏估計量的偏差,例如,男性層使用單重遺漏估計量的偏差為13人,組合式遺漏估計量的偏差為10人,與其未包括雙重遺漏估計量有直接關(guān)系,導(dǎo)致估計結(jié)果與真實值差距較大;三是基于均方誤差計算的變異系數(shù)來看,單重遺漏估計量估計的男性層、女性層及總體的變異系數(shù)分別為0.198、0.321、0.209,而組合式遺漏估計量的變異系數(shù)分別為0.191、0.317、0.200,可見組合式遺漏估計量的變異系數(shù)均小于單重遺漏估計量,表明采用均方誤差表示抽樣估計精度,組合式遺漏估計量仍然優(yōu)于單重遺漏估計量。

    五、結(jié) 語

    雖然政府統(tǒng)計部門設(shè)法在人口普查中登記完全,但普查遺漏無法避免,當(dāng)遺漏的人口數(shù)比重復(fù)登記的人口數(shù)多時,表現(xiàn)為凈遺漏,當(dāng)在普查中重復(fù)登記的人口數(shù)多于遺漏的人口數(shù)時,表現(xiàn)為凈多報。遺漏是凈遺漏的主要構(gòu)成部分,由于重復(fù)登記可能抵消遺漏,所以遺漏比起凈遺漏更能反映普查數(shù)據(jù)的質(zhì)量。凈遺漏為零,可能意味著沒有一個人遺漏,也沒有一個人重復(fù)登記,或者是遺漏和重復(fù)登記相互抵消,也就是說,從凈遺漏中無法區(qū)分出遺漏和重復(fù)登記數(shù)量。政府統(tǒng)計部門組織人口普查質(zhì)量評估工作的目標(biāo)主要是通過估計的遺漏和重復(fù)登記數(shù)目來發(fā)現(xiàn)人口普查登記工作中的漏洞,從而更好地完成下次人口普查登記工作任務(wù)。

    使用組合式遺漏估計量替代單重遺漏估計量有其必然性,然而,由于它需要對總體人口等概率分層,所以完成這種替代需要政府統(tǒng)計部門根據(jù)本國人口特點和覆蓋調(diào)查樣本規(guī)模設(shè)計相適應(yīng)的分層方案,以減少其異質(zhì)性偏差。

    組合式遺漏估計量優(yōu)勢明顯。首先,除了包括登記在覆蓋調(diào)查而未登記在普查中的單重遺漏人口外,還包括同時遺漏于這兩項調(diào)查的雙重遺漏人口,因而估計值更接近真實值;其次,它不只能夠查明單重遺漏人口的特征、居住位置及其遺漏程度,還能查明雙重遺漏人口的這些情況,這對未來普查操作方案的改進(jìn)具有重要意義。單重遺漏人口信息登記在覆蓋調(diào)查表,檢查該調(diào)查表可以知悉其姓名、性別、年齡、受教育程度、與戶主關(guān)系、普查時點居住地、覆蓋調(diào)查時點居住地、戶籍所在地。通過研究區(qū)域內(nèi)的負(fù)責(zé)人組織的入戶調(diào)查等手段可以查找到雙重遺漏人口,并獲悉其人口統(tǒng)計特征變量值和居住位置變量值。

    組合式遺漏估計量由單重遺漏估計量和覆蓋調(diào)查與普查相互獨立情況下的雙重遺漏估計量組成。單重遺漏估計量構(gòu)造的關(guān)鍵是這兩項調(diào)查名單的比對質(zhì)量。如果比對程序不合理,或者兩項調(diào)查名單登記的人口信息不完整,就可能影響比對效果,錯誤地將匹配人口當(dāng)作未匹配人口,或者將未匹配人口當(dāng)作匹配人口,從而產(chǎn)生比對誤差。雙重遺漏估計量建立的前提是覆蓋調(diào)查與普查獨立,否則產(chǎn)生交互作用偏差。交互作用偏差源于普查與覆蓋調(diào)查之間的因果相關(guān)性以及等概率人口層中的人口在普查與覆蓋調(diào)查中登記概率的異質(zhì)性。因此,在使用組合式遺漏估計量之前,要考慮到比對誤差和交互作用偏差是否存在,如果存在,則要采取恰當(dāng)方法將其列入,否則將低估或高估總體普查遺漏人口數(shù)。

    猜你喜歡
    估計量人口數(shù)遺漏
    來自動物星球的挑戰(zhàn)(二)小五狼遺漏的線索
    確定代表名額按戶籍人口數(shù)還是常住人口數(shù)
    人大研究(2022年3期)2022-04-13 00:47:04
    遺漏的光陰
    鴨綠江(2021年17期)2021-11-11 13:03:41
    基于BP神經(jīng)網(wǎng)絡(luò)的吉林省GDP預(yù)測研究
    市場周刊(2018年1期)2018-08-15 01:42:44
    淺談估計量的優(yōu)良性標(biāo)準(zhǔn)
    應(yīng)用品管圈降低腹腔鏡抗反流手術(shù)術(shù)前準(zhǔn)備遺漏率的實踐
    基于配網(wǎng)先驗信息的諧波狀態(tài)估計量測點最優(yōu)配置
    電測與儀表(2015年6期)2015-04-09 12:00:50
    你會收集數(shù)據(jù)嗎
    負(fù)極值指標(biāo)估計量的漸近性質(zhì)
    使用變異系數(shù)和Kurtosis系數(shù)的雙輔助變量的比估計
    免费在线观看日本一区| 欧美性猛交黑人性爽| 日本黄大片高清| 免费无遮挡裸体视频| 久久久久久久久中文| 99久国产av精品| 亚洲av免费高清在线观看| 欧美性感艳星| 哪里可以看免费的av片| 欧美在线黄色| 午夜福利高清视频| 中文字幕av在线有码专区| 中出人妻视频一区二区| www.www免费av| 男女之事视频高清在线观看| 国产三级黄色录像| 99热只有精品国产| 黄色女人牲交| 久久精品久久久久久噜噜老黄 | 精华霜和精华液先用哪个| 欧美精品啪啪一区二区三区| 亚洲成av人片在线播放无| 黄色丝袜av网址大全| 99精品久久久久人妻精品| 毛片女人毛片| 69av精品久久久久久| 一级毛片久久久久久久久女| 亚洲自偷自拍三级| 九九久久精品国产亚洲av麻豆| 亚洲成人精品中文字幕电影| 1000部很黄的大片| 成年版毛片免费区| 中文字幕av在线有码专区| 又黄又爽又刺激的免费视频.| 直男gayav资源| 婷婷丁香在线五月| 欧美色欧美亚洲另类二区| 一本久久中文字幕| 久久精品综合一区二区三区| 午夜福利欧美成人| 久久精品国产亚洲av香蕉五月| 欧美午夜高清在线| 久久精品综合一区二区三区| 乱人视频在线观看| 超碰av人人做人人爽久久| 男女视频在线观看网站免费| 亚洲国产精品成人综合色| 乱人视频在线观看| 久久久国产成人免费| 日本三级黄在线观看| 免费在线观看日本一区| 欧美激情国产日韩精品一区| 免费av不卡在线播放| 男人舔奶头视频| 黄色女人牲交| 久久久国产成人精品二区| 别揉我奶头 嗯啊视频| 日本在线视频免费播放| 精品99又大又爽又粗少妇毛片 | 国产高清视频在线观看网站| 午夜福利欧美成人| 亚州av有码| 国产黄a三级三级三级人| 成人特级av手机在线观看| 69av精品久久久久久| xxxwww97欧美| 国产在线精品亚洲第一网站| 日韩免费av在线播放| 国产黄a三级三级三级人| 成人av在线播放网站| 国产高清三级在线| 嫩草影视91久久| 久久久久久大精品| 成人亚洲精品av一区二区| 伊人久久精品亚洲午夜| 国产伦精品一区二区三区视频9| 亚洲国产精品成人综合色| 亚洲国产精品久久男人天堂| 亚洲精品一卡2卡三卡4卡5卡| 成人三级黄色视频| 国产精品久久视频播放| 亚洲在线观看片| 夜夜看夜夜爽夜夜摸| 国产精品,欧美在线| 桃红色精品国产亚洲av| 看免费av毛片| 亚洲美女黄片视频| 欧美成人a在线观看| 丰满人妻熟妇乱又伦精品不卡| 中文亚洲av片在线观看爽| 久久中文看片网| 禁无遮挡网站| 人人妻人人看人人澡| 人人妻人人澡欧美一区二区| 日日干狠狠操夜夜爽| 欧美极品一区二区三区四区| 午夜福利在线观看免费完整高清在 | 在线观看美女被高潮喷水网站 | 久久热精品热| 51午夜福利影视在线观看| 中出人妻视频一区二区| 黄色日韩在线| 床上黄色一级片| 亚洲最大成人手机在线| 赤兔流量卡办理| 亚洲精品在线观看二区| 性插视频无遮挡在线免费观看| 精品久久久久久久久久久久久| 欧美区成人在线视频| www.色视频.com| 禁无遮挡网站| 一本综合久久免费| 久久久久久久久大av| 18禁裸乳无遮挡免费网站照片| 国产成人a区在线观看| 99久久久亚洲精品蜜臀av| 亚洲av成人av| 久久99热6这里只有精品| 搡老熟女国产l中国老女人| 成人三级黄色视频| 色播亚洲综合网| 国产精品永久免费网站| 成年女人看的毛片在线观看| 亚洲18禁久久av| 99在线视频只有这里精品首页| 国产色婷婷99| 亚洲欧美日韩无卡精品| 日韩av在线大香蕉| 国产探花在线观看一区二区| 91av网一区二区| 国产精品久久久久久人妻精品电影| 99热这里只有精品一区| 变态另类成人亚洲欧美熟女| 淫秽高清视频在线观看| 直男gayav资源| 九九在线视频观看精品| 久久这里只有精品中国| 性色av乱码一区二区三区2| aaaaa片日本免费| 欧美潮喷喷水| 日本熟妇午夜| 神马国产精品三级电影在线观看| 亚洲激情在线av| 听说在线观看完整版免费高清| 女生性感内裤真人,穿戴方法视频| 两个人视频免费观看高清| 亚洲人与动物交配视频| 又爽又黄a免费视频| 中文字幕av在线有码专区| 亚洲 欧美 日韩 在线 免费| 国产午夜精品久久久久久一区二区三区 | 男女视频在线观看网站免费| 人人妻人人看人人澡| 亚洲精品成人久久久久久| 一边摸一边抽搐一进一小说| 91在线观看av| www.色视频.com| 久99久视频精品免费| 国产午夜精品论理片| 少妇丰满av| 男插女下体视频免费在线播放| 亚洲中文日韩欧美视频| 久久精品国产清高在天天线| xxxwww97欧美| 国产一级毛片七仙女欲春2| 午夜福利免费观看在线| 99久久成人亚洲精品观看| 我的老师免费观看完整版| 亚洲狠狠婷婷综合久久图片| 免费黄网站久久成人精品 | 久久久久性生活片| 九九在线视频观看精品| 天美传媒精品一区二区| 亚洲在线观看片| 亚洲精品一区av在线观看| 女人十人毛片免费观看3o分钟| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 夜夜爽天天搞| 中文字幕高清在线视频| 757午夜福利合集在线观看| 日韩欧美三级三区| 狂野欧美白嫩少妇大欣赏| 一个人看视频在线观看www免费| 午夜a级毛片| 久久久色成人| 伦理电影大哥的女人| 色吧在线观看| 在线观看午夜福利视频| 色综合欧美亚洲国产小说| 麻豆一二三区av精品| 亚洲av二区三区四区| 一级a爱片免费观看的视频| or卡值多少钱| 亚洲不卡免费看| 麻豆久久精品国产亚洲av| 深爱激情五月婷婷| 欧美最黄视频在线播放免费| 久久精品夜夜夜夜夜久久蜜豆| 精品欧美国产一区二区三| 欧美+日韩+精品| 久久精品国产清高在天天线| 如何舔出高潮| 亚洲精品在线观看二区| 在现免费观看毛片| 亚洲经典国产精华液单 | 老司机深夜福利视频在线观看| 国产欧美日韩一区二区三| 国产麻豆成人av免费视频| 国产毛片a区久久久久| 免费高清视频大片| 又爽又黄a免费视频| 精品人妻1区二区| 亚洲五月天丁香| a级一级毛片免费在线观看| eeuss影院久久| 最近在线观看免费完整版| 日日摸夜夜添夜夜添av毛片 | 日韩欧美免费精品| 真实男女啪啪啪动态图| 97人妻精品一区二区三区麻豆| 国产亚洲精品久久久com| 特级一级黄色大片| 亚洲精品色激情综合| 精品日产1卡2卡| 成人毛片a级毛片在线播放| 免费人成视频x8x8入口观看| 国产欧美日韩精品一区二区| 免费av毛片视频| 麻豆av噜噜一区二区三区| 国产亚洲精品久久久久久毛片| 久久久精品欧美日韩精品| www.熟女人妻精品国产| eeuss影院久久| 国产成人aa在线观看| 国产精品1区2区在线观看.| 琪琪午夜伦伦电影理论片6080| 内地一区二区视频在线| 99国产极品粉嫩在线观看| 97热精品久久久久久| a级一级毛片免费在线观看| 热99在线观看视频| 性欧美人与动物交配| 成人av一区二区三区在线看| 中亚洲国语对白在线视频| 99精品在免费线老司机午夜| 我的女老师完整版在线观看| 色综合欧美亚洲国产小说| 国产高潮美女av| 欧美zozozo另类| 免费人成在线观看视频色| 美女大奶头视频| 成年人黄色毛片网站| 欧美绝顶高潮抽搐喷水| 99热6这里只有精品| 少妇高潮的动态图| 露出奶头的视频| 亚洲国产精品久久男人天堂| 国产真实伦视频高清在线观看 | 国产精品美女特级片免费视频播放器| 欧美最新免费一区二区三区 | 国产v大片淫在线免费观看| 1000部很黄的大片| 18禁裸乳无遮挡免费网站照片| 一本综合久久免费| 成年女人永久免费观看视频| 久久国产精品人妻蜜桃| eeuss影院久久| 久99久视频精品免费| 免费av不卡在线播放| 最近视频中文字幕2019在线8| 亚洲天堂国产精品一区在线| 日韩中字成人| 国产视频内射| 99久久久亚洲精品蜜臀av| 中亚洲国语对白在线视频| 在线观看舔阴道视频| 美女黄网站色视频| 欧美另类亚洲清纯唯美| av在线观看视频网站免费| 日韩欧美精品免费久久 | 亚洲美女黄片视频| 狂野欧美白嫩少妇大欣赏| 少妇的逼水好多| av在线观看视频网站免费| 男人和女人高潮做爰伦理| 激情在线观看视频在线高清| or卡值多少钱| 欧美最黄视频在线播放免费| 亚州av有码| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品亚洲一区二区| 午夜精品在线福利| 黄色日韩在线| 一个人免费在线观看电影| 免费观看精品视频网站| 亚州av有码| 精品久久久久久成人av| 深夜a级毛片| 18禁黄网站禁片午夜丰满| 午夜两性在线视频| 久久亚洲真实| 伊人久久精品亚洲午夜| 国产精品女同一区二区软件 | 国语自产精品视频在线第100页| 国产精品日韩av在线免费观看| 国产av一区在线观看免费| 免费看a级黄色片| 久久久色成人| 亚洲第一区二区三区不卡| 国产精品久久久久久人妻精品电影| 国产精品人妻久久久久久| 精品一区二区免费观看| 欧美国产日韩亚洲一区| 一a级毛片在线观看| 99久久精品国产亚洲精品| 国产午夜福利久久久久久| 欧美成人a在线观看| 熟妇人妻久久中文字幕3abv| 国产精品自产拍在线观看55亚洲| 麻豆一二三区av精品| 久久久久久久久久黄片| 丰满人妻熟妇乱又伦精品不卡| 久久国产乱子伦精品免费另类| av女优亚洲男人天堂| 在线a可以看的网站| 日本a在线网址| 国产熟女xx| 90打野战视频偷拍视频| 99精品在免费线老司机午夜| 女生性感内裤真人,穿戴方法视频| 精品一区二区三区av网在线观看| 香蕉av资源在线| 欧美另类亚洲清纯唯美| 少妇人妻一区二区三区视频| 成年免费大片在线观看| 99久久精品一区二区三区| 久久久久久久久久黄片| 久久人妻av系列| 九九久久精品国产亚洲av麻豆| 蜜桃亚洲精品一区二区三区| 午夜福利在线在线| 99热这里只有是精品在线观看 | 亚洲精品456在线播放app | 乱人视频在线观看| 一本精品99久久精品77| 免费看a级黄色片| 性色av乱码一区二区三区2| 中文字幕免费在线视频6| 人妻夜夜爽99麻豆av| 国产精品久久视频播放| 国产亚洲精品久久久久久毛片| 亚洲精品影视一区二区三区av| 国产精品不卡视频一区二区 | 757午夜福利合集在线观看| 美女xxoo啪啪120秒动态图 | 一边摸一边抽搐一进一小说| 精品不卡国产一区二区三区| 99riav亚洲国产免费| 国产精品98久久久久久宅男小说| 久久久久国内视频| 精品不卡国产一区二区三区| 青草久久国产| 亚洲一区二区三区色噜噜| 久久九九热精品免费| 变态另类成人亚洲欧美熟女| 午夜福利视频1000在线观看| 婷婷亚洲欧美| 国产av不卡久久| 亚洲欧美清纯卡通| 亚洲国产欧洲综合997久久,| 露出奶头的视频| 青草久久国产| 成人欧美大片| 亚洲精品一区av在线观看| 中国美女看黄片| 日本一本二区三区精品| 日本与韩国留学比较| 嫁个100分男人电影在线观看| 精品一区二区三区av网在线观看| 亚洲真实伦在线观看| 亚洲av熟女| 十八禁国产超污无遮挡网站| 91字幕亚洲| 男女视频在线观看网站免费| 国产精品野战在线观看| 精品久久久久久久久亚洲 | 狠狠狠狠99中文字幕| 淫妇啪啪啪对白视频| 亚洲国产精品sss在线观看| 国产视频内射| 日韩中字成人| 久久伊人香网站| 国产欧美日韩一区二区三| 精品久久久久久久久久久久久| 99久久精品热视频| 国产大屁股一区二区在线视频| 99久久精品热视频| 久久国产精品人妻蜜桃| 色哟哟·www| 99精品在免费线老司机午夜| 成人特级av手机在线观看| 亚洲午夜理论影院| 观看美女的网站| 日韩国内少妇激情av| 欧美在线一区亚洲| 在线a可以看的网站| 桃红色精品国产亚洲av| 两个人的视频大全免费| 久久久久性生活片| 久久国产乱子免费精品| 丰满人妻熟妇乱又伦精品不卡| 他把我摸到了高潮在线观看| 国产 一区 欧美 日韩| 欧美精品国产亚洲| 少妇的逼水好多| 听说在线观看完整版免费高清| 国产免费男女视频| 久久亚洲真实| 全区人妻精品视频| 成人性生交大片免费视频hd| 人人妻,人人澡人人爽秒播| 国产精品亚洲av一区麻豆| 亚洲成人中文字幕在线播放| 日本在线视频免费播放| 欧美黄色片欧美黄色片| 99热精品在线国产| 欧美一区二区亚洲| 99久久九九国产精品国产免费| 成人无遮挡网站| 99视频精品全部免费 在线| 一本久久中文字幕| 每晚都被弄得嗷嗷叫到高潮| 国内精品久久久久久久电影| 无遮挡黄片免费观看| 亚洲人成网站在线播| 久久99热这里只有精品18| 乱码一卡2卡4卡精品| 日韩欧美三级三区| 日韩欧美在线乱码| 亚洲七黄色美女视频| av天堂中文字幕网| 麻豆国产av国片精品| 欧美一区二区精品小视频在线| 午夜激情欧美在线| 亚洲成人精品中文字幕电影| 直男gayav资源| 亚州av有码| 窝窝影院91人妻| 一个人免费在线观看的高清视频| 亚洲一区二区三区不卡视频| 日韩av在线大香蕉| 动漫黄色视频在线观看| 国产成人aa在线观看| 亚洲精品一区av在线观看| 午夜福利成人在线免费观看| 一区福利在线观看| 亚洲三级黄色毛片| 俄罗斯特黄特色一大片| 成年版毛片免费区| 免费观看的影片在线观看| 老熟妇仑乱视频hdxx| 在线观看66精品国产| 老鸭窝网址在线观看| 中文亚洲av片在线观看爽| 最近在线观看免费完整版| 国产av在哪里看| 亚洲精品456在线播放app | www.www免费av| 亚洲精品在线美女| av女优亚洲男人天堂| 欧美日本亚洲视频在线播放| 欧洲精品卡2卡3卡4卡5卡区| 变态另类丝袜制服| 欧美日韩国产亚洲二区| 久久久久国内视频| 中文字幕av成人在线电影| 亚洲精品456在线播放app | 青草久久国产| 国产精品久久久久久亚洲av鲁大| 亚洲人成网站在线播| 香蕉av资源在线| 成人一区二区视频在线观看| 我要搜黄色片| 国产黄片美女视频| 可以在线观看的亚洲视频| 人人妻人人澡欧美一区二区| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久久久黄片| netflix在线观看网站| 欧美zozozo另类| 日韩中文字幕欧美一区二区| 非洲黑人性xxxx精品又粗又长| 美女大奶头视频| 天堂影院成人在线观看| 亚洲狠狠婷婷综合久久图片| 国产乱人伦免费视频| 午夜精品久久久久久毛片777| 日本与韩国留学比较| 精品人妻一区二区三区麻豆 | av中文乱码字幕在线| 亚洲人成网站在线播| 久久精品影院6| 国产亚洲精品综合一区在线观看| 成人午夜高清在线视频| 人妻夜夜爽99麻豆av| 很黄的视频免费| av欧美777| 日韩大尺度精品在线看网址| 免费观看人在逋| 亚洲中文字幕日韩| 最新在线观看一区二区三区| 国产欧美日韩一区二区三| 大型黄色视频在线免费观看| 国产精品伦人一区二区| 日韩欧美三级三区| 一本综合久久免费| 欧美在线一区亚洲| ponron亚洲| 成人鲁丝片一二三区免费| 国产精品久久久久久精品电影| 啦啦啦观看免费观看视频高清| 欧美一区二区亚洲| 午夜两性在线视频| 在线看三级毛片| 日韩欧美一区二区三区在线观看| 黄色丝袜av网址大全| 美女黄网站色视频| 久久人妻av系列| 免费观看精品视频网站| 精品乱码久久久久久99久播| 久久精品国产亚洲av天美| 中文字幕av成人在线电影| 国产精品美女特级片免费视频播放器| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人免费| 国产一区二区三区视频了| 久久性视频一级片| 白带黄色成豆腐渣| 国产一区二区三区在线臀色熟女| 天堂av国产一区二区熟女人妻| 精品乱码久久久久久99久播| 高清在线国产一区| 精品久久久久久久久亚洲 | 日韩有码中文字幕| 不卡一级毛片| 久久久久国内视频| 国产精品98久久久久久宅男小说| 亚洲av熟女| 国产乱人伦免费视频| 好男人在线观看高清免费视频| 自拍偷自拍亚洲精品老妇| 一本精品99久久精品77| 一个人免费在线观看的高清视频| 成年版毛片免费区| 亚洲欧美激情综合另类| 精品一区二区免费观看| 好看av亚洲va欧美ⅴa在| 欧美成人一区二区免费高清观看| 国产乱人视频| 天堂影院成人在线观看| 国产精品自产拍在线观看55亚洲| 成人特级黄色片久久久久久久| 亚洲人成伊人成综合网2020| 国产精品免费一区二区三区在线| 小蜜桃在线观看免费完整版高清| 性色avwww在线观看| 最好的美女福利视频网| 级片在线观看| 男女之事视频高清在线观看| 欧美三级亚洲精品| 免费看日本二区| 国产成年人精品一区二区| 亚洲国产精品合色在线| 麻豆国产av国片精品| 搡女人真爽免费视频火全软件 | 国产成人av教育| 精品国产三级普通话版| 91在线观看av| 国产精华一区二区三区| 欧美bdsm另类| 国产精品人妻久久久久久| 性色avwww在线观看| 亚洲精品日韩av片在线观看| 亚洲在线观看片| 国产精品久久久久久人妻精品电影| 亚洲av成人不卡在线观看播放网| 国产在视频线在精品| 精品免费久久久久久久清纯| 国产高潮美女av| 国产av一区在线观看免费| 99国产综合亚洲精品| 我的老师免费观看完整版| 亚洲男人的天堂狠狠| 国产精品自产拍在线观看55亚洲| 一区二区三区激情视频| 老司机午夜福利在线观看视频| 床上黄色一级片| 国产老妇女一区| 又爽又黄无遮挡网站| 精品久久国产蜜桃| 国产乱人伦免费视频| 国产69精品久久久久777片| 亚洲熟妇中文字幕五十中出| h日本视频在线播放| 很黄的视频免费| 十八禁网站免费在线| 久久久久国产精品人妻aⅴ院| 亚洲人与动物交配视频| 亚洲av.av天堂| 亚洲一区二区三区不卡视频| 简卡轻食公司| 黄片小视频在线播放| 亚洲国产精品合色在线| x7x7x7水蜜桃| 一二三四社区在线视频社区8|