• <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ì)
    久久久久视频综合| 亚洲综合精品二区| 性色avwww在线观看| 亚洲综合精品二区| videos熟女内射| 亚洲国产日韩一区二区| 免费高清在线观看视频在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美少妇被猛烈插入视频| 精品少妇内射三级| 精品亚洲乱码少妇综合久久| 欧美日韩综合久久久久久| 日本免费在线观看一区| 亚洲精品自拍成人| av福利片在线| 日韩大片免费观看网站| 亚洲av福利一区| 成人国产麻豆网| 免费观看a级毛片全部| 看非洲黑人一级黄片| 亚洲国产av影院在线观看| 欧美日韩在线观看h| 亚洲欧美成人综合另类久久久| freevideosex欧美| 一级二级三级毛片免费看| 成人无遮挡网站| 卡戴珊不雅视频在线播放| 亚洲精品乱码久久久久久按摩| 亚洲欧美一区二区三区黑人 | 国产精品免费大片| 99视频精品全部免费 在线| 亚洲第一区二区三区不卡| 伊人亚洲综合成人网| 日韩人妻高清精品专区| 国产成人精品无人区| 高清午夜精品一区二区三区| 成年美女黄网站色视频大全免费 | 日日爽夜夜爽网站| 99久国产av精品国产电影| 一区二区三区乱码不卡18| 久久99一区二区三区| 精品熟女少妇av免费看| 亚洲,欧美,日韩| 国产日韩欧美在线精品| 十八禁高潮呻吟视频| 国产精品国产三级国产专区5o| 国产免费视频播放在线视频| 精品卡一卡二卡四卡免费| 免费观看性生交大片5| 在现免费观看毛片| 国产深夜福利视频在线观看| 日日摸夜夜添夜夜添av毛片| 国产熟女午夜一区二区三区 | 久久婷婷青草| 久久久午夜欧美精品| 国产男人的电影天堂91| 日韩av在线免费看完整版不卡| 国产精品久久久久成人av| 免费av中文字幕在线| 亚洲情色 制服丝袜| 亚洲欧美日韩卡通动漫| 午夜福利影视在线免费观看| 在线亚洲精品国产二区图片欧美 | 久久99一区二区三区| 男女免费视频国产| 亚洲精品日韩av片在线观看| 国产日韩欧美在线精品| 飞空精品影院首页| 特大巨黑吊av在线直播| 免费黄频网站在线观看国产| 亚洲精品一二三| 波野结衣二区三区在线| 国产成人一区二区在线| 人妻夜夜爽99麻豆av| 亚洲国产av新网站| 成人漫画全彩无遮挡| 99九九在线精品视频| 国产精品国产三级国产专区5o| 免费黄网站久久成人精品| 99久久精品一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 久久99精品国语久久久| 久久人人爽人人爽人人片va| 久久国产亚洲av麻豆专区| 插阴视频在线观看视频| 狠狠精品人妻久久久久久综合| 一个人免费看片子| 国产成人一区二区在线| 91精品三级在线观看| 国产精品久久久久久久电影| 我的老师免费观看完整版| 国产午夜精品久久久久久一区二区三区| 亚洲精品自拍成人| 一区二区av电影网| 日本欧美国产在线视频| 人妻 亚洲 视频| 亚洲激情五月婷婷啪啪| 日日爽夜夜爽网站| tube8黄色片| 成年美女黄网站色视频大全免费 | 久久久a久久爽久久v久久| 日本黄大片高清| 少妇的逼好多水| 最近2019中文字幕mv第一页| 我的老师免费观看完整版| 亚洲成人av在线免费| av在线观看视频网站免费| 中文精品一卡2卡3卡4更新| 好男人视频免费观看在线| 亚洲成人一二三区av| 久久久a久久爽久久v久久| videosex国产| 哪个播放器可以免费观看大片| 免费av中文字幕在线| av不卡在线播放| 成人二区视频| 国产成人一区二区在线| av网站免费在线观看视频| 精品国产一区二区久久| 日韩av免费高清视频| 亚洲成人手机| 国内精品宾馆在线| 精品人妻熟女av久视频| 亚洲av在线观看美女高潮| 欧美变态另类bdsm刘玥| 人妻夜夜爽99麻豆av| 日本av免费视频播放| 中文精品一卡2卡3卡4更新| 制服人妻中文乱码| 成年人午夜在线观看视频| 亚洲精品456在线播放app| 一级爰片在线观看| 麻豆成人av视频| 91午夜精品亚洲一区二区三区| 成人黄色视频免费在线看| 美女国产视频在线观看| 黑人欧美特级aaaaaa片| 日本黄色片子视频| 日韩欧美一区视频在线观看| 免费观看av网站的网址| 国产成人免费观看mmmm| 9色porny在线观看| 在现免费观看毛片| 亚洲美女搞黄在线观看| 赤兔流量卡办理| 大陆偷拍与自拍| 伊人亚洲综合成人网| 日韩三级伦理在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲美女视频黄频| 国产成人av激情在线播放 | 熟女人妻精品中文字幕| 中文字幕亚洲精品专区| 亚洲人成网站在线观看播放| 国产精品.久久久| 亚洲av日韩在线播放| 亚洲精品日韩av片在线观看| a级毛片免费高清观看在线播放| 国产欧美日韩综合在线一区二区| 美女国产视频在线观看| 制服人妻中文乱码| 一级二级三级毛片免费看| 亚洲人成77777在线视频| 日韩精品免费视频一区二区三区 | 成人午夜精彩视频在线观看| 人妻一区二区av| 亚洲欧美精品自产自拍| 人人妻人人澡人人看| 伦理电影免费视频| 亚洲国产色片| 在线精品无人区一区二区三| 观看av在线不卡| 搡老乐熟女国产| 国产成人aa在线观看| av黄色大香蕉| 亚洲无线观看免费| 成人国语在线视频| 97在线视频观看| 大片电影免费在线观看免费| 久久久欧美国产精品| 亚洲丝袜综合中文字幕| 一级毛片黄色毛片免费观看视频| 中文字幕人妻丝袜制服| 黑人高潮一二区| 最近2019中文字幕mv第一页| 免费看不卡的av| 大话2 男鬼变身卡| 七月丁香在线播放| 一本色道久久久久久精品综合| 午夜激情久久久久久久| 亚洲欧美日韩卡通动漫| 久久99热6这里只有精品| 人妻 亚洲 视频| 满18在线观看网站| 久久久精品94久久精品| 精品国产国语对白av| 亚洲精品成人av观看孕妇| 国产深夜福利视频在线观看| 色5月婷婷丁香| 亚洲国产日韩一区二区| 性色avwww在线观看| 啦啦啦啦在线视频资源| 精品一区二区三卡| 欧美+日韩+精品| 精品熟女少妇av免费看| 最黄视频免费看| a级毛片免费高清观看在线播放| 校园人妻丝袜中文字幕| 在线 av 中文字幕| 又黄又爽又刺激的免费视频.| 国产一区二区三区综合在线观看 | 久久久久久久久久人人人人人人| 日韩av免费高清视频| 久热久热在线精品观看| 母亲3免费完整高清在线观看 | 国产黄色免费在线视频| 丝袜美足系列| 午夜免费观看性视频| 99热国产这里只有精品6| 国产av国产精品国产| 美女中出高潮动态图| 最黄视频免费看| 亚洲欧美中文字幕日韩二区| 免费高清在线观看日韩| 国产黄片视频在线免费观看| 少妇被粗大的猛进出69影院 | 一本—道久久a久久精品蜜桃钙片| 欧美精品亚洲一区二区| 久久人妻熟女aⅴ| 日本av手机在线免费观看| 少妇精品久久久久久久| 亚洲美女视频黄频| 免费看av在线观看网站| 亚洲欧美一区二区三区黑人 | 少妇的逼好多水| 久久久久久久大尺度免费视频| 成年美女黄网站色视频大全免费 | 国产成人精品福利久久| 成人漫画全彩无遮挡| 国产精品麻豆人妻色哟哟久久| 亚洲激情五月婷婷啪啪| av不卡在线播放| 欧美日本中文国产一区发布| 狂野欧美激情性xxxx在线观看| 蜜桃久久精品国产亚洲av| 视频在线观看一区二区三区| 大香蕉97超碰在线| 亚洲成人av在线免费| 久久综合国产亚洲精品| 一区二区日韩欧美中文字幕 | 中文精品一卡2卡3卡4更新| 丝袜在线中文字幕| 国产精品国产三级专区第一集| 在线观看免费高清a一片| 午夜免费男女啪啪视频观看| 久久久久视频综合| 久久精品国产亚洲av天美| 亚洲精品久久久久久婷婷小说| 91精品伊人久久大香线蕉| 另类亚洲欧美激情| 又粗又硬又长又爽又黄的视频| 国产69精品久久久久777片| 日本黄大片高清| 91精品国产国语对白视频| 考比视频在线观看| 2018国产大陆天天弄谢| 亚洲欧美色中文字幕在线| 欧美bdsm另类| 91精品国产九色| 精品人妻偷拍中文字幕| av不卡在线播放| av有码第一页| 99九九线精品视频在线观看视频| 国产成人精品一,二区| 26uuu在线亚洲综合色| 在线免费观看不下载黄p国产| 成人亚洲精品一区在线观看| 91在线精品国自产拍蜜月| 考比视频在线观看| 精品久久国产蜜桃| 欧美97在线视频| 国产一区二区在线观看av| 性色av一级| 成人毛片60女人毛片免费| 精品国产国语对白av| 国产探花极品一区二区| 人人澡人人妻人| 午夜免费观看性视频| 桃花免费在线播放| 国产在线一区二区三区精| 亚洲天堂av无毛| 国产精品久久久久久精品电影小说| 国产成人freesex在线| 男女高潮啪啪啪动态图| 久久久精品免费免费高清| 久久久久视频综合| av在线老鸭窝| 色婷婷av一区二区三区视频| 99久久人妻综合| 亚洲av国产av综合av卡| 国产片特级美女逼逼视频| 精品国产一区二区三区久久久樱花| 黑人猛操日本美女一级片| 欧美丝袜亚洲另类| 涩涩av久久男人的天堂| 亚洲国产最新在线播放| 一本久久精品| 国产高清三级在线| 久久精品国产亚洲av涩爱| 精品国产一区二区久久| 亚洲精品国产av成人精品| 久久久国产欧美日韩av| 校园人妻丝袜中文字幕| 国产精品嫩草影院av在线观看| 亚洲人成网站在线播| 精品酒店卫生间| 久久久久久久久大av| 国产亚洲精品久久久com| tube8黄色片| 午夜福利网站1000一区二区三区| 亚洲精品av麻豆狂野| 久久av网站| 免费av中文字幕在线| 黄色欧美视频在线观看| 女人精品久久久久毛片| 老女人水多毛片| 九九久久精品国产亚洲av麻豆| 啦啦啦啦在线视频资源| 黄片无遮挡物在线观看| 亚洲精品中文字幕在线视频| 亚洲美女黄色视频免费看| 免费黄色在线免费观看| freevideosex欧美| 国产 精品1| 寂寞人妻少妇视频99o| 嫩草影院入口| 久久国产精品大桥未久av| 中国国产av一级| 少妇 在线观看| 制服人妻中文乱码| 久久精品国产亚洲网站| 97超碰精品成人国产| 亚洲精品自拍成人| 嫩草影院入口| 一区二区日韩欧美中文字幕 | 狂野欧美白嫩少妇大欣赏| 久久久精品免费免费高清| 美女大奶头黄色视频| 国产精品欧美亚洲77777| 热re99久久精品国产66热6| 欧美最新免费一区二区三区| 18禁裸乳无遮挡动漫免费视频| 男人爽女人下面视频在线观看| 精品久久久噜噜| 国产一区二区三区av在线| 好男人视频免费观看在线| 考比视频在线观看| 26uuu在线亚洲综合色| 春色校园在线视频观看| 久久狼人影院| 亚洲精品日韩在线中文字幕| 9色porny在线观看| 亚洲精品乱码久久久v下载方式| 夫妻性生交免费视频一级片| 久久久精品94久久精品| 自线自在国产av| 国产精品 国内视频| 久久人人爽av亚洲精品天堂| 欧美成人午夜免费资源| 91久久精品电影网| 午夜免费鲁丝| 一级爰片在线观看| 69精品国产乱码久久久| 美女cb高潮喷水在线观看| 午夜免费男女啪啪视频观看| 欧美bdsm另类| a 毛片基地| 看免费成人av毛片| 18禁裸乳无遮挡动漫免费视频| 国产黄色视频一区二区在线观看| 国产女主播在线喷水免费视频网站| 亚洲欧美清纯卡通| 免费观看在线日韩| 国产精品一区www在线观看| av网站免费在线观看视频| 欧美精品亚洲一区二区| 久久免费观看电影| 丰满饥渴人妻一区二区三| 中文欧美无线码| 新久久久久国产一级毛片| 蜜臀久久99精品久久宅男| 男男h啪啪无遮挡| 男女啪啪激烈高潮av片| 日本黄色日本黄色录像| 五月伊人婷婷丁香| 狂野欧美白嫩少妇大欣赏| 精品熟女少妇av免费看| 午夜老司机福利剧场| 国产探花极品一区二区| 制服人妻中文乱码| 97在线人人人人妻| 3wmmmm亚洲av在线观看| 国产 精品1| 亚洲成人av在线免费| 啦啦啦中文免费视频观看日本| 18禁观看日本| 街头女战士在线观看网站| 亚洲精品第二区| 欧美精品一区二区大全| 国内精品宾馆在线| 久久国产精品大桥未久av| 91在线精品国自产拍蜜月| 蜜桃在线观看..| 精品久久国产蜜桃| 一个人免费看片子| 九九久久精品国产亚洲av麻豆| 久久精品久久久久久久性| 亚洲精品久久久久久婷婷小说| 美女cb高潮喷水在线观看| 久久久国产精品麻豆| 在线观看美女被高潮喷水网站| 久热这里只有精品99| 精品一区二区免费观看| 丝袜在线中文字幕| 制服诱惑二区| 啦啦啦啦在线视频资源| 日韩一本色道免费dvd| 男人爽女人下面视频在线观看| 亚洲精品美女久久av网站| 亚洲精品国产色婷婷电影| 国产精品麻豆人妻色哟哟久久| 欧美日韩精品成人综合77777| 国产成人aa在线观看| 亚洲精品aⅴ在线观看| 成人无遮挡网站| 看免费成人av毛片| 欧美少妇被猛烈插入视频| 免费黄频网站在线观看国产| 五月伊人婷婷丁香| 少妇人妻精品综合一区二区| 亚洲第一区二区三区不卡| 日韩大片免费观看网站| 久久久久国产精品人妻一区二区| 一个人免费看片子| 亚洲成人av在线免费| 亚洲精品av麻豆狂野| 另类亚洲欧美激情| 午夜免费鲁丝| 国产男女超爽视频在线观看| 欧美最新免费一区二区三区| 能在线免费看毛片的网站| 永久网站在线| 久久久久国产精品人妻一区二区| 91国产中文字幕| 久久久久久人妻| 亚州av有码| 日韩一区二区视频免费看| 久久ye,这里只有精品| 欧美日本中文国产一区发布| 一区二区三区免费毛片| 成人毛片a级毛片在线播放| 韩国高清视频一区二区三区| 久久精品国产自在天天线| 久久久a久久爽久久v久久| 夫妻午夜视频| 国产在线视频一区二区| 能在线免费看毛片的网站| 在线看a的网站| 亚洲欧美日韩另类电影网站| av免费观看日本| 涩涩av久久男人的天堂| 2018国产大陆天天弄谢| 午夜福利,免费看| 91精品伊人久久大香线蕉| 久久人妻熟女aⅴ| 最近2019中文字幕mv第一页| av福利片在线| 内地一区二区视频在线| 国产成人免费观看mmmm| 精品一品国产午夜福利视频| 视频中文字幕在线观看| 久久久久久久久久久免费av| 免费日韩欧美在线观看| 成人手机av| 亚洲精品一区蜜桃| 夫妻午夜视频| 日本黄色日本黄色录像| 能在线免费看毛片的网站| 久久久久久久久久久久大奶| 香蕉精品网在线| 亚洲精品乱久久久久久| a级片在线免费高清观看视频| 国产日韩一区二区三区精品不卡 | 汤姆久久久久久久影院中文字幕| 亚洲av综合色区一区| 国产成人91sexporn| 久久这里有精品视频免费| 欧美性感艳星| 色94色欧美一区二区| 亚洲精品日本国产第一区| 国产精品成人在线| 十八禁高潮呻吟视频| 国国产精品蜜臀av免费| 国产在线视频一区二区| 日本欧美国产在线视频| 如日韩欧美国产精品一区二区三区 | 老司机亚洲免费影院| 寂寞人妻少妇视频99o| 亚洲av免费高清在线观看| 狂野欧美激情性bbbbbb| 看免费成人av毛片| 久久人人爽av亚洲精品天堂| 免费看不卡的av| 欧美另类一区| 国产又色又爽无遮挡免| 亚洲欧洲日产国产| 久久这里有精品视频免费| 国产男人的电影天堂91| 国产高清三级在线| 3wmmmm亚洲av在线观看| 亚洲欧美日韩另类电影网站| 日本黄大片高清| 久久精品国产亚洲av涩爱| 欧美性感艳星| 欧美人与善性xxx| 久久久久人妻精品一区果冻| 亚洲人成网站在线观看播放| 伦理电影免费视频| 亚洲精品一二三| 9色porny在线观看| 精品少妇黑人巨大在线播放| 亚洲精品久久午夜乱码| 国产成人freesex在线| 热99国产精品久久久久久7| 欧美精品一区二区免费开放| 久久久久久久亚洲中文字幕| 国产国拍精品亚洲av在线观看| 国产无遮挡羞羞视频在线观看| 色婷婷久久久亚洲欧美| 女性生殖器流出的白浆| 少妇猛男粗大的猛烈进出视频| 欧美日本中文国产一区发布| 国产免费视频播放在线视频| 大香蕉久久网| 日日爽夜夜爽网站| 看免费成人av毛片| 午夜久久久在线观看| tube8黄色片| 国产伦精品一区二区三区视频9| 精品一区二区三区视频在线| 夜夜看夜夜爽夜夜摸| 99久国产av精品国产电影| 插逼视频在线观看| 一本久久精品| 欧美三级亚洲精品| 国产老妇伦熟女老妇高清| 亚洲av欧美aⅴ国产| 一区二区三区乱码不卡18| 亚洲欧美一区二区三区黑人 | 哪个播放器可以免费观看大片| 久久99热这里只频精品6学生| 午夜福利视频在线观看免费| 精品久久久久久电影网| 麻豆精品久久久久久蜜桃| av网站免费在线观看视频| 最近2019中文字幕mv第一页| 亚洲精品视频女| 久久精品久久久久久久性| 看非洲黑人一级黄片| 精品亚洲成a人片在线观看| 国产一区亚洲一区在线观看| 在线观看免费视频网站a站| 狂野欧美激情性xxxx在线观看| 免费看不卡的av| 色视频在线一区二区三区| 久久久久国产网址| 国产精品嫩草影院av在线观看| 欧美亚洲 丝袜 人妻 在线| 国产高清三级在线| 91在线精品国自产拍蜜月| 22中文网久久字幕| av在线app专区| 国产成人av激情在线播放 | 亚洲国产欧美在线一区| 九草在线视频观看| 97精品久久久久久久久久精品| 少妇人妻 视频| av网站免费在线观看视频| 男人操女人黄网站| 男男h啪啪无遮挡| 国产高清三级在线| 成人二区视频| 久久久亚洲精品成人影院| 最近最新中文字幕免费大全7| 国产成人精品久久久久久| 久久韩国三级中文字幕| 18+在线观看网站| 99久久中文字幕三级久久日本| 一级毛片我不卡| 五月开心婷婷网| 制服人妻中文乱码| 亚洲欧美色中文字幕在线| 日本wwww免费看| 嘟嘟电影网在线观看| 永久免费av网站大全| 一本大道久久a久久精品| 久久国产精品大桥未久av| 日日摸夜夜添夜夜爱| 在线精品无人区一区二区三| 搡老乐熟女国产| 欧美老熟妇乱子伦牲交| 又粗又硬又长又爽又黄的视频| 青春草国产在线视频| 久久久亚洲精品成人影院| 男人爽女人下面视频在线观看|