田沛茹 ，張斯蒙 ，倪萍 ，艾聰 ，劉曉冰 ，孫秋紅 ，時(shí)景璞 *

慢性阻塞性肺疾?。–OPD）是一種常見(jiàn)的以持續(xù)氣流受限為特征的可以預(yù)防和治療的疾病，其全球發(fā)病率、病死率及致殘率均較高，全球疾病負(fù)擔(dān)研究顯示2016年COPD居于全球疾病死因順位的第4位［1］，預(yù)計(jì)至2020年COPD將成為全球死亡原因的第3位和全球疾病負(fù)擔(dān)的第5位［2］。有關(guān)研究表明，中國(guó)北方農(nóng)村地區(qū)COPD的死亡人數(shù)較多［3］，而目前關(guān)于我國(guó)北方農(nóng)村地區(qū)COPD死亡及過(guò)早死亡的綜合報(bào)道較少。為此，本研究對(duì)遼寧省農(nóng)村地區(qū)2009—2016年COPD死亡及過(guò)早死亡的變化趨勢(shì)進(jìn)行分析及預(yù)測(cè)，以了解我國(guó)北方農(nóng)村地區(qū)COPD死亡特點(diǎn)，為北方農(nóng)村地區(qū)的COPD防控提供更多科學(xué)依據(jù)。

1 資料與方法

1.1 資料來(lái)源 由彰武縣疾病預(yù)防控制中心和鳳城市疾病預(yù)防控制中心負(fù)責(zé)人員提供2009—2016年兩地的全死因監(jiān)測(cè)數(shù)據(jù)，包括地區(qū)、性別、民族、出生日期、死亡日期、最高診斷單位、最高診斷依據(jù)和根本死亡原因等資料。人口資料來(lái)自?xún)傻氐墓膊块T(mén)。根據(jù)國(guó)際疾病分類(lèi)第十版（ICD-10）［4］的標(biāo)準(zhǔn)進(jìn)行疾病分類(lèi)，研究編碼為J44的COPD死亡人群，標(biāo)準(zhǔn)人口采用2010年全國(guó)標(biāo)準(zhǔn)人口［5］。

1.2 研究方法 COPD過(guò)早死亡概率是根據(jù)WHO推薦的壽命表法［6］計(jì)算的30～69歲死于COPD的概率，30～69歲COPD無(wú)條件死亡概率計(jì)算公式為：69q30=1-Π69x=30（1-qx）。

根據(jù)WHO對(duì)慢性非傳染性疾病早死控制要求，2030年重大慢性病過(guò)早死亡概率比2015年降低30%［7］，2030年COPD過(guò)早死亡概率控制目標(biāo)=2015年COPD過(guò)早死亡概率×（1-30%）。

1.3 統(tǒng)計(jì)學(xué)方法 采用Excel錄入數(shù)據(jù)并進(jìn)行數(shù)據(jù)的整理，計(jì)算COPD死亡率、標(biāo)化死亡率、總死亡率、標(biāo)化總死亡率、過(guò)早死亡率、標(biāo)化過(guò)早死亡率、過(guò)早死亡概率、總過(guò)早死亡率、標(biāo)化總過(guò)早死亡率、總過(guò)早死亡概率；采用SPSS 20.0軟件進(jìn)行統(tǒng)計(jì)分析；計(jì)數(shù)資料采用相對(duì)數(shù)表示，組間比較采用χ2檢驗(yàn)；采用Joinpoint Regression Program 4.4.0軟件進(jìn)行Joinpoint回歸計(jì)算年度變化百分比（APC）和其95%CI，并進(jìn)行趨勢(shì)分析［8-9］；采用Matlab軟件建立GM（1，1）灰色模型進(jìn)行死亡率、標(biāo)化死亡率、過(guò)早死亡率、標(biāo)化過(guò)早死亡率的預(yù)測(cè)［10-11］，根據(jù)后驗(yàn)差比值（C）和小誤差概率（P）判斷模型是否合格（C＜0.350或P＞0.950為很好；0.350≤C＜0.500或0.800＜P≤0.950為好；0.500≤C＜0.650或0.700＜P≤0.800為合格；C≥0.650或P≤0.700為不合格），模型很好、好、合格可以用于預(yù)測(cè)，模型不合格則依次減少原始數(shù)據(jù)個(gè)數(shù)進(jìn)行建模，若減少數(shù)據(jù)后模型合格可進(jìn)行預(yù)測(cè)，若依舊不合格則不能進(jìn)行預(yù)測(cè)，預(yù)測(cè)精度用平均絕對(duì)百分誤差（MAPE）進(jìn)行評(píng)價(jià)（本研究中以MAPE ≤10.00%為好）。以P＜0.05為差異有統(tǒng)計(jì)學(xué)意義。

2 結(jié)果

2.1 遼寧省農(nóng)村地區(qū)2009—2016年COPD死亡情況及變化趨勢(shì)

2.1.1 遼寧省農(nóng)村地區(qū)2009—2016年COPD死亡情況 遼寧省農(nóng)村地區(qū)2009—2016年COPD死亡人數(shù)為1 912人，COPD總死亡率、標(biāo)化總死亡率分別為24.14/10萬(wàn)、21.24/10萬(wàn)，其中女性COPD死亡率、標(biāo)化死亡率大于男性，差異均有統(tǒng)計(jì)學(xué)意義（χ2=30.80、9 999.27，P值均＜0.001）；男性COPD死亡率、總死亡率均在2013年達(dá)到最大值，分別為28.16/10萬(wàn)和29.16/10萬(wàn)，而女性COPD死亡率在2011年達(dá)到最大值，為35.47/10萬(wàn)（見(jiàn)表1）。

2.1.2 遼寧省農(nóng)村地區(qū)2009—2016年不同性別COPD死亡變化趨勢(shì) 當(dāng)Joinpoint=0時(shí)，2009—2016年女性COPD標(biāo)化死亡率和COPD標(biāo)化總死亡率均隨年份增加呈下降趨勢(shì)（P＜0.05），年均下降速度分別為14.3%、11.3%；當(dāng)Joinpoint=1時(shí)，2011—2016年女性COPD標(biāo)化死亡率和COPD標(biāo)化總死亡率均隨年份增加呈下降趨勢(shì)（P＜0.05），年均下降速度分別為18.0%、21.4%（見(jiàn)表2）。

2.2 遼寧省農(nóng)村地區(qū)2009—2016年COPD過(guò)早死亡情況及變化趨勢(shì)

2.2.1 遼寧省農(nóng)村地區(qū)2009—2016年COPD過(guò)早死亡情況遼寧省農(nóng)村地區(qū)2009—2016年COPD過(guò)早死亡人數(shù)為389人，占COPD死亡人數(shù)的20.35%，COPD總過(guò)早死亡率和標(biāo)化總過(guò)早死亡率分別為8.33/10萬(wàn)和7.87/10萬(wàn)。不同性別COPD過(guò)早死亡率、標(biāo)化過(guò)早死亡率比較，差異均無(wú)統(tǒng)計(jì)學(xué)意義（χ2=0.47，P=0.490；χ2=0.08，P=0.784）；COPD總過(guò)早死亡概率從2009年的0.83%降到2016年的0.77%；以2016年為基準(zhǔn)，男性COPD過(guò)早死亡概率距離2030年控制目標(biāo)（1.00%）尚有0.58%的下降空間，女性COPD過(guò)早死亡概率距離2030年控制目標(biāo)（1.12%）尚有0.28%的下降空間（見(jiàn)表 3）［7］。

2.2.2 遼寧省農(nóng)村地區(qū)2009—2016年不同性別人群COPD過(guò)早死亡變化趨勢(shì) 當(dāng)Joinpoint=0時(shí)，2009—2016年男性COPD標(biāo)化過(guò)早死亡率、COPD標(biāo)化總過(guò)早死亡率均隨年份增加呈下降趨勢(shì)（P＜0.05），年均下降速度分別為9.5%、12.1%；當(dāng)Joinpoint=1時(shí)，2011—2016年女性COPD標(biāo)化過(guò)早死亡率、COPD總過(guò)早死亡率、標(biāo)化總過(guò)早死亡率均隨年份增加呈下降趨勢(shì)（P＜0.05），年均下降速度為22.8%（見(jiàn)表4）。2.3 遼寧省農(nóng)村地區(qū)2017—2020年COPD死亡和過(guò)早死亡情況預(yù)測(cè) 根據(jù)2009—2016年死亡數(shù)據(jù)建立GM（1，1）模型，女性COPD標(biāo)化死亡率和COPD標(biāo)化總死亡率模型合格可進(jìn)行預(yù)測(cè)，其均隨年份增加呈下降趨勢(shì)（P＜0.001）。而男性COPD死亡率、標(biāo)化死亡率、過(guò)早死亡率、標(biāo)化過(guò)早死亡率、女性COPD死亡率、過(guò)早死亡率、標(biāo)化過(guò)早死亡率、COPD總死亡率、總過(guò)早死亡率和標(biāo)化總過(guò)早死亡率模型均不合格，采用減少原始數(shù)據(jù)個(gè)數(shù)的方法進(jìn)行建模預(yù)測(cè)，重新建模后顯示，男性COPD死亡率、過(guò)早死亡率、COPD總死亡率隨年份的增加呈上升趨勢(shì)（P＜0.001），男性COPD標(biāo)化死亡率、標(biāo)化過(guò)早死亡率、女性COPD死亡率、過(guò)早死亡率、標(biāo)化過(guò)早死亡率、COPD總過(guò)早死亡率、標(biāo)化總過(guò)早死亡率均隨年份增加呈下降趨勢(shì)（P＜0.001），各模型的MAPE為1.18%～8.76%（見(jiàn)表5）。

表 2 遼寧省農(nóng)村地區(qū)2009—2016年不同性別COPD死亡變化趨勢(shì)（%）Table 2 The trends of COPD death in different genders in rural areas of Liaoning province from 2009 to 2016

表 3 遼寧省農(nóng)村地區(qū)2009—2016年COPD過(guò)早死亡情況Table 3 COPD premature death in rural areas of Liaoning province from 2009 to 2016

表 1 遼寧省農(nóng)村地區(qū)2009—2016年COPD死亡情況（1/10萬(wàn)）Table 1 COPD death in rural areas of Liaoning province from 2009 to 2016

3 討論

COPD是世界上嚴(yán)重影響人們健康的疾病之一，每天有數(shù)億人承受COPD帶來(lái)的痛苦和危害，世界各國(guó)COPD死亡率均有差別［12］。本研究結(jié)果顯示，遼寧省農(nóng)村地區(qū)2009—2016年COPD死亡率為24.14/10萬(wàn)，高于我國(guó)部分地區(qū)［13-14］，而低于美國(guó)等發(fā)達(dá)國(guó)家［15-16］，可能由于地理分布、人口構(gòu)成、自然環(huán)境和社會(huì)環(huán)境的不同使得接觸COPD危險(xiǎn)因素的機(jī)會(huì)不同。

表4 遼寧省農(nóng)村地區(qū)2009—2016年不同性別人群COPD過(guò)早死亡變化趨勢(shì)（%）Table 4 The trends of COPD premature death in different genders in rural areas of Liaoning province from 2009 to 2016

表 5 遼寧省農(nóng)村地區(qū)2017—2020年COPD死亡率和過(guò)早死亡率GM（1，1）灰色預(yù)測(cè)模型預(yù)測(cè)結(jié)果Table 5 The mortality and premature mortality GM（1，1） grey prediction results of COPD in rural areas of Liaoning province from 2017 to 2020

本研究采用Joinpoint模型對(duì)2009—2016年COPD死亡情況進(jìn)行趨勢(shì)分析，通過(guò)模型進(jìn)行數(shù)據(jù)擬合，把趨勢(shì)變化分成若干有統(tǒng)計(jì)學(xué)意義的趨勢(shì)區(qū)段，較傳統(tǒng)的人為分段更合理和科學(xué)，特別適用于對(duì)慢性非傳染性疾病發(fā)展的趨勢(shì)研究［17-18］。本研究結(jié)果顯示，2009—2016年女性COPD標(biāo)化死亡率和COPD標(biāo)化總死亡率均隨年份增加呈下降趨勢(shì)，年均下降水平分別為14.3%和11.3%，2011—2016年女性COPD標(biāo)化死亡率和COPD標(biāo)化總死亡率均隨年份增加呈下降趨勢(shì)，年均下降水平分別為18.0%和21.4%。國(guó)內(nèi)外學(xué)者調(diào)查也表明COPD死亡率有下降趨勢(shì)［19-20］，可能是因?yàn)榻陙?lái)國(guó)家經(jīng)濟(jì)實(shí)力增強(qiáng)，醫(yī)療水平和社會(huì)經(jīng)濟(jì)地位的提高［21］，減少了COPD的死亡危險(xiǎn)。

本研究結(jié)果顯示，2009—2016年男性COPD標(biāo)化過(guò)早死亡率和COPD標(biāo)化總過(guò)早死亡率隨年份的增加均呈下降趨勢(shì)，年均下降速度分別為9.5%和12.1%，2011年以后，女性COPD標(biāo)化過(guò)早死亡率也隨年份的增加呈明顯下降趨勢(shì)，年均下降速度為22.8%，與全國(guó)報(bào)道結(jié)果相接近［22-23］。另外在本地區(qū)8年間COPD過(guò)早死亡人數(shù)遠(yuǎn)小于總死亡人數(shù)，COPD過(guò)早死亡人數(shù)占COPD死亡人數(shù)的20.35%，可能是COPD發(fā)病、死亡年齡均較晚，＞69歲高齡死亡人數(shù)多于≤69歲的死亡人數(shù)，而高齡是COPD的危險(xiǎn)因素之一，應(yīng)重點(diǎn)關(guān)注。

本研究結(jié)果顯示，2009—2016年COPD女性死亡率和標(biāo)化死亡率大于男性，而不同性別COPD過(guò)早死亡率和標(biāo)化過(guò)早死亡率無(wú)差異，可能是因?yàn)镃OPD對(duì)本地區(qū)高齡群體中女性危害較大。室內(nèi)空氣污染是農(nóng)村女性COPD死亡的主要危險(xiǎn)因素［24］，農(nóng)村多采用釋放有害煙霧的生物燃料作為能源，且女性多操持家務(wù)，接觸引起室內(nèi)空氣污染的生物燃料機(jī)會(huì)多于男性。此外，2030年COPD過(guò)早死亡概率控制目標(biāo)顯示男性相比于女性需要下降的空間更大，提示未來(lái)需加強(qiáng)對(duì)男性COPD過(guò)早死亡危險(xiǎn)因素的監(jiān)測(cè)，減少男性COPD過(guò)早死亡的危險(xiǎn)因素。因此，應(yīng)減少高齡女性接觸生物燃料的機(jī)會(huì)，對(duì)30～69歲青壯年男性進(jìn)行控?zé)熃錈煹脑缙陬A(yù)防［25］，這對(duì)于實(shí)現(xiàn)2030年COPD過(guò)早死亡概率目標(biāo)有積極推動(dòng)作用。

GM（1，1）模型是基于已有的原始數(shù)據(jù)，預(yù)測(cè)未來(lái)的情況，理論可靠，對(duì)原始數(shù)據(jù)要求不高，對(duì)小樣本量和概率分布不明顯的數(shù)據(jù)預(yù)測(cè)效果較好，在疾病監(jiān)測(cè)中有重要意義［26-27］。本研究結(jié)果顯示，各模型的MAPE為1.18%～8.76%，預(yù)測(cè)精度較高，2017—2020年男性COPD標(biāo)化死亡率、男性COPD標(biāo)化過(guò)早死亡率、女性COPD死亡率、標(biāo)化死亡率、過(guò)早死亡率、標(biāo)化過(guò)早死亡率、COPD標(biāo)化總死亡率、總過(guò)早死亡率、標(biāo)化總過(guò)早死亡率將呈下降趨勢(shì)，但是男性COPD死亡率、過(guò)早死亡率、COPD總過(guò)早死亡率有增長(zhǎng)趨勢(shì)，應(yīng)引起重視，及時(shí)進(jìn)行疾病監(jiān)測(cè)。雖然COPD對(duì)農(nóng)村居民的總體死亡影響呈下降趨勢(shì)，但是受人口老齡化的影響［28］，COPD所導(dǎo)致的疾病危害仍居高不下，COPD仍是居民健康不容忽視的威脅。

綜上所述，雖然遼寧省農(nóng)村地區(qū)2009—2016年COPD死亡率和COPD過(guò)早死亡率總體呈下降趨勢(shì)，對(duì)未來(lái)預(yù)測(cè)顯示也將呈下降趨勢(shì)，但是COPD死亡形勢(shì)依然十分嚴(yán)峻，仍是農(nóng)村地區(qū)不容忽視的危險(xiǎn)疾病，應(yīng)分別對(duì)本地區(qū)高齡女性和青壯年男性采取有針對(duì)性的預(yù)防措施，提倡健康生活方式，減少居民接觸危險(xiǎn)因素機(jī)會(huì)［29］。早期預(yù)防COPD，開(kāi)展農(nóng)村人群健康教育，有利于從根本上減少COPD造成的危害。本研究所涉及數(shù)據(jù)僅是遼寧省部分農(nóng)村地區(qū)，尚不能完全涵蓋整個(gè)遼寧農(nóng)村地區(qū)，日后應(yīng)收集更多相關(guān)數(shù)據(jù)做出更具有說(shuō)服力的研究結(jié)果。本研究趨勢(shì)分析時(shí)間跨度較小，后續(xù)將不斷加入新數(shù)據(jù)，建立長(zhǎng)期的COPD死亡監(jiān)測(cè)系統(tǒng)，及時(shí)進(jìn)行動(dòng)態(tài)模型預(yù)測(cè)，提高預(yù)測(cè)精度。此外，模型預(yù)測(cè)反映的是數(shù)據(jù)的規(guī)律性，尚不能完全反映社會(huì)因素等非規(guī)律性因素對(duì)數(shù)據(jù)的影響，有待結(jié)合實(shí)際情況進(jìn)一步改進(jìn)和研究。

作者貢獻(xiàn)：田沛茹、時(shí)景璞進(jìn)行文章的構(gòu)思與設(shè)計(jì)、研究的實(shí)施與可行性分析，負(fù)責(zé)文章的質(zhì)量控制及審校，對(duì)文章整體負(fù)責(zé)，監(jiān)督管理；田沛茹進(jìn)行數(shù)據(jù)收集、結(jié)果的分析與解釋?zhuān)珜?xiě)論文；田沛茹、張斯蒙、倪萍、艾聰、劉曉冰、孫秋紅進(jìn)行數(shù)據(jù)整理、統(tǒng)計(jì)學(xué)處理、論文的修訂。

本文無(wú)利益沖突。