• <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ù)的雙輔助變量的比估計
    宅男免费午夜| 舔av片在线| 久久精品91无色码中文字幕| 成年免费大片在线观看| 可以在线观看毛片的网站| 国产黄片美女视频| www日本黄色视频网| 舔av片在线| 女警被强在线播放| 中文字幕最新亚洲高清| 精品国产亚洲在线| 波多野结衣高清作品| 久久久久国内视频| 美女 人体艺术 gogo| 久久热在线av| 日韩中文字幕欧美一区二区| 亚洲激情在线av| 国产精品久久久久久人妻精品电影| 国产视频内射| 在线观看舔阴道视频| 国产精品综合久久久久久久免费| 亚洲国产欧美人成| 好男人在线观看高清免费视频| 正在播放国产对白刺激| 黑人巨大精品欧美一区二区mp4| av欧美777| 久久精品国产亚洲av高清一级| 97超级碰碰碰精品色视频在线观看| 免费一级毛片在线播放高清视频| 久久欧美精品欧美久久欧美| 精品国产乱码久久久久久男人| 午夜老司机福利片| 真人做人爱边吃奶动态| 此物有八面人人有两片| 麻豆国产97在线/欧美 | 久久伊人香网站| 亚洲最大成人中文| 最好的美女福利视频网| 亚洲av第一区精品v没综合| 操出白浆在线播放| 国产亚洲精品一区二区www| 午夜福利在线在线| 国产亚洲精品久久久久5区| 久久久久性生活片| 欧美 亚洲 国产 日韩一| 人人妻人人澡欧美一区二区| 亚洲av熟女| 精品午夜福利视频在线观看一区| 777久久人妻少妇嫩草av网站| 99re在线观看精品视频| 最近最新免费中文字幕在线| 国产av一区在线观看免费| 亚洲精品一卡2卡三卡4卡5卡| av在线播放免费不卡| 国产高清有码在线观看视频 | 999久久久国产精品视频| 日本精品一区二区三区蜜桃| bbb黄色大片| 一级毛片高清免费大全| 精品久久久久久,| 免费观看人在逋| 亚洲一区高清亚洲精品| 亚洲最大成人中文| 亚洲成人国产一区在线观看| 一级作爱视频免费观看| 欧美日本视频| 99re在线观看精品视频| 男人的好看免费观看在线视频 | 精品久久久久久久毛片微露脸| 久久久国产成人免费| 精品国内亚洲2022精品成人| 99国产极品粉嫩在线观看| 精品久久久久久久毛片微露脸| 国产日本99.免费观看| 亚洲午夜理论影院| 人人妻人人澡欧美一区二区| 久久久精品欧美日韩精品| 日本撒尿小便嘘嘘汇集6| 搡老熟女国产l中国老女人| 亚洲精品一区av在线观看| 国产精品av久久久久免费| 国产亚洲欧美98| 久久精品亚洲精品国产色婷小说| 男女做爰动态图高潮gif福利片| 久久亚洲真实| 在线播放国产精品三级| 熟女电影av网| 亚洲精品粉嫩美女一区| 欧美日韩一级在线毛片| 搡老熟女国产l中国老女人| 日韩欧美精品v在线| 欧美日韩中文字幕国产精品一区二区三区| 男女视频在线观看网站免费 | 神马国产精品三级电影在线观看 | 色综合欧美亚洲国产小说| 两性午夜刺激爽爽歪歪视频在线观看 | 人人妻,人人澡人人爽秒播| or卡值多少钱| 美女扒开内裤让男人捅视频| 麻豆久久精品国产亚洲av| 亚洲精品粉嫩美女一区| 欧美人与性动交α欧美精品济南到| 熟女少妇亚洲综合色aaa.| 亚洲人成网站高清观看| 久久精品亚洲精品国产色婷小说| 久久国产精品影院| 大型黄色视频在线免费观看| 国产不卡一卡二| 19禁男女啪啪无遮挡网站| 一级片免费观看大全| 欧美人与性动交α欧美精品济南到| 91麻豆av在线| 一级毛片精品| 欧美日韩精品网址| 精品电影一区二区在线| 一区福利在线观看| or卡值多少钱| 两性夫妻黄色片| 久久久久国产一级毛片高清牌| 成人一区二区视频在线观看| av天堂在线播放| 麻豆国产97在线/欧美 | 久久久久久国产a免费观看| 亚洲18禁久久av| 久久久久久国产a免费观看| a级毛片在线看网站| 日本撒尿小便嘘嘘汇集6| 欧美最黄视频在线播放免费| 在线观看免费视频日本深夜| 亚洲成人精品中文字幕电影| 亚洲精品色激情综合| 岛国在线免费视频观看| 欧美另类亚洲清纯唯美| 国产黄片美女视频| www.自偷自拍.com| 日韩精品免费视频一区二区三区| 听说在线观看完整版免费高清| 18美女黄网站色大片免费观看| 国产精品久久电影中文字幕| 国产探花在线观看一区二区| 变态另类成人亚洲欧美熟女| 国内揄拍国产精品人妻在线| 午夜福利免费观看在线| 男女之事视频高清在线观看| tocl精华| 国产亚洲精品一区二区www| 一区福利在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 香蕉久久夜色| 欧美黄色片欧美黄色片| 中出人妻视频一区二区| 成人av一区二区三区在线看| 亚洲九九香蕉| 久久久久久久久免费视频了| 国产在线观看jvid| 别揉我奶头~嗯~啊~动态视频| netflix在线观看网站| 啦啦啦免费观看视频1| 欧美成人一区二区免费高清观看 | 国产黄片美女视频| 露出奶头的视频| 色av中文字幕| 亚洲全国av大片| а√天堂www在线а√下载| 国内久久婷婷六月综合欲色啪| 不卡一级毛片| 成人午夜高清在线视频| 看片在线看免费视频| 在线视频色国产色| 亚洲人成电影免费在线| av免费在线观看网站| 天堂√8在线中文| 国产97色在线日韩免费| 国产av麻豆久久久久久久| 欧美不卡视频在线免费观看 | 1024视频免费在线观看| 女人被狂操c到高潮| 精品一区二区三区av网在线观看| 不卡av一区二区三区| 成人特级黄色片久久久久久久| 国产精品免费视频内射| 老司机午夜福利在线观看视频| 亚洲人成网站高清观看| 成人高潮视频无遮挡免费网站| 999精品在线视频| 久久亚洲精品不卡| 亚洲欧美一区二区三区黑人| 国内精品久久久久久久电影| 观看免费一级毛片| 国产精品一区二区精品视频观看| 国产成人系列免费观看| 我要搜黄色片| 99热6这里只有精品| 脱女人内裤的视频| 悠悠久久av| 长腿黑丝高跟| 亚洲av第一区精品v没综合| 精品欧美国产一区二区三| 成人18禁在线播放| 亚洲va日本ⅴa欧美va伊人久久| 又粗又爽又猛毛片免费看| 哪里可以看免费的av片| 麻豆国产av国片精品| 91老司机精品| 99re在线观看精品视频| 亚洲人成伊人成综合网2020| 一边摸一边做爽爽视频免费| 大型av网站在线播放| 看免费av毛片| 国产免费av片在线观看野外av| 18禁黄网站禁片午夜丰满| 欧美精品亚洲一区二区| 精品国产乱码久久久久久男人| 久久人妻av系列| 午夜福利18| 欧美成人一区二区免费高清观看 | 亚洲专区国产一区二区| 欧美日韩乱码在线| 一级毛片精品| 欧美黑人欧美精品刺激| 国产精品永久免费网站| 神马国产精品三级电影在线观看 | 97碰自拍视频| 后天国语完整版免费观看| 欧美日本视频| 19禁男女啪啪无遮挡网站| 亚洲国产中文字幕在线视频| 最好的美女福利视频网| 免费在线观看黄色视频的| 国产精品98久久久久久宅男小说| 午夜亚洲福利在线播放| 看片在线看免费视频| 丰满的人妻完整版| 欧美一级毛片孕妇| 国产男靠女视频免费网站| 黄片大片在线免费观看| 国产精品99久久99久久久不卡| 露出奶头的视频| 国产aⅴ精品一区二区三区波| av中文乱码字幕在线| 俺也久久电影网| 亚洲专区字幕在线| 午夜福利免费观看在线| 欧美最黄视频在线播放免费| 免费在线观看完整版高清| 国产69精品久久久久777片 | 视频区欧美日本亚洲| 色尼玛亚洲综合影院| 人人妻,人人澡人人爽秒播| 亚洲五月婷婷丁香| 一级黄色大片毛片| 岛国在线免费视频观看| 一本大道久久a久久精品| 国产精品一区二区三区四区久久| 国产av一区在线观看免费| 男女午夜视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美三级亚洲精品| 精品久久久久久久末码| 不卡av一区二区三区| 亚洲中文av在线| 日韩精品免费视频一区二区三区| 国产精品亚洲av一区麻豆| 国产一区二区三区视频了| 亚洲一码二码三码区别大吗| 欧美乱色亚洲激情| 欧美性猛交╳xxx乱大交人| 国产成人av激情在线播放| 国产又色又爽无遮挡免费看| 亚洲免费av在线视频| 日韩欧美 国产精品| 一级片免费观看大全| 亚洲一区二区三区不卡视频| 露出奶头的视频| 精品国产乱子伦一区二区三区| 在线观看一区二区三区| 热99re8久久精品国产| 婷婷精品国产亚洲av| 99热只有精品国产| 1024手机看黄色片| 欧美最黄视频在线播放免费| 大型黄色视频在线免费观看| 极品教师在线免费播放| 久久久久国产一级毛片高清牌| 婷婷亚洲欧美| 床上黄色一级片| 久久久久精品国产欧美久久久| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲五月婷婷丁香| 丰满人妻熟妇乱又伦精品不卡| 两性夫妻黄色片| 日日爽夜夜爽网站| 国产一区在线观看成人免费| 制服诱惑二区| 99国产精品一区二区三区| 精品乱码久久久久久99久播| 久久久久久亚洲精品国产蜜桃av| 一本大道久久a久久精品| 日韩 欧美 亚洲 中文字幕| 丝袜美腿诱惑在线| 欧美3d第一页| 无限看片的www在线观看| 精华霜和精华液先用哪个| 两个人视频免费观看高清| 一本精品99久久精品77| 久久精品人妻少妇| 日韩av在线大香蕉| 黄色a级毛片大全视频| 日韩精品青青久久久久久| 最近最新免费中文字幕在线| 国产精品野战在线观看| 免费看a级黄色片| 18禁裸乳无遮挡免费网站照片| 国产成人精品久久二区二区免费| 97碰自拍视频| 嫩草影院精品99| 日本三级黄在线观看| 亚洲中文字幕一区二区三区有码在线看 | 色老头精品视频在线观看| av在线播放免费不卡| tocl精华| 国产亚洲精品一区二区www| 美女黄网站色视频| e午夜精品久久久久久久| 777久久人妻少妇嫩草av网站| 99久久精品热视频| 麻豆av在线久日| 日日爽夜夜爽网站| 婷婷精品国产亚洲av在线| 男女之事视频高清在线观看| 日韩欧美在线二视频| 一级黄色大片毛片| 99精品欧美一区二区三区四区| 国产欧美日韩一区二区精品| 国产精品亚洲美女久久久| 久久天堂一区二区三区四区| 欧美日韩一级在线毛片| 小说图片视频综合网站| 国产高清有码在线观看视频 | av国产免费在线观看| 99精品久久久久人妻精品| 黄色成人免费大全| 精品国产乱码久久久久久男人| 日日摸夜夜添夜夜添小说| 一进一出好大好爽视频| 老司机午夜十八禁免费视频| 国产免费男女视频| 色在线成人网| 91国产中文字幕| 曰老女人黄片| 变态另类丝袜制服| 最近最新中文字幕大全电影3| 美女大奶头视频| 午夜老司机福利片| 999精品在线视频| 久久久久久久精品吃奶| 在线观看舔阴道视频| 亚洲真实伦在线观看| 久久午夜亚洲精品久久| 天堂av国产一区二区熟女人妻 | 麻豆成人午夜福利视频| 一级黄色大片毛片| 99国产精品一区二区蜜桃av| 一区二区三区国产精品乱码| 97人妻精品一区二区三区麻豆| 亚洲精品久久国产高清桃花| 欧美日韩国产亚洲二区| 欧美黑人欧美精品刺激| 精品午夜福利视频在线观看一区| 悠悠久久av| 99国产精品一区二区蜜桃av| 香蕉久久夜色| 亚洲精品国产精品久久久不卡| 久久精品国产99精品国产亚洲性色| 亚洲国产精品合色在线| 99国产综合亚洲精品| 黄色毛片三级朝国网站| 精品乱码久久久久久99久播| 亚洲第一电影网av| 精华霜和精华液先用哪个| 国产视频内射| 夜夜看夜夜爽夜夜摸| 天天添夜夜摸| 国产精品一区二区免费欧美| 波多野结衣高清无吗| 1024视频免费在线观看| 亚洲av成人精品一区久久| www.自偷自拍.com| 俺也久久电影网| 久久久久免费精品人妻一区二区| 黄片小视频在线播放| 首页视频小说图片口味搜索| netflix在线观看网站| 国产成人系列免费观看| 一本一本综合久久| 波多野结衣巨乳人妻| av国产免费在线观看| 久99久视频精品免费| 欧美乱色亚洲激情| 亚洲熟妇中文字幕五十中出| 国产单亲对白刺激| 老汉色∧v一级毛片| 91av网站免费观看| 国产激情欧美一区二区| 男女做爰动态图高潮gif福利片| 18禁黄网站禁片午夜丰满| 天天躁夜夜躁狠狠躁躁| av片东京热男人的天堂| 国产不卡一卡二| 成人18禁在线播放| 亚洲熟妇中文字幕五十中出| 日日爽夜夜爽网站| or卡值多少钱| www日本在线高清视频| 日韩av在线大香蕉| 一级黄色大片毛片| 久久久国产成人免费| 久久精品人妻少妇| 国产成年人精品一区二区| 精品国产乱子伦一区二区三区| 久久精品影院6| 国产av麻豆久久久久久久| 国产精品爽爽va在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 制服丝袜大香蕉在线| www日本在线高清视频| 亚洲精华国产精华精| 亚洲电影在线观看av| 真人一进一出gif抽搐免费| 日本免费一区二区三区高清不卡| 97超级碰碰碰精品色视频在线观看| 国产在线观看jvid| 正在播放国产对白刺激| 国产黄a三级三级三级人| 久久久久久久久中文| 精品国产乱码久久久久久男人| 色老头精品视频在线观看| 久久精品成人免费网站| 麻豆久久精品国产亚洲av| 热99re8久久精品国产| 日韩成人在线观看一区二区三区| 我的老师免费观看完整版| 欧美在线一区亚洲| 在线观看www视频免费| 国产精品美女特级片免费视频播放器 | 熟女电影av网| 免费在线观看黄色视频的| 午夜影院日韩av| 精品午夜福利视频在线观看一区| 日韩欧美一区二区三区在线观看| 50天的宝宝边吃奶边哭怎么回事| 中文字幕熟女人妻在线| 久久久久国产精品人妻aⅴ院| 这个男人来自地球电影免费观看| 久久久久久九九精品二区国产 | 他把我摸到了高潮在线观看| 国产激情欧美一区二区| 一个人免费在线观看电影 | 一进一出好大好爽视频| 日日爽夜夜爽网站| 国产黄色小视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 午夜精品在线福利| av中文乱码字幕在线| 亚洲精品国产一区二区精华液| 国产精品,欧美在线| 国产1区2区3区精品| 精品久久久久久久毛片微露脸| 757午夜福利合集在线观看| 高清毛片免费观看视频网站| 日韩欧美国产一区二区入口| 色av中文字幕| 欧美黑人欧美精品刺激| 国产午夜精品论理片| 一级片免费观看大全| 欧美日韩黄片免| 欧美色欧美亚洲另类二区| 国内揄拍国产精品人妻在线| 免费在线观看日本一区| 国产黄色小视频在线观看| 黑人巨大精品欧美一区二区mp4| 一级作爱视频免费观看| 国产成人系列免费观看| 欧美久久黑人一区二区| 精品久久久久久久毛片微露脸| 欧美不卡视频在线免费观看 | 色在线成人网| 91大片在线观看| 久久久久久国产a免费观看| 免费在线观看视频国产中文字幕亚洲| 18美女黄网站色大片免费观看| 久久久久久人人人人人| 可以免费在线观看a视频的电影网站| 欧美午夜高清在线| 色综合婷婷激情| 欧美在线一区亚洲| 国产高清视频在线播放一区| 91麻豆av在线| 免费高清视频大片| 亚洲va日本ⅴa欧美va伊人久久| 嫁个100分男人电影在线观看| 伊人久久大香线蕉亚洲五| 成人亚洲精品av一区二区| 波多野结衣高清作品| 香蕉久久夜色| 午夜免费激情av| 人人妻,人人澡人人爽秒播| 日韩欧美在线二视频| 国产黄a三级三级三级人| 国产精品乱码一区二三区的特点| 国产亚洲av嫩草精品影院| 男女午夜视频在线观看| 国产aⅴ精品一区二区三区波| a级毛片a级免费在线| 成人国产综合亚洲| 欧美色视频一区免费| 国产精品久久久av美女十八| 夜夜躁狠狠躁天天躁| 黄频高清免费视频| 精品不卡国产一区二区三区| 很黄的视频免费| 丰满的人妻完整版| netflix在线观看网站| 18美女黄网站色大片免费观看| 又粗又爽又猛毛片免费看| 国产精品久久电影中文字幕| 亚洲中文av在线| 老司机在亚洲福利影院| 亚洲人与动物交配视频| 国产高清激情床上av| 天堂av国产一区二区熟女人妻 | 亚洲专区字幕在线| 成人手机av| 美女大奶头视频| 午夜精品一区二区三区免费看| 极品教师在线免费播放| 亚洲国产高清在线一区二区三| 国内精品久久久久久久电影| 亚洲激情在线av| 国产亚洲精品av在线| 高清毛片免费观看视频网站| 无遮挡黄片免费观看| xxx96com| а√天堂www在线а√下载| 人妻夜夜爽99麻豆av| 老司机福利观看| 老汉色av国产亚洲站长工具| 757午夜福利合集在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| av欧美777| 男男h啪啪无遮挡| 一级片免费观看大全| 国产三级中文精品| 色尼玛亚洲综合影院| 国产一区二区三区在线臀色熟女| 亚洲人成网站在线播放欧美日韩| 欧美3d第一页| 国产欧美日韩精品亚洲av| 国产精品国产高清国产av| 一级作爱视频免费观看| 婷婷六月久久综合丁香| 国产成人精品久久二区二区免费| 丝袜人妻中文字幕| 在线视频色国产色| 99精品在免费线老司机午夜| 免费搜索国产男女视频| 午夜福利在线在线| 日本一本二区三区精品| 两个人视频免费观看高清| 日韩免费av在线播放| 丰满人妻熟妇乱又伦精品不卡| 成在线人永久免费视频| 妹子高潮喷水视频| 婷婷丁香在线五月| 欧美一区二区国产精品久久精品 | 岛国在线免费视频观看| 男人舔女人的私密视频| 中文字幕最新亚洲高清| 一个人免费在线观看电影 | 久久久久久国产a免费观看| av福利片在线观看| e午夜精品久久久久久久| 天堂动漫精品| 欧美日韩精品网址| 两人在一起打扑克的视频| 深夜精品福利| 成年版毛片免费区| 中文字幕精品亚洲无线码一区| 他把我摸到了高潮在线观看| 无限看片的www在线观看| 成人av在线播放网站| 国产一区二区在线观看日韩 | av片东京热男人的天堂| 中文字幕人妻丝袜一区二区| 老熟妇乱子伦视频在线观看| 在线看三级毛片| 亚洲精品美女久久久久99蜜臀| 成年版毛片免费区| 桃色一区二区三区在线观看| 国产av不卡久久| 日本成人三级电影网站| 国产69精品久久久久777片 | 国产麻豆成人av免费视频| 免费电影在线观看免费观看| 岛国在线观看网站| aaaaa片日本免费| 岛国视频午夜一区免费看| 婷婷亚洲欧美| xxxwww97欧美| 99在线人妻在线中文字幕| 国产精品98久久久久久宅男小说| 最新在线观看一区二区三区| 亚洲av成人av| 久久久久国产精品人妻aⅴ院| 国内精品久久久久久久电影|