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

    四種要素與經(jīng)濟(jì)增長的空間溢出效應(yīng)分析

    2016-09-16 12:19:20曹曉俊
    關(guān)鍵詞:效應(yīng)模型

    曹曉俊

    (安徽財(cái)經(jīng)大學(xué) 統(tǒng)計(jì)與應(yīng)用數(shù)學(xué)學(xué)院,安徽 蚌埠 233030)

    四種要素與經(jīng)濟(jì)增長的空間溢出效應(yīng)分析

    曹曉俊

    (安徽財(cái)經(jīng)大學(xué)統(tǒng)計(jì)與應(yīng)用數(shù)學(xué)學(xué)院,安徽蚌埠233030)

    通過對(duì)1992-2013年我國30個(gè)省區(qū)面板數(shù)據(jù)模型的空間相關(guān)性檢驗(yàn),表明經(jīng)濟(jì)增長具有空間關(guān)聯(lián)。進(jìn)一步通過模型檢驗(yàn),最后確定使用空間隨機(jī)效應(yīng)和時(shí)間固定效應(yīng)的空間杜賓模型來分析物資資本、人力資本、勞動(dòng)力和全要素生產(chǎn)率對(duì)經(jīng)濟(jì)增長是否存在空間溢出效應(yīng)。最終模型的間接效應(yīng)值表明,物質(zhì)資本的空間溢出效應(yīng)較少小,勞動(dòng)力和全要素生產(chǎn)率的空間溢出效應(yīng)較為顯著,人力資本的空間溢出效應(yīng)不顯著。

    物質(zhì)資本;人力資本;空間相關(guān)性;空間計(jì)量經(jīng)濟(jì)模型;空間溢出

    一、引言

    物質(zhì)資本、勞動(dòng)力、人力資本和全要素生產(chǎn)率對(duì)經(jīng)濟(jì)增長作用的實(shí)證已經(jīng)有許多文獻(xiàn)。但是大多數(shù)文獻(xiàn)在構(gòu)建模型時(shí),以各個(gè)經(jīng)濟(jì)體為相互獨(dú)立、不存在相互影響的個(gè)體作為前提,忽視了經(jīng)濟(jì)體之間的空間相關(guān)性和空間依賴性,從而忽略了物質(zhì)資本、勞動(dòng)力、人力資本和全要素生產(chǎn)率的空間溢出效應(yīng),這會(huì)導(dǎo)致模型結(jié)果與經(jīng)濟(jì)理論存在不一致性、模型缺乏穩(wěn)健性等問題。

    空間計(jì)量經(jīng)濟(jì)學(xué)的產(chǎn)生與發(fā)展為分析物資資本、勞動(dòng)力、人力資本和全要素生產(chǎn)率對(duì)經(jīng)濟(jì)增長是否存在空間溢出效應(yīng)提供了一種新的方法。肖志勇①肖志勇:《人力資本、空間溢出與經(jīng)濟(jì)增長——基于空間面板數(shù)據(jù)模型的經(jīng)驗(yàn)分析》,《財(cái)經(jīng)科學(xué)》2010年第3期,第61-68頁。在分析人力資本對(duì)經(jīng)濟(jì)增長的空間溢出效應(yīng)時(shí),使用教育年限法來測算人力資本指標(biāo),回歸結(jié)果表明人力資本對(duì)經(jīng)濟(jì)增長存在空間溢出效應(yīng)。魏下海②魏下海:《人力資本、空間溢出與省際全要素生產(chǎn)率增長——基于三種空間權(quán)重測度的實(shí)證檢驗(yàn)》,《財(cái)經(jīng)研究》2010年第12期,第94-104頁。、王文靜等③王文靜,劉彤,李盛基:《人力資本對(duì)我國全要素生產(chǎn)率增長作用的空間計(jì)量研究》,《經(jīng)濟(jì)與管理》2014年第2期,第22-28頁。人在分析人力資本對(duì)全要素生產(chǎn)率增長作用時(shí),也使用了教育年限法來測算人力資本指標(biāo),模型結(jié)果表明人力資本對(duì)全要素生產(chǎn)率增長具有空間溢出效應(yīng)。本文使用羅植和趙安平④朱平芳,徐大豐:《中國城市人力資本的估算》,《經(jīng)濟(jì)研究》2007年第9期,第84-95頁。提出的對(duì)朱平芳、徐大豐⑤郭慶旺,趙志耘,賈俊雪:《中國省份經(jīng)濟(jì)的全要素生產(chǎn)率分析》,《世界經(jīng)濟(jì)》2005年第5期,第46-53頁。的收入法進(jìn)行修正的方法來核算人力資本指標(biāo)。通過對(duì)面板數(shù)據(jù)模型的空間相關(guān)性進(jìn)行檢驗(yàn),得出所建立的模型存在空間相關(guān)性的結(jié)論。之后,建立空間計(jì)量經(jīng)濟(jì)模型,對(duì)其進(jìn)行檢驗(yàn),得出適合的空間面板數(shù)據(jù)計(jì)量模型。最后,利用得到的最終模型進(jìn)行物質(zhì)資本、勞動(dòng)力、人力資本和全要素生產(chǎn)率的空間溢出效應(yīng)分析。

    二、人力資本理論模型

    新增長理論突破了傳統(tǒng)增長理論的局限性,強(qiáng)調(diào)經(jīng)濟(jì)增長是經(jīng)濟(jì)系統(tǒng)的內(nèi)部力量運(yùn)行的結(jié)果,尤其是內(nèi)生技術(shù)進(jìn)步的產(chǎn)物,也重視知識(shí)或?qū)I(yè)化人力資本的遞增收益、外溢效應(yīng)及邊干邊學(xué)效應(yīng)。將人力資本引入到增長模型中以后,投資于資本積累中的資源變動(dòng)對(duì)產(chǎn)出的效應(yīng)將會(huì)得到明顯提高,使模型具備潛力,可用于解釋國家之間出現(xiàn)收入差異的原因。人力資本積累模型是物質(zhì)資本和人力資本共同積累的經(jīng)濟(jì)增長簡明模型。在人力資本積累模型的基礎(chǔ)之上進(jìn)行變動(dòng),將勞動(dòng)的有效性用全要素生產(chǎn)率進(jìn)行替代,模型兩邊取對(duì)數(shù),得到如下的模型:

    其中:1nYit表示第i個(gè)省區(qū)第t年的地區(qū)生產(chǎn)總值的對(duì)數(shù)值;1nKit表示第i個(gè)省區(qū)第t年的物質(zhì)資本存量的對(duì)數(shù)值;1nLit表示第i個(gè)省區(qū)第t年的就業(yè)人員數(shù)量的對(duì)數(shù)值;1nHit表示第i個(gè)省區(qū)第t年的人力資本存量的對(duì)數(shù)值;1nTFPit表示第i個(gè)省區(qū)第t年的全要素生產(chǎn)率增長率累積變動(dòng)率的對(duì)數(shù)值,代表技術(shù)進(jìn)步;εit表示第i個(gè)省區(qū)第t年的隨機(jī)誤差項(xiàng)。上述理論模型是沒有考慮各變量之間空間相關(guān)性的一種理論模型。

    三、數(shù)據(jù)說明

    本文的所有數(shù)據(jù)均來源于《新中國六十年統(tǒng)計(jì)資料匯編》《中國統(tǒng)計(jì)年鑒》、各省市往年統(tǒng)計(jì)年鑒以及統(tǒng)計(jì)局?jǐn)?shù)據(jù)庫中的歷年數(shù)據(jù)。本文使用了1991-2013年間30各省區(qū)的面板數(shù)據(jù),不包括西藏、臺(tái)灣、香港和澳門的樣本數(shù)據(jù)。本文選取的變量說明如下:

    (一)地區(qū)生產(chǎn)總值

    收集到的各省區(qū)的地區(qū)生產(chǎn)總值數(shù)據(jù)為當(dāng)年價(jià)的數(shù)據(jù),通過計(jì)算1991為100的各省區(qū)的地區(qū)生產(chǎn)總值指數(shù),將1991年各省區(qū)的當(dāng)年價(jià)的地區(qū)生產(chǎn)總值依次乘以1991-2013年的相應(yīng)的(1991= 100)地區(qū)生產(chǎn)總值指數(shù),即可得到1991-2013年各省區(qū)以1991年價(jià)格為基礎(chǔ)的不變價(jià)地區(qū)生產(chǎn)總值。

    (二)物質(zhì)資本存量

    對(duì)于物質(zhì)資本存量的計(jì)算,大多數(shù)學(xué)者使用的是永續(xù)盤存法。根據(jù)如下公式計(jì)算:

    其中,Kit表示第i個(gè)省區(qū)第t年的物質(zhì)資本存量;Kit-1表示第i個(gè)省區(qū)第t-1年的物質(zhì)資本存量;δit表示第i個(gè)省區(qū)第t年的物質(zhì)資本存量的折舊率;Iit表示第i個(gè)省區(qū)第t年的名義總投資;PIFit表示第i個(gè)省區(qū)第t年的固定資產(chǎn)投資價(jià)格指數(shù)(1991=100),則Iit/PIFit表示第i個(gè)省區(qū)第t年的實(shí)際總投資。

    在實(shí)際中,用永續(xù)盤存法計(jì)算物質(zhì)資本存量時(shí),需要考慮五個(gè)方面的問題:一,基期物質(zhì)資本存量的確定,即1991年的物質(zhì)資本存量的確定;二,折舊率的確定;三、名義總投資的替代指標(biāo)的確定;四、各省區(qū)固定資產(chǎn)投資價(jià)格指數(shù)的數(shù)據(jù);五、對(duì)缺失數(shù)據(jù)的處理。

    本文使用如下公式計(jì)算1991年各省區(qū)的物質(zhì)資本存量,即:

    其中,Ki1991表示第i個(gè)省區(qū)1991年的物質(zhì)資本存量;Ii1992表示第i個(gè)省區(qū)1992年的(1991年價(jià)格)實(shí)際總投資;δi1991表示第i個(gè)省區(qū)1991年的物質(zhì)資本存量的折舊率;δi1991表示第i個(gè)省區(qū)1991年的地區(qū)生產(chǎn)總值的增長率,將1991年各省區(qū)的(上年=100)地區(qū)生產(chǎn)總值指數(shù)減去100即可得到。

    對(duì)于物質(zhì)資本存量的折舊率,許多學(xué)者沒有統(tǒng)一的意見。本文使用10%做為1991-2013年各省區(qū)的物質(zhì)資本存量的折舊率。

    名義總投資的替代指標(biāo)主要有全社會(huì)固定資產(chǎn)投資、資本形成總額和固定資本形成總額,本文使用歷年的固定資本形成總額作為當(dāng)年名義總投資的數(shù)據(jù)。

    固定資產(chǎn)投資價(jià)格指數(shù)可以將名義總投資轉(zhuǎn)換為實(shí)際總投資,以便剔除價(jià)格因素的影響。我國從1991年開始公布官方的固定資產(chǎn)投資價(jià)格指數(shù)數(shù)據(jù),對(duì)于缺失固定資產(chǎn)投資價(jià)格指數(shù)的年份,本文使用商品零售價(jià)格指數(shù)來替代。收集到的固定資產(chǎn)投資價(jià)格指數(shù)數(shù)據(jù)為上一年為100的指數(shù),需要將其轉(zhuǎn)化為以1991年為100的固定資產(chǎn)投資價(jià)格指數(shù)。將名義總投資除以1991為100的固定資產(chǎn)投資價(jià)格指數(shù)即可得到相應(yīng)年份的實(shí)際總投資。

    對(duì)于缺失數(shù)據(jù)的處理。1991-1995年缺失重慶的固定資本形成總額數(shù)據(jù),由于缺失年份較少,本文使用1991-1995年重慶的全社會(huì)固定資產(chǎn)投資額來替代;天津缺失1991年的固定資產(chǎn)投資價(jià)格指數(shù),浙江缺失1991-1992年的固定資產(chǎn)投資價(jià)格指數(shù),廣東缺失1991-2000年的固定資產(chǎn)投資價(jià)格指數(shù),海南缺失1991-1999年的固定資產(chǎn)投資價(jià)格指數(shù),對(duì)這些省區(qū)缺失的固定資產(chǎn)投資價(jià)格指數(shù)數(shù)據(jù)使用相應(yīng)年份的商品零售價(jià)格指數(shù)的數(shù)據(jù)代替。

    (三)就業(yè)人員數(shù)量

    對(duì)于就業(yè)人員的數(shù)量,可以查找相關(guān)資料直接得到。本文直接采用得到的各省區(qū)的就業(yè)人員數(shù)量。

    (四)人力資本存量

    對(duì)于人力資本存量的估算,可以根據(jù)出發(fā)點(diǎn)的角度不同,主要有三類:第一類為投入法;第二類為教育年限法;第三類是收入法。這三種方法有各自的優(yōu)缺點(diǎn)。本文使用羅植和趙安平(2014)提出的改進(jìn)的收入法來測算各省區(qū)的人力資本存量。

    測算各省區(qū)的人力資本存量時(shí)涉及到的主要公式如下:

    其中:Wit(1)*表示第i個(gè)省區(qū)第t年的單位人力資本的效率工資;Kit表示第i個(gè)省區(qū)第t年的物質(zhì)資本存量;Lit表示第i個(gè)省區(qū)第t年的勞動(dòng)力數(shù)量;hit表示第i個(gè)省區(qū)第t年的人均人力資本存量;Wit(hit)表示第i個(gè)省區(qū)第t年的勞動(dòng)者的平均工資。

    在計(jì)算單位人力資本的效率工資時(shí),假設(shè)a'= 0.60,β'=0.35,a'與β'的和不為零,主要是因?yàn)榱_植和趙安平(2014)所使用的生產(chǎn)函數(shù)中包括了難以度量的綜合制度變量。具體人均人力資本存量的推導(dǎo)過程可以參考羅植和趙安平的文章。在計(jì)算人均人力資本存量時(shí),由于各省區(qū)的全社會(huì)就業(yè)人員的平均工資的數(shù)據(jù)難以獲取,所以選擇各省區(qū)的在崗職工的平均工資,與之對(duì)應(yīng)的勞動(dòng)力數(shù)量指標(biāo)選擇各省區(qū)的在崗職工數(shù)量。各省區(qū)的在崗職工人數(shù)和在崗職工的平均工資可以通過搜索相關(guān)資料得到。其中,各省區(qū)的在崗職工的平均工資為當(dāng)年價(jià)格的工資水平,包含價(jià)格因素的影響,使用各省區(qū)1991 為100的地區(qū)生產(chǎn)總值縮減指數(shù)來計(jì)算各省區(qū)在崗職工的實(shí)際平均工資,結(jié)果剔除了價(jià)格因素的影響。將各省的人均人力資本存量乘以其在崗職工人數(shù),即得到各省的人力資本存量。

    (五)全要素生產(chǎn)率增長率累積變動(dòng)率

    全要素生產(chǎn)率增長率的測算方法較多,大致可分為兩類:第一類是增長會(huì)計(jì)法;第二類是經(jīng)濟(jì)計(jì)量法。其中,增長會(huì)計(jì)法的計(jì)算過程的假設(shè)條件較為苛刻,計(jì)算結(jié)果準(zhǔn)確性較差。本文使用經(jīng)濟(jì)計(jì)量法中的非參數(shù)數(shù)據(jù)包絡(luò)分析方法中的Malmquist指數(shù)法來計(jì)算全要素生產(chǎn)率增長率。

    與增長會(huì)計(jì)法相比而言,Malmquist指數(shù)來估算全要素生產(chǎn)率增長率避免了較為苛刻的理論假設(shè),而且可以將全要素生產(chǎn)率增長率分解為效率變化和技術(shù)進(jìn)步率兩個(gè)組成部分,得到的全要素生產(chǎn)率增長率的信息較為充分,可以做更為詳細(xì)的現(xiàn)實(shí)情況分析。

    在利用Malmquist指數(shù)法時(shí),選擇的投入變量為物質(zhì)資本存量和就業(yè)人員數(shù)量,選擇的產(chǎn)出變量為地區(qū)生產(chǎn)總值。本文使用的是30個(gè)省區(qū)的1991-2013年的面板數(shù)據(jù),Malmquist指數(shù)法得出的結(jié)果只有30個(gè)省區(qū)的1992-2013年的全要素生產(chǎn)率增長率數(shù)據(jù)。將得到的30個(gè)省區(qū)的1992-2013年全要素生產(chǎn)率增長率進(jìn)行累乘,得到30個(gè)省區(qū)1991-2013年共22年全要素生產(chǎn)率增長率累積變動(dòng)率。本文把全要素生產(chǎn)率增長率累積變動(dòng)率作為全要素生產(chǎn)率的替代指標(biāo)。

    本文為了保證數(shù)據(jù)時(shí)期的一致性,選擇1992-2013年30個(gè)省區(qū)的面板數(shù)據(jù)進(jìn)行建模分析,1992-2013年主要變量的描述統(tǒng)計(jì)見表1。

    表1 1992-2013年主要變量的描述統(tǒng)計(jì)

    三、模型介紹和實(shí)證分析

    (一)空間計(jì)量經(jīng)濟(jì)模型

    在區(qū)域經(jīng)濟(jì)分析中,由于空間依賴的存在違背了大多數(shù)古典統(tǒng)計(jì)和計(jì)量經(jīng)濟(jì)學(xué)分析中的樣本相互獨(dú)立的基本假設(shè)前提,因此使用古典計(jì)量經(jīng)濟(jì)學(xué)方法并不能獲取數(shù)據(jù)的空間依賴性,導(dǎo)致模型建立產(chǎn)生各種問題。使用空間計(jì)量經(jīng)濟(jì)學(xué)的方法來處理空間數(shù)據(jù),具有古典計(jì)量經(jīng)濟(jì)學(xué)所無法比擬的優(yōu)點(diǎn)。空間計(jì)量經(jīng)濟(jì)學(xué)是以各種空間依賴和空間相互作用現(xiàn)象作為主要的研究對(duì)象,截面數(shù)據(jù)和面板數(shù)據(jù)的空間計(jì)量模型構(gòu)成了空間計(jì)量經(jīng)濟(jì)學(xué)中組織各種模擬方法的框架。比較常用的空間計(jì)量經(jīng)濟(jì)模型主要有:空間滯后模型、空間誤差模型和空間杜賓模型。

    空間滯后模型(SLM)又稱空間自回歸模型:

    其中:Y是N×1階被解釋變量向量,W是N×N階空間權(quán)重矩陣,LN是N×1階單位向量,X是N×K階解釋變量矩陣,ε是N×1階干擾項(xiàng)向量,α是待估的截距項(xiàng)參數(shù),β是待估的解釋變量系數(shù)參數(shù),θ 是K×1階待估的參數(shù)向量,δ是空間自回歸系數(shù),λ是空間自相關(guān)系數(shù)。WY反映被解釋變量之間存在的內(nèi)生交互效應(yīng),即如果有兩個(gè)鄰近的個(gè)體,分別為A和B,個(gè)體A的被解釋變量變化會(huì)對(duì)鄰近個(gè)體B的被解釋變量產(chǎn)生影響,同時(shí)個(gè)體B的被解釋變量變化也會(huì)對(duì)鄰近個(gè)體A的被解釋變量產(chǎn)生影響。WX反映解釋變量之間存在的外生交互效應(yīng),即個(gè)體A的外生解釋變量變化對(duì)鄰近個(gè)體B的外生解釋變量產(chǎn)生影響,同時(shí)個(gè)體B的外生解釋變量變化會(huì)對(duì)鄰近個(gè)體A的外生解釋變量產(chǎn)生影響。Wu反映干擾項(xiàng)之間存在的交互效應(yīng),個(gè)體A的遺漏或不可觀測因素的變化對(duì)鄰近個(gè)體B的遺漏或不可觀測因素造成影響,同時(shí)個(gè)體B的遺漏或不可觀測因素的變化對(duì)鄰近個(gè)體A的遺漏或不可觀測因素產(chǎn)生影響。

    當(dāng)θ等于零向量時(shí),空間杜賓模型就變成空間滯后模型;當(dāng)θ=-δβ時(shí),空間杜賓模型轉(zhuǎn)化成空間誤差模型,此時(shí)λ=δ。由此可見,空間滯后模型和空間誤差模型是空間杜賓模型的特例,空間杜賓模型更具有普遍意義的空間計(jì)量經(jīng)濟(jì)模型。

    對(duì)于空間權(quán)重矩陣的構(gòu)造,主要有三種常用方法:第一種是0-1空間權(quán)重矩陣;第二種是基于地理距離的空間權(quán)重矩陣;第三種是基于經(jīng)濟(jì)距離的空間權(quán)重矩陣。本文使用的空間權(quán)重矩陣為0-1空間權(quán)重矩陣,即當(dāng)省區(qū)i與省區(qū)j相鄰,則Wij為1;否則,Wij為0。其中,在構(gòu)造0-1空間權(quán)重矩陣時(shí),將海南與廣東設(shè)為相鄰,取值為1。

    (二)空間相關(guān)性檢驗(yàn)

    檢驗(yàn)變量是否存在空間相關(guān)性的較為常用的方法主要有Moran I指數(shù),其計(jì)算公式如下:

    其中:S2表示方差;Yi表示第i個(gè)省區(qū)的觀測值;n表示省區(qū)總數(shù);Wij表示空間權(quán)重矩陣。Moran I的計(jì)算結(jié)果取值范圍為[-1,1]。

    上述Moran I指數(shù)適合于檢驗(yàn)截面數(shù)據(jù)模型的空間相關(guān)性,并不適用于面板數(shù)據(jù)模型的相關(guān)性檢驗(yàn)。但是,可以根據(jù)截面數(shù)據(jù)Moran I指數(shù)檢驗(yàn)的原理,將Moran I指數(shù)用于面板數(shù)據(jù)模型的空間相關(guān)性檢驗(yàn),公式如下:

    其中:WNT=IT○×W,W是空間權(quán)重矩陣;e是普通面板數(shù)據(jù)模型(不考慮空間相關(guān)性)的OLS估計(jì)殘差。

    在建立空間面板數(shù)據(jù)模型后,需要為所建立的模型進(jìn)行檢驗(yàn)以便確定使用空間滯后模型還是空間誤差模型或者空間杜賓模型較為合適。運(yùn)用LM檢驗(yàn)時(shí),一方面可以使用LM-lag檢驗(yàn)?zāi)P褪遣淮嬖诳臻g滯后被解釋變量的原假設(shè),另一方面可用LM-erro檢驗(yàn)?zāi)P褪遣淮嬖诳臻g自相關(guān)誤差項(xiàng)的原假設(shè)。在一定的置信水平上,若LM-lag或LM-erro檢驗(yàn)統(tǒng)計(jì)量值的結(jié)果說明拒絕原假設(shè),則表明采用空間滯后模型或空間誤差模型較為適合;否則,說明空間滯后被解釋變量沒有空間相關(guān)性或誤差項(xiàng)不存在空間相關(guān)性,使用普通面板數(shù)據(jù)模型即可。

    (三)空間相關(guān)性檢驗(yàn)結(jié)果

    本文先通過建立無空間交互效應(yīng)的面板數(shù)據(jù)模型,然后對(duì)其空間相關(guān)性進(jìn)行檢驗(yàn)。模型參數(shù)估計(jì)結(jié)果見表2。

    表2 無空間交互效應(yīng)的面板數(shù)據(jù)模型的參數(shù)估計(jì)結(jié)果

    根據(jù)以上的四個(gè)面板數(shù)據(jù)模型回歸結(jié)果,可以看出四個(gè)模型的參數(shù)估計(jì)值都在1%的顯著水平下顯著,并且除截距項(xiàng)外其他參數(shù)符號(hào)都為正,符合經(jīng)濟(jì)理論。在考慮模型存在空間固定效應(yīng)時(shí),模型的明顯增大,說明存在空間固定效應(yīng)的面板數(shù)據(jù)模型的擬合效果較好。Moran I指數(shù)的值都表明模型的變量應(yīng)該存在顯著的空間相關(guān)性。

    具體的空間相關(guān)性的形式可以依據(jù)LM-lag檢驗(yàn)和穩(wěn)健LM-lag檢驗(yàn)、LM-error檢驗(yàn)和穩(wěn)健LM-error檢驗(yàn)的結(jié)果來判斷。四個(gè)模型的LM-lag檢驗(yàn)在1%顯著性水平下都顯著拒絕模型不存在空間滯后被解釋變量的原假設(shè),穩(wěn)健LM-lag檢驗(yàn)只有存在空間和時(shí)間固定效應(yīng)的模型四沒有拒絕模型不存在空間滯后被解釋變量的原假設(shè),其他三個(gè)模型在1%顯著性水平下都顯著拒絕模型不存在空間滯后被解釋變量的原假設(shè);LM-error檢驗(yàn)在1%顯著性水平下都顯著拒絕模型不存在空間誤差項(xiàng)的原假設(shè),穩(wěn)健LM-error檢驗(yàn)只有存在空間固定效應(yīng)的模型二沒有拒絕模型不存在空間誤差項(xiàng)的原假設(shè),其他三個(gè)模型在1%顯著性水平下都拒絕了模型不存在空間誤差項(xiàng)的原假設(shè)。

    所以,空間相關(guān)性檢驗(yàn)結(jié)果說明模型存在空間相關(guān)性,可以使用空間計(jì)量經(jīng)濟(jì)模型來擬合模型,這樣有助使模型改善,以便更好地進(jìn)行經(jīng)濟(jì)理論檢驗(yàn)和實(shí)際情況分析。

    (四)空間面板數(shù)據(jù)模型的選擇

    由上述空間相關(guān)性形式的檢驗(yàn)結(jié)果可知,存在空間固定效應(yīng)的模型適合用空間滯后模型,存在空間和時(shí)間固定效應(yīng)的模型適合使用空間滯后模型。為了確定所建立的空間計(jì)量經(jīng)濟(jì)模型是否存在空間固定效應(yīng)的聯(lián)合顯著和時(shí)間固定效應(yīng)的聯(lián)合顯著,可以進(jìn)行似然比LR檢驗(yàn)。LR檢驗(yàn)的結(jié)果見表3。

    表3 LR檢驗(yàn)空間固定效應(yīng)的聯(lián)合顯著性和時(shí)間固定效應(yīng)的聯(lián)合顯著性結(jié)果

    通過將模型三的極大似然值與模型四的極大似然值相減再乘以-2,即得空間固定效應(yīng)的聯(lián)合顯著性檢驗(yàn)的LR值,結(jié)果表明在1%的顯著性水平下顯著拒絕空間固定效應(yīng)的聯(lián)合不顯著的原假設(shè)。同樣,將模型二的極大似然值減去模型四的極大似然值后再乘以-2,可得時(shí)間固定效應(yīng)的聯(lián)合顯著性檢驗(yàn)的LR值,檢驗(yàn)結(jié)果表明在1%的顯著性水平下顯著拒絕時(shí)間固定效應(yīng)的聯(lián)合不顯著的原假設(shè)。結(jié)合這兩個(gè)檢驗(yàn)結(jié)果,說明可以使用具有空間固定效應(yīng)和時(shí)間固定效應(yīng)的雙向固定效應(yīng)模型。

    本文先分別建立面板數(shù)據(jù)的空間和時(shí)間固定效應(yīng)的空間杜賓模型和空間隨機(jī)效應(yīng)且時(shí)間固定效應(yīng)的空間杜賓模型,然后對(duì)模型進(jìn)行空間滯后檢驗(yàn)和空間誤差檢驗(yàn),進(jìn)一步判斷模型的具體形式。面板數(shù)據(jù)的空間杜賓模型的估計(jì)結(jié)果見表4。

    根據(jù)模型的回歸結(jié)果可知,在都有時(shí)間固定效應(yīng)的情況下,有空間隨機(jī)效應(yīng)的空間杜賓模型的調(diào)整的可決系數(shù)要比有空間固定效應(yīng)的空間杜賓模型的調(diào)整的可決系數(shù)大的多,說明空間隨機(jī)效應(yīng)的空間杜賓模型擬合的效果較好。兩個(gè)模型的lnK、lnL、lnH和lnTFP的參數(shù)估計(jì)值都為正,只有l(wèi)nH的參數(shù)估計(jì)值不顯著,其余的參數(shù)估計(jì)值在1%的顯著性水平下顯著,說明目前人力資本對(duì)經(jīng)濟(jì)的增長影響效果不是很顯著,物質(zhì)資本、勞動(dòng)力及全要素生產(chǎn)率對(duì)經(jīng)濟(jì)增長的影響比較顯著。由此可見,當(dāng)前,我國人力資本水平對(duì)推動(dòng)經(jīng)濟(jì)增長的作用表現(xiàn)還不明顯,而資本、勞動(dòng)力和技術(shù)進(jìn)步對(duì)促進(jìn)經(jīng)濟(jì)增長的作用表現(xiàn)比較明顯。這從側(cè)面說明我國經(jīng)濟(jì)發(fā)展仍然依靠大量的資金投入和勞動(dòng)力投入,還沒有充分發(fā)揮人力資本的推動(dòng)作用,這樣會(huì)造成資源的浪費(fèi)和使用效率低下等問題。

    兩個(gè)模型的空間滯后解釋變量中只有資本存量的空間滯后變量系數(shù)的參數(shù)估計(jì)值在1%的顯著性水平下顯著不為零,說明資本存量存在外生交互效應(yīng),即一個(gè)省區(qū)的資本存量變化會(huì)影響鄰近省區(qū)的資本存量,鄰近省區(qū)資本存量的變動(dòng)也會(huì)對(duì)該省區(qū)的資本存量產(chǎn)生影響,而其他的空間滯后解釋變量系數(shù)的參數(shù)估計(jì)值不顯著,則說明勞動(dòng)力、人力資本存量和全要素生產(chǎn)率的外生交互影響較小。此外,兩個(gè)模型的空間滯后被解釋變量即地區(qū)生產(chǎn)總值的空間滯后變量,其系數(shù)估計(jì)值在1%顯著性水平下顯著不為零,表明地區(qū)生產(chǎn)總值存在內(nèi)生交互效應(yīng),即一個(gè)省區(qū)的地區(qū)生產(chǎn)總值發(fā)生變化會(huì)對(duì)鄰近省區(qū)的地區(qū)生產(chǎn)總值產(chǎn)生影響,鄰近省區(qū)的地區(qū)生產(chǎn)總值的變化同樣也會(huì)對(duì)該省區(qū)的地區(qū)生產(chǎn)總值產(chǎn)生變動(dòng)。

    為了進(jìn)一步確定模型的具體形式,對(duì)上述模型進(jìn)行空間滯后和空間誤差的Wald檢驗(yàn)。檢驗(yàn)結(jié)果見表5。

    表5 空間滯后和空間誤差的Wald檢驗(yàn)

    由空間滯后和空間誤差的Wald檢驗(yàn)結(jié)果可知,這兩個(gè)空間杜賓模型在1%的顯著性水平下都顯著拒絕了模型可以簡化為空間滯后模型和空間誤差模型的原假設(shè)。由此可見,可以拒絕空間滯后模型和空間誤差模型而使用空間杜賓模型。

    到此為止,可以確定使用空間杜賓模型來擬合模型,但還需要進(jìn)行Hausman檢驗(yàn)以確定使用空間固定效應(yīng)和時(shí)間固定效應(yīng)模型還是空間隨機(jī)效應(yīng)和時(shí)間固定效應(yīng)模型。Hausman檢驗(yàn)的結(jié)果見表6。

    表6Hausman檢驗(yàn)結(jié)果

    由Hausman檢驗(yàn)的結(jié)果看出,不能拒絕模型存在隨機(jī)效應(yīng)的原假設(shè),所以最終模型為空間隨機(jī)效應(yīng)和時(shí)間固定效應(yīng)的空間杜賓模型。

    (五)空間溢出效應(yīng)分析。

    表7 空間隨機(jī)效應(yīng)和時(shí)間固定效應(yīng)的空間杜賓模型的效應(yīng)估計(jì)值

    為了進(jìn)行空間溢出效應(yīng)分析,估計(jì)出了空間隨機(jī)效應(yīng)和時(shí)間固定效應(yīng)的空間杜賓模型的解釋變量的直接效應(yīng)值、間接效應(yīng)值和總效應(yīng)值以及對(duì)應(yīng)的檢驗(yàn)統(tǒng)計(jì)量和p值,估計(jì)結(jié)果見表7。

    其中,解釋變量的總效應(yīng)值等于其直接效應(yīng)值與間接效應(yīng)值之和,由于四舍五入的原因,造成計(jì)算結(jié)果存在一定的偏差。解釋變量的直接效應(yīng)值表示本省區(qū)的解釋變量變動(dòng)1%,使本省區(qū)的被解釋變量發(fā)生的變化等于該直接效應(yīng)值,解釋變量的直接效應(yīng)就相當(dāng)于解釋變量的彈性。解釋變量的間接效應(yīng)值表示本省區(qū)的解釋變量變動(dòng)1%,使其鄰近省區(qū)的被解釋變量發(fā)生的變化等于該間接效應(yīng)值,間接效應(yīng)值即空間溢出效應(yīng)值。

    檢驗(yàn)解釋變量是否存在空間溢出效應(yīng)主要根據(jù)解釋變量的間接效應(yīng)值來判斷,如果間接效應(yīng)值在一定的顯著性水平下顯著不為零,那么可以認(rèn)為該解釋變量存在空間溢出效應(yīng)。各個(gè)解釋變量的間接效應(yīng)值的顯著性檢驗(yàn)結(jié)果如下:

    資本存量的間接效應(yīng)值小于零,其p值比10%稍大一點(diǎn),可以認(rèn)為資本存量存在一定的空間溢出效應(yīng),但不是很顯著。一個(gè)省區(qū)的物質(zhì)資本存量增加1%,會(huì)使本省區(qū)的地區(qū)生產(chǎn)總值增加30.8%,但也會(huì)使鄰近省區(qū)的物質(zhì)資本存量減少4.2%。假定其他要素投入相同,在物質(zhì)資本存量總量一定的情況下,某個(gè)省區(qū)的物質(zhì)資本存量使用量增加會(huì)使該省區(qū)的地區(qū)生產(chǎn)總值提高,但是這會(huì)使其他鄰近省區(qū)的物質(zhì)資本存量得不到增加,造成鄰近地區(qū)經(jīng)濟(jì)增長得不到充足的物質(zhì)資本存量保證,結(jié)果就是這些省區(qū)的地區(qū)生產(chǎn)總值下降,經(jīng)濟(jì)下滑。

    勞動(dòng)力的間接效應(yīng)值在10%的顯著性水平下顯著不為零,說明勞動(dòng)力存在空間溢出效應(yīng),那么一個(gè)省區(qū)的勞動(dòng)力數(shù)量增加1%,會(huì)使本省區(qū)的地區(qū)生產(chǎn)總值增加17.5%,也會(huì)使鄰近省區(qū)的地區(qū)生產(chǎn)總值提高11.5%,本省區(qū)的地區(qū)生產(chǎn)總值增加幅度比鄰近省區(qū)的地區(qū)生產(chǎn)總值的增加幅度要大。勞動(dòng)力是地區(qū)經(jīng)濟(jì)發(fā)展不可或缺的要素,經(jīng)濟(jì)發(fā)達(dá)地區(qū)會(huì)吸引較多的勞動(dòng)力涌入,經(jīng)濟(jì)發(fā)達(dá)地區(qū)經(jīng)濟(jì)增長的同時(shí)也會(huì)帶動(dòng)鄰近地區(qū)相關(guān)產(chǎn)業(yè),如建筑業(yè)、餐飲、住宿等第二三產(chǎn)業(yè)的發(fā)展,使得鄰近地區(qū)的經(jīng)濟(jì)也會(huì)增長,但是經(jīng)濟(jì)發(fā)達(dá)地區(qū)的經(jīng)濟(jì)增長要比鄰近地區(qū)要多大得多,鄰近地區(qū)之間仍然存在較大的發(fā)展差距。

    人力資本的間接效應(yīng)值比較小且不顯著,說明人力資本不存在空間溢出效應(yīng)。這與肖志勇得出的人力資本存在空間溢出效應(yīng)的結(jié)論相反,主要原因是人力資本的測算方法存在差異,肖志勇使用的是教育年限法,對(duì)不同的教育階段賦權(quán),計(jì)算出平均受教育年限,來表示人力資本,而本文使用的是收入法,從收入角度來衡量人力資本。我國人力資本各地區(qū)分布不均,經(jīng)濟(jì)發(fā)達(dá)地區(qū)的居民有相對(duì)較高的收入,有更多的進(jìn)行教育、培訓(xùn)等繼續(xù)深造的機(jī)會(huì),對(duì)人力資本的投入較多,而經(jīng)濟(jì)相對(duì)落后地區(qū)的居民收入水平較低,而且無法利用這些機(jī)會(huì),對(duì)人力資本的投入較少??傮w而言,我國的人力資本水平仍然還需要進(jìn)一步提高,這樣才能使人力資本得到充分利用,從而帶動(dòng)經(jīng)濟(jì)增長。

    全要素生產(chǎn)率的間接效應(yīng)值大于零且在5%的顯著性水平下顯著,可見全要素生產(chǎn)率存在空間溢出效應(yīng)。那么一個(gè)省區(qū)的全要素生產(chǎn)率提高1%,本省區(qū)的地區(qū)生產(chǎn)總值會(huì)增加31.4%,與此同時(shí),鄰近省區(qū)的地區(qū)生產(chǎn)總值會(huì)提高7.7%,本省區(qū)的地區(qū)生產(chǎn)總值增加比重是鄰近省區(qū)增加比重的四倍多。地區(qū)技術(shù)發(fā)展程度與該地區(qū)的經(jīng)濟(jì)發(fā)展程度密切相關(guān)。我國地區(qū)之間的經(jīng)濟(jì)發(fā)展水平仍然不均,經(jīng)濟(jì)發(fā)達(dá)地區(qū)有進(jìn)行技術(shù)研發(fā)的資金支持和智力保障,鄰近的經(jīng)濟(jì)落后地區(qū)一方面缺乏資金,另一方面也缺乏技術(shù)人才,這種雙方面的資源缺乏造成經(jīng)濟(jì)落后地區(qū)的技術(shù)進(jìn)步較為緩慢,技術(shù)發(fā)展較為艱難的局面。但是,技術(shù)落后地區(qū)可以通過引進(jìn)或?qū)W習(xí)鄰近技術(shù)發(fā)達(dá)地區(qū)的技術(shù)等途徑獲得先進(jìn)的技術(shù)來發(fā)展本區(qū)域的經(jīng)濟(jì),帶動(dòng)自身經(jīng)濟(jì)增長。區(qū)域間的經(jīng)濟(jì)發(fā)展是相互聯(lián)動(dòng)的,一方面雙方彼此競爭,另一方面雙方彼此協(xié)助。為了達(dá)到共贏的目的,雙方可以密切加強(qiáng)技術(shù)研發(fā)合作與交流,以合作為主,競爭為輔,在合作與競爭中謀求共同發(fā)展。

    五、結(jié)論與啟示

    本文通過對(duì)面板數(shù)據(jù)模型的空間相關(guān)性進(jìn)行檢驗(yàn),確定面板數(shù)據(jù)模型應(yīng)該存在空間相關(guān)性。由此通過使用空間計(jì)量經(jīng)濟(jì)模型來進(jìn)行實(shí)證分析。通過模型檢驗(yàn),最終選擇了空間隨機(jī)效應(yīng)和時(shí)間固定效應(yīng)的空間杜賓模型。利用該模型進(jìn)行空間溢出分析,得出如下結(jié)論:物質(zhì)資本存量存在一定為負(fù)的空間溢出效應(yīng)但較小,勞動(dòng)力和全要素生產(chǎn)率存在較為顯著的空間溢出效應(yīng),勞動(dòng)力的空間溢出效應(yīng)大于全要素生產(chǎn)率的空間溢出效應(yīng),人力資本存量的空間溢出效應(yīng)不顯著。

    通過以上要素的空間溢出效應(yīng)分析,說明目前我國經(jīng)濟(jì)增長中勞動(dòng)力和技術(shù)進(jìn)步具有外部經(jīng)濟(jì)效應(yīng),本區(qū)域內(nèi)的勞動(dòng)力的增加和技術(shù)進(jìn)步,不僅有助于本區(qū)域內(nèi)的經(jīng)濟(jì)增長,而且會(huì)給鄰近區(qū)域的經(jīng)濟(jì)增長帶來促進(jìn)作用;而物質(zhì)資本存量具有外部不經(jīng)濟(jì)效應(yīng),本區(qū)域內(nèi)的物質(zhì)資本存量增加,只會(huì)有助于本區(qū)域的經(jīng)濟(jì)增長,而對(duì)鄰近區(qū)域的經(jīng)濟(jì)增長帶來負(fù)面效應(yīng),使鄰近區(qū)域經(jīng)濟(jì)發(fā)展減退;人力資本存量的外部經(jīng)濟(jì)效應(yīng)不是很明顯,表明人力資本對(duì)我國區(qū)域之間的經(jīng)濟(jì)增長的促進(jìn)效應(yīng)不足,我國人力資本水平還有待提高。我們應(yīng)該努力提高勞動(dòng)力、人力資本和全要素生產(chǎn)率的正向空間外溢效應(yīng),增加經(jīng)濟(jì)增長的有利因素,降低物質(zhì)資本的負(fù)向空間外溢效應(yīng),減少經(jīng)濟(jì)增長的不利因素。這樣才能進(jìn)一步推動(dòng)我國經(jīng)濟(jì)持續(xù)穩(wěn)定地增長。

    另外,區(qū)域之間的發(fā)展離不開競爭與合作。一個(gè)區(qū)域?yàn)榱吮緟^(qū)域的利益,會(huì)侵害其他鄰近區(qū)域的利益,造成區(qū)域之間發(fā)展的不均衡。但是,區(qū)域發(fā)展需要彼此之間的合作,通過合作,增加要素流動(dòng)性,提高要素的空間外溢程度,可以帶動(dòng)整體利益的提高,達(dá)到共同發(fā)展的目標(biāo)。

    Analysis of the spatial spillover effect of four factors and economic growth

    CAO Xiaojun

    This paper inspected the spatial correlation of the panel data model for 30 provinces and regions from 1992 to 2013,and found that economic growth had spatial association.Testing models were then used to decide the use of Durbin model that has spatial random effect and time fixed effect to analyze the spatial spillover effect among physical capital,human capital,labor and total factor productivity with economic growth.Finally,the indirect effect values of the model showed that the spatial spillover effect of physical capital was small,the spatial spillover effect of labor and total factor productivity was significant,and the spatial spillover effect of human capital was not significant.

    physical capital;human capital;spatial correlation;spatial econometric model;spatial spillover

    F064.1

    A

    1009-9530(2016)04-0045-07

    2016-03-20

    曹曉俊(1990-),男,安徽財(cái)經(jīng)大學(xué)統(tǒng)計(jì)學(xué)專業(yè)碩士研究生,導(dǎo)師葉安寧副教授;主要研究方向:經(jīng)濟(jì)統(tǒng)計(jì)。

    猜你喜歡
    效應(yīng)模型
    一半模型
    鈾對(duì)大型溞的急性毒性效應(yīng)
    懶馬效應(yīng)
    場景效應(yīng)
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    應(yīng)變效應(yīng)及其應(yīng)用
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    偶像效應(yīng)
    中文字幕免费在线视频6| 精品欧美国产一区二区三| 久久久久久久精品吃奶| 国产主播在线观看一区二区| 国产野战对白在线观看| 永久网站在线| 亚洲人成伊人成综合网2020| 一本一本综合久久| 真人一进一出gif抽搐免费| 亚洲精品影视一区二区三区av| 久久久精品欧美日韩精品| 男女下面进入的视频免费午夜| 1000部很黄的大片| 国产精品亚洲一级av第二区| 午夜福利18| 欧美在线黄色| 久久久久亚洲av毛片大全| 无遮挡黄片免费观看| 国产蜜桃级精品一区二区三区| 国产三级黄色录像| 国模一区二区三区四区视频| 久久久国产成人精品二区| 精品午夜福利视频在线观看一区| 精品一区二区三区av网在线观看| 看免费av毛片| 内射极品少妇av片p| 久久精品91蜜桃| 色在线成人网| 欧美一区二区国产精品久久精品| 亚洲精品影视一区二区三区av| 无人区码免费观看不卡| 日韩大尺度精品在线看网址| 三级男女做爰猛烈吃奶摸视频| 此物有八面人人有两片| 国产精品一区二区三区四区久久| 91久久精品国产一区二区成人| 精品免费久久久久久久清纯| 亚洲五月天丁香| 久久婷婷人人爽人人干人人爱| 首页视频小说图片口味搜索| 欧美黄色淫秽网站| 又爽又黄a免费视频| 久久精品国产99精品国产亚洲性色| 亚洲狠狠婷婷综合久久图片| 99久久久亚洲精品蜜臀av| 久久久久久久久久成人| 亚洲午夜理论影院| 久久精品国产亚洲av涩爱 | 黄色日韩在线| 免费在线观看亚洲国产| 天堂√8在线中文| 免费av毛片视频| 少妇高潮的动态图| 男女那种视频在线观看| 日本一二三区视频观看| 俺也久久电影网| 欧美在线一区亚洲| 变态另类成人亚洲欧美熟女| 亚洲国产精品久久男人天堂| 日韩欧美三级三区| 在线天堂最新版资源| 日韩欧美在线乱码| 两人在一起打扑克的视频| 露出奶头的视频| 最后的刺客免费高清国语| 中文字幕精品亚洲无线码一区| 亚洲精品乱码久久久v下载方式| 夜夜躁狠狠躁天天躁| 国产av一区在线观看免费| 波多野结衣巨乳人妻| 淫妇啪啪啪对白视频| 激情在线观看视频在线高清| av国产免费在线观看| 听说在线观看完整版免费高清| 真人一进一出gif抽搐免费| 欧洲精品卡2卡3卡4卡5卡区| 麻豆成人av在线观看| 男女下面进入的视频免费午夜| 国产精品久久久久久久电影| 国产精品免费一区二区三区在线| 久久久久国内视频| 日本三级黄在线观看| 一个人看视频在线观看www免费| 成人国产一区最新在线观看| 黄色视频,在线免费观看| 日韩 亚洲 欧美在线| 久久久精品大字幕| 欧美一区二区精品小视频在线| 丰满人妻一区二区三区视频av| 久久亚洲真实| av在线老鸭窝| 别揉我奶头 嗯啊视频| 成年版毛片免费区| 美女高潮喷水抽搐中文字幕| 俺也久久电影网| 久久精品国产亚洲av天美| 婷婷亚洲欧美| 国产精品一及| 日韩 亚洲 欧美在线| 国产在线男女| 欧美潮喷喷水| 91字幕亚洲| 国产一区二区亚洲精品在线观看| 熟女电影av网| 精品午夜福利在线看| 色精品久久人妻99蜜桃| 成人美女网站在线观看视频| 97碰自拍视频| 91在线精品国自产拍蜜月| 国产成人欧美在线观看| 一个人看的www免费观看视频| 一区福利在线观看| 国产亚洲精品久久久com| av在线老鸭窝| 亚洲精品在线观看二区| 岛国在线免费视频观看| 久久国产精品人妻蜜桃| 91久久精品电影网| 欧美性猛交╳xxx乱大交人| 免费一级毛片在线播放高清视频| 动漫黄色视频在线观看| 午夜福利免费观看在线| 精品日产1卡2卡| 国产精品国产高清国产av| 国内揄拍国产精品人妻在线| 国产亚洲av嫩草精品影院| 久久国产精品人妻蜜桃| 国产爱豆传媒在线观看| 听说在线观看完整版免费高清| 女人被狂操c到高潮| 99国产精品一区二区三区| 日韩中字成人| 嫩草影院入口| 国产探花在线观看一区二区| 国产野战对白在线观看| 国产精品久久久久久人妻精品电影| 亚洲人成网站在线播| 国产乱人伦免费视频| 亚洲综合色惰| 精品人妻熟女av久视频| 国产亚洲精品综合一区在线观看| 亚洲成人中文字幕在线播放| 俺也久久电影网| 国内精品久久久久精免费| 在线a可以看的网站| 高清日韩中文字幕在线| 最近在线观看免费完整版| 脱女人内裤的视频| 国产精品影院久久| 天美传媒精品一区二区| 老女人水多毛片| 日本黄大片高清| 国产精品日韩av在线免费观看| 非洲黑人性xxxx精品又粗又长| 亚洲国产日韩欧美精品在线观看| 亚洲欧美日韩无卡精品| 天美传媒精品一区二区| 国产精品一区二区免费欧美| 亚洲精品在线观看二区| 一区二区三区免费毛片| 成人欧美大片| 精品久久久久久久久亚洲 | 看黄色毛片网站| 在线观看av片永久免费下载| 亚洲性夜色夜夜综合| 中出人妻视频一区二区| 久久精品91蜜桃| 亚洲精品粉嫩美女一区| 又爽又黄无遮挡网站| 婷婷亚洲欧美| 国产亚洲欧美在线一区二区| 欧洲精品卡2卡3卡4卡5卡区| 好男人电影高清在线观看| 天美传媒精品一区二区| 高潮久久久久久久久久久不卡| 欧美+亚洲+日韩+国产| 日本免费a在线| 亚洲不卡免费看| 国产一区二区三区视频了| 免费高清视频大片| 亚洲精品456在线播放app | 午夜a级毛片| 嫁个100分男人电影在线观看| 免费人成在线观看视频色| 日本五十路高清| 女人十人毛片免费观看3o分钟| av黄色大香蕉| 午夜影院日韩av| 日韩av在线大香蕉| 精品一区二区三区av网在线观看| 午夜福利18| 国产精品久久久久久人妻精品电影| 日本一本二区三区精品| 亚洲 国产 在线| 午夜亚洲福利在线播放| 国产精品一区二区免费欧美| 亚洲成人中文字幕在线播放| 成人性生交大片免费视频hd| 欧美不卡视频在线免费观看| 久久人人爽人人爽人人片va | av在线蜜桃| 最近视频中文字幕2019在线8| 精品午夜福利视频在线观看一区| 人妻夜夜爽99麻豆av| 极品教师在线免费播放| 精品久久久久久久久久免费视频| 婷婷色综合大香蕉| 亚洲人成电影免费在线| h日本视频在线播放| 亚洲av成人av| 校园春色视频在线观看| 日本免费一区二区三区高清不卡| 婷婷精品国产亚洲av| 久久久久久久久大av| 在现免费观看毛片| 国产免费男女视频| 长腿黑丝高跟| 蜜桃久久精品国产亚洲av| 国内精品美女久久久久久| 欧美性猛交黑人性爽| 欧美中文日本在线观看视频| 国产高潮美女av| 国产高清有码在线观看视频| 看黄色毛片网站| 热99在线观看视频| 亚洲欧美激情综合另类| 国产精品久久久久久久电影| 亚洲欧美日韩卡通动漫| 赤兔流量卡办理| 日韩中字成人| 国产麻豆成人av免费视频| 又黄又爽又免费观看的视频| 欧美性猛交黑人性爽| 久久伊人香网站| 99久久久亚洲精品蜜臀av| 深爱激情五月婷婷| 国内精品久久久久精免费| 亚洲国产精品sss在线观看| 日韩亚洲欧美综合| 国内精品一区二区在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产又黄又爽又无遮挡在线| 欧美bdsm另类| 无人区码免费观看不卡| 欧美区成人在线视频| 99国产综合亚洲精品| 中文在线观看免费www的网站| 日本五十路高清| 美女cb高潮喷水在线观看| 老司机深夜福利视频在线观看| 美女黄网站色视频| 亚洲av电影不卡..在线观看| 神马国产精品三级电影在线观看| 精品欧美国产一区二区三| 国产精品亚洲美女久久久| 在线免费观看的www视频| 一级av片app| 一区二区三区高清视频在线| 亚洲av免费高清在线观看| 少妇熟女aⅴ在线视频| 少妇的逼好多水| 蜜桃久久精品国产亚洲av| av天堂在线播放| 免费观看的影片在线观看| 亚洲欧美日韩东京热| or卡值多少钱| 亚洲在线观看片| eeuss影院久久| 日韩 亚洲 欧美在线| 国产亚洲精品久久久久久毛片| 深爱激情五月婷婷| 18+在线观看网站| 国产精品免费一区二区三区在线| 丰满人妻熟妇乱又伦精品不卡| 欧美日本亚洲视频在线播放| 狂野欧美白嫩少妇大欣赏| 中文亚洲av片在线观看爽| 国产伦人伦偷精品视频| 国产亚洲精品综合一区在线观看| 精品久久久久久久久av| 国产色爽女视频免费观看| 18禁黄网站禁片免费观看直播| 无人区码免费观看不卡| ponron亚洲| 最后的刺客免费高清国语| 亚洲va日本ⅴa欧美va伊人久久| 女生性感内裤真人,穿戴方法视频| 中文在线观看免费www的网站| 午夜福利欧美成人| 国产69精品久久久久777片| 免费在线观看影片大全网站| 青草久久国产| 国产黄片美女视频| 亚洲黑人精品在线| 少妇熟女aⅴ在线视频| 我要看日韩黄色一级片| 国内精品久久久久精免费| 他把我摸到了高潮在线观看| 亚洲美女搞黄在线观看 | 国产黄片美女视频| 亚洲国产精品成人综合色| 欧美+日韩+精品| 久久久久久国产a免费观看| 亚洲av成人精品一区久久| 男女下面进入的视频免费午夜| 国产视频内射| 2021天堂中文幕一二区在线观| 别揉我奶头~嗯~啊~动态视频| 成人特级黄色片久久久久久久| 男女视频在线观看网站免费| 国产毛片a区久久久久| 国产精品自产拍在线观看55亚洲| 赤兔流量卡办理| 长腿黑丝高跟| 成人一区二区视频在线观看| 亚洲熟妇中文字幕五十中出| 午夜精品一区二区三区免费看| 国产精品久久久久久精品电影| 免费看日本二区| 999久久久精品免费观看国产| 色综合站精品国产| 欧洲精品卡2卡3卡4卡5卡区| av国产免费在线观看| x7x7x7水蜜桃| 国产乱人视频| 国产精品美女特级片免费视频播放器| 91久久精品电影网| 国产人妻一区二区三区在| 亚洲中文字幕日韩| 美女cb高潮喷水在线观看| av福利片在线观看| 欧美成人一区二区免费高清观看| 三级男女做爰猛烈吃奶摸视频| а√天堂www在线а√下载| 好看av亚洲va欧美ⅴa在| 久久久久久久午夜电影| 色综合婷婷激情| 国产一区二区亚洲精品在线观看| 精品久久久久久久人妻蜜臀av| 日本精品一区二区三区蜜桃| 日韩大尺度精品在线看网址| 国产成人福利小说| 美女被艹到高潮喷水动态| 草草在线视频免费看| 99热这里只有是精品在线观看 | 日本撒尿小便嘘嘘汇集6| 欧美+亚洲+日韩+国产| 精品熟女少妇八av免费久了| 午夜福利高清视频| 可以在线观看的亚洲视频| 国产免费一级a男人的天堂| 男人的好看免费观看在线视频| .国产精品久久| 毛片一级片免费看久久久久 | 999久久久精品免费观看国产| 在线观看午夜福利视频| 18+在线观看网站| 婷婷六月久久综合丁香| 国产久久久一区二区三区| 久久精品久久久久久噜噜老黄 | 亚洲国产精品合色在线| 亚洲黑人精品在线| 亚洲无线观看免费| 最近中文字幕高清免费大全6 | 成人特级黄色片久久久久久久| 18美女黄网站色大片免费观看| 少妇高潮的动态图| 老鸭窝网址在线观看| 欧美激情久久久久久爽电影| а√天堂www在线а√下载| 男女床上黄色一级片免费看| 国产欧美日韩一区二区精品| av在线蜜桃| 最近在线观看免费完整版| 久久精品91蜜桃| 给我免费播放毛片高清在线观看| 悠悠久久av| 亚洲成av人片免费观看| 黄色一级大片看看| a级毛片免费高清观看在线播放| 国产探花在线观看一区二区| 日韩欧美三级三区| 亚洲精品456在线播放app | 亚洲无线在线观看| 国产午夜精品论理片| 亚洲性夜色夜夜综合| 91狼人影院| av在线天堂中文字幕| 久久久久免费精品人妻一区二区| 国产伦人伦偷精品视频| 白带黄色成豆腐渣| 日韩欧美精品免费久久 | 亚洲国产色片| 久久精品人妻少妇| 国产亚洲欧美98| 国产真实伦视频高清在线观看 | 老司机深夜福利视频在线观看| 最近中文字幕高清免费大全6 | av在线蜜桃| 成人性生交大片免费视频hd| 亚洲精品亚洲一区二区| 最新在线观看一区二区三区| 最新中文字幕久久久久| 老女人水多毛片| 久久久久久久久久黄片| 精品国内亚洲2022精品成人| 噜噜噜噜噜久久久久久91| 亚洲国产日韩欧美精品在线观看| 日日干狠狠操夜夜爽| 国产白丝娇喘喷水9色精品| 天天躁日日操中文字幕| 全区人妻精品视频| 好男人在线观看高清免费视频| 久久久久久九九精品二区国产| 能在线免费观看的黄片| 午夜福利18| 人妻久久中文字幕网| 99国产精品一区二区三区| 色哟哟·www| 少妇的逼水好多| 综合色av麻豆| 国产探花极品一区二区| 国产精品久久视频播放| 国产免费一级a男人的天堂| bbb黄色大片| 亚洲 国产 在线| 国产亚洲欧美98| 在线观看av片永久免费下载| 欧美3d第一页| 亚洲片人在线观看| 老司机午夜福利在线观看视频| 欧美bdsm另类| 观看免费一级毛片| 国产精品三级大全| 内射极品少妇av片p| 色播亚洲综合网| 亚洲专区国产一区二区| 在现免费观看毛片| 午夜免费成人在线视频| 91麻豆av在线| 一进一出好大好爽视频| av专区在线播放| 舔av片在线| 熟妇人妻久久中文字幕3abv| 天天躁日日操中文字幕| av在线蜜桃| 国内精品久久久久久久电影| 18禁黄网站禁片免费观看直播| 直男gayav资源| 久久性视频一级片| 白带黄色成豆腐渣| 亚洲成av人片免费观看| 国产精品影院久久| 精品99又大又爽又粗少妇毛片 | 欧美黑人巨大hd| 九九在线视频观看精品| 麻豆国产av国片精品| 久久久久亚洲av毛片大全| 又爽又黄无遮挡网站| 热99在线观看视频| 脱女人内裤的视频| 午夜亚洲福利在线播放| 99视频精品全部免费 在线| 在现免费观看毛片| 国产野战对白在线观看| 久久久久精品国产欧美久久久| 国产欧美日韩精品亚洲av| 在线播放无遮挡| а√天堂www在线а√下载| 精品不卡国产一区二区三区| 别揉我奶头~嗯~啊~动态视频| 自拍偷自拍亚洲精品老妇| 人妻久久中文字幕网| 黄片小视频在线播放| 欧美潮喷喷水| 91在线精品国自产拍蜜月| 久久香蕉精品热| 国产淫片久久久久久久久 | 久99久视频精品免费| 欧美日本亚洲视频在线播放| 黄色丝袜av网址大全| 麻豆国产97在线/欧美| 特级一级黄色大片| 18禁裸乳无遮挡免费网站照片| 在线观看av片永久免费下载| 亚洲av不卡在线观看| 亚洲av美国av| 欧美高清成人免费视频www| 岛国在线免费视频观看| 国内精品美女久久久久久| 两人在一起打扑克的视频| 国产高清激情床上av| 高潮久久久久久久久久久不卡| 欧美性感艳星| 女同久久另类99精品国产91| 在线免费观看不下载黄p国产 | 中亚洲国语对白在线视频| 成人av一区二区三区在线看| 极品教师在线视频| 精品午夜福利在线看| 国产午夜精品论理片| 欧美xxxx黑人xx丫x性爽| 一进一出抽搐gif免费好疼| 成年女人毛片免费观看观看9| 色综合亚洲欧美另类图片| 天堂av国产一区二区熟女人妻| 亚洲av第一区精品v没综合| 日本黄色视频三级网站网址| 一级黄片播放器| avwww免费| 亚洲欧美日韩高清在线视频| 国产黄色小视频在线观看| 国产精品亚洲一级av第二区| 亚洲欧美激情综合另类| 国产欧美日韩一区二区三| 日日摸夜夜添夜夜添小说| 久久6这里有精品| 精品一区二区三区人妻视频| 欧美中文日本在线观看视频| 88av欧美| 国产探花极品一区二区| 久久午夜亚洲精品久久| 国产真实伦视频高清在线观看 | 三级国产精品欧美在线观看| 性欧美人与动物交配| 美女黄网站色视频| 国产亚洲精品综合一区在线观看| 国模一区二区三区四区视频| 久久精品人妻少妇| 丰满人妻熟妇乱又伦精品不卡| 亚洲av日韩精品久久久久久密| 精品福利观看| 国产成人aa在线观看| 欧美bdsm另类| 免费一级毛片在线播放高清视频| 99久久无色码亚洲精品果冻| 久久精品人妻少妇| 高清在线国产一区| 我要看日韩黄色一级片| 亚洲成人免费电影在线观看| xxxwww97欧美| 欧美黄色片欧美黄色片| 老司机深夜福利视频在线观看| 亚洲最大成人手机在线| 亚洲自拍偷在线| 精品久久久久久久久久免费视频| 一夜夜www| 亚洲一区二区三区色噜噜| 性欧美人与动物交配| 美女高潮喷水抽搐中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 国产精品嫩草影院av在线观看 | 美女高潮喷水抽搐中文字幕| 久久精品国产亚洲av涩爱 | 国产真实伦视频高清在线观看 | 一进一出抽搐gif免费好疼| 欧美日韩亚洲国产一区二区在线观看| 一二三四社区在线视频社区8| 人妻制服诱惑在线中文字幕| 国产欧美日韩一区二区精品| 男人舔女人下体高潮全视频| 男女那种视频在线观看| 亚洲七黄色美女视频| ponron亚洲| 国内久久婷婷六月综合欲色啪| 最近视频中文字幕2019在线8| 全区人妻精品视频| 久久精品国产99精品国产亚洲性色| 久久久久久久久久黄片| 国产成+人综合+亚洲专区| 毛片女人毛片| 国产不卡一卡二| 赤兔流量卡办理| 中亚洲国语对白在线视频| 国产高清视频在线播放一区| 亚洲av电影不卡..在线观看| 色综合站精品国产| 亚洲人成网站在线播放欧美日韩| 一级av片app| www.色视频.com| 国内精品久久久久久久电影| 久久中文看片网| 亚洲,欧美,日韩| 99在线视频只有这里精品首页| 中文字幕人成人乱码亚洲影| 中文字幕av在线有码专区| 桃红色精品国产亚洲av| 五月伊人婷婷丁香| 波多野结衣巨乳人妻| a在线观看视频网站| 一级a爱片免费观看的视频| 99久久精品一区二区三区| 午夜福利在线观看吧| av中文乱码字幕在线| 午夜精品一区二区三区免费看| 免费人成视频x8x8入口观看| 欧美一区二区精品小视频在线| 欧洲精品卡2卡3卡4卡5卡区| 99国产综合亚洲精品| 免费在线观看影片大全网站| 99久国产av精品| 九九热线精品视视频播放| 一区二区三区激情视频| 亚洲第一欧美日韩一区二区三区| 色综合站精品国产| 日本黄色视频三级网站网址| 欧美日韩亚洲国产一区二区在线观看| 噜噜噜噜噜久久久久久91| 亚洲一区二区三区色噜噜| 欧美国产日韩亚洲一区| 观看美女的网站| 老熟妇仑乱视频hdxx| 如何舔出高潮| 丰满乱子伦码专区| 久久精品91蜜桃|