李 原

（山西財經(jīng)大學(xué) 統(tǒng)計學(xué)院，太原 030006）

0 引言

上世紀(jì)80年代以來，我國人口流動規(guī)模越來越大。“六普”數(shù)據(jù)顯示，2010年中國大陸流動人口已超2.6億人，占總?cè)丝诘?6%以上，國家統(tǒng)計局調(diào)查數(shù)據(jù)顯示，2015年流動人口達(dá)2.92億。大規(guī)模的人口流動直接影響區(qū)域人口規(guī)模結(jié)構(gòu)和經(jīng)濟格局，進(jìn)而影響全國經(jīng)濟發(fā)展的結(jié)構(gòu)變遷、產(chǎn)業(yè)升級和人口城鎮(zhèn)化進(jìn)程，有關(guān)人口流動的研究在我國變得越來越重要。

圍繞我國人口流動模式問題，國內(nèi)很多學(xué)者進(jìn)行了相關(guān)探討，相關(guān)文獻(xiàn)的定量研究方法可分為三類：人口遷移指標(biāo)的描述統(tǒng)計分析法、多元回歸模型分析法及空間計量分析法。已有研究大都基于普查數(shù)據(jù)采用經(jīng)典統(tǒng)計方法分析我國人口流動問題。經(jīng)典統(tǒng)計視未知參數(shù)為固定常數(shù)，基于大樣本對總體進(jìn)行推斷，估計結(jié)果高度依賴樣本信息，遇到有偏樣本則會導(dǎo)致參數(shù)估計結(jié)果有偏,特別是對于人口流動空間數(shù)據(jù)而言，其同時具有小樣本和自相關(guān)的特性，基于大樣本推斷總體的傳統(tǒng)統(tǒng)計模型由于自身的局限受到更多挑戰(zhàn)。

而貝葉斯層次時空模型是貝葉斯層次模型和時空交互模型的結(jié)合，貝葉斯層次模型可在一定程度上克服小樣本缺陷，時空交互模型同時考慮空間相關(guān)和隨機效應(yīng)，貝葉斯層次時空模型為時空數(shù)據(jù)統(tǒng)計分析提供了全新的思路和方法，應(yīng)用該模型進(jìn)行人口流動分析，能解決經(jīng)典統(tǒng)計難以解決的問題，同時能克服以往人口流動分析對時間和空間交互作用的忽略。

實際上，貝葉斯層次時空分析法近年在國際上受到廣泛的重視，特別是在生態(tài)學(xué)、疾病傳播等方面的理論和方法探索越來越多，目前國內(nèi)也出現(xiàn)了相關(guān)的應(yīng)用，但該方法在人口流動特征分析中的應(yīng)用還未見。本文首次應(yīng)用貝葉斯層次時空模型，對我國2009—2014年間省際人口流入的時空演化特征進(jìn)行分析。

1 Bayesian層次時空模型

貝葉斯層次時空模型是貝葉斯層次模型和時空交互模型的結(jié)合，貝葉斯層次模型主要包括三個子模型：

θj為參數(shù)，φ是超參數(shù)，P(φ)是超先驗。P(θj,φ|yi)是參數(shù)和超參數(shù)的后驗分布，參數(shù)和超參數(shù)的后驗分布可通過似然函數(shù)、先驗分布和超先驗分布的乘積獲得：

后驗分布的計算涉及高維聯(lián)合概率密度函數(shù)估計，目前主要采用馬爾科夫鏈-蒙特卡洛法實現(xiàn)。為保證模型收斂性，可采用兩條MCMC鏈，收斂性通過Gelman-Rubin統(tǒng)計參數(shù)估計，其值越接近于1收斂性越好。

空間自相關(guān)性是空間數(shù)據(jù)的一大特點，不滿足獨立同分布特征，樣本數(shù)據(jù)有限時難以達(dá)到大樣本抽樣。而貝葉斯層次時空模型是貝葉斯層次模型與時空交互模型的結(jié)合，可以克服上述缺陷，該模型通過充分利用先驗信息，考慮時空交互項，同時考慮了數(shù)據(jù)的空間自相關(guān)性和隨機不確定性，以概率分布的形式描述時空特征，相應(yīng)的三個子模型為：

yit是時空觀測樣本值，θit是時空因變量，s(i)和α(t)分別代表穩(wěn)態(tài)空間相對風(fēng)險和總體變化趨勢，Ωit(i,t)為時空交互項，該項常用的表達(dá)形式為，代表研究時段內(nèi)各研究子區(qū)域的局部變化趨勢，bi度量從總體變化趨勢中分解出的局部變化程度，若bi＞0，說明其局部變化強度強于總體變化強度，bi＜0說明其局部變化強度弱于總體變化強度，t*=t-tm表示相對于研究期內(nèi)中間時點tm的時間跨度。εit是隨機噪聲。

計算所有參數(shù)后驗分布后，可根據(jù)各區(qū)域時空因變量的空間相對風(fēng)險s(i)的高低將其分為熱點、溫點和冷點區(qū)。該分類方法基于Richardson等（2004）提出的原則：若地區(qū)空間相對風(fēng)險大于1的后驗概率大于0.8,則屬于熱點區(qū)，若小于0.2則屬于冷點區(qū)，界于0.2和0.8之間為溫點區(qū)。與此類似，根據(jù)各區(qū)域觀測值局部變化趨勢bi可將其分為強變化、穩(wěn)定態(tài)和弱變化三種局部趨勢狀態(tài)。若局部趨勢參數(shù)bi＞0的后驗概率大于0.8，p(bi＞0)＞0.8，則為強變化區(qū)，若p(bi＞0)＜0.2則為弱變化區(qū)域，其余為穩(wěn)態(tài)區(qū)。

評價貝葉斯時空分解的穩(wěn)定組分對整個時空變化過程的解釋度，常用指標(biāo)是方差成分系數(shù)（VPC），VPC越大，則時空組分分解對時空變化過程的解釋度越高，表達(dá)式為：

2 實證分析

2.1 數(shù)據(jù)與指標(biāo)

本文研究區(qū)域范圍為中國大陸31個省級行政區(qū)域，由于普查數(shù)據(jù)滯后性太大，因此選用時效性較好的《中國人口和就業(yè)統(tǒng)計年鑒》數(shù)據(jù)，該年鑒可提供2009—2014年統(tǒng)計口徑一致的省際人口遷入率指標(biāo)，該指標(biāo)是某地區(qū)某時期的遷入人口占該地區(qū)同時期總?cè)丝诘谋戎?，遷入人口指住本鄉(xiāng)、鎮(zhèn)、街道，戶口在本鄉(xiāng)、鎮(zhèn)、街道，離開戶口登記地半年以上的人口。

2.2 描述統(tǒng)計分析

總的來看，2009—2014年我國各地區(qū)人口遷入率有大幅度上升，2009年，我國總平均遷入率為8.6%（表1），2011年達(dá)到19.5%后開始緩慢下降，2014年總平均遷入率達(dá)19%。2009年只有北京、上海、廣東、福建、浙江、天津、江蘇、內(nèi)蒙8個地區(qū)人口遷入率大于10%，2014年只有西藏、江西、甘肅、河南和河北5個地區(qū)人口遷入率低于10%。

表1 2009—2014年全國各地區(qū)人口遷入率統(tǒng)計值 （單位：%）

2009—2014年，我國各地區(qū)人口流入的格局大致穩(wěn)定，上海和北京一直是我國人口流入最高的兩大地區(qū)，具有對人口遷入的巨大吸引力。其次是浙江、福建和廣東，東部沿海地區(qū)的快速發(fā)展吸引了大量人口，三個地區(qū)的人口遷入率2010年后都達(dá)到30%以上，特別是浙江和福建，幾乎達(dá)到了40%。遷入率偏低的地區(qū)有西藏、甘肅、河北、河南、江西和云南，主要分布在中西部地區(qū)，這些地區(qū)2009年和2010年，人口遷入率都不到5%，2014年都不到10%，特別是西藏，一直是我國人口遷入率最低的地區(qū)。

2.3 中國省際人口流入貝葉斯時空模型設(shè)定與檢驗

遷入人口觀測數(shù)據(jù)屬于計數(shù)數(shù)據(jù)，且個人遷入概率不完全相同，存在過度散布情形，因此數(shù)據(jù)模型采用泊松-伽馬混合模型：

其中，yit表示i(i=1, 2,..., 31)地區(qū)t年時遷入人口數(shù)，λit表示i地區(qū)t年時的遷入人口均值，uit是i地區(qū)t年時的隨機效應(yīng)參數(shù)，rit是散度系數(shù)，相應(yīng)的概率密度函數(shù)為：

過程模型為：

pit為i地區(qū)t年遷入率，α為全國總體省際人口遷入基礎(chǔ)常數(shù)，si為i地區(qū)人口遷入的空間相對風(fēng)險參數(shù)，βt*+vt描述全國人口遷入率時間變化趨勢，bit*為各地區(qū)的局部變化趨勢，bi度量從總體變化趨勢中分解出的局部變化程度，若大于零說明其局部變化強度強于總體變化趨勢，若小于零說明局部變化強度弱于總體變化強度，εit是高斯噪聲隨機變量。

參數(shù)si和bi的先驗分布用Besag York Mollie（BYM）模型確定，BYM模型通過卷積運算實現(xiàn)空間結(jié)構(gòu)和非結(jié)構(gòu)隨機效應(yīng)的相互作用，數(shù)學(xué)表達(dá)式為：

l(y|θ,Θ)為樣本似然函數(shù)，Sit表示時空域，利用條件自回歸先驗分布表示空間結(jié)構(gòu)隨機效應(yīng)，空間鄰接矩陣W

采用一階“皇后”鄰接形式，其先驗分布形式為：St是時空隨機變量，rk和ρt分別代表時間和空間相關(guān)性參數(shù)，，I是單位矩陣，σt是t時期的總體方差，上述先驗分布對應(yīng)的概率密度函數(shù)為：

本文模型估計通過實現(xiàn)貝葉斯層次模型的專門軟件WinBUGS完成。為保證模型運行的收斂性，采用了兩條MCMC鏈，每條鏈的迭代次數(shù)均設(shè)置為20萬次，所有參數(shù)的Gelman-Rubin統(tǒng)計值都低于1.005，說明該模型收斂性較好。

2.4 中國省際人口流入貝葉斯時空模型統(tǒng)計分析結(jié)果

2.4.1 空間相對風(fēng)險和總體時間趨勢估計結(jié)果

根據(jù)各地區(qū)人口遷入率的空間相對風(fēng)險大于1的后驗概率把各地區(qū)分為熱點、溫點和冷點區(qū)域（見表2）。在研究期內(nèi)，人口遷入率高于全國總體水平的熱點區(qū)域有7個（按相對風(fēng)險中值估計值降序排序）：上海（3.15）、北京（2.77）、浙江（1.86）、福建（1.74）、天津（1.49）、廣東（1.48）、江蘇（1.34），主要集中在京津地區(qū)和東部沿海地區(qū)；溫點區(qū)域有11個（按相對風(fēng)險中值估計降序排序）：內(nèi)蒙（1.28）、寧夏（1.24）、遼寧（1.13）、海南（1.13）、山西（1.00）、黑龍江（1）、吉林（0.98）、新疆（0.97）、重慶（0.95）、青海（0.88）、河北（0.83），溫點區(qū)域主要集中在東北地區(qū)，包含華北地區(qū)和西北地區(qū)部分省份。冷點區(qū)域有13個（按相對風(fēng)險中值估計降序排序）：江西（0.97）、湖北（0.86）、安徽（0.84）、貴州（0.84）、陜西（0.83）、四川（0.82）、廣西（0.78）、云南（0.77）、山東（0.75）、甘肅（0.69）、湖南（0.68）、河南（0.64）、西藏（0.63），除山東省外，冷點區(qū)域都集中于我國中西部地區(qū)。冷、熱、溫點區(qū)域具體地理分布見下頁圖1。從圖1可見，我國人口遷入的冷、熱、溫點各自集中分布特征明顯，冷點區(qū)覆蓋了我國中西部的大部分區(qū)域。溫點地區(qū)除重慶市被冷點區(qū)環(huán)繞外，整個東北地區(qū)和華北地區(qū)的內(nèi)蒙、山西與河北以及寧夏接壤成片，面積廣袤的新疆和青海緊密相連，形成我國地區(qū)人口流入的溫點區(qū)。熱點區(qū)除京津兩地外大都集中于東南沿海地區(qū)，形成地區(qū)人口流入的環(huán)帶狀熱點區(qū)域群。

表2 全國各地區(qū)人口遷入總體相對風(fēng)險和局部趨勢估計結(jié)果

不可置否，地區(qū)人口遷入與地區(qū)經(jīng)濟發(fā)達(dá)程度密切相關(guān)，2014年人均GDP排在前9位的地區(qū)依次為天津、北京、上海、江蘇、浙江、內(nèi)蒙、遼寧、福建和廣東，除內(nèi)蒙和遼寧外，其余7?。ㄊ校┒紝儆谌丝谶w入的熱點區(qū)域，人口遷入的熱點區(qū)也是地區(qū)經(jīng)濟的發(fā)達(dá)區(qū)；13個人口遷入的冷點區(qū)域有9個區(qū)域（湖北、湖南、陜西、山東除外）人均GDP的排名位于最末10位，人口遷入的冷點區(qū)域也是地區(qū)經(jīng)濟欠發(fā)達(dá)區(qū)域；11個溫點區(qū)域其GDP排名除內(nèi)蒙和遼寧較靠前、山西較靠后外，其余9個溫點區(qū)域在全國的位次排在11位至21位之間。地區(qū)人口遷入的冷熱圖譜與地區(qū)經(jīng)濟發(fā)達(dá)程度圖譜高度一致。

圖1 各地區(qū)人口遷入熱點、溫點和冷點分布圖

圖2 人口遷入率總體時間變化曲線（95%置信區(qū)間）

全國省際人口遷入總體變化趨勢（圖2）表現(xiàn)為曲折上升，局部表現(xiàn)為先上升后下降，2009—2011年表現(xiàn)為較快的上升，2011—2014年呈現(xiàn)逐年緩慢下降的態(tài)勢。全國總體人口遷入率的局部走勢與我國整體經(jīng)濟走勢基本一致，2009—2014年，我國經(jīng)濟增長率依次為9.2%、10.6%、9.5%、7.7%、7.7%和7.3%，2009—2010年表現(xiàn)為上升趨勢，2010—2014表現(xiàn)出逐年下降趨勢，省際人口流入對經(jīng)濟放緩的響應(yīng)具有一定的時滯，所以人口流入下降趨勢的出現(xiàn)比經(jīng)濟增長的下降延遲了一年。

本文模型估算結(jié)果的方差成分系數(shù)VPC等于94%，說明該模型的時空分解組分對整個時空變化過程的解釋度比較高，但局部變化趨勢也在起作用，地區(qū)人口流入的局部變化趨勢影響著全國省際人口流入的空間相對關(guān)系。2.4.2 局部變化趨勢估計結(jié)果

考察各地區(qū)局部變化系數(shù)大于0的后驗概率估計（表2），根據(jù)Richardson原則把全國31個省區(qū)分為強變化區(qū)域、穩(wěn)定態(tài)區(qū)域和弱變化區(qū)域三個類別。我國大部分地區(qū)為人口局部變化的穩(wěn)態(tài)區(qū)，包括山西、內(nèi)蒙、遼寧、吉林、福建、江西、山東、海南、云南、西藏、陜西、青海和新疆13個地區(qū)，根據(jù)模型中局部趨勢的統(tǒng)計意義，這些地區(qū)的人口流入局部變化趨勢與全國總體變化趨勢相當(dāng)。穩(wěn)態(tài)區(qū)中只有福建是熱點區(qū)，說明福建人口流入程度高于全國平均水平但人口流入的局部變化與全國水平相當(dāng)。穩(wěn)態(tài)區(qū)中包含山西、內(nèi)蒙、遼寧、吉林、海南、青海、新疆7個溫點區(qū)，這7個區(qū)域人口流入水平和人口流入的局部變化都與全國平均水平相當(dāng)。穩(wěn)態(tài)區(qū)中江西、山東、云南、西藏、陜西為冷點區(qū)，這5個地區(qū)人口流入水平較全國低但人口流入局部變化與全國總體變化相當(dāng)。強變化區(qū)域有9個，包括安徽、河南、湖北、湖南、廣西、重慶、四川、貴州和甘肅，這9個地區(qū)的人口流入局部變化趨勢強于總體變化趨勢，又由于總體變化趨勢是上升的，所以這7個地區(qū)人口流入的增加速度高于總體平均增加速度。9個強變化區(qū)域中重慶為溫點區(qū)，其他8個強變化區(qū)為冷點區(qū)。弱變化區(qū)域有北京、天津、河北、黑龍江、上海、江蘇、浙江、廣東、寧夏這9個地區(qū)，弱變化區(qū)中大部分為熱點區(qū)，少部分為溫點區(qū)，不包含冷點區(qū)，包含的熱點區(qū)有6個：北京、天津、上海、江蘇、浙江、廣東，這些熱點區(qū)雖然人口流入率高于全國平均水平，但其人口流入率上升變化小于全國水平。弱變化區(qū)中的河北、黑龍江和寧夏為溫點區(qū)，3個溫點區(qū)人口流入水平與全國相當(dāng)，其局部變化趨勢弱于全國總體變化趨勢。

3 結(jié)論

本文首次應(yīng)用Bayesian層次時空模型，從總體和局部兩個層面分析了我國地區(qū)人口流入在空間和時間上的變化特征。目前我國各地區(qū)總體平均人口流入率已接近20%，各地區(qū)人口流入的空間相對風(fēng)險和局部變化趨勢存在著較大的差異，本文研究發(fā)現(xiàn)：第一，我國人口遷入的冷、熱、溫點各自集中分布特征明顯，地區(qū)人口遷入的冷熱圖譜與地區(qū)經(jīng)濟發(fā)展圖譜高度一致，近幾年，隨著我國經(jīng)濟增長的放緩，地區(qū)人口的遷入率在逐漸下降；第二，人口遷入的熱點區(qū)有7個，主要集中在京津地區(qū)和東部沿海地區(qū)，特別是上海和北京，人口遷入率一直是全國最高的。北京、天津、上海、江蘇、浙江、廣東這6個熱點區(qū)同時屬于局部變化的弱變化區(qū)，這些地區(qū)雖然人口流入率高于全國平均水平，但其人口流入率上升變化小于全國水平。熱點區(qū)中的福建屬于局部變化的穩(wěn)態(tài)區(qū)；第三，我國人口遷入的冷點區(qū)有13個，大都集中于中西部地區(qū)，冷點區(qū)中的安徽、河南、湖北、湖南、廣西、四川、貴州和甘肅這8個地區(qū)同時屬于局部強變化區(qū)，這些地區(qū)人口遷入率低于全國平均水平，但人口流入率的增長態(tài)勢高于全國。江西、山東、云南、西藏、陜西這5個冷點區(qū)同時屬于穩(wěn)態(tài)區(qū)，冷點區(qū)中無弱變化區(qū)；第四，溫點區(qū)域有11個，主要集中在東北地區(qū)，包含華北地區(qū)和西北地區(qū)部分省份。溫點區(qū)中的重慶屬于強變化區(qū)，說明重慶地區(qū)人口流入程度與全國平均水平相當(dāng)，但人口流入率的增加速度高于全國。河北、黑龍江和寧夏這3個溫點區(qū)屬于弱變化區(qū)，其他7個溫點區(qū)同時屬于局部變化的穩(wěn)態(tài)區(qū)。

貝葉斯層次時空模型是研究人口流動時空特征的有效方法?；诒疚难芯拷Y(jié)果，相關(guān)的政策制定需要充分認(rèn)識和掌握我國省際人口流入的空間差異和局部變化趨勢的異質(zhì)性。充分認(rèn)識北京、天津、上海、江蘇、浙江、廣東這6個熱點區(qū)的局部弱變化性，認(rèn)識安徽、河南、湖北、湖南、廣西、四川、貴州和甘肅這8個冷點區(qū)的局部強變化性，注重人口遷入冷熱與經(jīng)濟增長快慢之間以及經(jīng)濟發(fā)展區(qū)域格局和人口流入?yún)^(qū)域格局之間的關(guān)聯(lián)，通過協(xié)調(diào)區(qū)域經(jīng)濟發(fā)展引導(dǎo)人口流動趨于合理。

參考文獻(xiàn)：

[1]Gelman A,Rubin D B.Inference from Iterative Simulation Using Multiple Sequences.[J].Statist.Sci,1992(7).

[2]G Li,R Haining,S Richardson,N Best.Space-time Variability in Burglary Risk：A Bayesian Spatio-temporal Modeling Approach[J].Spatial Statistics,2014,(9).

[3]Richardson S,Thomson A,Best N,Elliot P.Interpreting Posterior Relative Risk Estimates in Disease-mapping Studies.[J]Environ.Health Perspect.,2004,112(9).

[4]Goldstein H,Browne W J,Rasbash J.Partitioning variation in multilevel models.[J],Underst.Statist,2002,(1).

[5]Besag J,York J,Mollie A.Bayesian Image Restoration,with two Applications in Spatial Statistics[J].Ann.Inst.Math.Stat,1991,(43).

[6]Gelman A.Prior Distribution for Variance Parameters in Hierarchical Models[J].Bayesian Anal,2006,(1).

[7]Lunn D.J.,Thomas A,Best N,Spiegelhalter D.WinBUGS-a Bayesian Modeling Framework：Concept,Structure,and Extensibility.[J]Stat.Comput,2000,(10).

[8]文雯,文小焱,胡珊等.貝葉斯層次模型在嵌套結(jié)構(gòu)調(diào)查數(shù)據(jù)中的應(yīng)用研究[J].中國衛(wèi)生統(tǒng)計,2015,32(2).