我國(guó)圍產(chǎn)兒死亡率時(shí)空變化及預(yù)測(cè)研究*

濱州醫(yī)學(xué)院公共衛(wèi)生與管理學(xué)院(264003) 趙星宇 王廣成 單海峰 楊寶順 丁素素 梁曉天 韓春蕾

【提 要】 目的 研究我國(guó)圍產(chǎn)兒死亡率的時(shí)空變化趨勢(shì)并進(jìn)行預(yù)測(cè)。方法 收集2005-2018年圍產(chǎn)兒死亡率相關(guān)數(shù)據(jù)，使用ArcGIS 10.2軟件繪制圍產(chǎn)兒死亡率分布地圖，用Stata 14.0軟件進(jìn)行空間相關(guān)分析，使用SPSS 19.0軟件進(jìn)行整合移動(dòng)平均自回歸模型(ARIMA)預(yù)測(cè)分析。結(jié)果 2005-2018年，圍產(chǎn)兒死亡率西部大于東部和中部，西藏和新疆是圍產(chǎn)兒死亡率的高發(fā)地區(qū)；圍產(chǎn)兒死亡率空間聚集性逐年上升(0.375到0.474)；預(yù)測(cè)得出2019-2022年圍產(chǎn)兒死亡率呈現(xiàn)下降趨勢(shì)(4.08‰到3.92‰)。結(jié)論 我國(guó)圍產(chǎn)兒死亡率呈下降趨勢(shì)，存在明顯的地理差異；圍產(chǎn)兒死亡率的降低需要臨近省份的共同努力。

婦嬰保健是建設(shè)健康中國(guó)的重要組成部分，是實(shí)現(xiàn)“健康中國(guó)2030”發(fā)展目標(biāo)的關(guān)鍵。對(duì)圍產(chǎn)兒死亡率的研究有助于婦幼保健的可持續(xù)發(fā)展[1]。胎兒體重達(dá)到1000克，或孕期滿28周至出生后7天以內(nèi)的死亡稱為圍產(chǎn)兒死亡[2]。國(guó)內(nèi)外學(xué)者對(duì)圍產(chǎn)兒死亡率相關(guān)問(wèn)題進(jìn)行了系列研究。如Xiang K、蔡桂舉、朱丹玲等對(duì)圍產(chǎn)兒死亡率現(xiàn)狀和分布的研究[3-5]，Ananth CV、杜兩省、莊艷艷等對(duì)圍產(chǎn)兒死亡率影響因素的研究[6-8]以及戴瓊對(duì)圍產(chǎn)兒死亡率的預(yù)測(cè)研究等[9]。

研究目前主要集中在對(duì)圍產(chǎn)兒死亡率的現(xiàn)狀、影響因素或預(yù)測(cè)的單一研究，缺少系統(tǒng)的研究。本文綜合使用空間計(jì)量分析與時(shí)間序列預(yù)測(cè)的方法，對(duì)圍產(chǎn)兒死亡率的時(shí)空分布現(xiàn)狀進(jìn)行研究，并預(yù)測(cè)未來(lái)我國(guó)圍產(chǎn)兒死亡率的發(fā)展趨勢(shì)，可為降低我國(guó)圍產(chǎn)兒死亡率、提高婦幼健康水平提供對(duì)策建議。

資料及方法

1.數(shù)據(jù)來(lái)源

收集我國(guó)2005-2018年31個(gè)省份的圍產(chǎn)兒死亡率數(shù)據(jù)，數(shù)據(jù)來(lái)自2006-2019年《中國(guó)衛(wèi)生統(tǒng)計(jì)年鑒》。

2.研究方法

空間相關(guān)性可以判斷圍產(chǎn)兒死亡率的分布是否存在統(tǒng)計(jì)上集聚或分散現(xiàn)象。在一定的顯著性水平下，Moran′s I 指數(shù)為正表示圍產(chǎn)兒死亡率存在空間集聚性；Moran′s I 指數(shù)為負(fù)則表示鄰近區(qū)域的圍產(chǎn)兒死亡率差異顯著，即存在空間分散性[10]。整合移動(dòng)平均自回歸模型(auto regressive integrated moving average model，ARIMA)即對(duì)差分平穩(wěn)的序列擬合自回歸移動(dòng)平均模型(auto regressive moving average model，ARMA)，是常用的時(shí)間序列預(yù)測(cè)方法。ARIMA(p，d，q)中，AR和MA分別表示自回歸和移動(dòng)平均，p、d、q分別為自回歸項(xiàng)數(shù)、平穩(wěn)序列所做的差分次數(shù)(階數(shù))和移動(dòng)平均項(xiàng)數(shù)[11]。

本文使用Stata 14.0軟件，利用Moran′s I檢驗(yàn)探究數(shù)據(jù)的空間相關(guān)性，使用SPSS19.0進(jìn)行ARIMA模型預(yù)測(cè)，使用ArcGIS10.2軟件繪制圍產(chǎn)兒死亡率分布地圖，分析其空間分布情況。以α=0.05作為顯著性標(biāo)準(zhǔn)。

3.指標(biāo)相關(guān)計(jì)算方法

圍產(chǎn)兒死亡率(‰)指使用每1000名嬰兒出生所對(duì)應(yīng)的圍產(chǎn)期死亡數(shù)。死亡率降幅(%)=(2018圍產(chǎn)兒死亡率-2005圍產(chǎn)兒死亡率)/2005圍產(chǎn)兒死亡率×100%。此外，相對(duì)誤差(%)=絕對(duì)誤差/真實(shí)值×100%。

結(jié) 果

1.我國(guó)圍產(chǎn)兒死亡率時(shí)間與空間分布

為研究我國(guó)不同地區(qū)圍產(chǎn)兒死亡率的走勢(shì)，將我國(guó)省份分為東部地區(qū)、中部地區(qū)與西部地區(qū)三部分[12]，繪制圍產(chǎn)兒死亡率走勢(shì)圖如圖1所示。

圖1 2005-2018年我國(guó)各地區(qū)總體圍產(chǎn)兒死亡率走勢(shì)圖

可見(jiàn)，我國(guó)圍產(chǎn)兒死亡率總體呈下降趨勢(shì)，由2005年的10.27‰下降到2018年的4.26‰。其中，東部地區(qū)圍產(chǎn)兒死亡率最低，西部地區(qū)最高，始終保持著“東—中—西”階梯式遞增態(tài)勢(shì)。三個(gè)地區(qū)中，中部圍產(chǎn)兒死亡率降幅最大(55.80%)，其次為西部地區(qū)(55.71%)，東部地區(qū)圍產(chǎn)兒死亡率降幅最少(52.67%)。

以2005和2018年為基準(zhǔn)，將31個(gè)省份劃分為圍產(chǎn)兒低死亡率(0～)、較低死亡率(5‰～)、較高死亡率(10‰～)和高死亡率(15‰～)四種類型。繪制圍產(chǎn)兒死亡率的空間分布圖，見(jiàn)圖2。

圖2 2005年(左)與2018年(右)我國(guó)圍產(chǎn)兒死亡率空間分布

2005年到2018年，我國(guó)各省份死亡率均呈現(xiàn)降低的趨勢(shì)。從空間分布圖來(lái)看，各省份圍產(chǎn)兒死亡率呈現(xiàn)明顯的空間差異，圍產(chǎn)兒死亡率由東部向西部呈現(xiàn)增高趨勢(shì)。經(jīng)濟(jì)落后的西部地區(qū)圍產(chǎn)兒死亡率較高，經(jīng)濟(jì)發(fā)達(dá)的東部地區(qū)圍產(chǎn)兒死亡率較低。2005年我國(guó)圍產(chǎn)兒死亡率為10.27‰，在發(fā)達(dá)的地區(qū)如北京、上海、天津等地，圍產(chǎn)兒死亡率均在7.00‰以下，其中上海市最低(2.52‰)。在中部地區(qū)如山西、河南、湖北、廣東等地，圍產(chǎn)兒死亡率介于7.00‰和13.00‰之間。在不發(fā)達(dá)的西部地區(qū)如新疆和西藏，圍產(chǎn)兒死亡率高達(dá)19.97‰和19.04‰。2018年我國(guó)圍產(chǎn)兒死亡率為4.26‰，各省份均有所改善。上海圍產(chǎn)兒死亡率依舊最低(2.38‰)，北京市(2.86‰)位居第二；西部地區(qū)死亡率均有顯著改善，新疆和西藏的圍產(chǎn)兒死亡率分別降為11.01‰和14.67‰。

2.圍產(chǎn)兒死亡率空間自相關(guān)分析

2005年和2018年我國(guó)圍產(chǎn)兒死亡率的全局莫蘭指數(shù)在0.375～0.495，均達(dá)到顯著水平(P<0.05)且具有上升趨勢(shì)。如圖3所示，圍產(chǎn)兒死亡率空間聚集程度逐漸上升。

圖3 2005-2018年圍產(chǎn)兒死亡率莫蘭指數(shù)走勢(shì)圖

3.我國(guó)圍產(chǎn)兒死亡率的預(yù)測(cè)

圍產(chǎn)兒死亡率數(shù)據(jù)經(jīng)過(guò)1次差分后(d=1)，自相關(guān)圖(ACF)與偏自相關(guān)圖(PACF)見(jiàn)圖4。

圖4 2005-2018年圍產(chǎn)兒死亡率數(shù)據(jù)經(jīng)一階差分后的ACF圖和PACF圖

自相關(guān)圖與偏自相關(guān)圖的結(jié)果顯示，置信區(qū)間在[-0.5，0.5]范圍內(nèi)，故在進(jìn)行一階差分之后的序列趨于平穩(wěn)。以圍產(chǎn)兒死亡率為因變量，采用SPSS 19.0擬合ARIMA模型。通過(guò)BIC指數(shù)與R2值對(duì)模型進(jìn)行篩選，最終得到ARIMA(1，1，2)模型。擬合R2=0.980，相對(duì)誤差均在10%以內(nèi)，說(shuō)明模型擬合良好(表1)。

表1 ARIMA模型對(duì)全國(guó)圍產(chǎn)兒死亡率擬合結(jié)果

2019-2022年我國(guó)圍產(chǎn)兒死亡率預(yù)測(cè)結(jié)果分別為4.08‰、4.06‰、3.86‰和3.92‰，呈整體下降趨勢(shì)(圖5)。

圖5 ARIMA模型對(duì)我國(guó)圍產(chǎn)兒死亡率的回代及預(yù)測(cè)情況

對(duì)31個(gè)省份圍產(chǎn)兒死亡率進(jìn)行ARIMA模型擬合，預(yù)測(cè)2019-2022年各省市圍產(chǎn)兒死亡率，擬合結(jié)果相對(duì)誤差均在10%以內(nèi)，擬合效果良好。將預(yù)測(cè)結(jié)果劃分為低死亡率(0～)、較低死亡率(3‰～)、較高死亡率(6‰～)和高死亡率(9‰～)四種類型，繪制空間分布圖如圖6所示。

圖6 2019-2022年我國(guó)圍產(chǎn)兒死亡率預(yù)測(cè)結(jié)果的空間分布

從空間分布圖來(lái)看，2019-2022年全國(guó)各省份圍產(chǎn)兒死亡率均呈現(xiàn)下降趨勢(shì)，我國(guó)圍產(chǎn)兒死亡率依然東低西高，但東西部地區(qū)差異具有減小的趨勢(shì)。其中上海的圍產(chǎn)兒死亡率由2019年的2.05‰降為2022年的1.63‰，西藏和新疆的死亡率分別由12.56‰和10.70‰降為5.42‰和8.43‰。

討 論

本文采用空間相關(guān)模型分析了圍產(chǎn)兒死亡率的時(shí)空變化趨勢(shì)，采用ARIMA模型對(duì)2005-2018年我國(guó)的省份數(shù)據(jù)進(jìn)行預(yù)測(cè)。從地區(qū)分布來(lái)看，我國(guó)圍產(chǎn)兒死亡率呈現(xiàn)東-中-西逐漸升高的趨勢(shì)。說(shuō)明圍產(chǎn)兒死亡率與地區(qū)經(jīng)濟(jì)狀況緊密相關(guān)，經(jīng)濟(jì)發(fā)達(dá)的東部地區(qū)相比相對(duì)落后的中西部地區(qū)圍產(chǎn)兒死亡率更低，這與前期研究結(jié)果一致[13]。明確圍產(chǎn)兒死亡率的地區(qū)差異，有針對(duì)性地對(duì)不同地區(qū)采取差異性的防治干預(yù)措施，有利于更好地做好圍產(chǎn)期保健工作[8]。中西部地區(qū)對(duì)進(jìn)一步降低全國(guó)的圍產(chǎn)兒死亡率至關(guān)重要。因此，需要加大對(duì)中西部地區(qū)尤其是西部地區(qū)的政策傾斜力度。以共建“一帶一路”為引領(lǐng)，加快西部大開(kāi)發(fā)戰(zhàn)略的實(shí)施力度，提高西部地區(qū)的經(jīng)濟(jì)和社會(huì)發(fā)展水平；同時(shí)加大對(duì)中西部地區(qū)的政府衛(wèi)生投入，出臺(tái)相應(yīng)政策留住更多的衛(wèi)生技術(shù)人員，以進(jìn)一步降低圍產(chǎn)兒死亡率[13]。

我國(guó)圍產(chǎn)兒死亡率歷年的莫蘭指數(shù)為正且逐漸上升，表明我國(guó)的圍產(chǎn)兒死亡率分布情況具有一定的空間集聚性且隨時(shí)間呈上升趨勢(shì)。即一個(gè)省份的圍產(chǎn)兒死亡率與周圍省份的圍產(chǎn)兒死亡率存在關(guān)聯(lián)，臨近省份的圍產(chǎn)兒死亡率存在潛在的相互依賴性，且這種依賴性隨著時(shí)間推移呈增大趨勢(shì)。由于空間聚集性的存在，要降低某一省份的圍產(chǎn)兒死亡率，需要綜合提升本省份以及周圍省份的經(jīng)濟(jì)水平、醫(yī)療衛(wèi)生水平、城市化水平及環(huán)境質(zhì)量等影響居民健康時(shí)空變化的重要因素[10]。

我國(guó)圍產(chǎn)兒死亡率的預(yù)測(cè)結(jié)果表明，2019-2022年我國(guó)各省份圍產(chǎn)兒死亡率將繼續(xù)呈現(xiàn)下降趨勢(shì)。說(shuō)明隨著時(shí)間推移，我國(guó)圍產(chǎn)醫(yī)學(xué)發(fā)展水平及婦幼衛(wèi)生工作質(zhì)量將逐步提高，兒童生存和兒童健康保健狀況將不斷改善[9]。利用ARIMA模型對(duì)圍產(chǎn)兒死亡率的預(yù)測(cè)具有較高的精度(10%以內(nèi))，可為母嬰健康政策及兒童早期發(fā)展政策的制定提供保障和科學(xué)依據(jù)。

值得一提的是，ARIMA模型的預(yù)測(cè)精度受樣本量的影響較大[14]，在變量個(gè)數(shù)一定的情況下，樣本量越多，模型的預(yù)測(cè)精度就會(huì)越高[11]。本文使用了2005-2018年14年的數(shù)據(jù)擬合ARIMA模型對(duì)圍產(chǎn)兒死亡率進(jìn)行預(yù)測(cè)，雖然預(yù)測(cè)的相對(duì)誤差均在10%以內(nèi)，但日后的研究可通過(guò)增加樣本量的方法進(jìn)一步提高預(yù)測(cè)精度。