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

    基于組合式三系統(tǒng)估計(jì)量的人口普查凈誤差估計(jì)

    2022-08-16 09:42:38胡桂華劉譽(yù)環(huán)
    統(tǒng)計(jì)與信息論壇 2022年8期
    關(guān)鍵詞:估計(jì)量人口數(shù)總體

    胡桂華,文 婷,劉譽(yù)環(huán)

    (重慶工商大學(xué) a.數(shù)學(xué)與統(tǒng)計(jì)學(xué)院;b.經(jīng)濟(jì)社會(huì)應(yīng)用統(tǒng)計(jì)重慶市重點(diǎn)實(shí)驗(yàn)室,重慶 400067)

    一、引 言

    在人口普查質(zhì)量評(píng)估領(lǐng)域,研究?jī)粽`差具有重要意義。凈誤差率是衡量一個(gè)國(guó)家或地區(qū)人口普查登記質(zhì)量的核心指標(biāo)。凈誤差率不足2%或超過5%表明普查登記質(zhì)量較高或登記質(zhì)量較低。凈誤差率也是修正普查登記人口數(shù)的主要依據(jù)。有些國(guó)家直接依據(jù)估計(jì)的凈誤差率修正普查登記人口數(shù),也有些國(guó)家基于凈誤差率估計(jì)值和其他數(shù)據(jù)來源對(duì)普查登記人口數(shù)進(jìn)行修正。因此,研究?jī)粽`差具有重要意義。

    凈誤差率計(jì)算的關(guān)鍵是估計(jì)總體實(shí)際人口數(shù)。目前大多數(shù)國(guó)家使用雙系統(tǒng)估計(jì)量來估計(jì)它。但是,雙系統(tǒng)估計(jì)量有一個(gè)內(nèi)含交互作用偏差的“軟肋”。研究表明,依據(jù)三次捕獲模型構(gòu)造的三系統(tǒng)估計(jì)量能夠較好地解決雙系統(tǒng)估計(jì)量的交互作用偏差問題。與雙系統(tǒng)估計(jì)量所不同的是,三系統(tǒng)估計(jì)量并不要求三份人口名單必須獨(dú)立,從而避免了人口名單之間不獨(dú)立引起的交互作用偏差。

    三系統(tǒng)估計(jì)量為不完整三維列聯(lián)表7個(gè)已知單元的人口數(shù)和1個(gè)未知單元(未登記在任何一份名單的人口構(gòu)成的子總體)的人口數(shù)估計(jì)量(稱為缺失單元估計(jì)量)之和。缺失單元估計(jì)量依據(jù)7個(gè)已知單元的人口數(shù)來估計(jì),其具體形式依據(jù)三份人口名單的統(tǒng)計(jì)關(guān)系確定。把這種模式下構(gòu)造的三系統(tǒng)估計(jì)量稱為普通三系統(tǒng)估計(jì)量,共有8種普通三系統(tǒng)估計(jì)量。

    普通三系統(tǒng)估計(jì)量有兩個(gè)局限性。一是在給出每種統(tǒng)計(jì)關(guān)系的普通三系統(tǒng)估計(jì)量的情況下,使用皮爾遜卡方統(tǒng)計(jì)量或?qū)?shù)似然比統(tǒng)計(jì)量,選擇最適合于既定樣本數(shù)據(jù)描述的普通三系統(tǒng)估計(jì)量。這兩個(gè)檢驗(yàn)統(tǒng)計(jì)量應(yīng)用的前提條件是,構(gòu)造每個(gè)已知單元的期望頻數(shù)估計(jì)量。這是一項(xiàng)困難的工作。二是普通三系統(tǒng)估計(jì)量要么未利用7個(gè)已知單元信息,要么不加區(qū)分地對(duì)待每個(gè)信息。抽樣推斷理論指出,在構(gòu)造總體參數(shù)估計(jì)量時(shí),既要盡可能利用已知信息,還要根據(jù)信息的重要性確定其使用的次數(shù)。在總體中,同時(shí)登記在三份人口名單的人口所占比例最大,是最重要的信息,其在估計(jì)量中使用的次數(shù)應(yīng)該最多;同時(shí)登記在兩份人口名單的人口所占比例相對(duì)較少,是次重要的信息,其在估計(jì)量中使用的次數(shù)應(yīng)相對(duì)少一些;而只登記在一份人口名單的人口更少,屬于不重要信息,其在估計(jì)量中使用的次數(shù)應(yīng)為最少。

    為避免上述缺陷,提出組合式三系統(tǒng)估計(jì)量。通過研究,為政府統(tǒng)計(jì)部門使用三系統(tǒng)估計(jì)量提供更多選擇。組合式三系統(tǒng)估計(jì)量的學(xué)術(shù)價(jià)值體現(xiàn)在兩個(gè)方面。第一,它為7個(gè)已知單元的人口數(shù)與1個(gè)缺失單元的人口數(shù)估計(jì)量之和。缺失單元估計(jì)量的分子為三份人口名單同一類統(tǒng)計(jì)關(guān)系下的3個(gè)普通三系統(tǒng)估計(jì)量的分子之和,分母為它們的分母之和。在計(jì)算總體的組合式三系統(tǒng)估計(jì)量的抽樣方差時(shí),除了采用分層刀切法計(jì)算每一等概率人口層的抽樣方差外,還要計(jì)算等概率人口層之間的協(xié)方差,否則低估或高估總體的組合式三系統(tǒng)估計(jì)量的抽樣方差。第二,它須在等概率人口層建立,否則產(chǎn)生異質(zhì)性偏差而影響其估計(jì)精度。

    組合式三系統(tǒng)估計(jì)量的應(yīng)用價(jià)值體現(xiàn)在兩個(gè)方面。一方面,它省去了使用對(duì)數(shù)似然比選擇最優(yōu)組合式三系統(tǒng)估計(jì)量這一繁瑣環(huán)節(jié),便于在實(shí)際部門推廣應(yīng)用。對(duì)政府統(tǒng)計(jì)部門的實(shí)際工作人員來說,他們希望得到的是便于操作應(yīng)用的、理論可靠的參數(shù)估計(jì)量。另一方面,它可以拓展到其他領(lǐng)域,解決其相關(guān)問題。例如,估計(jì)經(jīng)濟(jì)普查、農(nóng)業(yè)普查的凈誤差率、戶籍登記系統(tǒng)的覆蓋率、老年人口數(shù)目估計(jì)。

    創(chuàng)新之處體現(xiàn)在三個(gè)方面。第一,提出兩種抽樣登記的、人口移動(dòng)的組合式三系統(tǒng)估計(jì)量,并在后面的實(shí)證分析中進(jìn)行了抽樣估計(jì)精度比較。第二,將抽樣調(diào)查領(lǐng)域的分層二重抽樣方法移植到人口普查凈誤差估計(jì)領(lǐng)域,解決其事后計(jì)數(shù)調(diào)查樣本抽取及組合式三系統(tǒng)估計(jì)量各個(gè)單元估計(jì)量的構(gòu)造問題。第三,厘清了普通三系統(tǒng)估計(jì)量與組合式三系統(tǒng)估計(jì)量的內(nèi)在邏輯關(guān)系,這有助于讀者理解和應(yīng)用組合式三系統(tǒng)估計(jì)量。

    二、文獻(xiàn)綜述

    在捕獲—再捕獲模型基礎(chǔ)上,美國(guó)學(xué)者增加了一次捕獲,利用三次捕獲的結(jié)果構(gòu)造三次捕獲模型[1-5];將三次捕獲模型拓展到人口普查質(zhì)量評(píng)估領(lǐng)域構(gòu)造三份人口名單兩種統(tǒng)計(jì)關(guān)系下的全面登記的普通三系統(tǒng)估計(jì)量[6]。孟杰等研究了普通三系統(tǒng)估計(jì)量及其統(tǒng)計(jì)性質(zhì)[7]。胡桂華等研究了8種統(tǒng)計(jì)關(guān)系下的抽樣登記的、人口移動(dòng)的普通三系統(tǒng)估計(jì)量,為其應(yīng)用于人口普查凈誤差估計(jì)奠定了理論基礎(chǔ)[8-12]。

    美國(guó)普查局2010年人口普查質(zhì)量評(píng)估后,曾經(jīng)打算在未來使用普通三系統(tǒng)估計(jì)量。美國(guó)普查局專家Griffin撰寫了一篇這方面的研究論文,分析了雙系統(tǒng)估計(jì)量低估總體實(shí)際人口數(shù)的主要原因在于,未包括同時(shí)傾向于遺漏普查與事后計(jì)數(shù)調(diào)查的人口,并且使用不完整三維列聯(lián)表、最大似然估計(jì)和對(duì)數(shù)線性模型建立了三份人口名單各種統(tǒng)計(jì)關(guān)系下的普通三系統(tǒng)估計(jì)量[13]。后來由于未建立起覆蓋人口總體范圍廣的行政記錄系統(tǒng),以及未進(jìn)行大規(guī)模的試點(diǎn)研究,于是美國(guó)在2020年仍然使用雙系統(tǒng)估計(jì)量。

    中國(guó)原計(jì)劃在2020年使用普通三系統(tǒng)估計(jì)量,后來由于三個(gè)原因放棄。一是普查名單和事后計(jì)數(shù)調(diào)查名單與戶籍人口名單高度相似,即外表看起來是三份名單,其實(shí)是兩份人口名單,從而失去了構(gòu)建三系統(tǒng)估計(jì)量的意義[14]。二是尚未建立起科學(xué)完整的普查名單、事后計(jì)數(shù)調(diào)查名單與戶籍人口名單的比對(duì)方案。在2020年,中國(guó)仍然使用在等概率人口層建立的雙系統(tǒng)估計(jì)量[15]。三是中國(guó)并未建立起適合于普通三系統(tǒng)估計(jì)量所需要的復(fù)合人口行政記錄名單。由于每份人口行政記錄有其特定的用途和范圍,因此建立中國(guó)的復(fù)合人口行政記錄名單,不能只包括戶籍名單,除該名單外,還要把常住人口名單,出生和死亡醫(yī)學(xué)證明納入其中,尤其是要把這些名單中的人具體落實(shí)到每一個(gè)樣本普查小區(qū)。這項(xiàng)工作最好在普查標(biāo)準(zhǔn)時(shí)點(diǎn)前一個(gè)星期完成。復(fù)合人口行政記錄名單最理想的截止時(shí)間是普查標(biāo)準(zhǔn)時(shí)點(diǎn)。這在實(shí)際中難以做到。丹麥、挪威和冰島等北歐國(guó)家的行政記錄截止日在普查時(shí)點(diǎn)前一個(gè)星期。

    研究組合式三系統(tǒng)估計(jì)量的意義在于,將利用信息量不足的若干普通三系統(tǒng)估計(jì)量組合在一起,構(gòu)造出一個(gè)新的三系統(tǒng)估計(jì)量,從而間接利用了多個(gè)信息,避免了普通三系統(tǒng)估計(jì)量的缺陷。建立三系統(tǒng)估計(jì)量的初衷就是要盡可能多地使用多個(gè)信息來源。利用更多相關(guān)信息,有助于提高估計(jì)精度。從這個(gè)角度看,組合式三系統(tǒng)估計(jì)量具有研究?jī)r(jià)值,有助于提高凈誤差的抽樣估計(jì)精度。

    三、組合式三系統(tǒng)估計(jì)量及其相關(guān)理論

    (一)全面登記的雙系統(tǒng)估計(jì)量及普通三系統(tǒng)估計(jì)量

    我們首先依據(jù)三次捕獲模型要求,構(gòu)造事后計(jì)數(shù)調(diào)查對(duì)總體全面登記的普通三系統(tǒng)估計(jì)量。這么做的理由是,普通三系統(tǒng)估計(jì)量來源于三次捕獲模型。三次捕獲模型是在三次全面捕獲的情況下構(gòu)造的。

    普通三系統(tǒng)估計(jì)量由雙系統(tǒng)估計(jì)量拓展而來。為便于推導(dǎo)普通三系統(tǒng)估計(jì)量公式,簡(jiǎn)要介紹一下雙系統(tǒng)估計(jì)量是有必要的。構(gòu)造雙系統(tǒng)估計(jì)量的重要工具是不完整二維列聯(lián)表。在假設(shè)總體中的每一個(gè)人有同樣的概率登記在普查的情況下,總體中的人口在普查與事后計(jì)數(shù)調(diào)查中的登記情況列示在該表中,見表1。

    表1 不完整二維列聯(lián)表

    在表1中,右下腳標(biāo)1表示登記在普查或事后計(jì)數(shù)調(diào)查中,0表示未登記在普查或事后計(jì)數(shù)調(diào)查中。在假設(shè)總體人口有同樣的概率登記在普查或事后計(jì)數(shù)調(diào)查,以及事后計(jì)數(shù)調(diào)查對(duì)總體全面登記情況下,估計(jì)總體規(guī)模的雙系統(tǒng)估計(jì)量(Dual System Estimator,DSE)為:

    (1)

    從式(1)可以看出,x11,x10,x01是已知單元人口數(shù),而x00是未知單元人口數(shù),需要估計(jì)。未知單元的人口數(shù)可以通過式(2)來估計(jì):

    (2)

    為便于普通三系統(tǒng)估計(jì)量的建立,需要將總體中的人口在三份人口名單的登記情況列示在不完整三維列聯(lián)表,見表2。

    表2 不完整三維列聯(lián)表

    在表2中,xijk中的右下腳標(biāo)i,j,k為示性變量,均取值1或0,1表示在普查登記中,0表示不在普查登記中。除最后一個(gè)單元(ijk)=(000)的x000需要估計(jì)外,其他單元的人口數(shù)都是可以觀察到的。

    如果不存在比對(duì)誤差,兩項(xiàng)調(diào)查和人口行政記錄均是對(duì)同一總體的全面登記,總體人口有同樣概率登記在普查的情況下,那么估計(jì)總體規(guī)模的普通三系統(tǒng)估計(jì)量(Triple System Estimator,TSE)為:

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    (9)

    (二)全面登記的組合式三系統(tǒng)估計(jì)量

    全面登記的組合式三系統(tǒng)估計(jì)量為7個(gè)已知單元的人口數(shù)與組合式缺失單元估計(jì)量之和。

    1.全面登記的、無人口移動(dòng)的組合式三系統(tǒng)估計(jì)量

    不難看出,式(4)至式(6)屬于同類統(tǒng)計(jì)關(guān)系,而式(7)至式(9)屬于另一類統(tǒng)計(jì)關(guān)系。式(4)至式(6)存在的問題是,對(duì)7個(gè)已知單元的信息給予同樣的關(guān)注,賦予同樣的權(quán)數(shù)1。式(7)至式(9)只是使用了7個(gè)已知單元信息中的3個(gè),使用的信息量過少,必然影響估計(jì)精度。針對(duì)這種情況,可以構(gòu)造兩種組合式缺失單元估計(jì)量。一是將式(4)至式(6)的分子和分母分別相加,得到組合式三系統(tǒng)估計(jì)量式(10);二是將式(7)至式(9)的分子和分母分別相加,得到組合式缺失單元估計(jì)量式(11)。

    (10)

    (11)

    將式(10)和式(11)分別代入式(3),得到兩種組合式三系統(tǒng)估計(jì)量,見式(12)和式(15)。

    (12)

    (13)

    (14)

    (15)

    (16)

    (17)

    式(1)至式(17)未考慮普查標(biāo)準(zhǔn)時(shí)點(diǎn)與事后計(jì)數(shù)調(diào)查標(biāo)準(zhǔn)時(shí)點(diǎn)之間的人口移動(dòng)(不包括出生和死亡者)。但實(shí)際上,這兩個(gè)標(biāo)準(zhǔn)時(shí)點(diǎn)之間不可避免有人口移動(dòng)。有三種移動(dòng)人口:一是無移動(dòng)人口(Non-movers,縮寫為n),即普查標(biāo)準(zhǔn)時(shí)點(diǎn)和事后計(jì)數(shù)調(diào)查標(biāo)準(zhǔn)時(shí)點(diǎn)之間一直居住在本普查小區(qū)的人口。二是向內(nèi)移動(dòng)人口(In-movers,縮寫為i),即普查標(biāo)準(zhǔn)時(shí)點(diǎn)居住在其他普查小區(qū),在事后計(jì)數(shù)調(diào)查標(biāo)準(zhǔn)時(shí)點(diǎn)前來到本普查小區(qū)的人口。三是向外移動(dòng)人口(Out-movers,縮寫為o),即普查標(biāo)準(zhǔn)時(shí)點(diǎn)在本普查小區(qū),事后計(jì)數(shù)調(diào)查標(biāo)準(zhǔn)時(shí)點(diǎn)前離開本普查小區(qū)的人口。事后計(jì)數(shù)調(diào)查人口構(gòu)成有兩種方法:一是A構(gòu)成方法,由無移動(dòng)人口和向外移動(dòng)人口構(gòu)成。二是B構(gòu)成方法,由無移動(dòng)人口和向內(nèi)移動(dòng)人口構(gòu)成。

    對(duì)式(12)至式(17),均構(gòu)造兩種方法下的人口移動(dòng)的組合式三系統(tǒng)估計(jì)量。此時(shí),式(12)至式(17)中每個(gè)單元涉及事后計(jì)數(shù)調(diào)查的項(xiàng)一律變?yōu)閮身?xiàng)??紤]到向內(nèi)移動(dòng)人口不可能出現(xiàn)在本普查小區(qū)的行政記錄人口名單中,因?yàn)槠淦詹闃?biāo)準(zhǔn)時(shí)點(diǎn)在其他小區(qū),而不在本小區(qū),因此向內(nèi)移動(dòng)人口涉及到的人口行政記錄的項(xiàng)一律為零,例如x111i=0。另外,向外移動(dòng)人口不可能登記在本小區(qū)的事后計(jì)數(shù)調(diào)查人口名單中,因?yàn)槭潞笥?jì)數(shù)調(diào)查標(biāo)準(zhǔn)時(shí)點(diǎn)時(shí)不在本普查小區(qū),因此向外移動(dòng)人口涉及事后計(jì)數(shù)調(diào)查的項(xiàng)一律為零,例如x011o=0。

    2.式(12)至式(14)全面登記的、人口移動(dòng)的組合式三系統(tǒng)估計(jì)量

    如果采取A構(gòu)成法,那么式(12)至式(14)變?yōu)?

    (18)

    x1,A=x111n+x111o+x110n+x110o+x101n+x101o+x100n+x100o+x011n+x011o+x010n+x010o+x001n+x001o

    (19)

    由于x111o=x110o=x011o=x010o=0,因此式(19)變?yōu)槭?20):

    x1,A=x111n+x110n+x101n+x101o+x100n+x100o+x011n+x010n+x001n+x001o

    (20)

    (21)

    如果使用B構(gòu)成法,那么式(12)~(14)變?yōu)?

    (22)

    x1,B=x111n+x111i+x110n+x110i+x101n+x101i+x100n+x100i+x011n+x011i+x010n+x010i+x001n+x001i

    (23)

    由于x111i=x101i=x001i=x011i=0,因此式(23)變?yōu)?

    x1,B=x111n+x110n+x101n+x100n+x011n+x010n+x001n+x010i+x100i+x110i

    (24)

    (25)

    3.式(15)至式(17)全面登記的、人口移動(dòng)的組合式三系統(tǒng)估計(jì)量

    如果采取A構(gòu)成法,那么式(15)至式(17)變?yōu)?

    (26)

    x2,A=x111n+x111o+x110n+x110o+x101n+x101o+x100n+x100o+x011n+x011o+x010n+x010o+x001n+x001o

    (27)

    由于x111o=x110o=x011o=x010o=0,因此式(27)變?yōu)槭?28):

    x2,A=x111n+x110n+x101n+x101o+x100n+x100o+x011n+x010n+x001n+x001o

    (28)

    (29)

    如果使用B構(gòu)成法,那么式(15)至式(17)變?yōu)?

    (30)

    x2,B=x111n+x111i+x110n+x110i+x101n+x101i+x100n+x100i+x011n+x011i+x010n+x010i+x001n+x001i

    (31)

    由于x111i=x101i=x001i=x011i=0,因此式(31)變?yōu)?

    x2,B=x111n+x110n+x101n+x100n+x011n+x010n+x001n+x010i+x100i+x110i

    (32)

    (33)

    (三)抽樣登記的組合式三系統(tǒng)估計(jì)量

    事后計(jì)數(shù)調(diào)查是對(duì)樣本普查小區(qū)人口的抽樣登記[16-18]。這決定了只能建立抽樣登記的組合式三系統(tǒng)估計(jì)量。這項(xiàng)工作分為兩個(gè)步驟,第一步是在全面登記的組合式三系統(tǒng)估計(jì)量基礎(chǔ)上,寫出抽樣登記的組合式三系統(tǒng)估計(jì)量。第二步是采用有限總體概率樣本數(shù)據(jù),構(gòu)造抽樣登記的組合式三系統(tǒng)估計(jì)量每個(gè)單元人口數(shù)的估計(jì)量。

    1.式(12)至式(14)抽樣登記的、人口移動(dòng)的組合式三系統(tǒng)估計(jì)量

    如果采取A構(gòu)成法,那么:

    (34)

    (35)

    (36)

    如果使用B構(gòu)成法,那么:

    (37)

    (38)

    (39)

    2.式(15)至式(17)抽樣登記的、人口移動(dòng)的組合式三系統(tǒng)估計(jì)量

    如果使用A構(gòu)成法,那么:

    (40)

    (41)

    (42)

    如果使用B構(gòu)成法,那么:

    (43)

    (44)

    (45)

    (46)

    在式(46)中,yhgi表示層hg的最終進(jìn)入第二重樣本的普查小區(qū)i在三份人口名單登記的人口數(shù),例如同時(shí)在三份人口名單登記的無移動(dòng)人口數(shù)x111n,在普查中登記但未在事后計(jì)數(shù)調(diào)查和人口行政記錄登記的向外移動(dòng)人口數(shù)x100o;H為第一重抽樣前對(duì)總體普查小區(qū)劃分的總層數(shù),h為任意層;G是對(duì)第一重樣本普查小區(qū)進(jìn)一步劃分的總層數(shù),g是任意層;xhgi為層h的樣本普查小區(qū)是否在層g的示性函數(shù),如果在,xhgi=1,否則xhgi=0;Ihgi為層hg的第一重樣本普查小區(qū)是否進(jìn)入第二重樣本的示性函數(shù),如果進(jìn)入了,Ihgi=1,否則Ihgi=0;αhgi為第i樣本普查小區(qū)經(jīng)過二重抽樣后進(jìn)入到第二重樣本的最終抽樣權(quán)數(shù)。如果樣本答復(fù)率不足100%,則要根據(jù)樣本答復(fù)率調(diào)整初始抽樣權(quán)數(shù)。調(diào)整的辦法是,將初始抽樣權(quán)數(shù)除以樣本答復(fù)率。在后面的實(shí)際案例中,樣本答復(fù)率為100%。

    如果層h的普查小區(qū)總數(shù)為Nh,樣本普查小區(qū)數(shù)為nh,層hg的普查小區(qū)總數(shù)為Mhg,樣本普查小區(qū)數(shù)為mhg,并且每重抽樣采取簡(jiǎn)單隨機(jī)抽樣,每重抽樣的抽樣單位均為普查小區(qū),那么經(jīng)過兩重抽樣后,第i樣本普查小區(qū)的最終抽樣權(quán)數(shù)αhgi=(Nh/nh)(Mhg/mhg)。

    注意式(34)至式(45),需要在等概率人口層v建立。這里只是為了書寫便利,才在式(34)至式(45)中省去了v記號(hào)。在建立了等概率人口層v的組合式三系統(tǒng)估計(jì)量之后,還要建立總體的抽樣登記的、人口移動(dòng)的組合式三系統(tǒng)估計(jì)量,用P表示總體:

    (47)

    (48)

    (49)

    (50)

    (四)抽樣登記的、人口移動(dòng)的組合式三系統(tǒng)估計(jì)量的抽樣方差估計(jì)

    式(34)、式(37)、式(40)、式(43),以及式(47)至式(50)是較為復(fù)雜的估計(jì)量,其抽樣方差使用刀切法近似計(jì)算。刀切法實(shí)施的關(guān)鍵是復(fù)制權(quán)數(shù)的計(jì)算[19-20]。復(fù)制權(quán)數(shù)是指每剔除第一重樣本中的一個(gè)樣本普查小區(qū)后重新計(jì)算的進(jìn)入第二重樣本普查小區(qū)的抽樣權(quán)數(shù)。

    (51)

    (52)

    (53)

    (54)

    (55)

    (56)

    (57)

    (58)

    (59)

    (60)

    (61)

    等概率人口層抽樣登記的三系統(tǒng)估計(jì)量的刀切抽樣方差估計(jì)量為:

    (62)

    (63)

    (64)

    (65)

    總體的抽樣登記的三系統(tǒng)估計(jì)量的刀切抽樣方差估計(jì)量為:

    (66)

    (67)

    (68)

    (69)

    (五)總體的人口普查凈誤差估計(jì)

    (70)

    (71)

    (72)

    (73)

    在式(70)~(73)中,Cp是總體的、已知的最終普查登記人口數(shù),其抽樣方差為零。從式(70)~(73)可以看出,只要估計(jì)了總體的實(shí)際人口數(shù),將其減去普查登記人口數(shù),便可以很方便地估計(jì)人口普查凈誤差,并且兩者的抽樣方差相等。

    四、實(shí)證分析

    (一)基本信息及樣本數(shù)據(jù)

    為計(jì)算人口普查凈誤差,在實(shí)證分析中,應(yīng)該獲得實(shí)證范圍內(nèi)普查年樣本普查小區(qū)人口的微觀資料、各個(gè)等概率人口層的普查登記人口數(shù)。然而,我們沒有得到這樣的數(shù)據(jù),只是獲得了非普查年某行政區(qū)的相關(guān)資料。根據(jù)這些相關(guān)資料,可以估計(jì)總體實(shí)際人口數(shù),而總體實(shí)際人口數(shù)的估計(jì)是凈誤差計(jì)算的關(guān)鍵。從這個(gè)意義上來看,在凈誤差估計(jì)中,非普查年的微觀資料還是有使用價(jià)值的。

    實(shí)證資料來源于某行政區(qū)的20個(gè)樣本普查小區(qū)的三份人口名單。時(shí)間為2017年8月31日。研究采取分層二重抽樣方法,在第一重抽樣前,將其所有普查小區(qū)分在兩層:街道層共有1 970個(gè)普查小區(qū),記為h=1和N1=1 970;鎮(zhèn)層共有1 068個(gè)普查小區(qū),記為h=2和N2=1 068。這兩層總共有普查小區(qū)3 038個(gè),記為N=3 038。對(duì)h層樣本普查小區(qū),進(jìn)一步按照調(diào)查難度分為三層:調(diào)查難度大層,記為g=1;調(diào)查難度中等層,記為g=2;調(diào)查難度小層,記為g=3。在每重抽樣中,以普查小區(qū)為抽樣單位,簡(jiǎn)單隨機(jī)不重復(fù)抽取樣本。樣本形成過程見表3,樣本數(shù)據(jù)見表4至表7。

    總體人口劃分在四個(gè)等概率人口層:0~64歲男性;0~64歲女性;65歲及以上男性;65歲及以上女性。

    (二)基于組合式三系統(tǒng)估計(jì)量的估計(jì)結(jié)果

    利用表3至表7樣本數(shù)據(jù),以及使用式(34)至式(69),得到等概率人口層及總體實(shí)際人口數(shù)的估計(jì)值及其抽樣標(biāo)準(zhǔn)誤差估計(jì)值,見表8。從表8可以得到如下重要信息。

    表3 樣本形成、樣本量和抽樣權(quán)數(shù)

    第一,表8提供了4個(gè)等概率人口層及總體的實(shí)際人口數(shù)及其抽樣標(biāo)準(zhǔn)誤差估計(jì)值。各個(gè)等概率人口層的實(shí)際人口數(shù)使用組合式三系統(tǒng)估計(jì)量估計(jì)。組合式三系統(tǒng)估計(jì)量不能直接用來估計(jì)總體的實(shí)際人口數(shù)??傮w的抽樣標(biāo)準(zhǔn)誤差不是四個(gè)等概率人口層的抽樣標(biāo)準(zhǔn)誤差之和,而是四個(gè)等概率人口層的抽樣方差之和及它們之間協(xié)方差2倍之和總和的平方根。

    第二,使用兩種組合式三系統(tǒng)估計(jì)量估計(jì)的總體實(shí)際人口數(shù)在77萬人左右。按照中國(guó)平均每個(gè)普查小區(qū)250人常住人口來算,本行政區(qū)的3 038個(gè)普查小區(qū)的人口數(shù)應(yīng)為759 500人。可見,兩種組合式三系統(tǒng)估計(jì)量估計(jì)的結(jié)果與現(xiàn)實(shí)情況基本吻合,是有效的三系統(tǒng)估計(jì)量。如果把759 500人視作普查登記人口數(shù),那么基于A構(gòu)成法的第一個(gè)估計(jì)量估計(jì)的實(shí)際人數(shù)為769 339人,凈誤差為9 839人,凈誤差率為1.27%,這也屬于合理估計(jì)值。

    第三,對(duì)同一種組合式三系統(tǒng)估計(jì)量,使用A構(gòu)成法和B構(gòu)成法得到的估計(jì)結(jié)果不盡相同。如果采取第一種組合式三系統(tǒng)估計(jì)量,兩種構(gòu)成法提供的人口數(shù)估計(jì)值分別為769 339人和773 742人。采取第二種組合式三系統(tǒng)估計(jì)量,這兩種構(gòu)成法提供的人口數(shù)估計(jì)值分別為770 212人和775 640人。A構(gòu)成法直接提供的是普查標(biāo)準(zhǔn)時(shí)點(diǎn)的人口數(shù)估計(jì)值,與人口普查凈誤差對(duì)普查目標(biāo)總體實(shí)際人口數(shù)估計(jì)值的要求一致,但找到向外移動(dòng)人口有比較大的困難。B構(gòu)成法雖然也是提供普查標(biāo)準(zhǔn)時(shí)點(diǎn)的人口數(shù)估計(jì)值,但它是通過向內(nèi)移動(dòng)人口回憶普查標(biāo)準(zhǔn)時(shí)點(diǎn)的情況得到的,可能發(fā)生回憶誤差。不過B構(gòu)成法也有其優(yōu)勢(shì),那就是向內(nèi)移動(dòng)人口在本樣本普查小區(qū),不存在尋找人的問題,調(diào)查方便。因此,兩種構(gòu)成法各有利弊,需要進(jìn)行選擇。為了更好實(shí)現(xiàn)凈誤差估計(jì)目標(biāo),應(yīng)該選取A構(gòu)成法。為便于調(diào)查,則選取B構(gòu)成法。在實(shí)際應(yīng)用中,選取A構(gòu)成法是合適的。

    表4 0~64歲男性層各單元的常住人口數(shù)單位:人

    表5 0~64歲女性層各單元的常住人口數(shù)單位:人

    表6 65歲及以上男性層各單元的常住人口數(shù)單位:人

    表7 65歲及以上女性層各單元的常住人口數(shù)單位:人

    表8 基于組合式三系統(tǒng)估計(jì)量的人口數(shù)及其抽樣標(biāo)準(zhǔn)誤差估計(jì)值 單位:人

    (三)組合式三系統(tǒng)估計(jì)量與普通三系統(tǒng)估計(jì)量的抽樣估計(jì)精度比較

    表9 組合式三系統(tǒng)估計(jì)量與普通三系統(tǒng)估計(jì)量的比較 單位:人

    從表9可以看出:

    第一,兩個(gè)組合式三系統(tǒng)估計(jì)量估計(jì)的總體實(shí)際人口數(shù)分別為769 339人和773 742人,而普通三系統(tǒng)估計(jì)量估計(jì)的總體實(shí)際人口數(shù)為766 859人。如果普查登記人口數(shù)為759 500人,那么兩個(gè)組合式三系統(tǒng)估計(jì)量和普通三系統(tǒng)估計(jì)量估計(jì)的人口普查凈誤差分別為9 839人、14 242人和7 359人。這說明,由于利用的信息量不足,普通三系統(tǒng)估計(jì)量低估總體的凈誤差。組合式三系統(tǒng)估計(jì)量能避免低估總體凈誤差。

    第二,從總體實(shí)際人口數(shù)估計(jì)值的抽樣標(biāo)準(zhǔn)誤差來看,兩個(gè)組合式三系統(tǒng)估計(jì)量的抽樣標(biāo)準(zhǔn)誤差分別為5 018人和5 301人,而普通三系統(tǒng)估計(jì)量的抽樣標(biāo)準(zhǔn)誤差為5 390人。這意味著,組合式三系統(tǒng)估計(jì)量比起普通三系統(tǒng)估計(jì)量擁有更高的估計(jì)精度。

    第三,從各個(gè)等概率人口層實(shí)際人口數(shù)估計(jì)值的抽樣標(biāo)準(zhǔn)誤差來看,除0~64歲男性層外,兩個(gè)組合式三系統(tǒng)估計(jì)量的其他所有等概率人口層的抽樣標(biāo)準(zhǔn)誤差均小于普通三系統(tǒng)估計(jì)量。例如,0~64歲女性層,兩個(gè)組合式三系統(tǒng)估計(jì)量的抽樣標(biāo)準(zhǔn)誤差分別為1 910人和1 630人,而普通三系統(tǒng)估計(jì)量實(shí)際人口數(shù)估計(jì)值的抽樣標(biāo)準(zhǔn)誤差為1 932人。這表明兩個(gè)組合式三系統(tǒng)估計(jì)量更適合于等概率人口層實(shí)際人口數(shù)的估計(jì)。

    五、結(jié) 論

    第一,組合式三系統(tǒng)估計(jì)量是在普通三系統(tǒng)估計(jì)量的基礎(chǔ)上建立起來的,是一個(gè)全新的三系統(tǒng)估計(jì)量。與普通三系統(tǒng)估計(jì)量不同的是,它是把三份人口名單同一類統(tǒng)計(jì)關(guān)系下的三個(gè)普通三系統(tǒng)估計(jì)量的缺失單元估計(jì)量的分子和分母分別相加得到的。與普通三系統(tǒng)估計(jì)量相同的是,它的7個(gè)已知單元的人口數(shù)與普通三系統(tǒng)估計(jì)量的相同。

    第二,組合式三系統(tǒng)估計(jì)量須在等概率人口層建立。為此,需要選擇變量對(duì)總體分層。在選擇分層變量時(shí),要考慮到變量值能否獲得。個(gè)人性格雖然是影響人們是否參與人口普查的重要變量,但很難知悉一個(gè)人的性格,因此不應(yīng)該把性格作為對(duì)總體人口分層的重要變量。

    第三,刀切法適合于組合式三系統(tǒng)估計(jì)量的抽樣方差計(jì)算。組合式三系統(tǒng)估計(jì)量利用的是樣本資料,需要計(jì)算其抽樣方差。由于估計(jì)量結(jié)構(gòu)復(fù)雜,無法用代數(shù)方法直接計(jì)算其方差。然而,不能因此就放棄了方差的計(jì)算。用分層刀切法近似計(jì)算這個(gè)方差是可行的。

    第四,組合式三系統(tǒng)估計(jì)量對(duì)實(shí)際人口數(shù)的估計(jì)精度高于普通三系統(tǒng)估計(jì)量,希望被政府統(tǒng)計(jì)部門采納及推廣應(yīng)用。

    猜你喜歡
    估計(jì)量人口數(shù)總體
    確定代表名額按戶籍人口數(shù)還是常住人口數(shù)
    人大研究(2022年3期)2022-04-13 00:47:04
    用樣本估計(jì)總體復(fù)習(xí)點(diǎn)撥
    2020年秋糧收購(gòu)總體進(jìn)度快于上年
    外匯市場(chǎng)運(yùn)行有望延續(xù)總體平穩(wěn)發(fā)展趨勢(shì)
    基于BP神經(jīng)網(wǎng)絡(luò)的吉林省GDP預(yù)測(cè)研究
    直擊高考中的用樣本估計(jì)總體
    淺談估計(jì)量的優(yōu)良性標(biāo)準(zhǔn)
    基于配網(wǎng)先驗(yàn)信息的諧波狀態(tài)估計(jì)量測(cè)點(diǎn)最優(yōu)配置
    負(fù)極值指標(biāo)估計(jì)量的漸近性質(zhì)
    使用變異系數(shù)和Kurtosis系數(shù)的雙輔助變量的比估計(jì)
    少妇裸体淫交视频免费看高清| 夜夜躁狠狠躁天天躁| 嫁个100分男人电影在线观看| 久久这里只有精品中国| 国产成人精品久久二区二区免费| 欧美丝袜亚洲另类 | 人人妻人人看人人澡| 国产69精品久久久久777片 | 中文字幕高清在线视频| 免费在线观看成人毛片| 黄频高清免费视频| 宅男免费午夜| 嫩草影视91久久| 特级一级黄色大片| 天堂影院成人在线观看| bbb黄色大片| 啪啪无遮挡十八禁网站| 美女高潮喷水抽搐中文字幕| 午夜福利在线在线| 国产精品98久久久久久宅男小说| 小蜜桃在线观看免费完整版高清| 久久久久精品国产欧美久久久| 麻豆成人午夜福利视频| 搡老妇女老女人老熟妇| 在线免费观看的www视频| 日本与韩国留学比较| 白带黄色成豆腐渣| 成人精品一区二区免费| 淫妇啪啪啪对白视频| 成人亚洲精品av一区二区| 热99re8久久精品国产| 免费人成视频x8x8入口观看| 国产一区二区在线av高清观看| 国产私拍福利视频在线观看| av国产免费在线观看| 欧美一区二区精品小视频在线| 手机成人av网站| 一级毛片女人18水好多| 精品久久久久久成人av| 成人性生交大片免费视频hd| 色精品久久人妻99蜜桃| 女生性感内裤真人,穿戴方法视频| 久久精品夜夜夜夜夜久久蜜豆| 国内精品美女久久久久久| 特级一级黄色大片| 日本一本二区三区精品| 老司机午夜十八禁免费视频| a级毛片在线看网站| 极品教师在线免费播放| 日本撒尿小便嘘嘘汇集6| 性色avwww在线观看| 精品久久久久久久毛片微露脸| 国产精品久久久人人做人人爽| 国产精品久久久久久亚洲av鲁大| 国产精品一区二区精品视频观看| 日韩大尺度精品在线看网址| www国产在线视频色| 久久九九热精品免费| 琪琪午夜伦伦电影理论片6080| 国产精品久久久久久精品电影| 亚洲av片天天在线观看| 亚洲av电影不卡..在线观看| 亚洲成av人片免费观看| 色噜噜av男人的天堂激情| 99久久精品一区二区三区| 精品久久蜜臀av无| 1024手机看黄色片| 国产精品久久视频播放| 亚洲无线在线观看| 欧美午夜高清在线| 最新中文字幕久久久久 | 国产免费av片在线观看野外av| 免费看a级黄色片| 久久久久久九九精品二区国产| 国产精品1区2区在线观看.| 搡老熟女国产l中国老女人| 国产探花在线观看一区二区| 亚洲欧美精品综合久久99| 午夜激情福利司机影院| 午夜福利视频1000在线观看| 嫩草影院精品99| 欧美黑人欧美精品刺激| 久久婷婷人人爽人人干人人爱| 亚洲精品一区av在线观看| 999精品在线视频| 日韩中文字幕欧美一区二区| 国产精品亚洲美女久久久| 怎么达到女性高潮| 国产av一区在线观看免费| 国产69精品久久久久777片 | 亚洲午夜精品一区,二区,三区| 国产欧美日韩精品一区二区| 欧美国产日韩亚洲一区| 欧美中文日本在线观看视频| 久99久视频精品免费| 亚洲精品中文字幕一二三四区| АⅤ资源中文在线天堂| 国产精品一区二区三区四区久久| 亚洲美女视频黄频| 国产野战对白在线观看| 亚洲成人久久性| 日本黄色视频三级网站网址| 亚洲专区字幕在线| 国产精品野战在线观看| 国产又黄又爽又无遮挡在线| 欧美丝袜亚洲另类 | 免费高清视频大片| 国产成人精品无人区| 国产精品亚洲av一区麻豆| 无人区码免费观看不卡| 日韩高清综合在线| 欧美激情在线99| 19禁男女啪啪无遮挡网站| 美女黄网站色视频| 97超视频在线观看视频| 国产欧美日韩一区二区三| 在线永久观看黄色视频| 午夜a级毛片| 天天躁狠狠躁夜夜躁狠狠躁| 午夜免费成人在线视频| 国产一区二区三区视频了| 亚洲欧美精品综合一区二区三区| 搞女人的毛片| 欧美绝顶高潮抽搐喷水| 成人高潮视频无遮挡免费网站| 韩国av一区二区三区四区| 99国产精品99久久久久| 欧美一级a爱片免费观看看| 中文字幕最新亚洲高清| 精品日产1卡2卡| 亚洲欧美激情综合另类| 在线观看日韩欧美| 成人av在线播放网站| 一本一本综合久久| 在线观看午夜福利视频| 男女视频在线观看网站免费| 亚洲九九香蕉| av在线天堂中文字幕| 国产精品女同一区二区软件 | 老司机深夜福利视频在线观看| 两个人看的免费小视频| 亚洲av电影不卡..在线观看| 婷婷丁香在线五月| 又爽又黄无遮挡网站| av片东京热男人的天堂| 热99re8久久精品国产| 久99久视频精品免费| 欧美色视频一区免费| 一本久久中文字幕| 精品久久久久久久末码| 成人无遮挡网站| 久久亚洲真实| www.www免费av| 母亲3免费完整高清在线观看| 亚洲黑人精品在线| 国产高清有码在线观看视频| 亚洲av第一区精品v没综合| 国内毛片毛片毛片毛片毛片| 久久久国产精品麻豆| 亚洲精品美女久久久久99蜜臀| 婷婷六月久久综合丁香| 深夜精品福利| 欧美黑人欧美精品刺激| 亚洲狠狠婷婷综合久久图片| a级毛片a级免费在线| or卡值多少钱| 母亲3免费完整高清在线观看| 毛片女人毛片| 熟女人妻精品中文字幕| 国产单亲对白刺激| 老司机福利观看| 欧美中文综合在线视频| 蜜桃久久精品国产亚洲av| 欧美大码av| 九色国产91popny在线| 观看免费一级毛片| 一区二区三区国产精品乱码| 嫩草影视91久久| 男人舔奶头视频| 此物有八面人人有两片| 欧美成狂野欧美在线观看| 好男人在线观看高清免费视频| 香蕉国产在线看| 18禁黄网站禁片午夜丰满| 很黄的视频免费| 国产视频内射| 女人被狂操c到高潮| 精品欧美国产一区二区三| 亚洲专区国产一区二区| 听说在线观看完整版免费高清| 色播亚洲综合网| 欧美中文综合在线视频| 国产一级毛片七仙女欲春2| 国产私拍福利视频在线观看| 亚洲精品久久国产高清桃花| 日韩中文字幕欧美一区二区| 国产三级中文精品| 国产成人精品久久二区二区91| 久久久久国产一级毛片高清牌| 在线a可以看的网站| 精品乱码久久久久久99久播| 国内精品久久久久久久电影| 亚洲专区国产一区二区| 久久中文字幕人妻熟女| 欧美大码av| 黄色 视频免费看| 极品教师在线免费播放| 欧美大码av| 久久欧美精品欧美久久欧美| 免费看十八禁软件| 在线观看舔阴道视频| 草草在线视频免费看| 久久午夜亚洲精品久久| 国产毛片a区久久久久| 丰满人妻一区二区三区视频av | 欧美另类亚洲清纯唯美| 人妻夜夜爽99麻豆av| 久久中文字幕人妻熟女| 国产高清三级在线| 人人妻,人人澡人人爽秒播| 草草在线视频免费看| 性色avwww在线观看| 啦啦啦观看免费观看视频高清| 日本五十路高清| 亚洲精品色激情综合| 久久亚洲真实| 亚洲 欧美一区二区三区| av福利片在线观看| 国产精品一区二区三区四区久久| 欧美丝袜亚洲另类 | 男女视频在线观看网站免费| 成人欧美大片| 国产高清激情床上av| 久久午夜亚洲精品久久| 色老头精品视频在线观看| 亚洲av美国av| 久久草成人影院| 成人特级黄色片久久久久久久| 综合色av麻豆| 日韩欧美国产在线观看| 国产野战对白在线观看| 禁无遮挡网站| 99国产精品一区二区三区| 久久中文看片网| 免费在线观看亚洲国产| 亚洲成人免费电影在线观看| 国内毛片毛片毛片毛片毛片| 国产高清三级在线| 日韩 欧美 亚洲 中文字幕| 久久午夜亚洲精品久久| 激情在线观看视频在线高清| 欧美性猛交黑人性爽| 熟女少妇亚洲综合色aaa.| 日本黄大片高清| 叶爱在线成人免费视频播放| 国内精品久久久久精免费| 日韩成人在线观看一区二区三区| 真人一进一出gif抽搐免费| 特级一级黄色大片| 亚洲av成人av| 熟女少妇亚洲综合色aaa.| 午夜福利高清视频| 久久精品91无色码中文字幕| 香蕉久久夜色| 十八禁网站免费在线| 又爽又黄无遮挡网站| 午夜两性在线视频| 亚洲欧美日韩卡通动漫| 人人妻人人看人人澡| 久久婷婷人人爽人人干人人爱| 国产精品美女特级片免费视频播放器 | 免费av不卡在线播放| 亚洲国产欧美人成| 国产人伦9x9x在线观看| 天天躁狠狠躁夜夜躁狠狠躁| av女优亚洲男人天堂 | 毛片女人毛片| 欧美精品啪啪一区二区三区| 成年女人永久免费观看视频| 天堂影院成人在线观看| 亚洲欧美日韩无卡精品| 成人av在线播放网站| 手机成人av网站| 亚洲中文av在线| 国内久久婷婷六月综合欲色啪| 欧美+亚洲+日韩+国产| a级毛片在线看网站| 人人妻人人看人人澡| 国产精品一及| 成人三级做爰电影| 97碰自拍视频| 男人和女人高潮做爰伦理| 久久久久免费精品人妻一区二区| 国产精品一区二区三区四区久久| 国内久久婷婷六月综合欲色啪| 亚洲精品一区av在线观看| 美女cb高潮喷水在线观看 | 国产成人影院久久av| 欧美乱码精品一区二区三区| 亚洲精品美女久久久久99蜜臀| 中出人妻视频一区二区| 亚洲av成人不卡在线观看播放网| 青草久久国产| h日本视频在线播放| 一个人免费在线观看的高清视频| 精品一区二区三区av网在线观看| 草草在线视频免费看| 夜夜夜夜夜久久久久| 日本黄大片高清| 99久久无色码亚洲精品果冻| 国产成+人综合+亚洲专区| 亚洲精品美女久久久久99蜜臀| 亚洲国产中文字幕在线视频| 亚洲自拍偷在线| 午夜精品久久久久久毛片777| 国语自产精品视频在线第100页| 精品久久久久久久毛片微露脸| 两个人看的免费小视频| 精品国产三级普通话版| www国产在线视频色| 久久久久久人人人人人| 久久性视频一级片| 国产黄片美女视频| 观看美女的网站| 国产1区2区3区精品| 精品国产美女av久久久久小说| 欧美3d第一页| 久久欧美精品欧美久久欧美| 亚洲人成电影免费在线| 精品国产亚洲在线| 成人无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久| 国产三级黄色录像| 夜夜躁狠狠躁天天躁| 久久精品aⅴ一区二区三区四区| 亚洲欧美激情综合另类| 久久久久久人人人人人| 欧美色视频一区免费| 在线观看66精品国产| 日韩欧美免费精品| 99久久99久久久精品蜜桃| 哪里可以看免费的av片| 久久热在线av| 久久精品国产清高在天天线| www.精华液| 午夜福利18| 在线观看日韩欧美| 国产精品免费一区二区三区在线| 国产v大片淫在线免费观看| 国产三级中文精品| bbb黄色大片| 国产极品精品免费视频能看的| 色播亚洲综合网| 真人做人爱边吃奶动态| 岛国视频午夜一区免费看| 亚洲国产精品久久男人天堂| 一a级毛片在线观看| 久久这里只有精品19| 国产一区二区三区视频了| 精品久久久久久久末码| 手机成人av网站| 亚洲av第一区精品v没综合| 小蜜桃在线观看免费完整版高清| 亚洲国产欧美人成| 免费一级毛片在线播放高清视频| 亚洲国产欧洲综合997久久,| 久久人妻av系列| 偷拍熟女少妇极品色| 亚洲第一电影网av| 日韩欧美一区二区三区在线观看| 精品99又大又爽又粗少妇毛片 | 中文字幕人妻丝袜一区二区| 三级国产精品欧美在线观看 | 99国产精品一区二区蜜桃av| 午夜精品久久久久久毛片777| 精品99又大又爽又粗少妇毛片 | 久久精品综合一区二区三区| 国产黄色小视频在线观看| 一边摸一边抽搐一进一小说| 99在线视频只有这里精品首页| 久久草成人影院| 午夜成年电影在线免费观看| 少妇的逼水好多| 全区人妻精品视频| 国产精品亚洲一级av第二区| 亚洲国产精品999在线| 无遮挡黄片免费观看| 国产一区二区在线av高清观看| 亚洲精华国产精华精| 可以在线观看的亚洲视频| 欧美成人免费av一区二区三区| 色av中文字幕| 97超级碰碰碰精品色视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 在线观看一区二区三区| 美女午夜性视频免费| 久久精品aⅴ一区二区三区四区| 欧美av亚洲av综合av国产av| 日韩 欧美 亚洲 中文字幕| av国产免费在线观看| av国产免费在线观看| 在线免费观看的www视频| 久久香蕉国产精品| 久久亚洲真实| 在线看三级毛片| 成年人黄色毛片网站| 一二三四在线观看免费中文在| cao死你这个sao货| 激情在线观看视频在线高清| 天天躁日日操中文字幕| 婷婷精品国产亚洲av| 亚洲aⅴ乱码一区二区在线播放| 婷婷亚洲欧美| 日本 欧美在线| 久久国产精品影院| 十八禁人妻一区二区| 亚洲人成网站在线播放欧美日韩| a在线观看视频网站| 国产视频一区二区在线看| 男插女下体视频免费在线播放| 亚洲电影在线观看av| 成年人黄色毛片网站| 麻豆久久精品国产亚洲av| 久久久国产成人免费| 国产高清视频在线播放一区| 亚洲成人久久爱视频| 日韩欧美国产一区二区入口| 国产69精品久久久久777片 | 黄色丝袜av网址大全| 国产精品一及| 国内精品久久久久精免费| 变态另类丝袜制服| 成人精品一区二区免费| 国产精品精品国产色婷婷| 好看av亚洲va欧美ⅴa在| 99国产综合亚洲精品| 亚洲一区二区三区不卡视频| 亚洲av第一区精品v没综合| 色老头精品视频在线观看| 老鸭窝网址在线观看| 婷婷精品国产亚洲av| a在线观看视频网站| 不卡一级毛片| 欧美日韩福利视频一区二区| 波多野结衣高清无吗| av片东京热男人的天堂| 一级a爱片免费观看的视频| 午夜福利成人在线免费观看| 国产蜜桃级精品一区二区三区| 国产成人啪精品午夜网站| 国产亚洲精品久久久com| 亚洲人成网站在线播放欧美日韩| 欧美日韩一级在线毛片| 国产三级在线视频| 老司机午夜福利在线观看视频| 国产精品av视频在线免费观看| 午夜久久久久精精品| 亚洲精品色激情综合| 少妇的逼水好多| 亚洲av第一区精品v没综合| 三级男女做爰猛烈吃奶摸视频| 久久性视频一级片| 99国产极品粉嫩在线观看| 精品久久蜜臀av无| 亚洲精品一区av在线观看| 精品久久久久久成人av| 亚洲精品色激情综合| 中文字幕av在线有码专区| 亚洲九九香蕉| 国产乱人伦免费视频| 一本一本综合久久| 国产三级在线视频| 在线播放国产精品三级| 长腿黑丝高跟| 国产伦一二天堂av在线观看| 国产成人精品久久二区二区91| 老熟妇乱子伦视频在线观看| 18禁裸乳无遮挡免费网站照片| 精品国产超薄肉色丝袜足j| 国产精品美女特级片免费视频播放器 | 亚洲七黄色美女视频| 国产精品一区二区精品视频观看| 精品欧美国产一区二区三| 国产三级中文精品| 看片在线看免费视频| 五月玫瑰六月丁香| 最近最新中文字幕大全电影3| 国产又色又爽无遮挡免费看| 国产97色在线日韩免费| 日韩欧美国产一区二区入口| 亚洲成a人片在线一区二区| 国产精品香港三级国产av潘金莲| 日韩大尺度精品在线看网址| 国产精品女同一区二区软件 | 久久久久九九精品影院| 久久久久亚洲av毛片大全| 老司机福利观看| 日本精品一区二区三区蜜桃| 久久精品亚洲精品国产色婷小说| 十八禁人妻一区二区| 日韩欧美国产一区二区入口| 日韩欧美在线二视频| 一二三四社区在线视频社区8| 91在线观看av| 欧美成人一区二区免费高清观看 | 成人三级做爰电影| 亚洲av美国av| 无限看片的www在线观看| 俺也久久电影网| av天堂在线播放| www.熟女人妻精品国产| 90打野战视频偷拍视频| av天堂中文字幕网| av女优亚洲男人天堂 | 国内精品久久久久精免费| 在线观看美女被高潮喷水网站 | 人妻丰满熟妇av一区二区三区| 叶爱在线成人免费视频播放| 午夜两性在线视频| 国产97色在线日韩免费| 久久久久久久久久黄片| 9191精品国产免费久久| 亚洲精品久久国产高清桃花| 国产av在哪里看| 一本精品99久久精品77| 日韩三级视频一区二区三区| 99国产精品99久久久久| 嫩草影院精品99| 老司机在亚洲福利影院| 操出白浆在线播放| 国内少妇人妻偷人精品xxx网站 | 日本成人三级电影网站| 麻豆成人午夜福利视频| 精华霜和精华液先用哪个| 精品99又大又爽又粗少妇毛片 | 九色国产91popny在线| 欧美一区二区国产精品久久精品| 国产精品国产高清国产av| 欧美色视频一区免费| 18禁黄网站禁片免费观看直播| 亚洲av成人精品一区久久| 精品国产乱码久久久久久男人| 日本黄色视频三级网站网址| 国产成人欧美在线观看| 88av欧美| 岛国在线观看网站| 欧美日韩一级在线毛片| 久久亚洲精品不卡| 99久久精品一区二区三区| 欧美黑人巨大hd| 亚洲av成人不卡在线观看播放网| 级片在线观看| 我要搜黄色片| 999久久久精品免费观看国产| 美女高潮的动态| 精品乱码久久久久久99久播| 亚洲自偷自拍图片 自拍| 久久性视频一级片| 五月伊人婷婷丁香| 可以在线观看毛片的网站| cao死你这个sao货| 欧美+亚洲+日韩+国产| 亚洲欧洲精品一区二区精品久久久| 国产一区二区三区在线臀色熟女| 九九在线视频观看精品| 首页视频小说图片口味搜索| 国产欧美日韩精品亚洲av| 亚洲午夜精品一区,二区,三区| bbb黄色大片| 少妇丰满av| 精品一区二区三区av网在线观看| 性色av乱码一区二区三区2| 12—13女人毛片做爰片一| 97超级碰碰碰精品色视频在线观看| 香蕉av资源在线| 国产精品亚洲美女久久久| 成年人黄色毛片网站| 色噜噜av男人的天堂激情| 一级毛片精品| 久久久久国内视频| 无遮挡黄片免费观看| 亚洲 欧美一区二区三区| 亚洲精华国产精华精| 国产男靠女视频免费网站| 淫妇啪啪啪对白视频| 亚洲国产精品成人综合色| 999精品在线视频| 精品久久久久久成人av| 天天添夜夜摸| 久久精品91蜜桃| 色综合欧美亚洲国产小说| 久久久久久大精品| 成年女人看的毛片在线观看| 欧美另类亚洲清纯唯美| 亚洲av电影在线进入| 成人特级黄色片久久久久久久| 婷婷丁香在线五月| 国产三级黄色录像| 国产精品亚洲一级av第二区| 国产v大片淫在线免费观看| 欧美一区二区国产精品久久精品| 熟女电影av网| 欧美成人免费av一区二区三区| 欧美性猛交黑人性爽| 欧美三级亚洲精品| 色av中文字幕| 舔av片在线| 亚洲中文字幕日韩| 欧美极品一区二区三区四区| 久久久国产精品麻豆| 精品一区二区三区视频在线观看免费| 欧美另类亚洲清纯唯美| 亚洲性夜色夜夜综合| 12—13女人毛片做爰片一| av黄色大香蕉|