• <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ì)
    日本av免费视频播放| 国产精品亚洲av一区麻豆 | 在线观看免费日韩欧美大片| 在线天堂最新版资源| 久久女婷五月综合色啪小说| 又大又爽又粗| xxx大片免费视频| 卡戴珊不雅视频在线播放| www日本在线高清视频| 一区二区三区激情视频| 久久国产精品男人的天堂亚洲| 欧美黄色片欧美黄色片| 亚洲av欧美aⅴ国产| 欧美精品一区二区免费开放| 在线观看免费视频网站a站| 99久久综合免费| 一二三四中文在线观看免费高清| 最近中文字幕高清免费大全6| 丰满乱子伦码专区| 自拍欧美九色日韩亚洲蝌蚪91| 日韩精品有码人妻一区| 一区二区av电影网| 操美女的视频在线观看| 精品国产国语对白av| 欧美成人精品欧美一级黄| 国产精品一区二区在线观看99| 天堂中文最新版在线下载| 免费在线观看黄色视频的| 亚洲情色 制服丝袜| 国产精品久久久久成人av| 久久久久久久大尺度免费视频| av视频免费观看在线观看| 免费女性裸体啪啪无遮挡网站| 在线观看免费午夜福利视频| 男男h啪啪无遮挡| 性色av一级| h视频一区二区三区| 国产av码专区亚洲av| 日韩伦理黄色片| 久久97久久精品| 一级黄片播放器| 久久久精品免费免费高清| 亚洲七黄色美女视频| 各种免费的搞黄视频| 日韩欧美精品免费久久| 黑人巨大精品欧美一区二区蜜桃| 你懂的网址亚洲精品在线观看| 日本vs欧美在线观看视频| 午夜日韩欧美国产| 日日啪夜夜爽| 亚洲七黄色美女视频| 自拍欧美九色日韩亚洲蝌蚪91| 老司机亚洲免费影院| 国产成人系列免费观看| 欧美xxⅹ黑人| 亚洲国产中文字幕在线视频| 亚洲熟女毛片儿| 久久亚洲国产成人精品v| 亚洲精品aⅴ在线观看| 69精品国产乱码久久久| 久久久久久免费高清国产稀缺| 99久久人妻综合| 丝袜喷水一区| 午夜免费男女啪啪视频观看| 免费女性裸体啪啪无遮挡网站| 日韩制服骚丝袜av| 色视频在线一区二区三区| 又黄又粗又硬又大视频| 国产日韩欧美视频二区| 久久久久久久国产电影| 一本一本久久a久久精品综合妖精| 高清黄色对白视频在线免费看| 久久99一区二区三区| 国产极品天堂在线| 少妇被粗大猛烈的视频| xxx大片免费视频| 国产成人免费无遮挡视频| 欧美成人午夜精品| 女的被弄到高潮叫床怎么办| 涩涩av久久男人的天堂| 亚洲精品一区蜜桃| 黄片小视频在线播放| 精品福利永久在线观看| 午夜日本视频在线| 黄片无遮挡物在线观看| 欧美精品一区二区免费开放| 菩萨蛮人人尽说江南好唐韦庄| 成人18禁高潮啪啪吃奶动态图| 午夜91福利影院| 成人国产麻豆网| 最近中文字幕2019免费版| 51午夜福利影视在线观看| 夫妻性生交免费视频一级片| 亚洲第一区二区三区不卡| 这个男人来自地球电影免费观看 | av电影中文网址| 一级a爱视频在线免费观看| h视频一区二区三区| 亚洲精品视频女| 国产人伦9x9x在线观看| 成年av动漫网址| 久久国产精品大桥未久av| 久久久久视频综合| 亚洲欧美激情在线| 女人久久www免费人成看片| 亚洲欧洲国产日韩| 黄色视频在线播放观看不卡| 日韩,欧美,国产一区二区三区| 日韩不卡一区二区三区视频在线| 中文乱码字字幕精品一区二区三区| 男人添女人高潮全过程视频| 午夜av观看不卡| 天美传媒精品一区二区| 国产在线一区二区三区精| 黄频高清免费视频| 欧美另类一区| 免费看av在线观看网站| 精品免费久久久久久久清纯 | 曰老女人黄片| 亚洲视频免费观看视频| 亚洲国产av影院在线观看| 一本大道久久a久久精品| 一区福利在线观看| 看非洲黑人一级黄片| 天天躁日日躁夜夜躁夜夜| 亚洲国产看品久久| 99热国产这里只有精品6| netflix在线观看网站| 搡老岳熟女国产| 午夜日本视频在线| 亚洲色图综合在线观看| 巨乳人妻的诱惑在线观看| 久久天躁狠狠躁夜夜2o2o | 亚洲成人一二三区av| 女人久久www免费人成看片| 国产成人欧美在线观看 | 国产乱人偷精品视频| 极品人妻少妇av视频| 国产成人av激情在线播放| 高清视频免费观看一区二区| 男女免费视频国产| 高清av免费在线| 亚洲欧美清纯卡通| 国产亚洲精品第一综合不卡| 国产av精品麻豆| 中文字幕人妻熟女乱码| 日韩成人av中文字幕在线观看| 欧美精品av麻豆av| 热re99久久国产66热| xxx大片免费视频| 成人亚洲欧美一区二区av| 伊人亚洲综合成人网| 99香蕉大伊视频| 亚洲国产欧美在线一区| 中文字幕高清在线视频| 亚洲精品成人av观看孕妇| 国产免费又黄又爽又色| 在线看a的网站| 久久久久人妻精品一区果冻| 嫩草影视91久久| 国产高清不卡午夜福利| 国产精品秋霞免费鲁丝片| 成人漫画全彩无遮挡| 久热爱精品视频在线9| 国产成人午夜福利电影在线观看| 人成视频在线观看免费观看| 精品久久久久久电影网| 天美传媒精品一区二区| 午夜激情av网站| 女人被躁到高潮嗷嗷叫费观| 免费看不卡的av| 在线看a的网站| 亚洲av福利一区| 久久久久精品性色| 成人亚洲精品一区在线观看| 一级毛片电影观看| 亚洲专区中文字幕在线 | 亚洲av成人精品一二三区| 久久久久久久久免费视频了| 国语对白做爰xxxⅹ性视频网站| 女的被弄到高潮叫床怎么办| 尾随美女入室| 男的添女的下面高潮视频| 亚洲国产av新网站| 成人漫画全彩无遮挡| 欧美日韩精品网址| 97在线人人人人妻| 久久鲁丝午夜福利片| 日韩一卡2卡3卡4卡2021年| 国产 精品1| www.av在线官网国产| 国产精品亚洲av一区麻豆 | 日本午夜av视频| 国产精品久久久久成人av| 久久精品aⅴ一区二区三区四区| 亚洲av日韩在线播放| 精品一区二区三区四区五区乱码 | 男人操女人黄网站| √禁漫天堂资源中文www| 久久青草综合色| 欧美日韩精品网址| av在线app专区| 天天添夜夜摸| av网站在线播放免费| 妹子高潮喷水视频| 免费观看av网站的网址| 国产精品女同一区二区软件| 日韩 欧美 亚洲 中文字幕| 免费观看a级毛片全部| 一本大道久久a久久精品| 一本一本久久a久久精品综合妖精| av片东京热男人的天堂| 在线观看www视频免费| 一级毛片黄色毛片免费观看视频| 久久97久久精品| 91成人精品电影| 秋霞伦理黄片| 亚洲欧美激情在线| 人人妻人人添人人爽欧美一区卜| 久久久久网色| 在线观看一区二区三区激情| 国产男女超爽视频在线观看| 操美女的视频在线观看| 国产精品蜜桃在线观看| 日韩一本色道免费dvd| 美女福利国产在线| 国产午夜精品一二区理论片| 中文精品一卡2卡3卡4更新| 女的被弄到高潮叫床怎么办| 久久99精品国语久久久| 国产精品久久久久久精品电影小说| 19禁男女啪啪无遮挡网站| 国产男女内射视频| 波多野结衣av一区二区av| 91aial.com中文字幕在线观看| 不卡av一区二区三区| 十八禁网站网址无遮挡| 国产亚洲最大av| 亚洲国产精品国产精品| 人成视频在线观看免费观看| 国产麻豆69| 90打野战视频偷拍视频| 黄片无遮挡物在线观看| 久久综合国产亚洲精品| 热re99久久精品国产66热6| 两个人免费观看高清视频| 亚洲精品自拍成人| 一本—道久久a久久精品蜜桃钙片| 亚洲国产欧美网| 欧美日韩亚洲高清精品| 中文字幕人妻丝袜制服| 久久韩国三级中文字幕| 久久鲁丝午夜福利片| 免费av中文字幕在线| 久久久久久久国产电影| 91aial.com中文字幕在线观看| 一边摸一边抽搐一进一出视频| 女人久久www免费人成看片| 亚洲欧洲精品一区二区精品久久久 | 丝袜美足系列| 九草在线视频观看| 天堂8中文在线网| 成年动漫av网址| 精品一区二区免费观看| 国产精品蜜桃在线观看| 大陆偷拍与自拍| 人成视频在线观看免费观看| 久久久久视频综合| 美女福利国产在线| av线在线观看网站| 女人久久www免费人成看片| 欧美人与性动交α欧美精品济南到| 综合色丁香网| 久久久久久人人人人人| 国产精品.久久久| 成人黄色视频免费在线看| 久久午夜综合久久蜜桃| 亚洲国产欧美日韩在线播放| 高清视频免费观看一区二区| 国产精品成人在线| 看免费av毛片| 午夜激情av网站| 1024香蕉在线观看| 久久天躁狠狠躁夜夜2o2o | av一本久久久久| 国产深夜福利视频在线观看| 哪个播放器可以免费观看大片| 少妇 在线观看| 久久精品国产亚洲av高清一级| 久久午夜综合久久蜜桃| 久久久久久久大尺度免费视频| 精品人妻熟女毛片av久久网站| av福利片在线| 亚洲欧洲精品一区二区精品久久久 | 欧美国产精品va在线观看不卡| 国产精品久久久av美女十八| 日韩人妻精品一区2区三区| av在线播放精品| 宅男免费午夜| 久久久久人妻精品一区果冻| 亚洲视频免费观看视频| 亚洲精品国产色婷婷电影| 精品一区二区三卡| av网站在线播放免费| 91成人精品电影| 午夜老司机福利片| 亚洲图色成人| 成人国产麻豆网| 老司机靠b影院| 日韩大码丰满熟妇| 亚洲欧美日韩另类电影网站| 日韩人妻精品一区2区三区| 不卡av一区二区三区| 欧美最新免费一区二区三区| 丝袜美足系列| 免费高清在线观看日韩| 久久精品国产综合久久久| 国产一区二区激情短视频 | 欧美中文综合在线视频| av在线老鸭窝| 各种免费的搞黄视频| 十八禁高潮呻吟视频| 中文字幕人妻熟女乱码| 亚洲国产欧美日韩在线播放| 99精品久久久久人妻精品| 亚洲av综合色区一区| 婷婷色麻豆天堂久久| 国产片特级美女逼逼视频| 国产av码专区亚洲av| 精品卡一卡二卡四卡免费| 精品人妻熟女毛片av久久网站| 国产日韩一区二区三区精品不卡| 国产伦理片在线播放av一区| 亚洲国产精品一区三区| 亚洲,欧美,日韩| 亚洲av电影在线观看一区二区三区| 成人午夜精彩视频在线观看| 精品一区在线观看国产| 少妇被粗大的猛进出69影院| 久久久久人妻精品一区果冻| 国产毛片在线视频| 久久久久久久久免费视频了| 99re6热这里在线精品视频| 久久影院123| 久久久欧美国产精品| 国产精品 欧美亚洲| 亚洲国产看品久久| 男女边吃奶边做爰视频| 热99久久久久精品小说推荐| 成人18禁高潮啪啪吃奶动态图| 国产午夜精品一二区理论片| 中国国产av一级| 亚洲精品一二三| 亚洲精品国产色婷婷电影| 午夜福利一区二区在线看| √禁漫天堂资源中文www| 精品一品国产午夜福利视频| 精品一区在线观看国产| 亚洲国产精品999| 青春草亚洲视频在线观看| 国产午夜精品一二区理论片| 丰满迷人的少妇在线观看| 1024香蕉在线观看| 综合色丁香网| 精品酒店卫生间| 日日摸夜夜添夜夜爱| 国产欧美亚洲国产| 国产成人a∨麻豆精品| 亚洲中文av在线| 少妇猛男粗大的猛烈进出视频| 国产成人精品福利久久| 91精品伊人久久大香线蕉| videosex国产| 亚洲综合精品二区| 精品免费久久久久久久清纯 | 波多野结衣av一区二区av| 一级毛片 在线播放| 精品亚洲成国产av| 狂野欧美激情性xxxx| 亚洲第一av免费看| 美女扒开内裤让男人捅视频| 免费看不卡的av| 成人18禁高潮啪啪吃奶动态图| 亚洲,欧美精品.| 欧美日韩视频高清一区二区三区二| 亚洲久久久国产精品| 日本爱情动作片www.在线观看| 国产探花极品一区二区| 丁香六月欧美| 天天添夜夜摸| 在现免费观看毛片| 1024视频免费在线观看| 国产无遮挡羞羞视频在线观看| 日本av免费视频播放| 国产精品久久久久久久久免| 国产麻豆69| videos熟女内射| 亚洲一级一片aⅴ在线观看| 一级爰片在线观看| svipshipincom国产片| 久久久精品免费免费高清| 久久ye,这里只有精品| 国产伦理片在线播放av一区| 中文字幕av电影在线播放| 久久精品人人爽人人爽视色| 中文字幕精品免费在线观看视频| 大话2 男鬼变身卡| 日韩 亚洲 欧美在线| 1024视频免费在线观看| 亚洲精品国产区一区二| 亚洲欧美成人综合另类久久久| 国产乱人偷精品视频| 亚洲免费av在线视频| 国产一区有黄有色的免费视频| 色综合欧美亚洲国产小说| 午夜免费观看性视频| 男女国产视频网站| 国产在线视频一区二区| 久久精品国产综合久久久| 国产伦人伦偷精品视频| 狂野欧美激情性bbbbbb| 亚洲熟女精品中文字幕| av线在线观看网站| 啦啦啦在线观看免费高清www| 久热这里只有精品99| 亚洲精品国产色婷婷电影| 一本一本久久a久久精品综合妖精| 国产99久久九九免费精品| 久久精品人人爽人人爽视色| 欧美国产精品va在线观看不卡| 久久久久久久久久久久大奶| 妹子高潮喷水视频| 极品人妻少妇av视频| 一级爰片在线观看| 搡老乐熟女国产| 日本黄色日本黄色录像| 亚洲欧美一区二区三区久久| 国产成人精品久久二区二区91 | 黄频高清免费视频| 久热爱精品视频在线9| 亚洲美女黄色视频免费看| 午夜精品国产一区二区电影| 桃花免费在线播放| 精品一品国产午夜福利视频| 观看美女的网站| 久热这里只有精品99| 欧美人与善性xxx| 一级片'在线观看视频| 亚洲熟女毛片儿| 国产女主播在线喷水免费视频网站| 91精品伊人久久大香线蕉| 啦啦啦 在线观看视频| 国产视频首页在线观看| 欧美av亚洲av综合av国产av | 国产精品嫩草影院av在线观看| 国产精品一区二区在线观看99| 午夜精品国产一区二区电影| 不卡av一区二区三区| 90打野战视频偷拍视频| 国精品久久久久久国模美| 人人妻人人爽人人添夜夜欢视频| 国产成人啪精品午夜网站| 欧美变态另类bdsm刘玥| av在线播放精品| 精品国产露脸久久av麻豆| 国产老妇伦熟女老妇高清| 亚洲av中文av极速乱| 如何舔出高潮| 中文字幕色久视频| 热99久久久久精品小说推荐| 卡戴珊不雅视频在线播放| 成年美女黄网站色视频大全免费| 麻豆乱淫一区二区| 国产精品秋霞免费鲁丝片| 欧美av亚洲av综合av国产av | av线在线观看网站| 亚洲欧洲国产日韩| 999久久久国产精品视频| 69精品国产乱码久久久| 精品国产超薄肉色丝袜足j| 久久天躁狠狠躁夜夜2o2o | 最近的中文字幕免费完整| 天天躁夜夜躁狠狠久久av| 美女国产高潮福利片在线看| 大话2 男鬼变身卡| 一本久久精品| 人人妻人人澡人人爽人人夜夜| 亚洲精品av麻豆狂野| 日本黄色日本黄色录像| 黄色怎么调成土黄色| 亚洲男人天堂网一区| 精品第一国产精品| 午夜福利,免费看| 91精品国产国语对白视频| 欧美日韩一区二区视频在线观看视频在线| 秋霞在线观看毛片| 精品一区二区免费观看| 亚洲人成网站在线观看播放| 久久韩国三级中文字幕| 伊人亚洲综合成人网| 女的被弄到高潮叫床怎么办| 国产男人的电影天堂91| 99国产精品免费福利视频| 亚洲综合精品二区| 免费人妻精品一区二区三区视频| 精品人妻一区二区三区麻豆| 欧美精品av麻豆av| 久久精品人人爽人人爽视色| 亚洲欧美一区二区三区国产| 制服诱惑二区| 久久久国产欧美日韩av| 欧美久久黑人一区二区| 天天添夜夜摸| 韩国av在线不卡| 国产国语露脸激情在线看| 青春草视频在线免费观看| 观看美女的网站| 狠狠婷婷综合久久久久久88av| 精品少妇一区二区三区视频日本电影 | 黑人欧美特级aaaaaa片| 欧美日韩亚洲综合一区二区三区_| 一级黄片播放器| 菩萨蛮人人尽说江南好唐韦庄| 久久久精品免费免费高清| 卡戴珊不雅视频在线播放| 交换朋友夫妻互换小说| 亚洲专区中文字幕在线 | 不卡视频在线观看欧美| av不卡在线播放| 午夜av观看不卡| 免费少妇av软件| 久久青草综合色| 肉色欧美久久久久久久蜜桃| 天堂俺去俺来也www色官网| 最新在线观看一区二区三区 | 亚洲成人av在线免费| 秋霞在线观看毛片| 一级毛片 在线播放| 亚洲,欧美精品.| 最近手机中文字幕大全| 香蕉国产在线看| 国产日韩欧美在线精品| 久久99热这里只频精品6学生| 满18在线观看网站| 亚洲国产成人一精品久久久| 亚洲精品国产一区二区精华液| 国产成人精品福利久久| 国产精品一区二区精品视频观看| 国产精品秋霞免费鲁丝片| 久久久欧美国产精品| 亚洲国产最新在线播放| 久久精品亚洲av国产电影网| 天天躁狠狠躁夜夜躁狠狠躁| xxxhd国产人妻xxx| 满18在线观看网站| 天天躁夜夜躁狠狠躁躁| 伦理电影免费视频| 欧美日韩亚洲综合一区二区三区_| 激情视频va一区二区三区| 狂野欧美激情性xxxx| 久久久久精品人妻al黑| a 毛片基地| 日韩制服骚丝袜av| 欧美老熟妇乱子伦牲交| 超碰成人久久| av福利片在线| 国产精品久久久人人做人人爽| 午夜日韩欧美国产| 亚洲精品,欧美精品| 欧美av亚洲av综合av国产av | 制服丝袜香蕉在线| 三上悠亚av全集在线观看| 亚洲国产中文字幕在线视频| 欧美最新免费一区二区三区| 丰满饥渴人妻一区二区三| 亚洲成国产人片在线观看| 在线精品无人区一区二区三| av女优亚洲男人天堂| bbb黄色大片| 亚洲精品一二三| av不卡在线播放| 免费在线观看黄色视频的| 久久毛片免费看一区二区三区| av国产久精品久网站免费入址| av在线观看视频网站免费| 亚洲国产毛片av蜜桃av| 丁香六月欧美| 国产极品天堂在线| 日本欧美视频一区| 中文字幕高清在线视频| 国产精品偷伦视频观看了| 国产色婷婷99| 91精品伊人久久大香线蕉| av在线观看视频网站免费| 国产免费又黄又爽又色| 91精品国产国语对白视频| 亚洲人成电影观看| 国产色婷婷99| 午夜日韩欧美国产| 国产精品无大码| 久久精品熟女亚洲av麻豆精品| 王馨瑶露胸无遮挡在线观看| av在线观看视频网站免费| 欧美97在线视频| 91精品伊人久久大香线蕉| 国产黄频视频在线观看| 精品一区二区三卡| 99久久99久久久精品蜜桃| 久久精品亚洲av国产电影网| 男女无遮挡免费网站观看| netflix在线观看网站| 亚洲熟女精品中文字幕| av视频免费观看在线观看| 国产精品久久久久成人av| 视频区图区小说|