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

    人口老齡化與收入不平等的時(shí)變效應(yīng)

    2018-11-14 07:28:12呂康銀朱金霞
    稅務(wù)與經(jīng)濟(jì) 2018年6期
    關(guān)鍵詞:斷點(diǎn)協(xié)整人口老齡化

    孫 曄,呂康銀,朱金霞

    (1.東北師范大學(xué) 商學(xué)院,吉林 長春 130117; 2.魯東大學(xué) 商學(xué)院,山東 煙臺(tái) 264025)

    一、引 言

    隨著我國人口老齡化問題的凸顯,其與收入不平等之間的動(dòng)態(tài)關(guān)系已經(jīng)成為學(xué)術(shù)界研究的熱點(diǎn)之一,并且涌現(xiàn)出大量的相關(guān)研究文獻(xiàn)。作為我國經(jīng)濟(jì)社會(huì)發(fā)展中日益凸顯的兩大焦點(diǎn)問題,二者之間究竟存在著怎樣的關(guān)系,值得進(jìn)行深入的理論研究與探索。

    在實(shí)證研究方法上,可以將對該問題的研究分為基于線性模型的實(shí)證研究 (主要包括線性單方程回歸模型、協(xié)整分析、線性Granger因果檢驗(yàn)、線性向量自回歸模型和線性面板回歸模型等) 和基于非線性模型實(shí)證研究 (主要包括門限的單方程回歸、滾動(dòng)視窗Granger因果關(guān)系檢驗(yàn)、面板平滑遷移回歸模型等)。本文結(jié)合已有的研究思路,重點(diǎn)基于時(shí)域視角,使用含斷點(diǎn)的時(shí)變協(xié)整、線性Granger因果關(guān)系檢驗(yàn)、非線性Granger因果關(guān)系檢驗(yàn)法考察我國老齡化與收入不平等之間的時(shí)變的動(dòng)態(tài)關(guān)聯(lián)機(jī)制,并進(jìn)一步采用平滑遷移性向量自回歸模型 (STVAR) 檢驗(yàn)二者之間是否存在非對稱的影響及測度影響效果。

    二、數(shù)據(jù)選取與處理

    對于收入不平等的度量,已有研究文獻(xiàn)使用最多的指標(biāo)是基尼系數(shù) (董志強(qiáng)等,2012[1];劉金東等,2014[2];藍(lán)嘉俊等,2014[3];范洪敏和穆懷中,2018[4])。鑒于此,本文選擇我國基尼系數(shù)作為收入不平等的度量指標(biāo),記為Gini。對于老齡化的度量,存在大量的選擇空間,比如采用老年撫養(yǎng)比和少兒撫養(yǎng)比構(gòu)造度量老齡化的指標(biāo)*老齡化=老年撫養(yǎng)比/ (1+老年撫養(yǎng)比+少兒撫養(yǎng)比),具體參見藍(lán)嘉俊等 (2014)。和采用老年撫養(yǎng)比*65歲 (含65歲) 以上老年人口占15~64歲勞動(dòng)年齡人口的比重。度量老齡化。本文參考范洪敏和穆懷中 (2018)[4]的做法,選擇老年撫養(yǎng)比作為老齡化的代理變量,記為Aging。以上基尼系數(shù)原數(shù)據(jù)和老齡化原數(shù)據(jù)均來自于WIND數(shù)據(jù)庫,二者具體的實(shí)際序列如圖1、圖2所示。

    圖1 中國基尼系數(shù)序列 圖2 中國老齡化序列

    從圖1、圖2中繪制的原始數(shù)據(jù)序列可以看出我國基尼系數(shù)和老齡化之間在整體上存在緩慢的上升趨勢,即可能存在某一協(xié)整向量可以較好擬合二者間的長期相關(guān)關(guān)系;但在局部時(shí)間內(nèi)二者之間的協(xié)同性較差,具體表現(xiàn)在1990年至1999年間,基尼系數(shù)呈現(xiàn)出陡增驟降的特征,而此時(shí)的老齡化序列卻表現(xiàn)出具有較小波動(dòng)且逐漸上升的特點(diǎn);此外,在1999年之后基尼系數(shù)出現(xiàn)急劇上漲并持續(xù)保持較高的水平,相對應(yīng)的老齡化序列表現(xiàn)出穩(wěn)健上升的態(tài)勢。

    結(jié)合兩段時(shí)間內(nèi)每個(gè)時(shí)間序列變量均表現(xiàn)出獨(dú)自特點(diǎn)的現(xiàn)象,從統(tǒng)計(jì)學(xué)角度出發(fā),二者可能存在時(shí)變的協(xié)整關(guān)系,而此時(shí)如果采用常系數(shù)協(xié)整檢驗(yàn)方法可能有失準(zhǔn)確性和穩(wěn)健性。有鑒于此,本文主要的研究思路是首先檢驗(yàn)二者是否具有時(shí)變的協(xié)整關(guān)系;其次基于檢驗(yàn)結(jié)果通過構(gòu)建線性Granger因果關(guān)系檢驗(yàn)和非線性Granger因果關(guān)系檢驗(yàn)法進(jìn)一步確定二者之間的交互影響關(guān)系,即確定二者間的因果關(guān)系。

    三、我國收入不平等與老齡化的時(shí)變協(xié)整檢驗(yàn)

    基于圖1和圖2的描述,由于我國收入不平等與老齡化均在樣本期間表現(xiàn)出各自的特征,二者之間可能存在時(shí)變的協(xié)整關(guān)系,因此,需要構(gòu)造時(shí)變協(xié)整檢驗(yàn)?zāi)P蛠砜疾於咧g的長期關(guān)系是否穩(wěn)定。眾所周知,協(xié)整檢驗(yàn)是建立在變量均是同階單整序列基礎(chǔ)上的,因而在構(gòu)建時(shí)變協(xié)整模型前應(yīng)對時(shí)間序列數(shù)據(jù)進(jìn)行平穩(wěn)性檢驗(yàn)。

    (一)數(shù)據(jù)平穩(wěn)性檢驗(yàn)

    為了與數(shù)據(jù)選取中所標(biāo)記的符號一致,在數(shù)據(jù)平穩(wěn)性檢驗(yàn)中我們將基尼系數(shù)序列(Gini) 和老齡化序列(Aging)的一階差分分別標(biāo)記為△Gini和△Aging。同時(shí)為使數(shù)據(jù)的單位根檢驗(yàn)結(jié)果具有較好的穩(wěn)健性,分別采用ADF檢驗(yàn)法、PP檢驗(yàn)法和KPSS檢驗(yàn)法進(jìn)行數(shù)據(jù)的平穩(wěn)性檢驗(yàn)。數(shù)據(jù)平穩(wěn)性檢驗(yàn)結(jié)果如表1所示。

    表1 基尼系數(shù)和老齡化序列的單位根檢驗(yàn)結(jié)果

    注:“**”和“***”分別表示在5%、1%顯著概率水平下顯著。

    由表1的檢驗(yàn)結(jié)果可知,基于ADF檢驗(yàn)時(shí)基尼系數(shù)序列(Gini)和人口老齡化序列(Aging)均是非平穩(wěn)時(shí)間序列,其對應(yīng)的檢驗(yàn)統(tǒng)計(jì)值分別為-2.2705和-2.2930,均大于ADF檢驗(yàn)10%所對應(yīng)的臨界值-2.6422和-2.6299,即接受序列包含單位根的原假設(shè);而其差分序列的檢驗(yàn)統(tǒng)計(jì)值分別為-3.1061和-2.4622,均小于ADF檢驗(yàn)5%所對應(yīng)的臨界值-3.0124和-2.4622,即拒絕序列包含單位根的原假設(shè)。

    基于PP檢驗(yàn)時(shí)基尼系數(shù)序列(Gini)和人口老齡化序列(Aging)均是非平穩(wěn)時(shí)間序列,其對應(yīng)的檢驗(yàn)統(tǒng)計(jì)值分別為-2.4080和2.2930,均大于PP檢驗(yàn)10%所對應(yīng)的臨界值-2.6422和-2.6299,即接受序列包含單位根的原假設(shè)。而其差分序列的檢驗(yàn)統(tǒng)計(jì)值分別為-6.7653和-3.8872,均小于PP檢驗(yàn)1%所對應(yīng)的臨界值-3.7241,即拒絕序列包含單位根的原假設(shè)。同理,基于KPSS檢驗(yàn)時(shí),由于該檢驗(yàn)方法的原假設(shè)與ADF檢驗(yàn)和PP檢驗(yàn)的原假設(shè)相反,但檢驗(yàn)結(jié)果同樣表明基尼系數(shù)序列和人口老齡化序列均為非平穩(wěn)序列,而其一階差分序列為平穩(wěn)的時(shí)間序列。

    結(jié)合上述分析,可以得到基尼系數(shù)序列和人口老齡化序列均為一階單整序列,因此可以通過構(gòu)建時(shí)變協(xié)整檢驗(yàn)法捕捉二者之間是否具有時(shí)變的長期均衡關(guān)系,而不是僅僅存在一個(gè)或多個(gè)恒定的協(xié)整向量。

    (二)時(shí)變協(xié)整模型的構(gòu)建及檢驗(yàn)

    時(shí)變協(xié)整模型的構(gòu)建主要是識別出收入不平等和人口老齡化之間是否存在結(jié)構(gòu)斷點(diǎn),進(jìn)而導(dǎo)致常系數(shù)的協(xié)整檢驗(yàn)失效。通常情況下,進(jìn)行斷點(diǎn)識別使用最多的方法是Chow檢驗(yàn),但是該檢驗(yàn)方法需要外生設(shè)定斷點(diǎn),斷點(diǎn)的選擇往往依據(jù)的是重大經(jīng)濟(jì)事件,這一約束使得斷點(diǎn)獨(dú)立于樣本數(shù)據(jù)進(jìn)而導(dǎo)致傳統(tǒng)檢驗(yàn)統(tǒng)計(jì)量失效。考慮到上述方法的缺點(diǎn),參考Hansen (1992)[5]給出的一階單整時(shí)間序列的回歸模型參數(shù)是否具有時(shí)變性的檢驗(yàn)方法。該方法與Chow檢驗(yàn)最大的區(qū)別在于能夠內(nèi)生識別出結(jié)構(gòu)斷點(diǎn),并且可以進(jìn)一步判別協(xié)整方程的回歸系數(shù)是否具有非穩(wěn)定性。

    根據(jù)以上對基尼系數(shù)和人口老齡化原始數(shù)據(jù)進(jìn)行平穩(wěn)性檢驗(yàn)的結(jié)果,確定二者均為一階單整序列,可以通過構(gòu)建帶有變截距項(xiàng)的協(xié)整方程:

    yt=αt+βtxt+εt

    (1)

    其中yt代表基尼系數(shù)序列,xt代表人口老齡化序列,αt和βt分別代表時(shí)變截距項(xiàng)和時(shí)變系數(shù)項(xiàng)。εt為服從白噪聲過程的隨機(jī)誤差項(xiàng)。令Xt=(1,xt)′,Bt=(αt,βt)′,則式 (1) 可表示為:

    (2)

    如果式(2)中的協(xié)整向量不存在變異,則時(shí)變系數(shù)向量Bt為以常數(shù)。否則,如果上述協(xié)整向量存在變異,可以通過設(shè)定如下假設(shè)對上述結(jié)構(gòu)性斷點(diǎn)進(jìn)行檢驗(yàn)。假設(shè)t時(shí)刻為內(nèi)生結(jié)構(gòu)斷點(diǎn)(1

    (3)

    對于以上構(gòu)建的檢驗(yàn)協(xié)整方程是否具有時(shí)變特性的統(tǒng)計(jì)量,其中Fnt檢驗(yàn)統(tǒng)計(jì)量的計(jì)算最為便捷,并且當(dāng)樣本容量較大時(shí)該檢驗(yàn)統(tǒng)計(jì)量還有標(biāo)準(zhǔn)的漸進(jìn)分布。但如果結(jié)構(gòu)斷點(diǎn)并不是外生時(shí),即斷點(diǎn)的選擇和樣本數(shù)據(jù)具有顯著的相關(guān)關(guān)系,則Fnt檢驗(yàn)統(tǒng)計(jì)量對應(yīng)的檢驗(yàn)功效便會(huì)極大下降,導(dǎo)致檢驗(yàn)結(jié)果失準(zhǔn)。有鑒于此,Hansen(1992)[5]給出了SupF和MeanF檢驗(yàn)統(tǒng)計(jì)量,且二者能夠克服上述缺陷;但出現(xiàn)的另一個(gè)問題是二者計(jì)算相對復(fù)雜,沒有標(biāo)準(zhǔn)的漸進(jìn)分布。Hansen(1992)[5]利用蒙特卡洛模擬方法給出了SupF和MeanF檢驗(yàn)統(tǒng)計(jì)量1%、5%和10%顯著水平的臨界值,表2為基于SupF檢驗(yàn)統(tǒng)計(jì)量和MeanF檢驗(yàn)統(tǒng)計(jì)量的收入不平等和人口老齡化協(xié)整方程斷點(diǎn)檢驗(yàn)結(jié)果。

    表2 基于協(xié)整方程的結(jié)構(gòu)斷點(diǎn)檢驗(yàn)結(jié)果

    注:SupF和MeanF檢驗(yàn)統(tǒng)計(jì)量的1%、5%和10%的臨界值均來自Hansen (1992)。

    從表2的檢驗(yàn)結(jié)果可以看出,無論是SupF檢驗(yàn)統(tǒng)計(jì)量還是MeanF檢驗(yàn)統(tǒng)計(jì)量均能在10%的顯著水平下拒絕收入不平等和人口老齡化構(gòu)成的協(xié)整方程中協(xié)整向量不存在結(jié)構(gòu)斷點(diǎn)的原假設(shè);具體而言,SupF檢驗(yàn)統(tǒng)計(jì)值(11.40)大于10%臨界值(11.20),MeanF檢驗(yàn)統(tǒng)計(jì)值(7.91)大于1%臨界值(6.78)。這一檢驗(yàn)結(jié)果充分表明,在研究的樣本期間內(nèi),我國收入不平等和人口老齡化之間的長期均衡關(guān)系已經(jīng)發(fā)生了明顯的結(jié)構(gòu)變化。

    為了準(zhǔn)確識別出結(jié)構(gòu)斷點(diǎn)的具體時(shí)點(diǎn)位置,我們繪制了基于已知結(jié)構(gòu)斷點(diǎn)的F檢驗(yàn)統(tǒng)計(jì)量的時(shí)序圖,如圖3所示。從圖3可以看出,基于已知結(jié)構(gòu)斷點(diǎn)的F檢驗(yàn)統(tǒng)計(jì)量在1999年取得了最大值(對應(yīng)為11.40),這表明1999年前后我國收入不平等和人口老齡化之間的長期均衡關(guān)系發(fā)生了顯著的突變;而在1999年前后,對應(yīng)的F序列均低于10%顯著水平的臨界值,進(jìn)一步表明在統(tǒng)計(jì)學(xué)角度上,1999年我國收入不平等和人口老齡化之間的長期均衡機(jī)制發(fā)生了變異。

    圖3基尼系數(shù)與人口老齡化時(shí)變協(xié)整檢驗(yàn)

    我國收入不平等和人口老齡化之間的長期均衡關(guān)系在1999年發(fā)生了顯著突變的這一結(jié)果可以從原始數(shù)據(jù)中找到解釋?;趫D1和圖2中所繪制的原始數(shù)據(jù)時(shí)序圖,我們可以清晰地觀察到,1999年之前我國收入不平等呈現(xiàn)波幅大、增速快的特點(diǎn),而在1999年后則表現(xiàn)出緩慢上升、波幅減小的新特征;與之對應(yīng)的人口老齡化數(shù)據(jù)則在斷點(diǎn)前表現(xiàn)出增速較緩的特點(diǎn),而在斷點(diǎn)后呈現(xiàn)出增速加快且無絲毫下降趨勢的新特征。

    當(dāng)前全球幾乎所有的國家都面臨人口老齡化問題,只是程度不同而已。發(fā)達(dá)國家已經(jīng)進(jìn)入老齡化社會(huì),發(fā)展中國家也正在邁入老齡化社會(huì)。相較于歷史上大部分時(shí)期的人口結(jié)構(gòu),當(dāng)下我國人民生活水平大幅提高、社會(huì)保障制度不斷完善、醫(yī)療技術(shù)創(chuàng)新日新月異等因素帶來了人口平均預(yù)期壽命的增加,我國正在進(jìn)入持續(xù)40年的高速老齡化時(shí)期。出現(xiàn)老齡化這一新的人口現(xiàn)象,勢必會(huì)影響收入不平等和人口結(jié)構(gòu)的長期均衡關(guān)系,表現(xiàn)出顯著的時(shí)變特征。

    四、我國收入不平等與老齡化的因果關(guān)系檢驗(yàn)

    通常的情況下,檢驗(yàn)變量之間的因果關(guān)系多數(shù)基于線性Granger因果檢驗(yàn)方法。但當(dāng)變量之間存在時(shí)變的協(xié)整關(guān)系時(shí),或者說變量之間的長期均衡在某一時(shí)點(diǎn)發(fā)生了顯著變化,則基于線性Granger因果檢驗(yàn)方法得到的檢驗(yàn)結(jié)果往往會(huì)缺乏準(zhǔn)確性。結(jié)合前面的分析,我們知道我國收入不平等和人口老齡化之間存在時(shí)變的協(xié)整關(guān)系,因此,需要采用Péguin-Feissolle等(2008)[6]非線性Granger因果檢驗(yàn)法來識別出二者之間的非線性傳導(dǎo)機(jī)制。

    令△Gini和△Aging分別表示收入不平等的差分序列和人口老齡化的差分序列,且二者均為弱平穩(wěn)且遍歷的時(shí)間序列,則有二者構(gòu)成的二元隨機(jī)系統(tǒng)可以是:

    △Ginit=fG(△Ginit-1,△Ginit-2,…,△Ginit-p,△Agingt-1,△Agingt-2,…,△Agingt-p;φG)+εGt

    安徽六國化工是我國最早的磷酸二銨生產(chǎn)企業(yè),多年來,在磷資源高效利用上積累了豐富的技術(shù)和經(jīng)驗(yàn)。銅陵化學(xué)工業(yè)集團(tuán)黨委書記、董事長陳嘉生介紹:“作為與長江一路之隔的沿江企業(yè),六國化工十分重視環(huán)境保護(hù),不斷加大資金投入,環(huán)境整治成績顯著,實(shí)現(xiàn)了‘三廢’達(dá)標(biāo)排放,磷石膏綜合利用率達(dá)100%,在行業(yè)名列前茅?!彼硎?,平臺(tái)啟動(dòng)后,六國化工將繼續(xù)貫徹生態(tài)優(yōu)先、綠色發(fā)展的戰(zhàn)略決策,立足本業(yè),結(jié)合現(xiàn)代農(nóng)業(yè)發(fā)展需求,積極引進(jìn)國際國內(nèi)頂尖智力資源,不斷提高磷資源利用率和技術(shù)、服務(wù)水平,為長江經(jīng)濟(jì)帶區(qū)域綠色發(fā)展作出更大的貢獻(xiàn)。

    △Agingt=fA(△Ginit-1,△Ginit-2,…,△Ginit-p,△Agingt-1,△Agingt-2,…,△Agingt-p;φA)+εAt

    (4)

    基于二階泰勒近似和在二階泰勒近似基礎(chǔ)上提取主成分的刻畫△Ginit和△Agingt短期交互影響的動(dòng)態(tài)關(guān)系,均可采用單方程和方程系統(tǒng)的估計(jì)方法,構(gòu)建拉格朗日乘子形式的約束檢驗(yàn)統(tǒng)計(jì)量進(jìn)行Granger因果關(guān)系檢驗(yàn)。單方程和方程系統(tǒng)的拉格朗日乘子檢驗(yàn)統(tǒng)計(jì)量分別為:

    (5)

    (6)

    式(5)和式(6)分別表示基于單方程和方程系統(tǒng)估計(jì)方法得到的非線性Granger因果關(guān)系檢驗(yàn)拉格朗日乘子檢驗(yàn)統(tǒng)計(jì)量。其中RSS1和RSS0分別表示參數(shù)約束檢驗(yàn)中無約束回歸和約束回歸的殘差平方和,U1和U0分別為系統(tǒng)方程無約束和約束回歸對應(yīng)的T×2維的殘差矩陣。T表示樣本觀測值的個(gè)數(shù),m為系統(tǒng)中方程的個(gè)數(shù),tr(·)表示矩陣的跡。LMnonline-single和LMnonline-system檢驗(yàn)統(tǒng)計(jì)量均服從F分布。

    基于上述關(guān)于非線性Granger因果關(guān)系檢驗(yàn)方法,通過對單方程的拉格朗日乘子檢驗(yàn)統(tǒng)計(jì)值和方程系統(tǒng)的拉格朗日乘子檢驗(yàn)統(tǒng)計(jì)值的計(jì)算,我們得到收入不平等和人口老齡化之間的短期交互影響關(guān)系的檢驗(yàn)結(jié)果,如表3所示。此外,為充分捕捉到二者之間存在的非線性作用機(jī)制,我們同時(shí)采用了線性Granger因果檢驗(yàn)方法和非線性Granger因果檢驗(yàn)方法進(jìn)行實(shí)證研究。

    表3 變量間Granger因果檢驗(yàn)結(jié)果

    注:LMsingle和LMsystem分別基于單方程和方程系統(tǒng)構(gòu)造的檢驗(yàn)統(tǒng)計(jì)量,其中每一個(gè)統(tǒng)計(jì)量對應(yīng)的第一列為統(tǒng)計(jì)量檢驗(yàn)值,第二列為顯著性概率,表中加粗字體表示均在10%顯著概率下顯著。

    從表3的Granger因果檢驗(yàn)結(jié)果中可以清楚地看到,當(dāng)采用線性Granger因果關(guān)系檢驗(yàn)方法時(shí),得到的僅僅是存在收入不平等單向影響人口老齡化的結(jié)論。而當(dāng)采用非線性Granger因果關(guān)系檢驗(yàn)時(shí),則檢驗(yàn)結(jié)果表明收入不平等和人口老齡化在研究的樣本區(qū)間內(nèi)是互為因果關(guān)系的。

    基于非線性Granger因果檢驗(yàn)方法得到的結(jié)論更合理地捕捉到我國人口老齡化會(huì)影響收入不平等,而僅僅采用線性Granger因果檢驗(yàn)方法則未發(fā)現(xiàn)上述現(xiàn)象。這同時(shí)也表明人口老齡化作用于收入不平等存在非線性的傳導(dǎo)機(jī)制,而不是簡單地意味著人口結(jié)構(gòu)不同程度的改變均能帶來同等程度的收入不平等。

    對于收入不平等會(huì)影響到老齡化這一檢驗(yàn)結(jié)果,可以從我國少兒撫養(yǎng)比占總?cè)丝诘谋戎氐慕嵌冗M(jìn)一步理解。自20世紀(jì)90年代以來,我國少兒撫養(yǎng)比占總?cè)丝诘谋戎爻掷m(xù)下降,直到2010年左右才基本維持在近17%左右。[注]數(shù)據(jù)來源于WIND統(tǒng)計(jì)數(shù)據(jù)庫。其原因可能是由于收入不平等導(dǎo)致多數(shù)人撫養(yǎng)幼兒的成本負(fù)擔(dān)增加。人口老齡化既有來自于人口預(yù)期壽命增加引起的頂部老齡化,又有來自于生育率下降造成的底部老齡化[3];而底部老齡化問題同樣不容忽視。

    五、我國收入不平等與老齡化的脈沖響應(yīng)分析

    基于前面的計(jì)量模型,對應(yīng)的檢驗(yàn)結(jié)果表明收入不平等和老齡化之間存在顯著的非線性協(xié)整關(guān)系和非線性Granger因果關(guān)系。此前的計(jì)量模型主要是檢驗(yàn)二者在方程構(gòu)建方面是否存在時(shí)變系數(shù),但并未能形象地觀測出變量間隨時(shí)間的遞進(jìn)會(huì)表現(xiàn)出何種形式。為此,考慮到脈沖曲線分析法可以更直觀地觀測所考察的經(jīng)濟(jì)變量間的動(dòng)態(tài)影響路徑,本文通過構(gòu)建平滑遷移向量自回歸模型 (STVAR) 分析收入不平等和老齡化間的動(dòng)態(tài)影響效果。

    (一) 平滑遷移向量自回歸模型 (STVAR) 的構(gòu)建

    自Sims (1980)[7]提出向量自回歸模型后,該模型在分析研究復(fù)雜的經(jīng)濟(jì)問題上發(fā)揮了至關(guān)重要的作用。但隨著時(shí)間的推移,計(jì)量經(jīng)濟(jì)學(xué)家發(fā)現(xiàn)現(xiàn)實(shí)中的經(jīng)濟(jì)變量之間并不與線性向量自回歸模型中的假設(shè)一致;換言之,現(xiàn)實(shí)經(jīng)濟(jì)變量存在更復(fù)雜的傳導(dǎo)機(jī)制,并非是簡單的線性關(guān)系??紤]到上述線性向量自回歸模型的局限,參考Weise (1999)[8]構(gòu)建的邏輯平滑向量自回歸模型 (LSTVAR),分析收入不平等和老齡化之間的短期動(dòng)態(tài)非線性傳導(dǎo)機(jī)制。

    包含基尼系數(shù)一階差分序列、人口老齡化一階差分序列和實(shí)際產(chǎn)出同比增速序列作為內(nèi)生變量的p階滯后向量自回歸模型可以表示為:

    (7)

    其中,Xt=(x1t,x2t,x3t)′,x1t,x2t和x3t分別代表基尼系數(shù)一階差分序列、老齡化一階差分序列和實(shí)際產(chǎn)出同比增速序列。Γ0=(γ10,γ20,γ30)′為常數(shù)項(xiàng)對應(yīng)的系數(shù)列向量;Γi=(γi1,γi2,γi3)′,γij=(γij,1,γij,2,γij,3)′為第j個(gè)方程i階滯后內(nèi)生變量(Xt-i) 對應(yīng)的回歸系數(shù)向量;μt為隨機(jī)誤差列向量。由于內(nèi)生變量可能存在非線性數(shù)據(jù)生產(chǎn)機(jī)制,可以將式 (7) 表示的線性向量自回歸模型拓展為非線性平滑遷移向量自回歸模型的形式,其中非線性平滑向量自回歸模型中第j個(gè)方程可表示為:

    (8)

    其中,φ0j對應(yīng)的是非線性部分中常數(shù)項(xiàng)的回歸系數(shù),φij為非線性部分中內(nèi)生變量對應(yīng)的回歸系數(shù)列向量。Fj(·)表示轉(zhuǎn)移函數(shù),具體分為邏輯型轉(zhuǎn)移函數(shù)和指數(shù)型轉(zhuǎn)移函數(shù),其中邏輯型轉(zhuǎn)移函數(shù)形式為:

    F(zt;λ,c)={1+exp[-λ(zt-c)]}-1-1/2

    (9)

    其中,zt表示轉(zhuǎn)移變量,可以為內(nèi)生變量滯后項(xiàng)、外生變量和內(nèi)生變量的函數(shù);λ和c分別表示轉(zhuǎn)移斜率和位置參數(shù),轉(zhuǎn)移斜率越大,表示經(jīng)濟(jì)變量在不同區(qū)制的轉(zhuǎn)移速度越快,反之則反;特別地,當(dāng)轉(zhuǎn)移斜率趨于無窮時(shí),則轉(zhuǎn)移函數(shù)在位置參數(shù)附近表現(xiàn)為跳躍式的門檻模型,當(dāng)轉(zhuǎn)移斜率趨于零時(shí),則轉(zhuǎn)移函數(shù)值也同樣趨于零,此時(shí)的非線性模型退化為線性向量自回歸模型。繪制不同轉(zhuǎn)移斜率下轉(zhuǎn)移函數(shù)圖像如圖4所示。為簡單起見,將位置參數(shù)設(shè)置為0。換言之,線性向量自回歸模型和門檻向量自回歸模型是平滑遷移向量自回歸模型的特例。

    圖4不同轉(zhuǎn)移斜率對應(yīng)的轉(zhuǎn)移函數(shù)圖像

    從圖4中可以直觀地觀察到轉(zhuǎn)移斜率的變化對轉(zhuǎn)移函數(shù)圖像的影響,當(dāng)轉(zhuǎn)移斜率逐漸增大時(shí),其對應(yīng)的轉(zhuǎn)移函數(shù)在位置參數(shù)附近的轉(zhuǎn)變速度也變大。具體地,當(dāng)轉(zhuǎn)移斜率為0.5時(shí),轉(zhuǎn)移函數(shù)的函數(shù)值在0至1之間的變化速度最慢,當(dāng)轉(zhuǎn)移斜率增大到5時(shí),轉(zhuǎn)移函數(shù)對應(yīng)的函數(shù)值由0迅速變化到1。

    (二)非線性檢驗(yàn)

    對于平滑遷移向量自回歸模型的線性檢驗(yàn),由上面的分析可知,到轉(zhuǎn)移斜率λ=0時(shí)轉(zhuǎn)移函數(shù)對應(yīng)的函數(shù)值為0,此時(shí)平滑遷移向量自回歸模型退化為線性的向量自回歸模型,故可以將原假設(shè)設(shè)定為H0:λ=0,當(dāng)原假設(shè)成立時(shí),對應(yīng)的模型為線性向量自回歸,變量間具有線性關(guān)系;否則,拒絕原假設(shè)表示變量間具有非線性關(guān)系。我們注意到,當(dāng)原假設(shè)成立時(shí),式(8)中包含冗余的參數(shù)(即轉(zhuǎn)移斜率參數(shù)和位置參數(shù)在原假設(shè)成立的條件下是不可識別的),此時(shí)傳統(tǒng)的參數(shù)約束檢驗(yàn)統(tǒng)計(jì)量沒有漸進(jìn)的統(tǒng)計(jì)分別,故檢驗(yàn)方法失效。

    參考Weise(1999)[8]的做法,將式 (8) 在λ=0處進(jìn)行泰勒展開,得到線性檢驗(yàn)的輔助回歸方程:

    (10)

    (11)

    F分布類型的檢驗(yàn)統(tǒng)計(jì)量:

    (12)

    LR=T(log|Ω0|-log|Ω1|)~χ2[3m(mp+1)]

    (13)

    式 (13) 中,Ω0和Ω1分別表示平滑遷移向量自回歸模型對應(yīng)的輔助回歸方程的約束和無約束方差協(xié)方差矩陣。結(jié)合上述基于單方程的線性檢驗(yàn)和基于方程系統(tǒng)的線性檢驗(yàn)方法,包含測度收入不平等、人口老齡化和實(shí)際產(chǎn)出的三元平滑遷移向量自回歸模型的線性檢驗(yàn)結(jié)果如表4所示。

    表4 STVAR模型線性檢驗(yàn)結(jié)果

    注:方程1、方程2和方程3分別代表收入不平等方程、老齡化方程和實(shí)際產(chǎn)出方程;對于每一個(gè)轉(zhuǎn)移變量,第一行為檢驗(yàn)統(tǒng)計(jì)量,第二行為依據(jù)漸進(jìn)分布計(jì)算的顯著概率,第三行為利用自舉法 (bootstrap) 模擬的顯著性概率,抽樣次數(shù)設(shè)定為10 000次。

    由表4的模型線性檢驗(yàn)結(jié)果可以看出,當(dāng)老齡化一階差分滯后項(xiàng)和二階差分滯后項(xiàng)分別作為轉(zhuǎn)移變量時(shí),無論是基于單方程檢驗(yàn)還是方程系統(tǒng)檢驗(yàn),抑或是基于卡方檢驗(yàn)統(tǒng)計(jì)量和F檢驗(yàn)統(tǒng)計(jì)量,方程1、方程2和方程3均在1%顯著概率下拒絕模型為線性的原假設(shè),即方程中所包含的經(jīng)濟(jì)變量間存在顯著的非線性關(guān)聯(lián)機(jī)制。

    此外,由表4還可以看到,當(dāng)轉(zhuǎn)移變量為基尼系數(shù)一階滯后和二階滯后、實(shí)際產(chǎn)出同比增速的一階滯后或者二階滯后時(shí),模型同樣表現(xiàn)出非線性的特征。由于是基于老齡化對收入不平等的影響視角來分析變量間的短期動(dòng)態(tài)影響機(jī)制,因此選擇老齡化作為轉(zhuǎn)移變量,并依據(jù)老齡化的不同程度將研究樣本分為老齡化高區(qū)制和老齡化低區(qū)制,進(jìn)一步研究其對收入不平等是否具有顯著的非對稱影響,或者不同的影響程度。同時(shí)考慮到樣本數(shù)據(jù)是年度數(shù)據(jù)且樣本容量較少,故選擇老齡化一階差分滯后項(xiàng)作為轉(zhuǎn)移變量。后面的脈沖曲線分析同樣也將轉(zhuǎn)移變量設(shè)定為老齡化一階差分滯后項(xiàng),不再贅述。

    (三) 變量間脈沖響應(yīng)分析

    基于平滑遷移向量自回歸模型計(jì)算的脈沖響應(yīng)曲線與線性向量自回歸模型的脈沖曲線不同,對于線性模型而言,模型中的回歸系數(shù)在整個(gè)樣本區(qū)間為固定數(shù)值,而在非線性模型中,解釋變量前的回歸系數(shù)是非恒定在,如在不同的樣本區(qū)制內(nèi),同一解釋變量前的回歸系數(shù)大小不同或者符號迥異。此時(shí)采用線性模型計(jì)算脈沖相應(yīng)曲線有失偏頗,因此本文參考Koop (1996)[9]提出的計(jì)算非線性向量自回歸模型脈沖曲線的方法,通過計(jì)算估計(jì)出轉(zhuǎn)移函數(shù)中的轉(zhuǎn)移斜率和位置參數(shù),可以進(jìn)一步計(jì)算出在位置參數(shù)兩側(cè)的脈沖曲線函數(shù)。

    六、結(jié)論和政策建議

    人口老齡化是我國經(jīng)濟(jì)社會(huì)發(fā)展中必須直面的現(xiàn)實(shí)問題。從各國人口老齡化的歷史進(jìn)程上看,人口年齡結(jié)構(gòu)逐漸老化會(huì)給一個(gè)國家或地區(qū)的經(jīng)濟(jì)社會(huì)發(fā)展帶來一系列消極的影響。與此同時(shí),收入不平等會(huì)影響經(jīng)濟(jì)發(fā)展的驅(qū)動(dòng)力,并引發(fā)諸多社會(huì)問題,影響經(jīng)濟(jì)發(fā)展和社會(huì)穩(wěn)定。我國已經(jīng)進(jìn)入老齡化社會(huì),政府在制定政策改善收入不平等的同時(shí),應(yīng)該注意人口結(jié)構(gòu)的改變對收入不平等所產(chǎn)生的影響,加快完善我國老年人口的社會(huì)保障體系。

    猜你喜歡
    斷點(diǎn)協(xié)整人口老齡化
    一類無限可能問題的解法
    世界人口老齡化之住房問題
    英語文摘(2019年11期)2019-05-21 03:03:28
    技術(shù)創(chuàng)新視角下人口老齡化對經(jīng)濟(jì)增長的影響
    外商直接投資對我國進(jìn)出口貿(mào)易影響的協(xié)整分析
    河南金融發(fā)展和城鄉(xiāng)居民收入差距的協(xié)整分析
    主導(dǎo)電回路發(fā)生斷點(diǎn)故障判斷方法探討
    人口老齡化背景下的財(cái)政支出與經(jīng)濟(jì)增長
    內(nèi)蒙古地區(qū)人口老齡化問題研究
    中國商論(2016年33期)2016-03-01 01:59:44
    中國居民消費(fèi)與經(jīng)濟(jì)增長的協(xié)整關(guān)系檢驗(yàn)
    TKScope仿真調(diào)試Cortex-M3內(nèi)核的高級手段
    久久久久国内视频| 中文亚洲av片在线观看爽| 国产极品粉嫩免费观看在线| а√天堂www在线а√下载| 中文字幕人成人乱码亚洲影| 国产一区二区三区视频了| 欧美精品啪啪一区二区三区| 精品久久久久久久久久久久久 | a级毛片在线看网站| av有码第一页| 村上凉子中文字幕在线| 欧美三级亚洲精品| 搞女人的毛片| 中文字幕av电影在线播放| 正在播放国产对白刺激| 午夜老司机福利片| 久久精品91无色码中文字幕| 亚洲精品色激情综合| 亚洲国产精品sss在线观看| 欧美成人午夜精品| 看黄色毛片网站| 亚洲男人的天堂狠狠| 成年人黄色毛片网站| 天天添夜夜摸| 少妇被粗大的猛进出69影院| 午夜视频精品福利| 日韩精品中文字幕看吧| 午夜成年电影在线免费观看| tocl精华| 一夜夜www| 亚洲,欧美精品.| 精品人妻1区二区| 精品少妇一区二区三区视频日本电影| 九色国产91popny在线| 欧美色欧美亚洲另类二区| a在线观看视频网站| 免费高清在线观看日韩| 日本 av在线| 久久久久亚洲av毛片大全| 久热爱精品视频在线9| 中文资源天堂在线| 人人澡人人妻人| 国产精品野战在线观看| 99在线视频只有这里精品首页| 青草久久国产| 成人手机av| 国产欧美日韩一区二区精品| 午夜福利视频1000在线观看| 欧美乱妇无乱码| 成人午夜高清在线视频 | 性欧美人与动物交配| 国产真实乱freesex| tocl精华| 国产精华一区二区三区| 最近最新中文字幕大全免费视频| 欧美日韩中文字幕国产精品一区二区三区| 久久精品国产清高在天天线| 国产三级在线视频| 欧美又色又爽又黄视频| 少妇 在线观看| 亚洲va日本ⅴa欧美va伊人久久| 在线国产一区二区在线| 久久久久久久久中文| 日日爽夜夜爽网站| 啦啦啦观看免费观看视频高清| 国产亚洲av嫩草精品影院| 美女高潮到喷水免费观看| 日本精品一区二区三区蜜桃| av有码第一页| 美女高潮到喷水免费观看| 波多野结衣高清作品| 久久亚洲精品不卡| www日本黄色视频网| 久久精品夜夜夜夜夜久久蜜豆 | 国产高清videossex| 欧美一级毛片孕妇| 久久久久久亚洲精品国产蜜桃av| 丝袜人妻中文字幕| 麻豆av在线久日| 久久精品国产亚洲av高清一级| 国内毛片毛片毛片毛片毛片| 日韩视频一区二区在线观看| 欧美性猛交╳xxx乱大交人| 少妇裸体淫交视频免费看高清 | 夜夜看夜夜爽夜夜摸| 白带黄色成豆腐渣| 欧美激情久久久久久爽电影| 757午夜福利合集在线观看| 不卡av一区二区三区| 变态另类成人亚洲欧美熟女| 黑人巨大精品欧美一区二区mp4| 国产一级毛片七仙女欲春2 | 欧美激情高清一区二区三区| 变态另类丝袜制服| 精品国产美女av久久久久小说| 亚洲片人在线观看| 男人舔女人下体高潮全视频| 日韩 欧美 亚洲 中文字幕| 精品国产一区二区三区四区第35| 91成人精品电影| 亚洲自拍偷在线| 午夜福利高清视频| 一区二区三区高清视频在线| 大型黄色视频在线免费观看| 欧美又色又爽又黄视频| 国产一卡二卡三卡精品| 国产免费av片在线观看野外av| 国产av在哪里看| 午夜福利18| 欧美+亚洲+日韩+国产| 欧美+亚洲+日韩+国产| 久久狼人影院| 一区福利在线观看| 天堂√8在线中文| 国产精品av久久久久免费| 日本a在线网址| 成人午夜高清在线视频 | 黄片播放在线免费| 天堂影院成人在线观看| 亚洲av中文字字幕乱码综合 | 久久久久国产精品人妻aⅴ院| 亚洲国产日韩欧美精品在线观看 | 免费一级毛片在线播放高清视频| 日日摸夜夜添夜夜添小说| 成人精品一区二区免费| 人妻丰满熟妇av一区二区三区| 午夜久久久久精精品| 91大片在线观看| 久久热在线av| 高潮久久久久久久久久久不卡| 中文字幕最新亚洲高清| 亚洲国产日韩欧美精品在线观看 | 欧美国产精品va在线观看不卡| 亚洲一区二区三区不卡视频| 欧美国产精品va在线观看不卡| 免费无遮挡裸体视频| 国产视频一区二区在线看| 国产一级毛片七仙女欲春2 | 国内精品久久久久久久电影| 校园春色视频在线观看| 欧美另类亚洲清纯唯美| 亚洲av电影不卡..在线观看| 亚洲国产精品久久男人天堂| 亚洲国产精品久久男人天堂| 欧美日韩亚洲国产一区二区在线观看| 国产精品 欧美亚洲| www.精华液| 免费在线观看日本一区| 色综合婷婷激情| 亚洲国产日韩欧美精品在线观看 | 精品高清国产在线一区| 一级片免费观看大全| 18禁美女被吸乳视频| 日本 欧美在线| 91国产中文字幕| 国产成人一区二区三区免费视频网站| 免费在线观看亚洲国产| 亚洲成人久久爱视频| 91国产中文字幕| 亚洲五月婷婷丁香| 亚洲人成电影免费在线| 18美女黄网站色大片免费观看| 欧美日韩中文字幕国产精品一区二区三区| 久久久水蜜桃国产精品网| 老汉色av国产亚洲站长工具| 在线免费观看的www视频| 亚洲精品色激情综合| 久久久久久久午夜电影| 日韩视频一区二区在线观看| 免费观看人在逋| 91av网站免费观看| 精品免费久久久久久久清纯| 国产高清有码在线观看视频 | 精品不卡国产一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲av成人不卡在线观看播放网| 国产野战对白在线观看| 国内揄拍国产精品人妻在线 | 一个人观看的视频www高清免费观看 | 最近最新中文字幕大全免费视频| 成人午夜高清在线视频 | 婷婷亚洲欧美| 亚洲中文字幕一区二区三区有码在线看 | 亚洲 欧美 日韩 在线 免费| 国产爱豆传媒在线观看 | 中文亚洲av片在线观看爽| 亚洲熟女毛片儿| 国产精品香港三级国产av潘金莲| 男女做爰动态图高潮gif福利片| 在线视频色国产色| 午夜福利一区二区在线看| 99热只有精品国产| 国产黄a三级三级三级人| 精品国产乱码久久久久久男人| 亚洲国产精品成人综合色| 免费在线观看成人毛片| 淫妇啪啪啪对白视频| 老汉色∧v一级毛片| 国产精品一区二区三区四区久久 | 一级作爱视频免费观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲在线自拍视频| 女性被躁到高潮视频| 一进一出抽搐gif免费好疼| 日韩成人在线观看一区二区三区| 宅男免费午夜| 亚洲七黄色美女视频| 成人手机av| 久久久久久亚洲精品国产蜜桃av| 国内久久婷婷六月综合欲色啪| 欧美久久黑人一区二区| 一级毛片高清免费大全| 激情在线观看视频在线高清| 男人舔女人的私密视频| 亚洲第一av免费看| 国产黄a三级三级三级人| 丝袜美腿诱惑在线| 女警被强在线播放| 国产精品一区二区免费欧美| 变态另类成人亚洲欧美熟女| 国产在线观看jvid| 中文字幕人妻丝袜一区二区| 免费人成视频x8x8入口观看| 午夜老司机福利片| 91大片在线观看| 国产精品99久久99久久久不卡| 身体一侧抽搐| 国产在线观看jvid| 亚洲国产欧洲综合997久久, | 精品久久久久久成人av| 国产成人av激情在线播放| 成人永久免费在线观看视频| 黄网站色视频无遮挡免费观看| 999久久久国产精品视频| 巨乳人妻的诱惑在线观看| 岛国视频午夜一区免费看| 国产成人精品无人区| 欧美中文日本在线观看视频| 欧美最黄视频在线播放免费| 侵犯人妻中文字幕一二三四区| 好男人电影高清在线观看| 黄片小视频在线播放| 性色av乱码一区二区三区2| 真人一进一出gif抽搐免费| 老司机福利观看| 黑人操中国人逼视频| 国产一级毛片七仙女欲春2 | 欧美日韩亚洲国产一区二区在线观看| 怎么达到女性高潮| 成人特级黄色片久久久久久久| 亚洲欧美日韩无卡精品| 男女之事视频高清在线观看| 国产国语露脸激情在线看| 50天的宝宝边吃奶边哭怎么回事| 成年版毛片免费区| 久久国产乱子伦精品免费另类| 一二三四社区在线视频社区8| cao死你这个sao货| 国产精品亚洲av一区麻豆| 亚洲国产日韩欧美精品在线观看 | 在线十欧美十亚洲十日本专区| 欧美日韩黄片免| 18禁黄网站禁片免费观看直播| 深夜精品福利| 18禁美女被吸乳视频| 欧美色视频一区免费| 老司机午夜福利在线观看视频| 久久热在线av| 欧美日韩瑟瑟在线播放| 亚洲国产欧美日韩在线播放| av在线播放免费不卡| 哪里可以看免费的av片| 成人午夜高清在线视频 | 精品少妇一区二区三区视频日本电影| 老司机靠b影院| 国产主播在线观看一区二区| 久久久久久久精品吃奶| 97超级碰碰碰精品色视频在线观看| 久久久久久久久免费视频了| 手机成人av网站| 在线视频色国产色| 日本免费a在线| 久9热在线精品视频| 观看免费一级毛片| 精品国产超薄肉色丝袜足j| 日韩欧美一区二区三区在线观看| 亚洲一码二码三码区别大吗| 成人国语在线视频| 在线观看舔阴道视频| 自线自在国产av| a级毛片a级免费在线| 99精品久久久久人妻精品| 精品久久久久久久久久免费视频| 日本成人三级电影网站| 久久亚洲精品不卡| 男人的好看免费观看在线视频 | 国产激情久久老熟女| 国产免费男女视频| 婷婷精品国产亚洲av在线| 成人亚洲精品一区在线观看| 午夜福利视频1000在线观看| 观看免费一级毛片| 精品国产亚洲在线| 国产av一区二区精品久久| 午夜免费激情av| 国产精品永久免费网站| 精品福利观看| 国产一卡二卡三卡精品| 国产精品久久视频播放| 丰满人妻熟妇乱又伦精品不卡| 欧美绝顶高潮抽搐喷水| 亚洲第一av免费看| 精品久久久久久久末码| 欧美成人免费av一区二区三区| 国产亚洲精品一区二区www| 亚洲真实伦在线观看| 成人亚洲精品一区在线观看| 国产真实乱freesex| 伦理电影免费视频| 老司机在亚洲福利影院| 欧美av亚洲av综合av国产av| 久久久精品欧美日韩精品| 在线观看日韩欧美| 两个人看的免费小视频| 午夜激情av网站| 一a级毛片在线观看| 亚洲无线在线观看| 1024手机看黄色片| 我的亚洲天堂| 黄网站色视频无遮挡免费观看| 国产精品久久久av美女十八| 亚洲性夜色夜夜综合| 热99re8久久精品国产| 怎么达到女性高潮| 久久精品91蜜桃| 国产一卡二卡三卡精品| a级毛片在线看网站| 熟女少妇亚洲综合色aaa.| 两性午夜刺激爽爽歪歪视频在线观看 | 人成视频在线观看免费观看| 精品欧美一区二区三区在线| 日本熟妇午夜| 亚洲自拍偷在线| 亚洲av第一区精品v没综合| www国产在线视频色| 国产精品乱码一区二三区的特点| 色在线成人网| 久久久久久久精品吃奶| 99热这里只有精品一区 | 一级a爱片免费观看的视频| 精品久久久久久久末码| 精品高清国产在线一区| 国产主播在线观看一区二区| 1024视频免费在线观看| 久久人妻av系列| 亚洲国产高清在线一区二区三 | 国产一区二区三区在线臀色熟女| 国产成人一区二区三区免费视频网站| 日日爽夜夜爽网站| 午夜福利高清视频| 亚洲av电影在线进入| 人人妻人人澡人人看| 欧美日本视频| 久久 成人 亚洲| 男人舔奶头视频| 国产精品日韩av在线免费观看| 可以在线观看的亚洲视频| 精品高清国产在线一区| 男人舔奶头视频| 亚洲av熟女| 欧美国产日韩亚洲一区| 亚洲人成网站在线播放欧美日韩| 国产精品九九99| 可以免费在线观看a视频的电影网站| 熟妇人妻久久中文字幕3abv| 2021天堂中文幕一二区在线观 | 夜夜躁狠狠躁天天躁| ponron亚洲| 香蕉丝袜av| 在线播放国产精品三级| 色哟哟哟哟哟哟| 亚洲人成77777在线视频| 91麻豆av在线| 精品不卡国产一区二区三区| 国产午夜精品久久久久久| 老司机福利观看| 亚洲中文字幕日韩| 午夜福利免费观看在线| 中出人妻视频一区二区| 黄色成人免费大全| 岛国视频午夜一区免费看| 最近在线观看免费完整版| 亚洲av成人一区二区三| 18禁黄网站禁片免费观看直播| 精品熟女少妇八av免费久了| 精品不卡国产一区二区三区| 日本免费a在线| 亚洲av片天天在线观看| 国产精品,欧美在线| 国产单亲对白刺激| 一级a爱视频在线免费观看| 亚洲久久久国产精品| 欧美黄色淫秽网站| 亚洲专区国产一区二区| 亚洲精品粉嫩美女一区| 国产单亲对白刺激| 少妇 在线观看| 欧美性猛交黑人性爽| 精品久久蜜臀av无| 亚洲人成网站在线播放欧美日韩| or卡值多少钱| 国产亚洲欧美精品永久| 一区二区三区激情视频| 一个人免费在线观看的高清视频| svipshipincom国产片| 免费在线观看影片大全网站| 久久中文字幕一级| 一本精品99久久精品77| 国产av在哪里看| 大香蕉久久成人网| 热re99久久国产66热| 精品国产一区二区三区四区第35| 1024视频免费在线观看| 亚洲 欧美一区二区三区| 久久久久亚洲av毛片大全| 亚洲av熟女| 日本一本二区三区精品| 亚洲国产毛片av蜜桃av| 男人操女人黄网站| 一二三四在线观看免费中文在| 丰满的人妻完整版| 亚洲色图 男人天堂 中文字幕| 亚洲中文日韩欧美视频| 一级a爱片免费观看的视频| 亚洲熟妇熟女久久| 亚洲成人久久爱视频| 欧美乱色亚洲激情| 婷婷六月久久综合丁香| 青草久久国产| 9191精品国产免费久久| 深夜精品福利| 十分钟在线观看高清视频www| 老司机在亚洲福利影院| 免费高清视频大片| 成人亚洲精品av一区二区| 好男人电影高清在线观看| 桃红色精品国产亚洲av| 变态另类成人亚洲欧美熟女| av在线天堂中文字幕| 制服人妻中文乱码| 午夜福利在线在线| 99riav亚洲国产免费| 人妻丰满熟妇av一区二区三区| 亚洲美女黄片视频| 青草久久国产| 美女国产高潮福利片在线看| 中文字幕人妻丝袜一区二区| xxxwww97欧美| 亚洲熟女毛片儿| 亚洲人成77777在线视频| 国产av一区在线观看免费| 美女大奶头视频| 国产精品美女特级片免费视频播放器 | 国产久久久一区二区三区| 国产爱豆传媒在线观看 | 91av网站免费观看| 亚洲国产精品久久男人天堂| 嫩草影视91久久| 中国美女看黄片| 韩国av一区二区三区四区| 亚洲成av片中文字幕在线观看| 久久草成人影院| 精品电影一区二区在线| 日本 欧美在线| netflix在线观看网站| 久久精品国产清高在天天线| 精品午夜福利视频在线观看一区| 日本黄色视频三级网站网址| 日本免费a在线| 可以在线观看毛片的网站| 日本免费一区二区三区高清不卡| 日本免费a在线| 50天的宝宝边吃奶边哭怎么回事| 精品国产一区二区三区四区第35| 一级作爱视频免费观看| 午夜激情av网站| 久久久久久大精品| 观看免费一级毛片| 美女免费视频网站| 国产极品粉嫩免费观看在线| 午夜两性在线视频| 亚洲九九香蕉| 久久青草综合色| 在线视频色国产色| tocl精华| 侵犯人妻中文字幕一二三四区| 午夜福利免费观看在线| 91大片在线观看| 啦啦啦免费观看视频1| 女性生殖器流出的白浆| 国产成人啪精品午夜网站| 男女午夜视频在线观看| 亚洲人成电影免费在线| 一边摸一边做爽爽视频免费| 18禁美女被吸乳视频| 精品国产乱子伦一区二区三区| 亚洲精品久久成人aⅴ小说| 成年免费大片在线观看| 老熟妇仑乱视频hdxx| 日韩高清综合在线| 白带黄色成豆腐渣| 91av网站免费观看| 九色国产91popny在线| 午夜福利欧美成人| 精品无人区乱码1区二区| 999久久久精品免费观看国产| 亚洲av第一区精品v没综合| 91在线观看av| videosex国产| 黄色a级毛片大全视频| 在线观看一区二区三区| 亚洲成a人片在线一区二区| 身体一侧抽搐| 午夜福利在线观看吧| 国产亚洲精品一区二区www| 欧美激情极品国产一区二区三区| 亚洲国产中文字幕在线视频| 91成人精品电影| 国产精品美女特级片免费视频播放器 | 久久午夜综合久久蜜桃| 久久久国产欧美日韩av| tocl精华| 日本黄色视频三级网站网址| 亚洲欧美一区二区三区黑人| 色综合婷婷激情| 19禁男女啪啪无遮挡网站| 淫妇啪啪啪对白视频| 久久天堂一区二区三区四区| 精品久久蜜臀av无| av中文乱码字幕在线| 日本一区二区免费在线视频| 成人av一区二区三区在线看| tocl精华| 国产精品综合久久久久久久免费| 一区二区三区精品91| 黑人巨大精品欧美一区二区mp4| 欧美zozozo另类| 欧美激情久久久久久爽电影| 极品教师在线免费播放| 久久热在线av| 黄色女人牲交| 国内久久婷婷六月综合欲色啪| 国产又爽黄色视频| 久久久久久亚洲精品国产蜜桃av| ponron亚洲| 亚洲人成电影免费在线| 欧美成人午夜精品| 十八禁人妻一区二区| 亚洲欧美精品综合一区二区三区| 日本a在线网址| 757午夜福利合集在线观看| 黄频高清免费视频| 中文字幕精品免费在线观看视频| 91成年电影在线观看| 成人精品一区二区免费| 琪琪午夜伦伦电影理论片6080| 啦啦啦观看免费观看视频高清| 宅男免费午夜| 国产av一区在线观看免费| 久久亚洲精品不卡| 午夜福利成人在线免费观看| 美女大奶头视频| 狂野欧美激情性xxxx| 99re在线观看精品视频| 在线播放国产精品三级| 两个人免费观看高清视频| √禁漫天堂资源中文www| 亚洲第一青青草原| 久久午夜综合久久蜜桃| 久久99热这里只有精品18| 久久久久国产一级毛片高清牌| 久久欧美精品欧美久久欧美| 欧美丝袜亚洲另类 | 欧美不卡视频在线免费观看 | 变态另类丝袜制服| 中文字幕久久专区| 女生性感内裤真人,穿戴方法视频| 看免费av毛片| 国产精品九九99| 黄色毛片三级朝国网站| 亚洲熟妇熟女久久| 真人做人爱边吃奶动态| 亚洲性夜色夜夜综合| 熟妇人妻久久中文字幕3abv| 日本免费a在线| 曰老女人黄片| 18美女黄网站色大片免费观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精华国产精华精| 一级a爱视频在线免费观看| 中国美女看黄片| 亚洲 欧美 日韩 在线 免费| 制服诱惑二区| 成人av一区二区三区在线看| 在线观看免费午夜福利视频| 男人舔女人的私密视频| av片东京热男人的天堂| www日本黄色视频网| www.熟女人妻精品国产| 亚洲精品色激情综合| 99精品在免费线老司机午夜| 久久久久久国产a免费观看| 亚洲av成人一区二区三| 国产伦人伦偷精品视频| 免费无遮挡裸体视频|