采用時間序列泊松回歸分析天津市大氣污染物對心腦血管疾病死亡的急性影響*

王德征，江國虹，顧清，張輝，徐忠良、宋桂德，張穎、沈成鳳

采用時間序列泊松回歸分析天津市大氣污染物對心腦血管疾病死亡的急性影響*

王德征，江國虹**，顧清，張輝，徐忠良、宋桂德，張穎、沈成鳳

目的：了解天津市大氣污染物對居民心腦血管疾病死亡的急性影響，探索早期預(yù)警方法，為心腦血管疾病的預(yù)防控制提供依據(jù)。

方法：采用天津市疾病預(yù)防控制中心收集的居民全死因監(jiān)測數(shù)據(jù)， 氣象資料和大氣污染物資料來源于天津市氣象局和天津市環(huán)境監(jiān)測中心。采用時間序列的泊松回歸GAM進行天津市每日大氣污染物與居民心腦血管疾病死亡危險度分析， 同時控制氣象因素、長期趨勢、星期幾效應(yīng)以及人口數(shù)等混雜因素的影響，進行單污染物和多污染物分析。

結(jié)果：單因素GAM分析顯示大氣中SO2、NO2、PM10效應(yīng)在當日達到最大，風險RR值分別增加1.13% [95%CI（0.76～1.51） %]、 0.78% [95%CI（ 0.41～1.15）%]和 0.61% [95%CI（0.51～0.71）%]；SO2、NO2、PM10日均濃度每升高10 μg/m3，0～5日后心腦血管疾病死亡的風險RR值分別增加0.70% [95%CI（0.47～0.94）]%， 0.51% [95%CI（0.27～0.74）%]和 0.16% [95%CI（0.06～0.27）%]；多因素分析顯示，SO2、NO2、PM10日均濃度每增加10 μg/m3，心腦血管疾病死亡的風險RR值分別增加0.77% [95%CI（0.58～0.97）%]，0.41% [95%CI（0.05～0.78）%]和 0.38% [95%CI（0.12～0.64）%]。

結(jié)論：天津市大氣污染物能增加人群心腦血管疾病死亡風險，可能對心腦血管疾病急性發(fā)作起到誘導作用，有必要開展大氣污染物治理，建立預(yù)警系統(tǒng)，合理分配衛(wèi)生資源，降低心腦血管疾病死亡危害。

心腦血管疾病；死亡率；廣義可加模型; 大氣污染物; 急性影響

Methods: The Mortality for patients with cardio cerebral vascular disease from 2001-01 to 2009-12 was from Tianjin Centers for Disease Control and Prevention, the meteorological and air pollution data were from Tianjin meteorological bureau and Tianjin environmental monitoring station respectively. The Single and multiple generalized additive model (GAM) extended poisson regression analysis was performed to calculate the relationship between air pollution and cardio-cerebral vascular disease mortality by controlling the time trends, weather, the day of week and air population.

Results: Results: Single GAM analysis indicated that when the effect of SO2, NO2and PM10in the air reached the maximum in the day, the RR values for the mortality in patients with cardio-cerebral vascular disease increased 1.13%

[95%CI (0.76-1.51)%], 0.78% [95%CI ( 0.41-1.15)%] and 0.61% [95%CI ( 0.51-0.71)%] respectively; when the average concentration of SO2, NO2and PM10 increasing 10μg/m3per day, after 0-5 days, the RR values for the mortality elevated 0.70% [95%CI (0.47-0.94)%], 0.51% [95%CI (0.27-0.74)%] and 0.16% [95%CI (0.06-0.27)%] respectively. Multiple GAM analysis presented that when SO2, NO2and PM10increasing 10 μg/m3per day, the RR values for the mortality elevated 0.77% [95%CI (0.58-0.97)%], 0.41% [95%CI (0.05-0.78) %] and 0.38% [95%CI (0.12-0.64%)%] respectively. Conclusion: The air pollution could increase the mortality risk in patients with cardio-cerebral vascular disease, it is necessary to establish the prevention system in order to decrease the mortality risk in those patients.

(Chinese circulation Journal, 2014,29:453.)

在全球范圍內(nèi)， 60%的由心腦血管疾病導致的負擔都發(fā)生在發(fā)展中國家[1]。我國城市和農(nóng)村冠心病在全死因順位中居第三或第四，腦卒中在全死因順位中居第二或第三[2]，是導致我國居民死亡的主要原因之一。心腦血管疾病是一種多因素疾病，是遺傳因素和環(huán)境因素共同作用的結(jié)果[3]，特別是隨著工業(yè)化、城市化和人口快速增長，人類活動導致的環(huán)境變化特別是可吸入顆粒物（PM10）、二氧化硫（SO2）、二氧化氮（ NO2）、一氧化碳（CO） 、臭氧（O3）等的增加對心腦血管疾病發(fā)病、死亡的影響已引起人們廣泛關(guān)注。天津市是典型的北方工業(yè)城市, 探討大氣污染物對天津市居民心腦血管疾病死亡的急性影響，探索早期預(yù)警方法，為心腦血管疾病的預(yù)防控制提供依據(jù)。

1 材料和方法

材料：采用2001-01-01至2009-12-31天津市疾病預(yù)防控制中心所收集的天津市居民全死因監(jiān)測數(shù)據(jù)。心腦血管疾病診斷分類按照“疾病和有關(guān)健康問題的國際統(tǒng)計分類 （ICD） ”：ICD-9: 390-459；ICD-10: I00-I99。以實際死亡日期為分析時間點。采用天津市氣象局和天津市環(huán)境監(jiān)測中心收集的天津市同時段氣象指標（包括每日平均溫度、平均相對濕度、平均大氣壓、平均風速等，以及大氣污染物資料包括8個國控測點的可吸入顆粒物PM10、SO2及NO2的每日平均值，每日24小時24次監(jiān)測值的算術(shù)平均值）。對大氣污染物資料、氣象資料中的缺失值，采用前后各兩天數(shù)值的中位數(shù)為缺失值當天的數(shù)值，缺失值填補比例分別為0.06%和0.91%。

統(tǒng)計方法：統(tǒng)計軟件為SPSS11.5軟件和SAS8.1。對于總?cè)丝趤碚f，每日心腦血管疾病死亡為小概率事件，作為一種時間序列資料，其實際分布近似泊松分布。應(yīng)采用時間序列的泊松回歸廣義可加模型（GAM）。具體模型為:

[Yi ：觀察日；i 當天的死亡人數(shù)， E（ Yi ） ：觀察日i日死亡數(shù)的預(yù)期值；α：截距；X：對應(yīng)變量產(chǎn)生線性影響指示變量；β：通過回歸模型估計的指示變量系數(shù)： f ：非參數(shù)自然三次光滑樣條函數(shù)；Z ：對應(yīng)變量發(fā)生非線性影響的變量，∑：加和]。

通過擬合單因素GAM并進行偏差性檢驗，同時考慮氣象因素、人口結(jié)構(gòu)、星期幾效應(yīng)以及長期趨勢等混雜因素，分析大氣單污染物對心腦血管疾病死亡的影響，并用相對危險度（RR）評價危害的大小。考慮到污染物對人體危害可能的滯后效應(yīng)，選取遲滯日（滯后天數(shù)）1～5天，分析大氣污染物對心腦血管疾病死亡影響的變化情況。為檢驗?zāi)Ｐ偷姆€(wěn)定性,采用廣義交叉確認（ GCV） 法進行不同統(tǒng)計學模型的擬合優(yōu)度檢驗[4]。

對單因素分析有意義、相對危險度較大且模型擬合優(yōu)度較好的污染物、滯后天數(shù)等，并考慮可能的混雜因素，建立共同擬合多因素GAM，并進行偏差性檢驗和擬合優(yōu)度檢驗。

以2000年世界標準人口計算年齡、性別標化死亡率，采用趨勢卡方檢驗對心腦血管疾病的年死亡率進行趨勢分析。

采用Spearman等級相關(guān)分析大氣污染物、氣象因素指標與心腦血管疾病死亡的相關(guān)性。

2 結(jié)果

2001～2009年天津市每日心腦血管疾病死亡、氣象因素和大氣污染物濃度分布：心腦血管疾病死

亡、氣象因素和大氣污染物日間分布存在明顯的周期性（圓分布檢驗，P＜0.01），且均不符合正態(tài)分布（P＜0.05）。2001～2009年天津市每日中位心腦血管疾病死亡數(shù)為77例；氣象因素中，中位氣溫、濕度、氣壓、風速分別為16 ℃、69 %、1 015 hPa、10 Km/h；大氣污染物中的SO2、NO2和PM10中位濃度分別為4.90 10 μg/m3、 4.68 10 μg/m3和10.32 10 μg/m3。表1

2001～2009年天津市每日心腦血管疾病死亡數(shù)與氣象因素和大氣污染物因素Spearman相關(guān)分析：SO2、NO2每日濃度同心腦血管疾病每日死亡數(shù)間呈正相關(guān)（r 分別為0.479、0.250），差異有統(tǒng)計學意義（ P＜0.01）；PM10每日濃度同心腦血管疾病每日死亡數(shù)間呈負相關(guān)（r=-0.023），但差異無統(tǒng)計學意義（ P=0.182）。每日氣象因素中，每日溫度、每日濕度與心腦血管疾病每日死亡數(shù)呈負相關(guān)（r 分別為-0.572和-0.164），差異有統(tǒng)計學意義（ P＜0.01）；每日氣壓、每日風速與心腦血管疾病每日死亡數(shù)呈正相關(guān)（r 分別為0.479和0.063），差異有統(tǒng)計學意義（P＜0.01）。同時每日大氣污染物、氣象因素內(nèi)部及相互之間呈現(xiàn)不同程度的相關(guān)關(guān)系，具體結(jié)果見表2。

表1 2001～2009年天津市每日心腦血管疾病死亡例數(shù)、氣象因素和大氣污染物濃度分布

表2 2001～2009年天津市每日心腦血管疾病死亡數(shù)與氣象因素和大氣污染物因素Speaman相關(guān)分析

大氣污染物對不同遲滯日心腦血管疾病死亡的單因素GAM分析：經(jīng)過可能的混雜因素調(diào)整后大氣污染物日均濃度對當日及不同滯天數(shù)的心腦血管疾病死亡的泊松回歸GAM分析顯示: SO2、NO2、PM10日均濃度每升高10 μg/m3，對心腦血管疾病死亡的風險在當日達到最大，其RR值及95%可信區(qū)間分別增加 1.13% [95%CI（0.76～1.51 ） %]、0.78% [95%CI（ 0.41～1.15）%]和 0.61% [95%CI（0.51～0.71）%]； 5日后SO2、NO2、PM10日均濃度每升高10 μg/m3，心腦血管疾病死亡的風險RR值分別增加0.70% [95%CI（0.47～0.94 ） %], 0.51% [95%CI（0.27～0.74）%]和 0.16% [95%CI（0.06～0.27）%]，其滯后作用隨著時間延長而逐漸減小。表3

大氣污染物及相關(guān)因素對心腦血管疾病死亡的多因素GAM分析：由污染物對心腦血管疾病死亡影響的單因素分析，選取當日（Lag0）SO2、NO2、PM10及可能的混雜因素進行多因素泊松GAM分析，SO2、NO2、PM10日均濃度每增加10 μg/m3，心腦血管疾病死亡的風險RR值及95% IC 分別增加0.77%（0.58～0.97）%，0.41% （0.05～0.78）%和 0.38%（0.12～0.64）%。與此同時，低溫、低氣壓是當日（Lag0）心腦血管疾病死亡的危險因素（P＜0.05）; 隨著每年人口數(shù)的增加，心腦血管疾病死亡數(shù)也在增加，P=0.0000。結(jié)果見表4。

表3 大氣污染物對不同遲滯日心腦血管疾病死亡的單因素GAM分析

表4 大氣污染物及相關(guān)因素對心腦血管疾病死亡的多因素GAM分析

3 討論

本研究結(jié)果表明，天津市2001至2009年SO2、NO2和PM10平均水平，高于國家環(huán)境空氣質(zhì)量標準[5]二級水平，其中NO2平均水平低于北京和上海，SO2平均水平高于北京和上海，PM10平均水平在上海和北京之間[6]。天津市心腦血管疾病死亡、大氣污染物和氣象因素日間分布存在明顯的相關(guān)性，低氣溫能增加心腦血管疾病死亡的風險。Scarborough 等[7]對英格蘭1999～2004年冠心病死亡數(shù)據(jù)的分析表明，人群冠心病的死亡率隨著溫度的降低和日照時間的縮短而增加。此前關(guān)于天津市心腦血管病死亡周期性的時間序列分析[8]也表明，冠心病、腦卒中死亡率在冬季，特別是1月份呈現(xiàn)高水平。天津市四季分明, 采暖期污染物的水平明顯高于非采暖期污染物水平。Choi等[9]研究認為大氣污染物與心腦血管疾病的關(guān)系, 必須要考慮氣象因素的影響，其中對血壓產(chǎn)生影響的大氣污染物因季節(jié)不同而不同, 在夏季主要是PM10和NO2, 而在冬季則是SO2和O3。

本研究大氣污染物單因素分析與許多研究結(jié)果相似，多因素分析模型同時引入各污染物，危險度分析結(jié)果穩(wěn)健[10,11]。

許多研究認為大氣污染物對人體健康的急性影響存在滯后效應(yīng) ，Analitis等[12]對歐洲29個城市的Meta分析結(jié)果說明大氣污染物對心腦血管系統(tǒng)疾病的滯后效應(yīng)可能比對呼吸系統(tǒng)疾病明顯。本研究選取了大氣污染物當日及滯后1～5天的心腦血管疾病數(shù)據(jù)進行分析，結(jié)果顯示心腦血管疾病死亡的滯后作用隨著時間延長而逐漸減小。

目前，大氣污染物對人體心血管系統(tǒng)的影響和機制研究是初步的和很不深入的。有研究表明，大氣污染物對心血管系統(tǒng)健康的影響中, 氧化應(yīng)激和炎性反應(yīng)是至關(guān)重要的, 可以引起動脈粥樣斑塊、血管內(nèi)皮、血小板和心律的變化[13]，大氣污染物暴露會損傷微血管功能，進而引起NO介導的交感神經(jīng)興奮[14]。提示大氣污染物可能對心血管病急性發(fā)病起到誘導作用。

綜上所述，本研究發(fā)現(xiàn)SO2、NO2和PM10日平均濃度短期升高是人群心腦血管疾病死亡的危險因素，有必要開展大氣污染物治理，合理分配衛(wèi)生資源，降低心腦血管死亡危害。

[1]Kim AS, Johnston SC. Global variation in the relative burden of stroke and ischemic heart disease. Circulation, 2011, 124: 314-323.

[2]中華人民共和國衛(wèi)生部. 2010中國衛(wèi)生統(tǒng)計年鑒. 北京: 中國協(xié)和醫(yī)科大學出版社, 2010.

[3]高閱春, 何繼強, 姜騰勇, 等. 冠心病患者冠狀動脈病變嚴重程度與冠心病危險因素的相關(guān)分析. 中國循環(huán)雜志, 2012, 27: 178-181.

[4]馮國雙, 陳景武, 張國英. 用GAM 程序擬合光滑樣條非參數(shù)回歸.數(shù)理醫(yī)藥學雜志, 2005, 18: 403-406.

[5]中華人民共和國環(huán)境保護部. 《環(huán)境空氣質(zhì)量標準》（GB3095-2012）. 北京：中國環(huán)境科學出版社, 2012.

[6]宋桂香, 江莉莉, 陳國海, 等. 上海市大氣氣態(tài)污染物與居民每日死亡關(guān)系的時間序列研究. 環(huán)境與健康雜志, 2006, 23: 390-393.

[7]Scarborough P, Allender S, Rayner M, et al. Contribution of climate and air pollution to variation in coronary heart disease mortality rates in England. PLoS One, 2012, 7: e32787. Epub 2012 Mar 12.

[8]王德征, 江國虹, 宋桂德, 等. 求和自回歸滑動平均模型結(jié)合圓分布法分析腦卒中死亡率動態(tài)規(guī)律. 中華流行病學雜志, 2009, 30：82-84.

[9]Choi JH, Xu QS, Park SY, et al . Seasonal variation of eff ect of air pollution on blood pressure. J Epidemiol Community Health, 2007, 61: 314-318.

[10]Yorifuji T, Kawachi I, Sakamoto T, et al. Associations of outdoor air pollution with hemorrhagic stroke mortality. J Occup Environ Med, 2011, 53: 124-126.

[11]Scaife A, Barclay J, Hillis GS, et al. Lack of effect of nitrogen dioxide exposure on heart rate variability in patients with stable coronary heart disease and impaired left ventricular systolic function. Occup Environ Med, 2012, 69: 587-91.

[12]Analitis A, Katsouyanni K, Dimakopoulou K, et al. Short-term effects of ambient particles on cardiovascular and respiratory mortality. Epidemiology, 2006, 17: 230-233.

[13]黃新金, 李菊香, 顏素娟, 等. 冠狀微血管痙攣性心絞痛患者血內(nèi)源性一氧化氮合酶抑制物的變化. 中國循環(huán)雜志, 2007, 22: 100-102.

[14]Knuckles TL, Stapleton PA, Minarchick VC, et al. Air pollution particulate matter collected from an Appalachian mountaintop mining site induces microvascular dysfunction. Microcirculation, 2013, 20: 158-169.

A Time-series Study for Acute Effect of Air Pollution on Mortality in Patients With Cardio-cerebral Vascular Disease in Tianjin City

WANG De-zheng, JIANG Guo-hong, GU Qing, ZHANG Hui, XU Zhong-liang, SONG Gui-de, ZHANG Ying, SHEN Cheng-feng. Tianjin Center for Diseases Control and Prevention, Tianjin (300011), China

Objective: To explore the acute effect of air pollution on mortality for patients with cardio cerebral vascular disease and to provide the basis for disease prevention and control.

Cardio-cerebral vascular disease; Mortality; Generalized additive model; Air pollution; Acute effect

2013-09-16）

（助理編輯：曹洪紅）

國家973計劃課題（2011CB503801）

300011 天津市，天津市疾病預(yù)防控制中心

王德征 主治醫(yī)師 碩士 主要從事慢性病流行病學 Email：wangdezheng1001@163.com**為共同第一作者 通訊作者：顧清

Email：guqing315@126.com

R541

A

1000-3614（2014）06-0453-05

10.3969/j.issn.1000-3614.2014.06.015